NTT模板（无讲解）

 #include<bits/stdc++.h>//只是在虚数部分改了一下
 using namespace std;
 typedef long long int ll;
 const ll maxn=1E6+;
 const ll mod=;
 const ll G=;
 const ll Gi=;
 ll n,m,limit,r[maxn*],len,f[maxn],g[maxn];
 ll qpow(ll x,ll y)
 {
     ll ans=,base=x;
     while(y)
     {
         if(y&)ans=ans*base%mod;
         base=base*base%mod;
         y>>=;
     }
     return ans;
 }
 int re(int x)
 {
     int sum=;
     for(int i=;i<len;++i)sum=sum*+((x&(<<i))>);
     return sum;
 }
 void FFT(ll*A,int g)
 {
     for(int i=;i<limit;++i)
         if(i<r[i])swap(A[i],A[r[i]]);
     for(int i=;i<=limit;i*=)
     {
         ll w;
         if(g==)w=qpow(G,(mod-)/i);
         else w=qpow(Gi,(mod-)/i);
         for(int j=;j<limit/i;++j)
         {
             ll d=;
             for(int k=;k<i/;++k)
             {
                 ll a=A[i*j+k],b=d*A[i*j+i/+k]%mod;
                 A[i*j+k]=(a+b)%mod;
                 A[i*j+i/+k]=(a-b+mod)%mod;
                 d=d*w%mod;
             }
         }
     }
 }
 int main()
 {
     ios::sync_with_stdio(false);
     cin>>n>>m;
     for(int i=;i<=n;++i)cin>>f[i];
     for(int i=;i<=m;++i)cin>>g[i];
     limit=;
     len=;
     while(limit<=n+m+)
     {
         limit*=;
         ++len;
     }
     for(int i=;i<limit;++i)r[i]=re(i);
     FFT(f,);
     FFT(g,);
     for(int i=;i<=limit;++i)f[i]=f[i]*g[i];
     FFT(f,-);
     ll g=qpow(limit,mod-);
     for(int i=;i<=n+m;++i)cout<<(f[i]*g%mod+mod)%mod<<' ';
     cout<<endl;
     return ;
 }