BZOJ3625

此处不作题解。。只有细节

记住两个多项式相乘的时候必须 在高位留出空，

也就是代码里的memset(b+n,0,n<<2);

然后 在%x^n的操作结束后一定要把>=n的部分清零，保持答案在%x^n意义下。(当然，如果能记得，之后用的时候再清零也行）

还有就是从本题学到两个卡常技巧：

一、

① Time:24732ms
void ntt(int a[],int n,int d=1){
    for (int i=0;i<n;++i) if (i<rev[i]) swap(a[i],a[rev[i]]);
    for (int m=2;m<=n;m<<=1){
        int k=m>>1; LL wm=po(G,d*(mo-1)/m);

        for (int i=1;i<k;++i) w[i]=w[i-1]*wm%mo;    //瓶颈在这？
 
        for (int i=0;i<n;i+=m){
            for (int j=i;j<i+k;++j){
                int u=a[j],v=1LL*w[j-i]*a[j+k]%mo;
                a[j]=u+v<mo?u+v:u+v-mo;
                a[j+k]=u<v?mo-v+u:u-v;
            }
        }
    }
    if (d==-1){
        LL x=po(n,mo-2);
        for (int i=0;i<n;++i) a[i]=x*a[i]%mo;
    }
}

② Time:31844ms
void ntt(int a[],int n,int d=1){
    for (int i=0;i<n;++i) if (i<rev[i]) swap(a[i],a[rev[i]]);
    for (int m=2;m<=n;m<<=1){
        int k=m>>1,wm=po(G,d*(mo-1)/m);
        for (int i=0;i<n;i+=m){
            LL w=1;
            for (int j=i;j<i+k;++j){
                int u=a[j],v=w*a[j+k]%mo;
                a[j]=u+v<mo?u+v:u+v-mo;
                a[j+k]=u<v?mo-v+u:u-v;
                w=w*wm%mo;
            }
        }
    }
    if (d==-1){
        LL x=po(n,mo-2);
        for (int i=0;i<n;++i) a[i]=x*a[i]%mo;
    }
}

二、

① Time:21844ms
int w[N]; 
int wm=po(G,d*(mo-1)/m);
for (int i=1;i<k;++i) w[i]=1LL*w[i-1]*wm%mo;

② Time:24732ms    所以瓶颈应该是LL的读写速度(吧。。。)
LL w[N];
LL wm=po(G,d*(mo-1)/m);
for (int i=1;i<k;++i) w[i]=w[i-1]*wm%mo;

③ Time:24696ms    这个是为了证明并不是LL*LL和LL*int的瓶颈
LL w[N];
int wm=po(G,d*(mo-1)/m);
for (int i=1;i<k;++i) w[i]=w[i-1]*wm%mo;

贴代码

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const int mo=998244353;
const int G=3;
const int N=300000;
int n,m,rev[N],x,b[N],c[N],t1[N],t2[N],t3[N],w[N];
int po(int x,int y){
    if (y<0) y=y%(mo-1)+mo-1; int z=1;
    for(;y;y>>=1,x=1LL*x*x%mo) if (y&1) z=1LL*z*x%mo;
    return z;
}
void ntt(int a[],int n,int d=1){
    for (int i=0;i<n;++i) if (i<rev[i]) swap(a[i],a[rev[i]]);
    for (int m=2;m<=n;m<<=1){
        int k=m>>1,wm=po(G,d*(mo-1)/m);
        for (int i=1;i<k;++i) w[i]=1LL*w[i-1]*wm%mo;
        for (int i=0;i<n;i+=m){
            for (int j=i;j<i+k;++j){
                int u=a[j],v=1LL*w[j-i]*a[j+k]%mo;
                a[j]=u+v<mo?u+v:u+v-mo;
                a[j+k]=u<v?mo-v+u:u-v;
            }
        }
    }
    if (d==-1){
        LL x=po(n,mo-2);
        for (int i=0;i<n;++i) a[i]=x*a[i]%mo;
    }
}
void inv(int *a,int n,int *b){
    if (n==1){
        b[0]=po(a[0],mo-2);
        b[1]=0; return;
    }
    inv(a,n>>1,b); int m=n<<1;
    for (int i=0;i<n;++i) t3[i]=a[i];
    memset(b+n,0,n<<2); memset(t3+n,0,n<<2);
    for (int i=0;i<m;++i) rev[i]=(rev[i>>1]>>1)+(i&1)*(m>>1);
    ntt(t3,m); ntt(b,m);
    for (int i=0;i<m;++i) b[i]=(2LL+mo-1LL*t3[i]*b[i]%mo)*b[i]%mo;
    ntt(b,m,-1);
    memset(b+n,0,n<<2);
}
void gen(int *a,int n,int *b){
    if (n==1){
        b[0]=1; b[1]=0; return;
    }
    gen(a,n>>1,b); int m=n<<1;
    for (int i=0;i<n;++i) t1[i]=a[i];            
    memset(t1+n,0,n<<2);
    inv(b,n,t2);
    for (int i=0;i<m;++i) rev[i]=(rev[i>>1]>>1)+(i&1)*(m>>1);
    ntt(t1,m); ntt(t2,m);
    for (int i=0;i<m;++i) t1[i]=1LL*t1[i]*t2[i]%mo;
    ntt(t1,m,-1); LL x=mo+1>>1;
    for (int i=0;i<n;++i) b[i]=x*(b[i]+t1[i])%mo;
    memset(b+n,0,n<<2);
}
int main(){
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;++i){
        scanf("%d",&x);
        if (x<=m) c[x]=mo-4;
    }
    n=1; w[0]=1;
    while (n<=m) n<<=1;
    c[0]=1; gen(c,n,b);
    if (b[0]) (++b[0])%=mo;
    else{
        for (int i=0;i<n;++i) b[i]=mo-b[i];
        b[0]?--b[0]:b[0]=mo-1;
    }
    inv(b,n,c);
    for (int i=1;i<=m;++i) printf("%d
",(c[i]<<1)%mo);
    return 0;
}