NTT模板（无讲解）

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