P3803 【模板】多项式乘法（NTT）

传送门

NTT好像是比FFT快了不少

然而感觉不是很看得懂……主要是点值转化为系数表示那里……

upd：大概已经搞明白是个什么玩意儿了……吧……

 //minamoto
 #include<bits/stdc++.h>
 #define R register
 #define fp(i,a,b) for(R int i=a,I=b+1;i<I;++i)
 #define fd(i,a,b) for(R int i=a,I=b-1;i>I;--i)
 #define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
 using namespace std;
 char buf[<<],*p1=buf,*p2=buf;
 inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,,<<,stdin),p1==p2)?EOF:*p1++;}
 int read(){
     R int res,f=;R char ch;
     while((ch=getc())>''||ch<'')(ch=='-')&&(f=-);
     for(res=ch-'';(ch=getc())>=''&&ch<='';res=res*+ch-'');
     return res*f;
 }
 char sr[<<],z[];int C=-,Z=;
 inline void Ot(){fwrite(sr,,C+,stdout),C=-;}
 void print(R int x){
     if(C><<)Ot();if(x<)sr[++C]='-',x=-x;
     while(z[++Z]=x%+,x/=);
     while(sr[++C]=z[Z],--Z);sr[++C]=' ';
 }
 const int N=3e6+,P=,Gi=;
 inline int add(R int x,R int y){return x+y>=P?x+y-P:x+y;}
 inline int dec(R int x,R int y){return x-y<?x-y+P:x-y;}
 inline int mul(R int x,R int y){return 1ll*x*y-1ll*x*y/P*P;}
 int ksm(int x,int y){
     R int res=;
     for(;y;y>>=,x=mul(x,x))if(y&)res=mul(res,x);
     return res;
 }
 int A[N],B[N],O[N],r[N],lim=,n,m,l;
 void NTT(int *A,int ty){
     fp(i,,lim-)if(i<r[i])swap(A[i],A[r[i]]);
     for(R int mid=;mid<lim;mid<<=){
         int D=(mid<<),Wn=ksm(ty==?:Gi,(P-)/D);O[]=;
         fp(i,,mid-)O[i]=mul(O[i-],Wn);
         for(R int j=;j<lim;j+=D){
             for(R int k=;k<mid;++k){
                 int x=A[j+k],y=mul(O[k],A[j+k+mid]);
                 A[j+k]=add(x,y),A[j+k+mid]=dec(x,y);
             }
         }
     }
     if(ty==-)for(R int i=,inv=ksm(lim,P-);i<lim;++i)A[i]=mul(A[i],inv);
 }
 int main(){
 //    freopen("testdata.in","r",stdin);
     n=read(),m=read();while(lim<=n+m)lim<<=,++l;
     fp(i,,lim-)r[i]=(r[i>>]>>)|((i&)<<(l-));
     fp(i,,n)A[i]=read();fp(i,,m)B[i]=read();
     NTT(A,),NTT(B,);
     fp(i,,lim-)A[i]=mul(A[i],B[i]);
     NTT(A,-);
     fp(i,,n+m)print(A[i]);return Ot(),;
 }