【基础算法】多项式三大运算 & C++实现

简介: 多项式三大运算 & C++实现

●多项式计算


一维多项式就是包含一个变量的多项式,一个一维多项式示例如下:

88d6ec418e0d6cde3e516b46e65af47c_60dfd07e2e144f05aa1dafa5f9e8f387.png

一维多项式求值就是对于上述多项式,计算在指定的x处的函数值。一个通用的计算多项式值的算法可以采用递推的方式,可以将上述多项式变为如下的等价形式(提公因式变形):

66947d3b29dcb8be986afd2fec5fb10b_6a5e1d37c7cb489b9c355a3cff8ebf4a.png

下面,我们按照这个算法来计算我们程序中的多项式。宏定义中已规定多项式系数的个数为7,并且在程序中我们输入系数值,则多项式如下所示:

7ce28132f7739975f314ad355a70a456_3263318075674be394f7614e2d0e1ede.png

输入的6个未知数x的值分别为,-2、-0.5、1、2、3.7、4 ,带入上式多项式则可知我们要求以下几个多项式的值:

6799e5ea335bcfd6061c787eb1a1628a_c510b3be89744b35be06a361b7e16686.png

最后在程序中用计算多项式值的算法对其进行计算,并求得结果。

#include<iostream>
using namespace std;
#define sizea 7  //多项式系数个数
#define sizex 6  //未知数x的个数
class calculator {
public:
  double polynomials_1(double x)
  {
  result = a[sizea - 1];
  for(int i=sizea-2;i>=0;i--)
  { 
    result = result * x + a[i];   
    //P(x)=A(n-1)*X^n-1+A(n-2)*X^n-2+...+A(1)*x+A(0)
    //P(x)=(...((A(n-1)X+A(n-2))X+A(n-3))X+...+A1)x+A(0)
  }
  return result;
  }
  double a[sizea];
  double x[sizex];
  double result;
};
void polynomials(calculator &c)
{
  for (int i = 0; i < sizex; i++)
  {
  c.polynomials_1(c.x[i]);
  cout << "x=" << c.x[i] << " " << "p(x)=" << c.result << endl;
  }
}
void text()
{
  calculator c;
  cout << "请输入系数:" << endl;
  for(int i=0;i<sizea;i++)
  { 
  cin >> c.a[i];
  }
  cout << "请输入x的值:" << endl;
  for(int j=0;j<sizex;j++)
  { 
  cin >> c.x[j];
  }
  polynomials(c);
}
int main()
{
  text();
}

cbf3fbcea8b2a6545012d5e974e41c83_75e7272ad9884c25b19e7b7c281690b1.png


●多项式乘法


多项式乘法就是将两个多项式进行相乘,最后得到一个新的多项式。如下所示:

d26cf2c7b22ec8accb48c5a149c4874d_d6801236c1094962ad998c4eea2c92fa.png

则这两个多项式相乘的结果如下:

437c387ea84b7f34ac63af0baa2a76c6_49d75e2be44c4a59bb19405c0f8a4bfe.png

在下面的程序中,宏定义已确定多项式①的系数个数为6个,多项式②的系数个数为4个,则两个多项式如下所示:

2365cf1c0f0439df252eff544b0ad8d7_e8420e8699854a19b888fd39d4e9aca1.png

最后在程序中运用两个多项式相乘的算法公式对其进行计算,并且求得结果。

#include<iostream>
using namespace std;
#define sizea 6  //多项式①的系数个数
#define sizeb 4  //多项式②的系数个数
#define sizec 9
class calculator {
public:
  void polynomialsmull()
  {
  for(int i=0;i<sizea;i++)
  { 
    for(int j=0;j<sizeb;j++)
    { 
    c[i + j] += a[i] * b[j];
    /*  
                    [0]=[0]*[0]/[1]=[0]*[1]/[2]=[0]*[2]/[3]=[0]*[3]
           [1]=[1]*[0]/[2]=[1]*[1]/[3]=[1]*[2]/[4]=[1]*[3]
           [2]=[2]*[0]/[3]=[2]*[1]/[4]=[2]*[2]/[5]=[2]*[3]
           ...
       */
    }
  }
  }
  void showresult()
  {
  cout << "两多项式相乘各项的系数如下:" << endl;
  for (int k = 0; k < sizec; k++)
  {
    cout << c[k] << " ";
  }
  }
  double a[sizea];
  double b[sizeb];
  double c[sizec] = {0};  //先将该数组每项置为0
};
void text()
{
  calculator c;
  cout << "请输入多项式①的系数a:" << endl;
  for(int i=0;i<sizea;i++)
  { 
  cin >> c.a[i];
  }
  cout << "请输入多项式②的系数b:" << endl;
  for(int j=0;j<sizeb;j++)
  { 
  cin >> c.b[j];
  }
  c.polynomialsmull();
  c.showresult();
}
int main()
{
  text();
}

4da78b529699c328311d22c9c13f3474_2a409cac9be9418eac6e6e00a279dfae.png


●多项式除法


多项式除法就是将两个多项式进行相除,最后得到一个商多项式和余多项式,如下所示得两个多项式:

e68df02e18ba1b1230b404f39d220e06_5cf43463bb3d4c4e905fe266de7fc44e.png

通过推算我们可以知道商多项式为 r(x)的最高次数为k=m-n,余多项式的最高次数为n-2.在数学上通过综合除法来计算商多项式 r(x)和余多项式为 l(x)中的各个系数,结合程序可进行进一步的细致理解。商多项式 l(x)的各个系数由下列递推算法来求得:

1c07cf5f47b09fd6eb4c21f56ba3a3bb_1a29397ecbb442c18c8b4b39eb545793.png

              在下面的程序中,宏定义已确定多项式①的系数个数为5个,多项式②的系数个数为3个,则两个多项式如下所示:

4e2153ea6e29f63843af248bff5d3728_8cf912898095469fac87f7f0bda9331b.png

 最后在程序中运用两个多项式相除的算法公式及其数学内容对其进行计算,并且求得结果。

#include<iostream>
using namespace std;
#define sizea 5   //多项式①的系数项个数
#define sizeb 3   //多项式②的系数项个数
#define sizer 3
#define sizel 2
class calculator {
public:
  void polynomials_div()
  {
  //计算商多项式系数和余多项式系数
  int m, n;
  n = sizea - 1;
  for(int i=sizer;i>0;i--)
  { 
    r[i - 1] = a[n] / b[sizeb - 1];    //商
    m = n;
    for (int j = 1; j <= sizeb - 1; j++)
    {
    a[m - 1] -= r[i - 1] * b[sizeb - j - 1];  //余
    m -= 1;
    }
    n -= 1;
  }
  for (int k = 0; k < sizel; k++)
  {
    l[k] = a[k];
  }
  }
  void showresult()
  {
  cout << "请输出商多项式:" << endl;
  for (int i = 0; i < sizer; i++)
  {
    cout << r[i]<<" ";
  }
  cout << endl;
  cout << "请输出余多项式:" << endl;
  for (int i = 0; i < sizel; i++)
  {
    cout << l[i] <<" ";
  }
  }
  double a[sizea];
  double b[sizeb];
  double r[sizer] = { 0 };
  double l[sizel] = { 0 };
};
void text()
{
  calculator c;
  cout << "请输入多项式①的系数a:" << endl;
  for(int i=0;i<sizea;i++)
  { 
  cin >> c.a[i];
  }
  cout << "请输入多项式②的系数b:" << endl;
  for(int j=0;j<sizeb;j++)
  { 
  cin >> c.b[j];
  }
  c.polynomials_div();
  c.showresult();
}
int main()
{
  text();
}

01ccb4a9636ff46c3bec254afb6ff83a_7847b48f0471492882e1ceae5a5bd065.png

目录
相关文章
|
2月前
|
算法 测试技术 C++
【动态规划算法】蓝桥杯填充问题(C/C++)
【动态规划算法】蓝桥杯填充问题(C/C++)
|
2月前
|
存储 算法 C++
高精度算法(加、减、乘、除,使用c++实现)
高精度算法(加、减、乘、除,使用c++实现)
570 0
高精度算法(加、减、乘、除,使用c++实现)
|
2月前
|
存储 算法 决策智能
【算法】博弈论(C/C++)
【算法】博弈论(C/C++)
|
2月前
|
存储 算法 C++
【算法】哈希映射(C/C++)
【算法】哈希映射(C/C++)
|
2月前
|
机器学习/深度学习 人工智能 算法
【算法】最长公共子序列(C/C++)
【算法】最长公共子序列(C/C++)
|
2月前
|
人工智能 算法 BI
一篇带你速通差分算法(C/C++)
一篇带你速通差分算法(C/C++)
|
2月前
|
人工智能 算法 C++
一篇带你速通前缀和算法(C/C++)
一篇带你速通前缀和算法(C/C++)
|
2月前
|
存储 算法 C++
弗洛伊德(Floyd)算法(C/C++)
弗洛伊德(Floyd)算法(C/C++)
|
2月前
|
存储 算法 程序员
迪杰斯特拉(Dijkstra)算法(C/C++)
迪杰斯特拉(Dijkstra)算法(C/C++)
|
2月前
|
算法 C++
【算法解题思想】动态规划+深度优先搜索(C/C++)
【算法解题思想】动态规划+深度优先搜索(C/C++)