Деление полиномов без класов и векторов

@matviiv · Регистрация: 09.02.2021

Студворк — интернет-сервис помощи студентам

Делаю калькулятор для полиномов в будущем пригодиться и вот не получається деление в чем ошибка и как исправить?

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
void division(int* arr1, int n, int* arr2, int m, int* result, int& nm) {
    nm = n + m;
    for (int i = 0; i <= nm; i++)
    {
        result[i] = 0;
    }
    for (int i = 0; i <= nm; i++)
    {
        for (int k = 0; k <= m; k++)
        {
            if (i - k >= 0 && i - k <= n)
                result[i] = arr1[i - k] / arr2[k];
        }
    }
}

@XLAT · 21.04.2021, 21:59

Сообщение от matviiv

так вот что такое asign

не стесняйтесь пользоваться гуглом:
https://en.cppreference.com/w/... tor/assign
http://espressocode.top/vector-assign-in-c-stl/

@matviiv · 21.04.2021, 22:04 **[ТС]**

XLAT, я вот переделивал под свою функцию ну и опять что-то такое

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
if (n >= m) {
        while (n >= m) {
            // d equals D shifted right
            d.assign(d.size(), 0);
 
            for (int i = 0; i <= m; i++)
                arr1[i + n - m] = arr2[i];
            int* d = new int[n];
 
            cout<< d;
            int* q = new int[n];
            q[n - m] = arr1[n / arr2[m]];
 
            cout<< q;
 
            // d equals d * q[dN-dD]
            for (int i = 0; i < n-m + 1; i++)
                d[i] = d[i] * q[n - m];
 
            cout<< d;
            int* N = new int[m];
            // N equals N - d
            for (int i = 0; i < n + 1; i++)
                N[i] = N[i] - d[i];
            n--;
 
            cout<< N;
            cout << "-----------------------" << endl << endl;
 
        }
    }

Добавлено через 3 минуты
XLAT, просто под мои функции для калькулятора не подходит

@XLAT · 21.04.2021, 22:08

Сообщение от matviiv

ну и опять что-то такое

мыслим логически:
получилось "что-то такое", как следствие вашей переделки.
Решение: не нужно переделывать!

Сообщение от matviiv

просто под мои функции для калькулятора не подходит

короче, ваш калькулятор колим на две части:
1. интерфейсную часть
2. чёрный ящик с алгоритмом. то бишь сию функцию.

вам нужно передать в функцию из демо два ИНИЦИАЛИЗИРОВАННЫХ вектора.
и вернуть результат то же вектор.
но это ж мелочь...

@matviiv · 21.04.2021, 22:16 **[ТС]**

XLAT, ну просто мне нужно без вектора,я хочу понять именно так,я даже думал через схему горнера но не получаеться

Добавлено через 1 минуту
XLAT,

C++
1
2
3
4
5
6
7
8
9
10
void division(int* arr1, int n, int* arr2, int m, int* result, int& nm) {
    for (int j = 0; j <= m; j++) {
        for (int i = 0; i <= n; i++)
        {
            result[i] = arr1[i] + arr2[i] * result[i - 1];
            cout << result[i] << "  ";
        }
        cout << endl;
    }
}

может как-то так можно?

Добавлено через 2 минуты
XLAT, просто я уже 2 день сижу,я уже хочу так доделать ето действие и суперпозицию что просто уже не могу

@XLAT · 21.04.2021, 22:27

divide

Сообщение от matviiv

ну просто мне нужно без вектора

а я обожаю векторы и без них даже на дискотеки не хожу.

вот ваша функция на языке С++:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
/// C++17
///----------------------------------------------------------------------------|
/// Есть 2 полинома например        х^3 - 12*х^2 - 42
///                        и        х   -  3
///                        и должно x^2 -  9*x - 27
#include <iostream>
#include <iterator>
#include <vector>
 
using namespace std;
typedef vector<double> Poly;
 
// does:  prints all members of vector
// input: c - const string with the name of the vector
//        A - reference to polynomial (vector)
void Print(const char* name, const Poly &A)
{   cout << name << "(" << A.size()-1 << ") = [ ";
    copy(A.begin(), A.end(), ostream_iterator<decltype(A[0])>(cout, " "));
    cout << "]\n";
}
 
///----------------------------------------------------------------------------|
/// Ваша функция.(идентификатор есть глагол, а не существительное)
///----------------------------------------------------------------------------:
Poly divide(Poly& N, Poly& D)
{   Poly d, q, r, err;        // vectors - N / D == q && N % D == r
    size_t dN = N.size() - 1,
           dD = D.size() - 1,
           dd, dq, dr; // degrees of vectors
    size_t i;                  // loop counter
 
    // setting the degrees of vectors
    cout << "Enter the degree of N: " << dN;
    cout << "Enter the degree of D: " << dD;
    dq = dN-dD;
    dr = dN-dD;
 
    if( dD < 1 || dN < 1 )
    {   cerr << "Error: degree of D and N must be positive.\n";
        return err;
    }
 
    // allocation and initialization of vectors
    N.resize(dN+1);
    cout << "Enter the coefficients of N:"<<endl;
    for ( i = 0; i <= dN; i++ )
    {   cout << "N[" << i << "]= " << N[i];
    }
 
    D.resize(dN+1);
    cout << "Enter the coefficients of D:"<<endl;
    for ( i = 0; i <= dD; i++ )
    {   cout << "D[" << i << "]= " << D[i];
    }
 
    d.resize(dN+1);
    q.resize(dq+1);
    r.resize(dr+1);
 
    cout << "-- Procedure --" << endl << endl;
    if( dN >= dD )
    {   while(dN >= dD)
        {   // d equals D shifted right
            d.assign(d.size(), 0);
 
            for( i = 0 ; i <= dD ; i++ )
                d[i+dN-dD] = D[i];
            dd = dN;
 
            Print( "d", d );
 
            // calculating one element of q
            q[dN-dD] = N[dN]/d[dd];
 
            Print( "q", q );
 
            // d equals d * q[dN-dD]
            for( i = 0 ; i < dq + 1 ; i++ )
                d[i] = d[i] * q[dN-dD];
 
            Print( "d", d );
 
            // N equals N - d
            for( i = 0 ; i < dN + 1 ; i++ )
                N[i] = N[i] - d[i];
            dN--;
 
            Print( "N", N );
            cout << "-----------------------" << endl << endl;
 
        }
    }
 
    // r equals N
    for( i = 0 ; i <= dN ; i++ )
        r[i] = N[i];
 
    cout << "=========================" << endl << endl;
    cout << "-- Result --" << endl << endl;
 
    Print( "q", q );
    Print( "r", r );
 
    return q;
}
 
///----------------------------------------------------------------------------|
/// Старт.
///----------------------------------------------------------------------------:
int main()
{
    Poly a = {-42, 0, -12, 1};
    Poly b = {-3, 1};
 
    Poly res = divide(a, b);
 
    Print( "res: ", res );
}

внутри осталось убрать демо-лог, шоп ничаво там не отсвечивало...

@matviiv · 21.04.2021, 23:48 **[ТС]**

XLAT, ну я не про етот код,но ладно

Добавлено через 33 минуты
XLAT,

C++
1
2
3
4
5
6
7
8
9
10
11
12
void division(int* arr1, int n, int* arr2, int m, int* result, int& nm) {
    for (int j = 0; j <= (n - m); j++)
    {
        result[j] = arr1[n - j] / arr2[m];
 
        for (int i = 0; i <= n; i++)
        {
            if (i >= (n - (m + j)))
                arr1[i] = (arr1[i] - (arr2[i + j - (n - m)] * result[j]));
        }
    }
}

еще была такая попытка но не получается, я уже не знаю что делать, хотел завтра утром показать преподавателю (мы с ним общаемся дополнительно, готовлюсь к олимпиаде)

Добавлено через 46 минут
XLAT, поможете?чесно я уже задолбался,как-то сумма ,разница,и умножение с схемой горнера легко а вот ети 2 подпункта уже так достали

