Суперпозиция полиномов без классов

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

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

Привет всем, захотел написать калькулятор для полиномов и нашел интересный раздел суперпозиции, поэтому написал вот такую функцию но не получилось, в чем ошибка? P.s. код пишу без классов(если нужно то ввод полиномов такой например:1+2х+3х^2 и тд.)

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
void s(int* mas1, int n, int* mas2, 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++){
                result[i] += mas1[i - k] * pow(mas2[k],i);
        }
    }
}

@Volga_ · 21.04.2021, 09:34

По математике:
$https://www.cyberforum.ru/cgi-bin/latex.cgi?P(x)=a_0+a_1x+...+a_nx^n$
$https://www.cyberforum.ru/cgi-bin/latex.cgi?Q(x)=b_0+b_1x+...+b_mx^m$
$https://www.cyberforum.ru/cgi-bin/latex.cgi?Result(x)=P(x)Q(x)=\displaystyle\sum_{k=0}^{n+m}\displaystyle\sum_{i+j=k}^{} a_ib_jx^k$
Поэтому получим код снизу (переделаю ваш код):

C++
1
2
3
4
5
6
7
8
9
10
11
12
void s(int* a, int n, int* b, int m, int* result, int& nm) 
{
    nm = n + m;
    for (int k = 0; k <= nm; k++)
    {
        result[k] = 0;
        for (int i = 0; i <= n; i++)
            for (int j = 0; j <= m; j++)
                if (i + j == k)
                    result[k] += a[i] * b[j];
    }
}

Попробуйте !

Или можно другой вариант с тестом:

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
#include <iostream>
using namespace std;
 
int* s(int* a, int n, int* b, int m, int& nm)
{
    nm = n + m;
    int* result = new int[nm + 1];
    if (result)
    {
        for (int k = 0; k <= nm; k++)
        {
            result[k] = 0;
            for (int i = 0; i <= n; i++)
                for (int j = 0; j <= m; j++)
                    if (i + j == k)
                        result[k] += a[i] * b[j];
        }
        return result;
    }
    else
        return nullptr;
}
 
int main()
{
    const int n = 3;
    int* a = new int[n + 1]{ 1, -3, 2, 4 };
    cout << "a:" << endl;
    for (int i = 0; i <= n; i++)
        cout << a[i] << " ";
    cout << endl;
 
    const int m = 2;
    int* b = new int[m + 1]{ -2, 3, 2 };
    cout << "b:" << endl;
    for (int i = 0; i <= m; i++)
        cout << b[i] << " ";
    cout << endl;
 
    int nm = 0;
    int* result = s(a, n, b, m, nm);
    if (result)
    {
        cout << "Result:" << endl;
        for (int i = 0; i <= nm; i++)
            cout << result[i] << " ";
        cout << endl;
        delete[] result;
    }
    else
        cout << "No result !" << endl;
    delete[] a;
    delete[] b;
 
    system("pause");
    return 0;
}

@matviiv · 21.04.2021, 10:12 **[ТС]**

Volga_, извини,я просто проверил то ответ неверний,тоесть я ввел два полинома 1+2х+3х^2 и вишло не 27x4+36x3+36x2+16x+6

@Volga_ · 21.04.2021, 10:38

matviiv, покажите ваш код полный, пожалуйста.

Добавлено через 15 минут

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
#include <iostream>
using namespace std;
 
int* s(int* a, int n, int* b, int m, int& nm)
{
    nm = n + m;
    int* result = new int[nm + 1];
    if (result)
    {
        for (int k = 0; k <= nm; k++)
        {
            result[k] = 0;
            for (int i = 0; i <= n; i++)
                for (int j = 0; j <= m; j++)
                    if (i + j == k)
                        result[k] += a[i] * b[j];
        }
        return result;
    }
    else
        return nullptr;
}
 
int main()
{
    const int n = 2;
    int* a = new int[n + 1]{ 1, 2, 3 };
    cout << "a:" << endl;
    for (int i = 0; i <= n; i++)
    {
        if (i == 0)
            cout << a[i] << "+";
        else if (i == 1)
            cout << a[i] << "x" << "+";
        else if (i == n)
            cout << a[i] << "x^" << i;
        else
            cout << a[i] << "x^" << i << "+";
    }
    cout << endl;
 
    const int m = 2;
    int* b = new int[m + 1]{ 1, 2, 3 };
    cout << "b:" << endl;
    for (int i = 0; i <= m; i++)
    {
        if (i == 0)
            cout << b[i] << "+";
        else if (i == 1)
            cout << b[i] << "x" << "+";
        else if (i == m)
            cout << b[i] << "x^" << i;
        else
            cout << b[i] << "x^" << i << "+";
    }
    cout << endl;
 
    int nm = 0;
    int* result = s(a, n, b, m, nm);
    if (result)
    {
        cout << "Result:" << endl;
        for (int i = 0; i <= nm; i++)
        {
            if (i == 0)
                cout << result[i] << "+";
            else if (i == 1)
                cout << result[i] << "x" << "+";
            else if (i == nm)
                cout << result[i] << "x^" << i;
            else
                cout << result[i] << "x^" << i << "+";
        }
        cout << endl;
        delete[] result;
    }
    else
        cout << "No result !" << endl;
    delete[] a;
    delete[] b;
 
    system("pause");
    return 0;
}

