大多数东西只有代码啦
主要是因为数学公式太烦了
快速傅里叶变换 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