@XLAT · 21.04.2021, 23:49

если чуть ужать:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Poly divide(Poly& N, Poly& D)
{   size_t dN = N.size(),
           dD = D.size();
 
    const size_t dq = dN-dD--;
 
    N.resize(dN);
    D.resize(dN);
 
    Poly d(dN--);
    Poly q(dq+1);
 
    for(size_t a = dq; dN >= dD; dN--, a--)
    {
        for(size_t i = 0; i <= dD; i++) d[i+a] = D[i];
        q[a] = N[dN]/d[dN];
        for(size_t i = 0; i <= dq; i++) d[i]   = d[i] * q[a];
        for(size_t i = 0; i <= dN; i++) N[i]   = N[i] - d[i];
    }
 
    /// N.resize(dN +1); 
    /// Print( "r", N );
    return q;
}

@matviiv · 22.04.2021, 01:14 **[ТС]**

XLAT,

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
int division(int* arr1, int n, int* arr2, int m, int* result, int& nm) {
    size_t dN = n,
        dD = m;
 
    const size_t dq = dN - dD--;
    int* d=new int[dN--];
    int* q=new int[dq + 1];
 
    for (size_t a = dq; dN >= dD; dN--, a--)
    {
        for (size_t i = 0; i <= dD; i++) d[i + a] = arr2[i];
        q[a] = arr1[dN] / arr2[dN];
        for (size_t i = 0; i <= dq; i++) d[i] = d[i] * q[a];
        for (size_t i = 0; i <= dN; i++) arr1[i] = arr1[i] - d[i];
    }
    return* q;
}

Почему то visual studio не пропускает

@XLAT · 22.04.2021, 07:38

Сообщение от matviiv

Почему то visual studio не пропускает

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
int* division(int* arr1, int n, int* arr2, int m)
{   size_t dN = n,
           dD = m;
 
    const size_t dq = dN - dD--;
    int* d=new int[dN--];
    int* q=new int[dq + 1];
 
    for (size_t a = dq; dN >= dD; dN--, a--)
    {   for (size_t i = 0; i <= dD; i++) d[i + a] = arr2[i];
        q[a] = arr1[dN] / arr2[dN];
        for (size_t i = 0; i <= dq; i++) d[i] = d[i] * q[a];
        for (size_t i = 0; i <= dN; i++) arr1[i] = arr1[i] - d[i];
    }
    return q;
}

@matviiv · 22.04.2021, 08:29 **[ТС]**

XLAT, опять не пропускает

@XLAT · 22.04.2021, 08:38

matviiv, ладно вот по схеме Горнера:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
/// C++17
///----------------------------------------------------------------------------|
/// Есть 2 полинома например        х^3 - 12*х^2 - 42
///                        и        х   -  3
///                        и должно x^2 -  9*x - 27
#include <iostream>
 
using namespace std;
 
void Print(const char* name, int* m, size_t n)
{   cout << name << "(" << n-1 << ") = [ ";
    for(size_t i = 0; i < n; ++i) std::cout << m[i] << ' ';
    cout << "]\n";
}
 
///----------------------------------------------------------------------------|
/// Деление полинома на полином первого порядка по схеме Горнера.
///----------------------------------------------------------------------------:
int* divide(int* arr1, int n1, int* arr2, int n2)
{
    int n3 = n1 - 1;
    int* res = new int[n3];
 
    res[0] = arr1[0];
    int a = -arr2[1];
 
    for(int i = 1; i < n1; ++i)
    {   res[i] = a * res[i-1] + arr1[i];
    }
    return res;
}
 
///----------------------------------------------------------------------------|
/// Старт.
///----------------------------------------------------------------------------:
int main()
{   int a[] = {1, -12, 0, -42};
    int b[] = {1, -3};
 
    Print("a: ", a, 4);
    Print("b: ", b, 2);
 
    int* res = divide(a, 4, b, 2);
 
    Print("res: ", res, 3);
}

Сообщение от matviiv

опять не пропускает

