#include #include using namespace std; using namespace atcoder; typedef long long ll; using mint=modint998244353; int main(){ int N,M; cin>>N>>M; vector> X(15,vector(0)); vector>> Y(15,vector>(0)); X[0].assign(2*M+1,1); for(int i=0;i<14;i++){ X[i+1]=convolution(X[i],X[i]); }for(int i=0;i<15;i++){ Y[i].assign(1<<(15-i),vector(1,1)); } Y[14][1]=X[14]; vector t(17,0); t[1]=1; for(int i=2;i<17;i++){ t[i]=t[i-1]*2; } for(int i=13;i>=0;i--){ for(int j=0;j<(int)(Y[i+1].size());j++){ Y[i][j<<1]=Y[i+1][j]; Y[i][(j<<1)+1]=convolution(Y[i+1][j],X[i]); vector Z((t[i+2]-1)*2*M+1); if(Y[i][j<<1].size()>=Z.size()){ for(int k=0;k<(int)(Z.size());k++){ Z[k]=Y[i][j<<1][(Y[i][j<<1].size()>>1)+k-(Z.size()>>1)]; } Y[i][j<<1]=Z; }if(Y[i][(j<<1)+1].size()>=Z.size()){ for(int k=0;k<(int)(Z.size());k++){ Z[k]=Y[i][(j<<1)+1][(Y[i][(j<<1)+1].size()>>1)+k-(Z.size()>>1)]; } Y[i][(j<<1)+1]=Z; } } } mint ANS=0; mint P=2*M+1; mint V=P.pow(N); mint W=P.inv(); mint U=1; for(int i=1;i<=N;i++){ V*=W; U*=P; ANS+=(U-Y[0][i][M])*V*(N-i+1); } cout<<(ANS/2).val()<<"\n"; }