Форум программистов, компьютерный форум, киберфорум
С++ для начинающих
Войти
Регистрация
Восстановить пароль
 
Рейтинг 5.00/7: Рейтинг темы: голосов - 7, средняя оценка - 5.00
0 / 0 / 0
Регистрация: 25.12.2016
Сообщений: 39
1

Деление двух многочленов

12.05.2017, 00:01. Показов 1335. Ответов 7
Метки нет (Все метки)

Привет всем! Было дано задание реализовать деление многочленов с комплексными коэффициентами через метод класса. Точнее сказать только с делением и возникли проблемы. Класс комплексных чисел выполнен верно(может не самым рациональным способом). При деление основывал свой алгоритм на этом Деление многочленов . В самом методе было использовано несколько дополнительных функций (без которых возможно можно было обойтись, но оптимизацией я решил заняться только после получения готового и правильного метода).
Проблема: считает неверно при некоторых значениях.
Надо: найти ошибку в алгоритме и/или исправить код (желательно без использования заумностей)
Ниже сам код(неполный, без ненужных данной проблеме строк)

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
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
#include <iostream> 
#include <stack> 
#include <cctype> 
#include <cmath> 
#include <cstdio> 
#include <iomanip> 
#include <map> 
#include <algorithm> 
#include <string> 
#include <cstring> 
#include <algorithm>
using namespace std;
#define druj cout << endl
double proovi(double a) {
    if (a == -0.00000000000000000)
        return 0.00000000000000000;
    else
        return a;
}
class Complex {
private:
    double Re;
    double Im;
public:
    //Конструкторы 
    Complex();
    Complex(double);
    Complex(double, double);
    //Операции 
 
    //кто не понял тот поймет
    Complex& SetParts(double, double);
 
    //Равенство 
    Complex& operator=(const Complex&);
 
    //Сложение 
    Complex operator +(Complex);
 
    //Унарное вычитание 
    Complex operator-();
 
    //Вычитание 
    Complex operator-(Complex);
 
    //Уможение 
    Complex operator*(Complex);
 
