BZOJ 3625 多项式求逆+多项式开根

思路：

RT

//By SiriusRen
#include <bits/stdc++.h>
using namespace std;
const int N=<<,mod=;
int A[N],C[N],invC[N],c[N],d[N],R[N],tmp[N],xx,len,sqrA[N],F[N];
typedef long long ll;
int power(ll x,int y){
    ll res=;
    while(y){
        if(y&)res=res*x%mod;
        x=x*x%mod,y>>=;
    }return res;
}
void NTT(int *a,int n,int f){
    int m=,L=;
    for(;m<n;m<<=)L++;
    for(int i=;i<n;i++)R[i]=(R[i>>]>>)|((i&)<<(L-));
    for(int i=;i<n;i++)if(i<R[i])swap(a[i],a[R[i]]);
    for(int l=;l<n;l<<=){
        int wn=power(,((mod-)/(l<<)*f+mod-)%(mod-));
        for(int j=;j<n;j+=(l<<)){
            int w=;
            for(int k=;k<l;k++,w=1ll*w*wn%mod){
                int x=a[j+k],y=1ll*a[j+k+l]*w%mod;
                a[j+k]=(x+y)%mod,a[j+k+l]=(x-y+mod)%mod;
            }
        }
    }
    if(f==-){
        int ni=power(n,mod-);
        for(int i=;i<n;i++)a[i]=1ll*a[i]*ni%mod;
    }
}
void get_inv(int *a,int *b,int n){
    if(n==){b[]=power(a[],mod-);return;}
    get_inv(a,b,n>>);
    memcpy(tmp,a,sizeof(int)*n),memset(tmp+n,,sizeof(int)*n);
    NTT(tmp,n<<,),NTT(b,n<<,);
    for(int i=;i<n<<;i++)b[i]=(1ll*b[i]*(-1ll*b[i]*tmp[i]%mod)%mod+mod)%mod;
    NTT(b,n<<,-),memset(b+n,,sizeof(int)*n);
}
void get_root(int *a,int *b,int n){
    if(n==){b[]=;return;}
    get_root(a,b,n>>),memset(d,,sizeof(int)**n);
    get_inv(b,d,n);
    memcpy(c,a,sizeof(int)*n),memset(c+n,,sizeof(int)*n);
    NTT(c,n<<,),NTT(b,n<<,),NTT(d,n<<,);
    for(int i=;i<n<<;i++)b[i]=(1ll*c[i]*d[i]%mod+b[i])%mod*%mod;
    NTT(b,n<<,-),memset(b+n,,sizeof(int)*n);
}
signed main(){
    int n,m;
    scanf("%d%d",&n,&m);
    for(len=;len<=m;len<<=);A[]=;
    for(int i=;i<=n;i++)scanf("%d",&xx),C[xx]=,A[xx]=mod-;
    get_root(A,sqrA,len),sqrA[]++,get_inv(sqrA,F,len);
    for(int i=;i<=m;i++)printf("%d\n",(F[i]<<)%mod);
}