#define _USE_MATH_DEFINES #include using namespace std; //template #define rep(i,a,b) for(int i=(int)(a);i<(int)(b);i++) #define ALL(v) (v).begin(),(v).end() typedef long long int ll; const int inf = 0x3fffffff; const ll INF = 0x1fffffffffffffff; const double eps=1e-12; templateinline bool chmax(T& a,T b){if(ainline bool chmin(T& a,T b){if(a>b){a=b;return 1;}return 0;} //end templatestruct fp { unsigned v; static unsigned get_mod(){return mod;} unsigned inv() const{ int tmp,a=v,b=mod,x=1,y=0; while(b)tmp=a/b,a-=tmp*b,swap(a,b),x-=tmp*y,swap(x,y); if(x<0){x+=mod;} return x; } fp():v(0){} fp(ll x):v(x>=0?x%mod:mod+(x%mod)){} fp operator-()const{return fp(-v);} fp pow(ll t){fp res=1,b=*this; while(t){if(t&1)res*=b;b*=b;t>>=1;} return res;} fp& operator+=(const fp& x){if((v+=x.v)>=mod)v-=mod; return *this;} fp& operator-=(const fp& x){if((v+=mod-x.v)>=mod)v-=mod; return *this;} fp& operator*=(const fp& x){v=ll(v)*x.v%mod; return *this;} fp& operator/=(const fp& x){v=ll(v)*x.inv()%mod; return *this;} fp operator+(const fp& x)const{return fp(*this)+=x;} fp operator-(const fp& x)const{return fp(*this)-=x;} fp operator*(const fp& x)const{return fp(*this)*=x;} fp operator/(const fp& x)const{return fp(*this)/=x;} bool operator==(const fp& x)const{return v==x.v;} bool operator!=(const fp& x)const{return v!=x.v;} }; using Fp=fp<>; templatestruct factorial { vector Fact,Finv,Inv; factorial(int maxx){ Fact.resize(maxx); Finv.resize(maxx); Inv.resize(maxx); Fact[0]=Fact[1]=Finv[0]=Finv[1]=Inv[1]=1; rep(i,2,maxx){Fact[i]=Fact[i-1]*i;} Finv[maxx-1]=Fact[maxx-1].inv(); for(int i=maxx-1;i>=2;i--){Finv[i-1]=Finv[i]*i; Inv[i]=Finv[i]*Fact[i-1];} } T fact(int n,bool inv=0){if(inv)return Finv[n];else return Fact[n];} T inv(int n){return Inv[n];} T nPr(int n,int r){if(n<0||nstruct NTT{ vector rt,irt; NTT(int lg=21){ const unsigned m=T(-1).v; T prt=p; rt.resize(1<>w),ig=g.inv(); for(int i=0;i& f,bool inv=0){ int n=f.size(); if(inv){ for(int i=1;i>1;i;i>>=1)for(int j=0;j conv(vector a,vector b,bool same){ if(a.empty() and b.empty())return vector(); int n=a.size()+b.size()-1,m=1; while(m ntt; inline vector multiply(vector a,vector b,bool same=0){return ntt.conv(a,b,same);}; factorial fact(1048576); templatestruct Poly{ vector f; Poly(){} Poly(int _n):f(_n){} Poly(vector _f){f=_f;} T& operator[](const int i){return f[i];} T eval(T x){T res,w=1; for(T v:f)res+=w*v,w*=x; return res;} int size()const{return f.size();} Poly resize(int n){Poly res=*this; res.f.resize(n); return res;} void shrink(){while(!f.empty() and f.back()==0)f.pop_back();} Poly inv()const{ assert(f[0]!=0); int n=f.size(); Poly res(1); res[0]=f[0].inv(); for(int k=1;kf.size())f.resize(g.size()); rep(i,0,g.size())f[i]+=g[i]; shrink(); return *this; } Poly& operator-=(Poly g){ if(g.size()>f.size())f.resize(g.size()); rep(i,0,g.size())f[i]-=g[i]; shrink(); return *this; } Poly& operator*=(Poly g){f=multiply(f,g.f); shrink(); return *this;} Poly& operator/=(Poly g){ if(g.size()>f.size())return *this=Poly(); reverse(ALL(f)); reverse(ALL(g.f)); int n=f.size()-g.size()+1; f.resize(n); g.f.resize(n); *this*=g.inv(); f.resize(n); reverse(ALL(f)); shrink(); return *this; } Poly& operator%=(Poly g){*this-=*this/g*g; shrink(); return *this;} Poly diff(){Poly res(f.size()-1); rep(i,0,res.size())res[i]=f[i+1]*(i+1); return res;} Poly inte(){Poly res(f.size()+1); for(int i=res.size()-1;i;i--)res[i]=f[i-1]*fact.inv(i); return res;} Poly log(){ assert(f[0]==1); int n=f.size(); Poly res=diff()*inv(); res=res.inte(); return res.resize(n); } Poly exp(){ assert(f[0]==0); int n=f.size(); Poly res(1),g(1); res[0]=g[0]=1; for(int k=1;kstruct Multipoint_evaluation{ int n; vector xs; vector> buf; Multipoint_evaluation(){} Multipoint_evaluation(const vector& _xs):n(_xs.size()),xs(_xs),buf(4*n){pre(0,n,1);} void pre(int l,int r,int k){ if(r-l==1){buf[k].f={xs[l]*-1,1}; return;} int m=(l+r)>>1; pre(l,m,k*2); pre(m,r,k*2+1); buf[k]=buf[k*2]*buf[k*2+1]; } vector run(const Poly& f){ vector res(n); function,int,int,int)> dfs=[&](Poly g,int l,int r,int k){ g%=buf[k]; if(r-l<=128){rep(i,l,r)res[i]=g.eval(xs[i]); return;} int m=(l+r)>>1; dfs(g,l,m,k*2); dfs(g,m,r,k*2+1); }; dfs(f,0,n,1); return res; } Poly build(const vector& ys){ Poly w=buf[1].diff(); auto vs=run(w); function(int,int,int)> dfs=[&](int l,int r,int k){ if(r-l==1){Poly res(1); res[0]=ys[l]/vs[l]; return res;} int m=(l+r)>>1; return buf[k*2+1]*dfs(l,m,k*2)+buf[k*2]*dfs(m,r,k*2+1); }; Poly res=dfs(0,n,1); res.resize(n); return res; } }; int main(){ int n,q; cin>>n>>q; vector a(n); rep(i,0,n){ ll x; cin>>x; a[i]=Fp(-x+1); } Multipoint_evaluation mul(a); auto res=mul.buf[1]; rep(_,0,q){ int x; cin>>x; cout<