излучение 5G наводит помехи на мой волшебный шар...
но видно, что есть вероятность неудачной копипасты...

@matviiv · 22.04.2021, 08:47 **[ТС]**

XLAT,

@XLAT · 22.04.2021, 10:48

matviiv, вот всё до кучи:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
/// C++17
///----------------------------------------------------------------------------|
/// Есть 2 полинома например        х^3 - 12*х^2 - 42
///                        и        х   -  3
///                        и должно x^2 -  9*x - 27
#include <iostream>
#include <vector>
 
using namespace std;
typedef vector<int> Poly;
 
void Print(const char* name, int* m, size_t n)
{   cout << name << "(" << n-1 << ") = [ ";
    for(size_t i = 0; i < n; ++i) std::cout << m[i] << ' ';
    cout << "]\n";
}
 
#define  l(v)  std::cout << #v << " = " << (v) << "\n";
 
///----------------------------------------------------------------------------|
/// Деление полинома на полином первого порядка по схеме Горнера.
///----------------------------------------------------------------------------:
int* divide(int* arr1, int n1, int* arr2, int n2)
{   int n3 = n1 - 1;
    int* res = new int[n3];
 
    res[0] = arr1[0];
    int a = -arr2[1];
 
    for(int i = 1; i < n1; ++i)
    {   res[i] = a * res[i-1] + arr1[i];
    }
    return res;
}
 
///----------------------------------------------------------------------------|
/// STL Euclidean division of polynomials.
///----------------------------------------------------------------------------:
Poly divide(Poly& N, Poly& D)
{   size_t dN = N.size(),
           dD = D.size();
 
    const size_t dq = dN-dD--;
 
    Poly d(dN--, 0);
    Poly q(dq+1, 0);
 
    for(size_t a = dq; dN >= dD; dN--, a--)
    {
        for(size_t i = 0; i <= dD; i++) d[i+a] = D[i];
        q[a] = N[dN] / d[dN];
        for(size_t i = 0; i <= dq; i++) d[i]   = d[i] * q[a];
        for(size_t i = 0; i <= dN; i++) N[i]   = N[i] - d[i];
    }
    return q;
}
 
///----------------------------------------------------------------------------|
/// Euclidean division of polynomials.
///----------------------------------------------------------------------------:
int* division(int* arr1, int n, int* arr2, int m)
{   size_t dN = n,
           dD = m;
 
    const size_t dq = dN - dD--;
    int* d=new int[dN--  ]; for(int i = n; i--;) d[i] = 0;
    int* q=new int[dq + 1];
 
    for (size_t a = dq; dN >= dD; dN--, a--)
    {   for (size_t i = 0; i <= dD; i++) d[i + a] = arr2[i];
        q[a] = arr1[dN] / d[dN];
        for (size_t i = 0; i <= dq; i++) d[i] = d[i] * q[a];
        for (size_t i = 0; i <= dN; i++) arr1[i] = arr1[i] - d[i];
    }
    return q;
}
 
///----------------------------------------------------------------------------|
/// Старт.
///----------------------------------------------------------------------------:
int main()
{
    ///--------------------------------|
    /// Gorner.                        |
    ///--------------------------------:
    {
        ///        x^3  x^2  x^1  x^0
        int a[] = {1,  -12,  0,  -42};
        int b[] = {1,  -3};
 
        Print("  a: ", a, 4);
        Print("  b: ", b, 2);
 
        int* Gor = divide(a, 4, b, 2);
        Print("Gorner        : ", Gor, 3);
 
        delete[] Gor;
    }
 
    ///--------------------------------|
    /// Euclidean.                     |
    /// Порядок коэффициентов обратный.| <---(!)
    ///--------------------------------:
    {
        int a[] = {-42, 0, -12, 1, 0};
        int b[] = {-3, 1, 0};
 
        Print("\n  a: ", a, 4);
        Print("  b: "  , b, 2);
 
        int* Euc = division(a, 4, b, 2);
        Print("Euclidean C :  ", Euc, 3);
 
        delete[] Euc;
    }
 
    ///--------------------------------|
    /// Euclidean C++.                 |
    /// Порядок коэффициентов обратный.| <---(!)
    ///--------------------------------:
    {
        Poly   a = {-42, 0, -12, 1};
        Poly   b = {-3, 1};
        Print( "\n  a: ", &a[0], a.size());
        Print( "  b: "  , &b[0], b.size());
 
        Poly res = divide(a, b);
 
        Print( "Euclidean STL: ", &res[0], res.size());
    }
}

