【NTT】【多项式】多项式快速插值（log^3）模板

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define SF scanf
#define PF printf
#define MAXN 3000010
#define MOD 998244353
using namespace std;
const int G=3;
int n,m,b[MAXN],WN[MAXN];
int buf[MAXN*4];
int *bgn[MAXN],*ncnt=buf;
int fsp(int x,int y){
   
    int res=1;
    while(y){
   
        if(y&1)
            res=1ll*res*x%MOD;
        x=1ll*x*x%MOD;
        y>>=1;
    }
    return res;
}
void ntt(int A[],int N,int flag){
   
    for(int i=1,j=0;i<N;i++){
   
        for(int d=N;j^=d>>=1,~j&d;);
        if(i<j)
            swap(A[i],A[j]);	
    }
    for(int i=1,id=0;i<N;i<<=1,id++){
   
//        int wn=fsp(G,(MOD-1)/(i<<1));
		int wn=WN[id];
        if(flag) wn=fsp(wn,MOD-2);
        for(int j=0;j<N;j+=i<<1){
   
            int w=1;
            for(int k=0;k<i;k++,w=1ll*w*wn%MOD){
   
                int x=A[j+k],y=1ll*w*A[i+j+k]%MOD;
                A[j+k]=(x+y)%MOD;
                A[i+j+k]=(x-y+MOD)%MOD;
            }
        }
    }
    if(flag) for(int i=0,invN=fsp(N,MOD-2);i<N;i++) A[i]=1ll*A[i]*invN%MOD;
}
void mul(int A[],int N,int B[],int M,int res[]){
   
    static int A1[MAXN],B1[MAXN],res1[