NTT多项式板子

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int input=1000000;
char InPut[1000005];
inline char GetChar()
{
    if(input>=1000000)fread(InPut,1,1000000,stdin),input=0;
    return InPut[input++];
}
#define GetChar() getchar()
template<typename Typename>
inline void Read(Typename &s)
{
    register bool f=false;register char c=GetChar();
    s=0.0;
    while(c>'9'||c<'0'){if(c=='-')f=true;c=GetChar();}
    while(!(c>'9'||c<'0')){s=s*10+c-48;c=GetChar();}
    if(f)s=-s;
}
const long long Mod=998244353;
const int Maxn=3e5;
template<typename Typename>
inline void Swap(Typename &x,Typename &y)
{register Typename swap=x;x=y;y=swap;}
inline long long GetPow(long long base,long long index)
{
    long long power=1;base%=Mod;
    while(index!=0)
    {
        if(index&1)power=power*base%Mod;
        index>>=1;base=base*base%Mod;
    }
    return power;
}
int Rev[Maxn];long long Gn[31],Gi[31];
inline void RevInit(int limit,int len)
{
    for(register int i=1;i<limit;i++)
        Rev[i]=(Rev[i>>1]>>1)|((i&1)<<(len-1));
}
inline void NTTInit()
{
    for(register int i=1,g=1;g<=30;g++,i<<=1)
    {
        Gn[g]=GetPow(3,(Mod-1)/(i<<1));
        Gi[g]=GetPow(Gn[g],Mod-2);
    }
}
inline void NTT(long long *Arr,int limit,bool flag=true)
{
    for(register int i=0;i<limit;i++)
        if(i<Rev[i])Swap(Arr[i],Arr[Rev[i]]);
    for(register int i=1,g=1;i<limit;i<<=1,g++)
    {
        register long long T=flag?Gn[g]:Gi[g];
        for(register int j=0;j<limit;j+=(i<<1))
        {
            register long long W=1;
            for(register int k=0;k<i;k++,W=(W*T)%Mod)
            {
                register long long x=Arr[j+k],y=W*Arr[j+k+i]%Mod;
                Arr[j+k]=(x+y)%Mod;Arr[j+k+i]=(x-y+Mod)%Mod;
            }
        }
    }
    if(flag)return;
    long long I=GetPow(limit,Mod-2);
    for(register int i=0;i<limit;i++)Arr[i]=Arr[i]*I%Mod;
}
long long Copy[Maxn];
inline void GetInv(long long *Arr,long long *Inv,int n)
{
    int now=1,limit=2,len=1;
    Inv[0]=GetPow(Arr[0],Mod-2);
    while(now<=n*2)
    {
        RevInit(limit,len);
        for(register int i=0;i<now;i++)Copy[i]=Arr[i];
        for(register int i=now;i<limit;i++)Copy[i]=0;
        NTT(Copy,limit);NTT(Inv,limit);
        for(register int i=0;i<limit;i++)
            Inv[i]=((2ll-Inv[i]*Copy[i]%Mod)+Mod)%Mod*Inv[i]%Mod;
        NTT(Inv,limit,false);
        for(register int i=now;i<limit;i++)Inv[i]=0;
        now<<=1;limit<<=1;len++;
    }
}
inline void GetDev(long long *Arr,long long *Inv,int n)
{
    for(register int i=1;i<n;i++)Inv[i-1]=i*Arr[i]%Mod;
    Inv[n-1]=0;
}
inline void GetInvDev(long long *Arr,long long *Inv,int n)
{
    for(register int i=1;i<n;i++)Inv[i]=Arr[i-1]*GetPow(i,Mod-2)%Mod;
    Inv[0]=0;
}
long long TmpAi[Maxn],TmpBi[Maxn];
inline void GetLn(long long *Arr,long long *Inv,int n)
{
    for(register int i=0;i<=(n<<1);i++)TmpAi[i]=0;
    for(register int i=0;i<=(n<<1);i++)TmpBi[i]=0;
    GetDev(Arr,TmpAi,n);
    GetInv(Arr,TmpBi,n);
    int limit=1,len=0;
    while(limit<(n<<1))limit<<=1,len++;
    RevInit(limit,len);
    NTT(TmpAi,limit);NTT(TmpBi,limit);
    for(register int i=0;i<limit;i++)TmpAi[i]=TmpAi[i]*TmpBi[i]%Mod;
    NTT(TmpAi,limit,false);
    GetInvDev(TmpAi,Inv,n);
}
long long TmpCi[Maxn],TmpDi[Maxn];
inline void GetExp(long long *Arr,long long *Inv,int n)
{
    int len=1,limit=2,now=1;
    Inv[0]=1;
    while(now<=n*2)
    {
        RevInit(limit,len);
        for(register int i=0;i<limit;i++)TmpCi[i]=0;
        for(register int i=0;i<now;i++)TmpDi[i]=Arr[i];
        for(register int i=now;i<limit;i++)TmpDi[i]=0;
        GetLn(Inv,TmpCi,now);
        NTT(Inv,limit);NTT(TmpCi,limit);NTT(TmpDi,limit);
        for(register int i=0;i<limit;i++)
            Inv[i]=(1ll-TmpCi[i]+TmpDi[i]+Mod)*Inv[i]%Mod;
        NTT(Inv,limit,false);
        for(register int i=now;i<limit;i++)Inv[i]=0;
        now<<=1;limit<<=1;len++;
    }
}
long long F[Maxn],G[Maxn],RG[Maxn],RF[Maxn],_G[Maxn],Q[Maxn],R[Maxn];
int n,m,limit,len;
int main()
{
    NTTInit();
    scanf("%d%d",&n,&m);
    for(register int i=0;i<=n;i++)scanf("%lld",&F[i]),RF[n-i]=F[i];
    for(register int i=0;i<=m;i++)scanf("%lld",&G[i]),RG[m-i]=G[i];
    for(register int i=n-m+1;i<=m;i++)RG[i]=0;
    for(register int i=n-m+1;i<=n;i++)RF[i]=RG[i]=0;
    GetInv(RG,_G,n-m);

    limit=1,len=0;
    while(limit<=n+n-m)limit<<=1,len++;
    RevInit(limit,len);
    NTT(RF,limit);NTT(_G,limit);
    for(register int i=0;i<limit;i++)RF[i]=RF[i]*_G[i]%Mod;
    NTT(RF,limit,false);
    for(register int i=0;i<=n-m;i++)Q[i]=RF[n-m-i],printf("%lld ",Q[i]);
    for(register int i=2*n-m+1;i<limit;i++)RF[i]=0;
    printf("\n");

    limit=1,len=0;
    while(limit<=n)limit<<=1,len++;
    RevInit(limit,len);
    NTT(G,limit);NTT(Q,limit);
    for(register int i=0;i<limit;i++)G[i]=G[i]*Q[i]%Mod;
    NTT(G,limit,false);
    for(register int i=0;i<m;i++)R[i]=(F[i]-G[i]+Mod)%Mod,printf("%lld ",R[i]);

    return 0;
}