Тест:

Code
1
2
3
4
5
6
7
a:
1+2x+3x^2
b:
1+2x+3x^2
Result:
1+4x+10x^2+12x^3+9x^4
Press any key to continue . . .

Ок или нет ?
Почему вы должны получить 27x4+36x3+36x2+16x+6 из два полинома 1+2х+3х^2 ? Не совсем понял по математике !

@matviiv · 21.04.2021, 10:57 **[ТС]**

Volga_, ну если в 1+2х+3х^2 подставить 1+2х+3х^2 то будет 1+2(1+2х+3х^2)+3(1+2х+3х^2)^2

@Volga_ · 21.04.2021, 15:34

Ага, теперь я понял ваш вопрос. Подумаю. Пожалуйста, снять "лучший ответ" для меня в посте #2.

Добавлено через 4 часа 36 минут
matviiv, честно я не смогу помочь вам если не использую vector (без массива). Это мой вариант:

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
#include <iostream>
#include <vector>
using namespace std;
 
vector<int> sum2Vec(vector<int> vec1, vector<int> vec2)
{
    vector<int> ret;
    int n1 = vec1.size();
    int n2 = vec2.size();
    if (n1 > n2)
    {
        for (int i = 0; i < n1; i++)
        {
            if (i < n2)
                ret.push_back(vec1[i] + vec2[i]);
            else
                ret.push_back(vec1[i]);
        }
    }
    else
    {
        for (int i = 0; i < n2; i++)
        {
            if (i < n1)
                ret.push_back(vec1[i] + vec2[i]);
            else
                ret.push_back(vec2[i]);
        }
    }
    return ret;
}
 
vector<int> mul2Vec(vector<int> vec1, vector<int> vec2)
{
    vector<int> ret;
    int n1 = vec1.size();
    int n2 = vec2.size();
    int n = n1 + n2 - 2;
    for (int k = 0; k < n + 1; k++)
    {
        int retTmp = 0;
        for (int i = 0; i < n1; i++)
            for (int j = 0; j < n2; j++)
                if (i + j == k)
                    retTmp += vec1[i] * vec2[j];
        ret.push_back(retTmp);
    }
    return ret;
}
 
vector<int> mulK_Vec(vector<int> vec, int k)
{
    vector<int> ret;
    int n = vec.size();
    for (int i = 0; i < n; i++)
        ret.push_back(k * vec[i]);
    return ret;
}
 
void PrintResult(vector<int> vec)
{
    int n = vec.size();
    for (int i = 0; i < n; i++)
    {
        if (i == 0)
            cout << vec[i];
        else if (i == 1)
            if (vec[i] < 0)
                cout << vec[i] << "x";
            else
                cout << "+" << vec[i] << "x";
        else
            if (vec[i] < 0)
                cout << vec[i] << "x^" << i;
            else
                cout << "+" << vec[i] << "x^" << i;
    }
    cout << endl;
}
 
int main()
{
    vector<int> a = { 1, 2, 3 };
    int na = a.size();
    cout << "a:" << endl;
    PrintResult(a);
 
    vector<int> b = { 1, 2, 3 };
    cout << "b:" << endl;
    PrintResult(b);
 
    vector<int> mul = { 1 };
    vector<int> result = { a[0] };
    if (na > 1)
        for (int i = 1; i < na; i++)
        {
            mul = mul2Vec(mul, b);
            result = sum2Vec(result, mulK_Vec(mul, a[i]));
        }
    
    cout << "result:" << endl;
    PrintResult(result);
 
    system("pause");
    return 0;
}

Результат:

Code
1
2
3
4
5
6
a:
1+2x+3x^2
b:
1+2x+3x^2
result:
6+16x+36x^2+36x^3+27x^4

Подумаю о вариантах с использованием массива потом.

@matviiv · 21.04.2021, 15:41 **[ТС]**

Volga_, понял,буду ждать

@matviiv · 21.04.2021, 21:47 **[ТС]**

Volga_,

C++
1
2
3
4
5
6
7
8
9
10
11
int superposition(int* arr1, int n, int* arr2, int m, int* result, int& nm) {
    int* h = new int[n*m];
    int k = n * m;
    for (int i = n*m- 2; i >= 0; i--)
    {
        *h = mull(h,n*m, arr2,m,h,k);
        h[0] += arr1[i];
    }
    return *h;
 
}

вот я думал через умножение но чет не получилосб

@Volga_ · 21.04.2021, 21:51

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

вот я думал через умножение но чет не получилосб

Это нелегко. Я не продумаю ничего !

@matviiv · 22.04.2021, 21:15 **[ТС]**

Ну как я понял то без функции умножение и схеми Горнера очень сложно,но с ними почему-то не могу

@Volga_ · 22.04.2021, 21:30

matviiv, при работы с вектором ты не боишься изменения размера ! Но когда с массивом размера массива result изменяется через каждой раз вызова функции. Удаление - выделение памяти массива очень сложно. Проблема у меня именно изменение размера массива ! Это недостатка массива.

@matviiv · 22.04.2021, 21:57 **[ТС]**

Volga_, да тоже столкнулся с етой проблемой

@matviiv · 23.04.2021, 09:30 **[ТС]**

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
void superposition(int* arr1, int n, int* arr2, int m, int* result, int& nm)
{
    int* tmp_arr = new int[1];
    tmp_arr[0] = 0;
    int tmp_n = 0;
 
    for (int i = 0; i <= n; i++)
    {
        int res_n;
        int* res_tmp=new int[maxN];
 
        /* DEBUG PRINT
        cout << "(";
        print(tmp_arr, tmp_n);
        cout << ") * (";
        print(arr2, m);
        cout << ") = ";
        */
        mull(tmp_arr, tmp_n, arr2, m, res_tmp, res_n);
        /*
        print(res_tmp, res_n);
        cout << endl;
        */
        res_tmp[0] += arr1[n - i];
 
        delete[] tmp_arr;
        tmp_arr = res_tmp;
        tmp_n = res_n;
    }
 
    result = tmp_arr;
    nm = tmp_n;
}

Вот,но только в чем-то ошибка,виводит только 2 полином

Добавлено через 12 минут
Исправил,но только виводит мусор,вот как визиваю

C++
1
2
3
4
5
cout << "Result:" << endl;
        int* res = new int[maxN];
        int nm = 0;
        superposition(mass, n, mass1, m, res, nm);
        print(res,nm);

Добавлено через 3 минуты
Все исправил

Новые блоги и статьи Все статьи Все блоги /
Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .	PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .
От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .

@matviiv 429 / 313 / 129 Регистрация: 09.02.2021 Сообщений: 1,096 Записей в блоге: 1
	21.04.2021, 10:12 [ТС]
	Volga_, извини,я просто проверил то ответ неверний,тоесть я ввел два полинома 1+2х+3х^2 и вишло не 27x4+36x3+36x2+16x+6 0

@matviiv 429 / 313 / 129 Регистрация: 09.02.2021 Сообщений: 1,096 Записей в блоге: 1
	21.04.2021, 10:57 [ТС]
	Volga_, ну если в 1+2х+3х^2 подставить 1+2х+3х^2 то будет 1+2(1+2х+3х^2)+3(1+2х+3х^2)^2 0

@matviiv 429 / 313 / 129 Регистрация: 09.02.2021 Сообщений: 1,096 Записей в блоге: 1
	21.04.2021, 15:41 [ТС]
	Volga_, понял,буду ждать 0

@matviiv 429 / 313 / 129 Регистрация: 09.02.2021 Сообщений: 1,096 Записей в блоге: 1
	22.04.2021, 21:15 [ТС]
	Ну как я понял то без функции умножение и схеми Горнера очень сложно,но с ними почему-то не могу 0

@Volga_ Модератор 5208 / 2925 / 1509 Регистрация: 14.12.2018 Сообщений: 5,266 Записей в блоге: 1
	22.04.2021, 21:30
	matviiv, при работы с вектором ты не боишься изменения размера ! Но когда с массивом размера массива result изменяется через каждой раз вызова функции. Удаление - выделение памяти массива очень сложно. Проблема у меня именно изменение размера массива ! Это недостатка массива. 1

@matviiv 429 / 313 / 129 Регистрация: 09.02.2021 Сообщений: 1,096 Записей в блоге: 1
	22.04.2021, 21:57 [ТС]
	Volga_, да тоже столкнулся с етой проблемой 1