думаю вопросов по redefinition multiple initialization у вас не должно уже быть...

@matviiv · 22.04.2021, 14:52 **[ТС]**

XLAT,

C++
1
2
3
4
5
6
7
8
9
10
11
void division(int* arr1, int n, int* arr2, int m, int* result) {
    for (int k = (n - m); k >= 0; k--)
    {
        result[k] = arr1[k + m] / arr2[m];
 
        for (int j = k + m - 1; j >= k; j--)
        {
            arr1[j] -= arr2[j - k] * result[k];
        }
    }
}

Рабочая функция получилась вот,спасибо что мучался со мной целий день

@SergioO · 23.04.2021, 13:58

Сообщение от matviiv

что такое asign

позволю ответить: обычно assign, в ЯП , это метод который добавляет нечто, не изменяя исходного объекта, просто добавдвляеет функциональность

Новые блоги и статьи Все статьи Все блоги /
Учебный мини "цифровой двойник" = (Сервер = Java + Flask_питон )+(3D = Гугловскаая 4T nvideo + google colab) +( WMS = Odoo на Docker) anaschu 07.07.2026 Увидел в требованиях к вакансии создателя имитационных моделей связку (REst + Api + WMS + Anylogic + питон + 3D Nvideo), и решил посмотреть, что это такое. Из этого всего я работал более менее. . .	Море в объективе kumehtar 06.07.2026 Хотел написать много проникновенных слов, но передумал. Оно само - говорит. Кто слышит тот слышит.	Doom для терминала без стрельбы и монстров. 3D Raycasting на ascii. dcc0 05.07.2026 Попросил нейронную сеть deepai. org написать рейкастинг 3D с библиотекой ncurses для Linux. Чтобы можно было ходить на стрелочки. Чтобы стены были отрисованы символами. Справилась. Первый вариант. . .	Установка статуса документа по условию Maks 05.07.2026 Алгоритм из решения ниже реализован на нетиповом документе "НарядПутевка" разработанного в КА2. Задача: в табличной части "Материалы" документа при записи автоматически устанавливать статус. . .
Сезонность и суточность закисления почв anaschu 04.07.2026 200 часов это все равно моловато. Есть ситуации, но нестандартные, когда смена происходит за 5 лет. Но обычно это 50 лет и более. Наверное, закисление почвы происходит сезонно в средней. . .	В чем ценность человеческого опыта в глобальном смысле? kumehtar 03.07.2026 Возможно, ценность человека не в том, что он однажды достигает мудрости, а в том, что он становится носителем карты пути. Он знает не только истину, но и последовательность внутренних изменений,. . .	интеграция AnyLogic с самописным REST API и переход на Odoo anaschu 03.07.2026 Успешная интеграция AnyLogic с самописным REST API и переход на промышленную Odoo WMS Сегодня проделал огромный путь от простой симуляции физических процессов до построения полноценной. . .	Поиск всех путей на ориентированном графе. Linux dcc0 02.07.2026 Переработка старого кода из моей статьи. Через несколько переработок от PHP кода к C89 (надеюсь, 89). Но довольно запутанно получилось. Код для Linux. Но если убрать time и то, что с ним. . .

@matviiv 429 / 313 / 129 Регистрация: 09.02.2021 Сообщений: 1,096 Записей в блоге: 1
	22.04.2021, 08:29 [ТС]
	XLAT, опять не пропускает 0

@matviiv 429 / 313 / 129 Регистрация: 09.02.2021 Сообщений: 1,096 Записей в блоге: 1
	22.04.2021, 08:47 [ТС]
	XLAT, 0