#include #include #include using namespace atcoder; using mint = modint; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf 1000000000000000000 template vector convolution_mint(vector a,vector b){ static constexpr unsigned long long M0 = 998244353; static constexpr unsigned long long M1 = 754974721; static constexpr unsigned long long M2 = 469762049; vector aa(a.size()),bb(b.size()); rep(i,a.size())aa[i] = a[i].val(); rep(i,b.size())bb[i] = b[i].val(); auto c0 = convolution(aa,bb); auto c1 = convolution(aa,bb); auto c2 = convolution(aa,bb); vector ret(c0.size()); vector m = {M0,M1,M2}; rep(i,c0.size()){ ret[i] += c0[i]; { long long cur = c0[i]; cur %= M1; cur = c1[i]-cur; if(cur<0)cur += M1; cur *= 416537774; cur %= M1; mint m = M0; ret[i] += m*cur; cur *= M0; cur += c0[i]; cur %= M2; cur = c2[i] - cur; if(cur<0)cur += M2; cur *= 429847628; cur %= M2; m *= M1; ret[i] += m * cur; } } return ret; } int main(){ mint::set_mod(1234567891); int n; cin>>n; long long m; cin>>m; vector a(n); rep(i,n)cin>>a[i]; vector dp(25005,0); dp[0] = 1; rep(i,n){ vector ndp(25005,0); rep(j,dp.size()){ if(dp[j]==0)continue; ndp[j] += dp[j]; ndp[j+a[i]] += dp[j]; } swap(dp,ndp); } vector ddp(25005,0); ddp[0] = 1; while(m!=0){ vector ndp = convolution_mint(dp,ddp); vector temp(25005,0); rep(i,ndp.size()){ if(i%2!=m%2)continue; temp[i>>1] = ndp[i]; } swap(temp,ddp); m>>=1; } cout<