luogu1919 A*BProblem升级版 (FFT)

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

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

 struct Cpx{
     double x,y;
     Cpx(double xx=,double yy=){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<''||c>'') c=getchar();
     int i=N-;
     while(c>=''&&c<='') A[i--].x=(int)(c-''),c=getchar();
 }

 void fft(Cpx *A,int opt){
     for(int i=;i<N;i++) if(i<rev[i]) swap(A[i],A[rev[i]]);
     for(int l=;l<N;l<<=){
         Cpx wn=Cpx(cos(Pi/l),opt*sin(Pi/l));int step=l<<;
         for(int i=;i<N;i+=step){
             Cpx w=Cpx(,);
             for(int k=;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*-;
     rd(X);rd(Y);
     for(i=,j=;i<M;i<<=,j++);N=i;
     for(i=;i<N;i++) rev[i]=(rev[i>>]>>)|((i&)<<(j-));
     fft(X,);fft(Y,);
     for(i=;i<N;i++) X[i]=X[i]*Y[i];fft(X,-);
     for(i=;i<M;i++)
         j=(int)(X[i].x/N+0.5),ans[i+]=(ans[i]+j)/,ans[i]=(ans[i]+j)%;
     for(i=M;i>=&&!ans[i];i--);
     for(;i>=;i--) printf("%d",ans[i]);
 }