FFT中的一个常见小问题
这里不细说FFT的内容,详细内容看这些就足以了解大概了
主要是对学习过程中一个容易困扰的小问题进行解释,以便于理解
- 用FFT将多项式的系数转换为点值时,原系数数组a最后存的是不同的点值,而不是只有第一个是点值
这一点最开始困扰了我很久
设 $A(x)=a_0+a_1x+a_2x^2+...+a_{n?1}x^{n?1}$
则可将其移项 $A(x)=(a_0+a_2x^2+...+a_{n?2}x^{n?2})+(a_1x+a_{3}x^3+...+a_{n?1}x^{n?1})$
a的下标为偶数的放在一起 $A_1(x)=a_0+a_2x+...+a_{n?2}xn^{2?1}$
a的下标为奇数的放在一起 $A_2(x)=a_1+a_3x+...+a_{n?1}xn^{2?1}$
则$A(x)=A_1(x^2)+xA_2(x^2)$
注意此处为$x^2$所以有
$A(-x)=A_1(x^2)-xA_2(x^2)$
由于单位根的特殊性质,有
性质一 $ω_n^{k+\frac{n}{2}}=-ω_n^k$
性质二 $ω_n^k=ω_{2n}^{2k}$
所以才有了代码中的那两行
cpp for (int i=0;i<=mid-1;++i){ buf[i]=a[i]+w*a[i+mid]; buf[i+mid]=a[i]-w*a[i+mid]; w=w*wn; }
也就是说,我们可以由一个答案进而算出另外一个答案,这里可以理解为递推
所以当我们的递归递到最下面一层后往上走时每次都是将目前答案个数扩大两倍,而且这些答案是由不同的x算出来的,而且由于性质一,我们在计算过程中所用到的不同的$ω^{x*k}$是没有问题的
原文地址:https://www.cnblogs.com/Morning-Glory/p/10283697.html
时间: 2024-11-05 21:00:16