阅读背景:

多项式的一些操作

来源:互联网 

大多数东西只有代码啦
主要是因为数学公式太烦了

快速傅里叶变换 FFT

void FFT(Complex *a,int tp){
    for(int i=0;i<Lim;i++)
    if(i<rev[i])swap(a[i],a[rev[i]]);
    for(int pos=1;pos<Lim;pos<<=1){
        Complex w; w.a=cos(pi/pos),w.b=tp*sin(pi/pos);
        for(int R=pos<<1,j=0;j<Lim;j+=R){
            Complex p; p.a=1,p.b=0;
            for(int k=j;k<pos+j;k++,p=p*w){
                Complex x=a[k],y=p*a[k+pos];
                a[k]=x+y,a[k+pos]=x-y;
            }
        }
    }
}void FF



你的当前访问异常,请进行认证后继续阅读剩余内容。

分享到: