多项式相关

本文介绍了多项式乘法的方法,包括从朴素暴力算法到利用点值表达的优化,最终引入快速傅里叶变换(FFT)以实现O(NlogN)的时间复杂度。内容涵盖单位复根的性质、消去引理、折半引理和求和引理,阐述了FFT如何通过指数奇偶性将问题规模减半,从而提高计算效率。
摘要由CSDN通过智能技术生成

引用

从多项式乘法到快速傅里叶变换
【快速傅里叶变换】【FFT】

准备知识

多项式表达

对于一个多项式我们可以写成系数表达和点值表达
对于多项式 f(x)=i=ni=0aixi
n+1,n

  1. 系数表达
    写成 n+1 维向量形式 a⃗ =(a0,a1,,an)
  2. 点值表达
    我们要选取任意 n+1 个点值 x0,,xn 求出它的 f(xi) ,得到 { (xi,f(xi)):0<=i<=n,iZ}
    n+1 个点值可以任意选取, FFT 的高效性也是因为这 n+1 个点的特殊选取
  3. 卷积运算
    • 点值表达有什么优越性?
    • 对于两个多项式,我们分别得到了他们的点值表达,我们就能将两个点值表达式相乘(不足的位也要算出来) O(N) 时间得到两个多项式乘积的多项式的点值表达式
      这个乘积运算就是卷积
      我们举个例子来看
      A(x)=x2+2x1     { (0,1),(1,2),(2,7),(3,14),(4,23)}
      B(x)=x2x+2       { (0,2),(1,2),(2,4),(3,8),(4,14)}
      C(x)=x4+x3x2+5x2    { (0,2),(1,4),(2,28),(3,112),(4,322)}
  4. 复数
    • n 次单位复根:满足
  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
。一元多项式的标准表达式可以用一个数组来表示,数组下标表示项的次数,数组元素表示该次数项的系数。例如,一个三次多项式f(x)=2x^3-3x^2+4x+1可以表示为: ```c float poly[4] = {1, 4, -3, 2}; ``` 其中,数组下标0对应x^0,即常数项系数为1,数组下标1对应x^1,即一次项系数为4,依此类推。可以通过循环遍历数组来输出多项式的标准表达式。 多项式相加的实现可以使用类似的循环遍历数组的方法,将同一次数的项的系数相加即可。例如,将两个三次多项式f(x)=2x^3-3x^2+4x+1和g(x)=x^3+2x^2+3x+5相加得到: ```c float poly1[4] = {1, 4, -3, 2}; float poly2[4] = {5, 3, 2, 1}; float result[4]; // 存放结果多项式的数组 int i; for(i=0; i<4; i++) { result[i] = poly1[i] + poly2[i]; } ``` 多项式相乘的实现比较复杂,可以使用嵌套循环遍历两个多项式的所有项,将相应次数项的系数相乘再累加得到结果多项式的相应项系数。例如,将f(x)和g(x)相乘得到: ```c float poly1[4] = {1, 4, -3, 2}; float poly2[4] = {5, 3, 2, 1}; float result[7]; // 存放结果多项式的数组,最高次为6 int i, j; for(i=0; i<7; i++) { result[i] = 0; for(j=0; j<=i; j++) { if(j<4 && (i-j)<4) { // 保证下标不越界 result[i] += poly1[j] * poly2[i-j]; } } } ``` 上述代码中,外层循环遍历结果多项式的所有项,内层循环遍历两个多项式的所有可能的项的乘积,并累加到结果多项式的相应项系数中。注意,结果多项式的最高次数是两个多项式最高次数之和减1。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值