    //Деление 
    Complex operator/(Complex);
    //Модуль 
    double Abs();
    //Вывод 
    friend ostream& operator<<(ostream&, const Complex&);
    //Ввод 
    friend istream& operator>>(istream&, Complex&);
    //Re
    double GetRe();
    //Im
    double GetIm();
    //возведение в степень
 
};
//Конструкторы 
Complex::Complex() {
    Re = 0;
    Im = 0;
}
Complex::Complex(double Vvod) {
    Re = Vvod;
    Im = 0;
}
Complex::Complex(double Chisl, double Znam) {
    Re = Chisl;
    Im = Znam;
}
Complex& Complex::operator= (const Complex &c) {
    Re = c.Re;
    Im = c.Im;
    return (*this);
}
//Сложение 
Complex Complex::operator+(Complex a) {
    Complex result;
    result.Re = Re + a.Re;
    result.Im = Im + a.Im;
    return result;
}
//Унарное вычитание 
Complex Complex::operator-() {
    return Complex(-Re, -Im);
}
//Вычитание 
Complex Complex::operator-(Complex a) {
    Complex result;
    result.Re = Re - a.Re;
    result.Im = Im - a.Im;
    return result;
}
//Умножение 
Complex Complex::operator*(Complex a) {
    Complex result;
    result.Re = Re*a.Re - Im*a.Im;//ac-bd 
    result.Im = Im*a.Re + Re*a.Im;//bc+ad 
    return result;
}
//Деление 
Complex Complex::operator/(Complex a) {
    Complex result;
    result.Re = (Re*a.Re + Im*a.Im) / (a.Re*a.Re + a.Im*a.Im);//(ac+bd)/(c*c+d*d) 
    result.Im = (Im*a.Re - Re*a.Im) / (a.Re*a.Re + a.Im*a.Im);//(bc-ad)/(c*c+d*d) 
    if (result.Im == -0.00000000000000000)
        result.Im = 0;
    if (result.Re == -0.00000000000000000)
        result.Re = 0;
 
    return result;
}
//Модуль 
double Complex::Abs() {
    return sqrt(Re*Re + Im*Im);
}
//Вывод 
ostream &operator<<(ostream &out, const Complex &c) {
    out << fixed << setprecision(4) << proovi(c.Re);
    if (c.Im >= 0)
        out << "+";
    out << proovi(c.Im) << fixed << setprecision(4) << "i";
 
    return out;
}
//Ввод 
istream& operator>>(istream &in, Complex& c) {
    in >> c.Re >> c.Im;
    return in;
}
//кто не понял тот поймет
Complex & Complex::SetParts(double a, double b)
{
    Re = a;
    Im = b;
    return (*this);
}
//
double Complex::GetRe()
{
    return Re;
}
double Complex::GetIm()
{
    return Im;
}
 
 
void reverse_(Complex *&a, int count) {
    //if (count % 2 == 1)
    //  count--;
    for (int i = 0; i < count / 2; i++)
        swap(a[i], a[count - i - 1]);
}
Complex last_nonull(Complex *a, int razm) {
    for (int i = razm; i >= 0; i--)
        if (a[i].GetIm() != 0 || a[i].GetRe() != 0)
            return a[i];
    Complex o;
    return o;
}
int new_degree(Complex *a, int old_degree) {
    for (int i = old_degree; i >= 0; i--)
        if (a[i].GetIm() != 0 || a[i].GetRe() != 0)
            return i;
    return 0;
}
Complex* sdvig(Complex *a, int NEW, int old_degree) {
    Complex o;
    Complex *res = new Complex[NEW + 1];
    int temp_k = 0;
    for (int i = NEW - old_degree; i <= NEW; i++) {
        res[i] = a[temp_k];
        temp_k++;
    }
    return res;
}
 
 
class Polinom {
private:
    int degree;
    Complex *koef;
public:
    //Конструктор по умолчанию
    Polinom();
    //Конструктор инициализации
    Polinom(int, Complex*);
    //Конструктор копирования
    Polinom(const Polinom&);
    //Оператор присваивания
    Polinom& operator=(const Polinom&);
    //
    int GetDegree();
    //Деструктор
    ~Polinom();
    //Вывод
    friend ostream& operator<<(ostream&, const Polinom&);
    //увеличение степени
    void uvelich(int);
    // + - / * %
    Polinom operator+(Polinom);
    Polinom operator-(Polinom);
    Polinom operator/(Polinom);
    Polinom operator*(Polinom);
    Polinom operator%(Polinom);
    Complex operator()(Complex);
    Polinom derived(int);
    void SetDegree(int);
    Complex GetKoef(int);
 
};
//Конструктор по умолчанию
Polinom::Polinom()
{
    degree = 0;
    koef = new Complex[1];
    Complex a(0);
    koef[0] = a;
}
//Конструктор инициализации
Polinom::Polinom(int a, Complex *x)
{
    degree = a;
    //if (koef)
    //  delete[] koef;
    koef = new Complex[a + 1];
    for (int i = 0; i <= a; i++) {
        koef[i] = x[i];
    }
}
//Конструктор копирования
///Чек аут
Polinom::Polinom(const Polinom &x)
{
    //if (koef)
    //  delete[] koef;
 
    degree = x.degree;
    koef = new Complex[degree + 1];
    for (int i = 0; i <= degree; i++)
        koef[i] = x.koef[i];
}
//Деструктор
Polinom::~Polinom()
{
    delete[] koef;
}
//Оператор присваивания 
///Чек Аут
Polinom&  Polinom::operator=(const Polinom& x)
{
    //if (koef)
    //  delete[] koef;
    degree = x.degree;
    koef = new Complex[degree + 1];
    for (int i = 0; i <= degree; i++)
        koef[i] = x.koef[i];
    return (*this);
}
//Получение степени многочлена    
int Polinom::GetDegree()
{
    return degree;
}
//Вывод
ostream& operator<<(ostream& out, const Polinom& x) {
    int temp = new_degree(x.koef, x.degree);
    for (int i = temp; i >= 0; i--) {
        if (i == 0)
            out << "(" << x.koef[i] << ")" << "x^" << i;
        else
            out << "(" << x.koef[i] << ")" << "x^" << i << " + ";;
    }
    return out;
}
//Перегрузка +
Polinom Polinom::operator+(Polinom x)
{
    int res_degree = 0;
    if (degree > x.degree)
    {
        res_degree = degree;
        //увеличение степени(выделение места для коэфов)
        x.uvelich(degree);
    }
    else
    {
        res_degree = x.degree;
        //
        uvelich(x.degree);
    }
    Complex *res_koef = new Complex[res_degree + 1];
    for (int i = 0; i <= res_degree; i++)
        res_koef[i] = koef[i] + x.koef[i];
    for (int i = res_degree; i >= 0; i--) {
        if (res_koef[i].GetIm() != 0 || res_koef[i].GetRe() != 0)
        {
            res_degree = i;
            Polinom res(res_degree, res_koef);
            /*
            cout << "res_degree=" << res_degree << endl;
            for (int j = 0; j <= res_degree; j++)
            cout << "res_koef[" << j << "]=" << res_koef[i] << endl;
            */
            return res;
        }
    }
    Polinom o;
    return o;
}
//Перегрузка -
Polinom Polinom::operator-(Polinom x)
{
    int res_degree = 0;
    if (degree > x.degree)
    {
        res_degree = degree;
        //увеличение степени(выделение места для коэфов)
        x.uvelich(degree);
    }
    else
    {
        res_degree = x.degree;
        //
        uvelich(x.degree);
    }
    Complex *res_koef = new Complex[res_degree + 1];
    for (int i = 0; i <= res_degree; i++)
        res_koef[i] = koef[i] - x.koef[i];
    for (int i = res_degree; i >= 0; i--) {
        if (res_koef[i].GetIm() != 0 || res_koef[i].GetRe() != 0)
        {
            res_degree = i;
            Polinom res(res_degree, res_koef);
            return res;
        }
    }
    Polinom o;
    return o;
}
//uvelich
//Не повышает степень самого полинома, но выделяет место для новых коэффициентов
void Polinom::uvelich(int new_degree)
{
    Complex *a = new Complex[new_degree + 1];
    for (int i = 0; i <= degree; i++)
        a[i] = koef[i];
    Complex o;
    for (int i = degree + 1; i <= new_degree; i++)
        a[i] = o;
    //Polinom res(new_degree, a);
    //res.SetDegree(degree);
    koef = new Complex[new_degree + 1];
    for (int i = 0; i <= new_degree; i++)
        koef[i] = a[i];
}
//Set
void Polinom::SetDegree(int a)
{
    degree = a;
}
//Get
Complex Polinom::GetKoef(int index)
{
    return koef[index];
}
//Перегрузка /
Polinom Polinom::operator/(Polinom P)
{
    Complex o;
    //0 Частные случаи
    if (P.degree > degree)
        return *this;
    if (P.degree == 0)
    {
        Complex *res = new Complex[degree + 1];
        for (int i = 0; i <= degree; i++)
            res[i] = koef[i] / P.koef[0];
        Polinom result(new_degree(res, degree), res);
        return result;
    }
    //1
 
 
 
    Complex *res = new Complex[degree + 1];
    //Создание массива делимого многочлена и его заполнение
    Complex *delim = new Complex[degree + 1];
    for (int i = 0; i <= degree; i++)
        delim[i] = koef[i];
    //Создание массива для делителя многочлена и его заполнение
    Complex *delit_old = new Complex[degree + 1];
    for (int i = 0; i <= P.degree; i++)
        delit_old[i] = P.koef[i];
    for (int i = P.degree + 1; i <= degree; i++)
        delit_old[i] = o;
 
    int schet_for_temp = 0;
    while (1) {
 
        //Повышение степени делителя
        int NEW = new_degree(delim, degree);//NEW-степень до которой надо повысить степень делителя
        int OLD = new_degree(delit_old, degree);//OLD-степень нашего делителя
        if (NEW + 1 >= OLD)
            delit_old = sdvig(delit_old, NEW, OLD);
        else
        {
            for (int i = 0; i <= OLD; i++)
                delit_old[i] = P.koef[i];
        }
        //2
        Complex temp = last_nonull(delim, new_degree(delim, degree - schet_for_temp)) / last_nonull(delit_old, new_degree(delit_old, degree - schet_for_temp));//последний ненулевой коэф делимого
        res[schet_for_temp] = temp;
        schet_for_temp++;
 
        //Новое делимое равно delim-delit*temp
        for (int i = 0; i <= new_degree(delim, degree); i++)
            delim[i] = delim[i] - delit_old[i] * temp;
        //переопределение delit
        delit_old = new Complex[degree + 1];
        for (int i = 0; i <= P.degree; i++)
            delit_old[i] = P.koef[i];
 
        if (new_degree(delim, degree) < P.degree)
            break;
    }
 
    reverse_(res, new_degree(res, degree) + 1);
    int NEW = new_degree(res, degree);
    Polinom result(NEW, res);
    return result;
}
void zapolnenie(Complex *&x, int deg) {
    string s;
    cin >> s;
    //if (x)
    //  delete[] x;
    x = new Complex[deg + 1];
    for (int i = 0; i <= deg; i++) {
 
        s.erase(0, 1);
        string temp1;
        while (s[0] != ',') {
            temp1 += s[0];
            s.erase(0, 1);
        }
        s.erase(0, 1);
        string temp2;
        while (s[0] != ')') {
            temp2 += s[0];
            s.erase(0, 1);
        }
        s.erase(0, 1);
        double tmp1, tmp2;
        tmp1 = atof(temp1.c_str());
        tmp2 = atof(temp2.c_str());
 
        Complex tmp(tmp1, tmp2);
        x[i] = tmp;
    }
}
Про функцию proovi, если возникли вопросы.
Кликните здесь для просмотра всего текста

