[模板]线性递推+BM

暴力版本：

 #include<bits/stdc++.h>
 #define mod 998244353
 using namespace std;
 typedef long long int ll;
 const int maxn=1E4+;
 const int inf=;

 ll A[maxn],B[maxn],tmp[maxn];
 vector<ll>wait[maxn];

 inline ll qpow(ll x,ll y)
 {
     ll ans=,base=x;
     while(y)
     {
         if(y&)
             ans=ans*base%mod;
         base=base*base%mod;
         y>>=;
     }
     return ans;
 }
 struct linear
 {
     vector<ll>val,a;
     int fail[maxn],cnt;
     ll delta[maxn];
     inline void mul(ll*A,ll*B)
     {
         int k=a.size();
         for(int i=;i<*k;++i)
             tmp[i]=;
         for(int i=;i<k;++i)
             for(int j=;j<k;++j)
                 tmp[i+j]=(tmp[i+j]+A[i]*B[j])%mod;
         for(int i=*k-;i>=k;--i)
             for(int j=;j<k;++j)
                 tmp[i-j-]=(tmp[i-j-]+tmp[i]*a[j])%mod;
         for(int i=;i<k;++i)
             A[i]=tmp[i];
     }
     void qpowA(ll*A,ll*B,ll k)
     {
         while(k)
         {
             if(k&)
                 mul(A,B);
             mul(B,B);
             k>>=;
         }
     }
     ll quickRecursion(ll n)
     {
         int k=a.size();
         if(n<k)
             return val[n];
         for(int i=;i<k;++i)
             A[i]=B[i]=;
         A[]=;
         B[]=;
         qpowA(A,B,n);
         ll ans=;
         for(int i=;i<k;++i)
             ans=(ans+A[i]*val[i])%mod;
         return ans;
     }
     void BM()
     {
         cnt=;
         for(int i=;i<(int)val.size();++i)
             wait[i].clear();
         for(int i=;i<(int)val.size();++i)
         {
             delta[i]=val[i];
             for(int j=;j<(int)wait[cnt].size();++j)
                 delta[i]=(delta[i]-wait[cnt][j]*val[i-j-]%mod+mod)%mod;
             if(delta[i]==)
                 continue;
             fail[cnt]=i;
             if(cnt==)
             {
                 wait[++cnt].resize(i+);
                 continue;
             }
             int len=i-fail[cnt-]+wait[cnt-].size(),pos=cnt-;
             for(int j=;j<cnt;++j)
                 if(i-fail[j]+(int)wait[j].size()<len)
                     len=i-fail[j]+wait[j].size(),pos=j;
             ll g=delta[i]*qpow(delta[fail[pos]],mod-)%mod;
             ++cnt;
             wait[cnt]=wait[cnt-];
             while((int)wait[cnt].size()<len)
                 wait[cnt].push_back();
             wait[cnt][i-fail[pos]-]=(wait[cnt][i-fail[pos]-]+g)%mod;
             for(int j=;j<(int)wait[pos].size();++j)
                 wait[cnt][i-fail[pos]+j]=(wait[cnt][i-fail[pos]+j]-g*wait[pos][j]%mod+mod)%mod;
         }
         a.clear();
         a=wait[cnt];
         for(int i=;i<(int)a.size();++i)
             cout<<a[i]<<" ";
         cout<<endl;
     }
 }T;
 int main()
 {
     ios::sync_with_stdio(false);
     ll n,m;
     cin>>n>>m;
     for(int i=;i<=n;++i)
     {
         ll x;
         cin>>x;
         T.val.push_back(x);
     }
     T.BM();
     cout<<T.quickRecursion(m)<<endl;
     return ;
 }