#include using namespace std; #define rep(i,a,b) for(int i=(int)(a);i<(int)(b);i++) #define ALL(v) (v).begin(),(v).end() using ll=long long int; 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;} templatestruct fp { int v; static int get_mod(){return mod;} int 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(ll x=0){init(x%mod+mod);} fp& init(int x){v=(x>=1;} return res;} fp& operator+=(const fp& x){return init(v+x.v);} fp& operator-=(const fp& x){return init(v+mod-x.v);} 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;} friend istream& operator>>(istream& is,fp& x){return is>>x.v;} friend ostream& operator<<(ostream& os,const fp& x){return os<; 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(n<0)return 0; return (inv?Finv[n]:Fact[n]);} T inv(int n){if(n<0)return 0; return Inv[n];} T nPr(int n,int r,bool inv=0){if(n<0||nstruct NTT{ vector rt,irt; NTT(int lg=21){ unsigned m=T::get_mod()-1; T prt=p; rt.resize(lg); irt.resize(lg); rep(k,0,lg){ rt[k]=-prt.pow(m>>(k+2)); irt[k]=rt[k].inv(); } } void ntt(vector& f,bool inv=0){ int n=f.size(); if(inv){ for(int m=1;m>=1;){ T w=1; for(int s=0,t=0;s mult(const vector& a,const vector& b,bool same=0){ if(a.empty() or b.empty())return vector(); int n=a.size()+b.size()-1,m=1<<__lg(n*2-1); vector res(m); rep(i,0,a.size()){res[i]=a[i];} ntt(res); if(same)rep(i,0,m)res[i]*=res[i]; else{ vector c(m); rep(i,0,b.size())c[i]=b[i]; ntt(c); rep(i,0,m)res[i]*=c[i]; } ntt(res,1); res.resize(n); return res; } }; NTT ntt; vector mult(const vector& a,const vector& b,bool same=0){ return ntt.mult(a,b,same); } factorial fact(2010101); templatestruct Poly:vector{ Poly(int n=0){this->assign(n,T());} Poly(const vector& f){this->assign(ALL(f));} T eval(const T& x){ T res; for(int i=this->size()-1;i>=0;i--)res*=x,res+=this->at(i); return res; } Poly rev()const{Poly res=*this; reverse(ALL(res)); return res;} void shrink(){while(!this->empty() and this->back()==0)this->pop_back();} Poly inv()const{ assert(this->front()!=0); const int n=this->size(); Poly res(1); res.front()=T(1)/this->front(); for(int k=1;kthis->size())this->resize(g.size()); rep(i,0,g.size()){(*this)[i]+=g[i];} shrink(); return *this; } Poly& operator-=(const Poly& g){ if(g.size()>this->size())this->resize(g.size()); rep(i,0,g.size()){(*this)[i]-=g[i];} shrink(); return *this; } Poly& operator*=(const Poly& g){ *this=mult(*this,g,0); shrink(); return *this; } Poly& operator/=(const Poly& g){ if(g.size()>this->size()){ this->clear(); return *this; } Poly g2=g; reverse(ALL(*this)); reverse(ALL(g2)); int n=this->size()-g2.size()+1; this->resize(n); g2.resize(n); *this*=g2.inv_fast(); this->resize(n); // reverse(ALL(*this)); shrink(); return *this; } Poly& operator%=(const Poly& g){*this-=*this/g*g; shrink(); return *this;} Poly diff()const{ Poly res(this->size()-1); rep(i,0,res.size())res[i]=(*this)[i+1]*(i+1); return res; } Poly inte()const{ Poly res(this->size()+1); for(int i=res.size()-1;i;i--)res[i]=(*this)[i-1]*fact.inv(i); return res; } Poly log()const{ assert(this->front()==1); const int n=this->size(); Poly res=diff()*inv_fast(); res=res.inte(); // res.resize(n); return res; } Poly exp()const{ assert(this->front()==0); const int n=this->size(); Poly res(1),g(1); res.front()=g.front()=1; for(int k=1;ksize(); Poly res=*this,g(n); g[0]=1; rep(i,1,n)g[i]=g[i-1]*c*fact.inv(i); rep(i,0,n){res[i]*=fact.fact(i);} res=res.rev(); res*=g; res.resize(n); res=res.rev(); rep(i,0,n){res[i]*=fact.fact(i,1);} return res; } Poly inv_fast()const{ const int n=this->size(); Poly res(1); res.front()=T(1)/this->front(); for(int k=1;ksize(); if(n==1)return Poly({T(1)}); Poly b(2),c(1),z1,z2(2); b[0]=c[0]=z2[0]=z2[1]=1; b[1]=(*this)[1]; for(int k=2;k>1)z[i]=0; ntt.ntt(z); rep(i,0,k)z[i]*=-z1[i]; ntt.ntt(z,1); c.insert(c.end(),z.begin()+(k>>1),z.end()); z2=c; z2.resize(k*2); ntt.ntt(z2); Poly x=*this; x.resize(k); x=x.diff(); x.resize(k); ntt.ntt(x); rep(i,0,k)x[i]*=y[i]; ntt.ntt(x,1); Poly bb=b.diff(); rep(i,0,k-1)x[i]-=bb[i]; x.resize(k*2); rep(i,0,k-1){x[k+i]=x[i]; x[i]=0;} ntt.ntt(x); rep(i,0,k*2)x[i]*=z2[i]; ntt.ntt(x,1); x.pop_back(); x=x.inte(); rep(i,k,min(n,k*2))x[i]+=(*this)[i]; rep(i,0,k)x[i]=0; ntt.ntt(x); rep(i,0,k*2)x[i]*=y[i]; ntt.ntt(x,1); b.insert(b.end(),x.begin()+k,x.end()); } b.resize(n); return b; } Poly pow(ll t){ int n=this->size(),k=0; while(k=n)return res; n-=t*k; Poly g(n); T c=(*this)[k],ic=T(1)/c; rep(i,0,n)g[i]=(*this)[i+k]*ic; g=g.log(); for(auto& x:g)x*=t; g=g.exp_fast(); // c=c.pow(t); rep(i,0,n)res[i+t*k]=g[i]*c; return res; } }; Fp nth(Poly p,Poly q,ll n){ while(n){ Poly base(q),np,nq; for(int i=1;i<(int)q.size();i+=2)base[i]=-base[i]; p*=base; q*=base; for(int i=n&1;i<(int)p.size();i+=2)np.emplace_back(p[i]); for(int i=0;i<(int)q.size();i+=2)nq.emplace_back(q[i]); swap(p,np); swap(q,nq); n>>=1; } return p[0]/q[0]; } vector chi; Poly MultivariateConvolution(const Poly& f,const Poly& g,int k){ int n=f.size(),m=1<<__lg(4*n-1); if(k==0)return Poly({f[0]*g[0]}); vector F(k,Poly(m)),G(k,Poly(m)); rep(i,0,n){ F[chi[i]][i]=f[i]; G[chi[i]][i]=g[i]; } for(auto& v:F)ntt.ntt(v); for(auto& v:G)ntt.ntt(v); rep(x,0,m){ Poly tmp(k); rep(i,0,k)rep(j,0,k){ tmp[(i+j)%k]+=F[i][x]*G[j][x]; } rep(i,0,k)F[i][x]=tmp[i]; } for(auto& v:F)ntt.ntt(v,1); Poly res(n); rep(i,0,n)res[i]=F[chi[i]][i]; return res; } const int ten[]={1,10,100,1000,10000,100000}; const Fp g=Fp(6).pow((Fp::get_mod()-1)/10); const Fp ig=g.inv(); const Fp inv10=Fp(10).inv(); void ntt10(vector& a){ Fp cur=1; vector b(10); rep(i,0,10){ for(int j=a.size()-1;j>=0;j--){ b[i]*=cur; b[i]+=a[j]; } cur*=g; } a=b; } void intt10(vector& a){ Fp cur=1; vector b(10); rep(i,0,10){ for(int j=a.size()-1;j>=0;j--){ b[i]*=cur; b[i]+=a[j]; } cur*=ig; } rep(i,0,10)b[i]*=inv10; a=b; } Fp nCr(ll n,int r){ Fp ret=1; for(ll x=n-r+1;x<=n;x++)ret*=x; ret*=fact.fact(r,1); return ret; } void Pow(Poly& f,ll t,int k){ Poly ret(f.size()),add(f.size()); // ret[0]=1; // while(t){ // if(t&1){ // ret=MultivariateConvolution(ret,f,a); // } // f=MultivariateConvolution(f,f,a); // t>>=1; // } Fp f0=f[0]; const Fp invf0=f0.inv(); Fp coeff=f0.pow(t); f[0]-=f0; add[0]=1; rep(d,0,50){ if(t>n>>k>>m>>t; { chi.resize(ten[t]); vector a(t,10); rep(x,0,ten[t]){ int y=x; rep(i,0,k-1){ y/=a[i]; chi[x]+=y; } chi[x]%=t; } } Poly a(ten[k]); rep(_,0,n){ int x; cin>>x; a[x]+=1; } rep(d,t,k){ rep(base,0,ten[k]){ if((base/ten[d])%10)continue; vector tmp(10); rep(x,0,10)tmp[x]=a[x*ten[d]+base]; ntt10(tmp); rep(x,0,10)a[x*ten[d]+base]=tmp[x]; } } rep(base,0,ten[k-t]){ Poly tmp(ten[t]); rep(i,0,ten[t])tmp[i]=a[i+base*ten[t]]; Pow(tmp,m,t); rep(i,0,ten[t])a[i+base*ten[t]]=tmp[i]; } rep(d,t,k){ rep(base,0,ten[k]){ if((base/ten[d])%10)continue; vector tmp(10); rep(x,0,10)tmp[x]=a[x*ten[d]+base]; intt10(tmp); rep(x,0,10)a[x*ten[d]+base]=tmp[x]; } } rep(i,0,ten[k])cout<