luogu1919 A*BProblem升级版 (FFT)

把一个n位数看做n-1次的多项式，每一项的系数是反过来的每一位
最后每一项系数进进位搞一搞就行了
（数组一定要开到2的次数..要不然极端数据会RE）

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=132000;
const double Pi=acos(-1);

struct Cpx{
    double x,y;
    Cpx(double xx=0,double yy=0){x=xx;y=yy;}
}X[maxn],Y[maxn];
Cpx operator *(Cpx a,Cpx b){return Cpx(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}
Cpx operator +(Cpx a,Cpx b){return Cpx(a.x+b.x,a.y+b.y);}
Cpx operator -(Cpx a,Cpx b){return Cpx(a.x-b.x,a.y-b.y);}
int N,M,rev[maxn],ans[maxn];

void rd(Cpx *A){
    char c=getchar();
    while(c<'0'||c>'9') c=getchar();
    int i=N-1;
    while(c>='0'&&c<='9') A[i--].x=(int)(c-'0'),c=getchar();
}

void fft(Cpx *A,int opt){
    for(int i=0;i<N;i++) if(i<rev[i]) swap(A[i],A[rev[i]]);
    for(int l=1;l<N;l<<=1){
        Cpx wn=Cpx(cos(Pi/l),opt*sin(Pi/l));int step=l<<1;
        for(int i=0;i<N;i+=step){
            Cpx w=Cpx(1,0);
            for(int k=0;k<l;k++,w=w*wn){
                Cpx a=A[i+k],b=A[i+k+l]*w;
                A[i+k]=a+b;A[i+k+l]=a-b;
            }
        }
    }
}

int main(){
    int i,j,k;
    scanf("%d",&N);M=N*2-1;
    rd(X);rd(Y);
    for(i=1,j=0;i<M;i<<=1,j++);N=i;
    for(i=0;i<N;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<(j-1));
    fft(X,1);fft(Y,1);
    for(i=0;i<N;i++) X[i]=X[i]*Y[i];fft(X,-1);
    for(i=0;i<M;i++)
        j=(int)(X[i].x/N+0.5),ans[i+1]=(ans[i]+j)/10,ans[i]=(ans[i]+j)%10;
    for(i=M;i>=0&&!ans[i];i--);
    for(;i>=0;i--) printf("%d",ans[i]);
}