Функция proovi нужна была, потому что в выводе периодически появлялись минус нули, которых быть не должно. Сам не знаю почему они появлялись, но они тоже появлялись именно при делении. Буду также благодарен за объяснение этого странного, на мой взгляд, феномена. Хоть это проблема была исправлена костылем, было бы интересно узнать ее причины. Но она все равно является дополнительной и не обязательной.
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
12.05.2017, 00:01
Ответы с готовыми решениями:

Деление многочленов от двух переменных
Есть многочлен от двух переменных, заданный следующей структурой: struct Monom { int...

Деление многочленов
Помогите,пожалуйста! Застрял,конкретно. Есть класс многочленов,представленный в виде двусвязного...

Деление многочленов
https://www.cyberforum.ru/cpp-beginners/thread328191.html Из этой темы был представлен алгоритм...

Деление многочленов 4 степени
Доброго времени суток! Необходимо было реализовать алгоритм, делящий два многочлена 4 степени с...

7
Форумчанин
Эксперт CЭксперт С++
8165 / 5013 / 1436
Регистрация: 29.11.2010
Сообщений: 13,455
12.05.2017, 00:29 2
Цитата Сообщение от ApocFel Посмотреть сообщение
Буду также благодарен за объяснение этого странного, на мой взгляд, феномена
Проблема в представлении чисел с плавающей точкой.
0
0 / 0 / 0
Регистрация: 25.12.2016
Сообщений: 39
12.05.2017, 00:35  [ТС] 3
Доп. квест сделан
0
21 / 21 / 10
Регистрация: 11.09.2015
Сообщений: 103
12.05.2017, 00:48 4
ApocFel, если ты случайно пропустил, может тебе поможет логика системы операций для действительных многочленов, что я опубликовал только что
Деление многочленов

Добавлено через 5 минут
Несколько странно изобретать свои представления и операции для комплексных чисел, имея под рукой STL с классом complex. Это как если бы на курсах трактористов и экскаваторщиков учили работать лопатой и киркой.

Добавлено через 6 минут
В любом умножение-деление чисел и многочленов - разные вещи.
0
0 / 0 / 0
Регистрация: 25.12.2016
Сообщений: 39
12.05.2017, 10:53  [ТС] 5
Хотелось бы именно исправление, а не новый алгоритм.
0
21 / 21 / 10
Регистрация: 11.09.2015
Сообщений: 103
12.05.2017, 13:49 6
ApocFel, разбирать чужой код труднее, чем писать свой. Я сделал лёгкую работу, а ты помоложе, справишься и с трудной
Кликните здесь для просмотра всего текста

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
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
#include <algorithm> // sort
#include <complex>
#include <fstream>
#include <iomanip>
#include <iostream>
#include <limits>
#include <list>
#include <sstream>
#include <string>
#include <utility> // pair
#include <vector>
#include <cstdio>
using namespace std;
 
using Item = pair<int,complex<double>>;
using Poly = list<Item>;
 
vector<complex<double>> AllCoef (const Poly* pol, int pwr) {
  vector<complex<double>> coef (pwr+1, 0.0);
  for (const Item& itm : *pol)
    coef[pwr-itm.first] = itm.second;
  return coef; // implicit move(), compiler optimization
} // AllCoef()
 
void NormPolynom (Poly* pol) {
  if (pol->size() == 0)
    return;
  // Привести подобные степени
  auto bgn = begin(*pol);
  while (bgn != end(*pol)) {
    auto p = *bgn;
    int pwr = p.first;
    complex<double> cft = p.second;
    auto nxt = bgn; ++nxt;
    while (nxt != end(*pol)) {
      auto q = *nxt;
      if (q.first != pwr)
        ++nxt;
      else {
        cft += q.second;
        nxt = pol->erase(nxt);
      }
    }
    if (cft == p.second)
      ++bgn;
    else {
      p.second = cft;
      auto del = bgn;
      bgn = pol->insert(bgn, p);
      pol->erase(del);;
    }
  }
  // Упорядочить по убыванию степеней
  pol->sort ([] (Item& p1, Item& p2)
    { return p1.first > p2.first; });
  // Удалить старшие нули
  bgn = begin(*pol);
  while (bgn != end(*pol) && real((*bgn).second) == 0 &&
      imag((*bgn).second) == 0) {
    pol->erase(bgn);
    bgn = begin(*pol);
  }
} // NormPolynom()
 
Poly* AddPolynom (Poly* pol1, Poly* pol2) {
  Poly* pol = new Poly;
  int pwr = max (pol1->front().first, pol2->front().first);
  vector<complex<double>> cf1 = AllCoef (pol1,  pwr);
  vector<complex<double>> cf2 = AllCoef (pol2,  pwr);
  for (int i = 0; i <= pwr; ++i) {
    Item itm = make_pair (pwr - i, cf1[i] + cf2[i]);
    if (real(itm.second) != 0 || imag(itm.second))
      pol->emplace_back(itm);
  }
  //NormPolynom (pol);
  return pol;
} // AddPolynom()
 
Poly* DivPolynom (Poly* pol1, Poly* pol2, Poly** rem) {
  const int pwr1 = pol1->front().first;
  const int pwr2 = pol2->front().first;
  if (pwr1 < pwr2)
    throw logic_error ("DivPolynom: pwr1 < pwr2");
  Poly* pol = new Poly;
  *rem = new Poly;
  vector<complex<double>> cf1 = AllCoef (pol1,  pwr1);
  vector<complex<double>> cf2 = AllCoef (pol2,  pwr2);
  for (int i = 0; i <= pwr1 - pwr2; ++i) {
    if (real(cf1[i]) == 0 && imag(cf1[i]) == 0) continue;
    // Вычисляем и размещаем коэффициент частного
    complex<double> c = cf1[i] / cf2[0];
    int k = pwr1 - i - pwr2;
    pol->emplace_back (make_pair (k, c));
    // Пересчитываем коэффициенты остатка
    cf1[i] = 0;
    for (unsigned j = 1; j < cf2.size(); ++j)
      cf1[i+j] -= cf2[j] * c;
  }
  // Размещаем коэффициенты остатка в полиноме
  double eps = numeric_limits<double>::epsilon();
  for (unsigned i = 1; i < cf1.size(); ++i) {
    if (abs (cf1[i]) > eps)
      (*rem)->emplace_back (make_pair (pwr1 - i, cf1[i]));
  }
  NormPolynom (pol);
  NormPolynom (*rem);
  return pol;
} // DivPolynom()
 
Item GetItem (string txt) { // только p r,i
  while (txt[0] == ' ' && txt.size() > 0)
    txt.erase (begin(txt));
  int pwr = stoi (txt);
  int blank = txt.find (" ");
  txt = txt.substr (blank+1);
  double real = stod(txt);
  int comma = txt.find (",");
  txt = txt.substr (comma+1);
  double imag = stod(txt);
  if (blank < 0 || comma < 0)
    throw runtime_error ("invalid Item");
  return make_pair (pwr, complex<double>(real, imag));
}
  
Poly* GetPolynom (istream& in) {
  Poly* pol = new Poly;
  string line, str1;
  getline (in, line);
  istringstream elm1 (line);
  while (getline (elm1, str1, ';')) {
    //cout << "\n<" << str1 << ">";
    Item itm = GetItem (str1);
    cout << "\nitm "
      << real(itm.second) << " "  << imag(itm.second);
    pol->emplace_back (make_pair (itm.first, itm.second));
  }
  return pol;
} // GetPolynom()
 
Poly* MultPolynom (Poly* pol1, Poly* pol2) {
  Poly* pol = new Poly;
  for (Item& p1 : *pol1) {
    for (Item& p2 : *pol2) {
      Item itm = make_pair (p1.first + p2.first
        ,p1.second * p2.second);
      if (real(itm.second) != 0 || imag(itm.second))
        pol->emplace_back (itm);
    }
  }
  NormPolynom (pol);
  return pol;
} // MultPolynom()
 
void PrintPolynom (string title, Poly* pol, ofstream& log) {
  cout << title << '\n';
  log << title << '\n';
  if (pol->size() == 0) {
    cout << 0;
    log << 0;
    return;
  }
  bool firstTime = true;
  for (Item& p : *pol) {
    if (! firstTime) {
      cout << "; ";
      log << "; ";
    }
    cout << p.first << " "
      << real(p.second) << "," << imag(p.second);
    log  << p.first << " "
      << real(p.second) << "," << imag(p.second);
    firstTime = false;
  }
} // PrintPolynom()
 
Poly* SubtPolynom (Poly* pol1, Poly* pol2) {
  Poly* pol = new Poly;
  int pwr = max (pol1->front().first, pol2->front().first);
  vector<complex<double>> cf1 = AllCoef (pol1,  pwr);
  vector<complex<double>> cf2 = AllCoef (pol2,  pwr);
  for (int i = 0; i <= pwr; ++i) {
    Item itm = make_pair (pwr - i, cf1[i] - cf2[i]);
    if (real(itm.second) != 0 || imag(itm.second))
      pol->emplace_back(itm);
  }
  //NormPolynom (pol);
  return pol;
} // SubtPolynom()
 
void main () { // При вводе с клавиатуры (cin)
  // конец полинома по клавише Enter
  string logName (__FILE__);
  logName = logName.replace (logName.find(".cpp"),4,".txt");
  istringstream sin1 ("0 1,2; 1 2,3; 2 1,2; 0 1,2");
  istringstream sin2 ("0 1,1"); // i1
  istringstream sin3 ("1 1,0; 0 1,0"); // x+1
  ofstream log (logName);
  Poly *pol1, *pol2, *pol3, *pol4, *pol5, *pol6, *rem;
  try {
    pol1 = GetPolynom (sin1);
    PrintPolynom ("pol1", pol1, log);
    NormPolynom (pol1);
    PrintPolynom ("\nNormalized pol1", pol1, log);
    pol2 = AddPolynom (pol1, pol1);
    PrintPolynom ("\npol2", pol2, log);
    pol3 = SubtPolynom (pol2, pol1);
    PrintPolynom ("\npol3", pol3, log);
    pol4 = GetPolynom (sin2);
    PrintPolynom ("\npol4", pol4, log);
    pol5 = MultPolynom (pol4, pol4);
    PrintPolynom ("\npol5", pol5, log);
    pol6 = DivPolynom (pol5, pol4, &rem);
    PrintPolynom ("\npol6", pol6, log);
    PrintPolynom ("\nrem", rem, log);
  } catch (exception err) {
    cout << string("\nBad input data: ") + err.what();
    exit(0);
  }
}


Добавлено через 14 минут
В GetPolynom() застряли строки 131,133,134 от отладочного вывода.

Добавлено через 1 час 18 минут
Две строки создания и размещения элемента в GetPolynom() можно заменить одной:
pol->emplace_back (GetItem (str1));
0
0 / 0 / 0
Регистрация: 25.12.2016
Сообщений: 39
12.05.2017, 16:40  [ТС] 7
Вы выполняете другую задачу.
Цитата Сообщение от ApocFel Посмотреть сообщение
Надо: найти ошибку в алгоритме и/или исправить код
Вы же предлагаете другой код.
0
4455 / 2074 / 263
Регистрация: 01.03.2013
Сообщений: 5,511
Записей в блоге: 22
12.05.2017, 16:42 8
Другой кот - это лучший вариант исправления существующего
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
12.05.2017, 16:42

Заказываю контрольные, курсовые, дипломные работы и диссертации здесь.

Деление многочленов(полиномов). доделать класс
Суть задания - сделать класс-полином со основными арифметическими операциями между многочленами....

Нужны советы как реализовать сложение, вычитание, умножение полиномов/многочленов и деление на число
Всем здравствуйте! надеюсь написать с вашей помощью программку для сложения, вычитания и умножения...

Сумма двух многочленов
#include&lt;iostream&gt; using namespace std; int main() { float* a, *b; float x, r, R, v, S = 0;...

Перемножение двух многочленов
CPolynom operator * (CPolynom &amp;A, CPolynom &amp;B) { int p = A.getPower() + B.getPower(); CPolynom...

Умножение двух многочленов
Здравствуйте! Помогите умножить два многочлена вида: (a0 + a1x + a2x2 + ... + anxn) * (b0 + b1x +...

Поместить в массив с[m] коэффициенты произведения двух многочленов.
В массивах a и b хранятся коэффициенты двух многочленов степеней k и l. Поместить в массив с...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
8
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.