大数乘法傅里叶c语言,FFT（快速傅里叶求大数乘法hdoj1402）

好久看过这么纠结的算法了，导论上的资料很神很好，作者果然是数学功底扎实啊，每一个函数都可以用数学语言描述出来,好神奇的卷积公式，看了好几遍才看懂，囧啊，果然是菜到没人鸟的鸟。

如果导论看懂了的话应该没，写程序应该没太大问题，主要是一些细节方面的问题了。如傅里叶变换过程中的运算用算数存取。

还好C++中有复数类complex,把它当成普通的数运算就行了。而傅里叶变换的原理理论上是没有误差的，计算机里用double型存储可精确到10^-16,也就是只要四舍五入就可以得到正确结果了。

可能是最近神经有点紧张，思维也变得迟钝，理解东西真是费解，现在终于可以松口气了！

#include #include #include #include #define N 50010// using namespace std; const double pi=acos(-1); const complex I(0,1); const double eps=1e-6; char aa[N],bb[N]; int ans[4*N];//char ans[4*N];!!!!!!! complex a[4*N],b[4*N];//an-1,an-2,……,a1,a0 int n; void initial() { int lenA=strlen(aa),lenB=strlen(bb); n=max(lenA,lenB); double t=log2(n); int tt=int(t); if(t-tt>0)tt++; n=1< * a) { for(int i=1; i>=1); if(j>i)swap(a[i],a[j]); } } void iterative_FFT(complex *a,int sig) { bitReverse(a); for(int m=2; m<=n; m<<=1)//是把a[]按m的长度分组， { int mh =m>>1; for(int i=0; i wi=exp(i*sig*pi/mh*I); for(int j=i; j u=a[j],t=wi*a[k];; a[j]=u+t; a[k]=u-t; } } } if(sig==-1)for(int i=0; i=0)cout<