#include<bits/stdc++.h>
using namespace std;
const int mod=998244353,ggg=3,invggg=(mod+1)/ggg;
int n,d;
struct matr{
	int n,m;
	vector<vector<int> >v;
	void init(int nn,int mm){
		n=nn;m=mm;
		vector<int>q;
		for(int i=0;i<m;i++)
			q.push_back(0);
		for(int i=0;i<n;i++)
			v.push_back(q);
	}
};
matr mul(matr a,matr b){
	matr ret;ret.init(a.n,b.m);
	for(int i=0;i<a.n;i++)
		for(int k=0;k<b.n;k++)
			for(int j=0;j<b.m;j++)
				ret.v[i][j]=(ret.v[i][j]+1ll*a.v[i][k]*b.v[k][j])%mod;
	return ret;
}
matr Mypow(matr a,int times){
	if(times==1)
		return a;
	matr ret=Mypow(a,times/2);
	ret=mul(ret,ret);
	if(times&1)
		ret=mul(ret,a);
	return ret;
}
const int maxn=262144;
int tmp[maxn],F[maxn],G[maxn];
long long mypow(long long x,int y){
	long long ret=1;
	while(y){
		if(y&1)ret=ret*x%mod;
		y>>=1;x=x*x%mod;
	}return ret;
}
int revs[maxn];
void initrev(int Len){
	for(int i=1;i<Len;++i){
		revs[i]=revs[i>>1]>>1;
		if(i&1)revs[i]|=(Len>>1);
	}
}
void DFT(int *f,int n,int rev){
	int g0=rev==1?ggg:invggg;
	initrev(n);
	for(register int i=0;i<n;++i)
		if(i<revs[i])
			f[i]^=f[revs[i]]^=f[i]^=f[revs[i]];
	for(register int h=2;h<=n;h<<=1){
		int gn=mypow(g0,(mod-1)/h);
		for(register int j=0;j<n;j+=h){
			int gk=1;
			for(register int k=j;k<j+(h>>1);++k){
				int e=f[k],o=1ll*gk*f[k+(h>>1)]%mod;
				f[k]=(e+o)%mod,f[k+(h>>1)]=(e+mod-o)%mod;
				gk=1ll*gk*gn%mod;
			}
		}
	}
	if(rev==-1){
		int invv=mypow(n,mod-2);
		for(register int i=0;i<n;++i)
			f[i]=1ll*f[i]*invv%mod; 
	}
}
long long dpn[2][100003];
vector<int>mat[2][2][2][2],dmat[2][2][2][2];
int f[2][200003],g[2][200003],h[2][200003],ff[262144],gg[262144],hh[262144];
void init(){
	long long nw=1,x=(n+1)/2,y=(n+1)/2;
	for(int i=0;i<=d;i++){
		f[0][i]=nw;
		if(i%2==0)
			g[0][i]=nw;
		else
			h[0][i]=nw;
		if(i<=n+1)
			if((n+1-i)%2==0)
				nw=nw*y%mod,
				y--,
				nw=nw*mypow(x-y,mod-2)%mod;
			else
				x++,
				nw=nw*x%mod,
				nw=nw*mypow(x-y,mod-2)%mod;
		else break;
	}nw=1;x=n/2;y=n/2;
	for(int i=0;i<=d;i++){
		f[1][i]=nw;
		if(i%2==0)
			g[1][i]=nw;
		else
			h[1][i]=nw;
		if(i<=n)
			if((n-i)%2==0)
				nw=nw*y%mod,
				y--,
				nw=nw*mypow(x-y,mod-2)%mod;
			else
				x++,
				nw=nw*x%mod,
				nw=nw*mypow(x-y,mod-2)%mod;
		else break;
	}
	for(int i=0;i<=d;i++)
		F[i]=f[0][i],
		G[i]=f[1][i];
	DFT(F,maxn,1);
	DFT(G,maxn,1);
	for(int i=0;i<maxn;i++)
		F[i]=1ll*F[i]*G[i]%mod;
	DFT(F,maxn,-1);
	for(int i=0;i<maxn;i++)
		ff[i]=F[i];
	for(int i=0;i<maxn;i++)F[i]=0,G[i]=0;
	for(int i=0;i<=d;i++)
		F[i]=g[0][i],
		G[i]=h[1][i];
	DFT(F,maxn,1);
	DFT(G,maxn,1);
	for(int i=0;i<maxn;i++)
		F[i]=1ll*F[i]*G[i]%mod;
	DFT(F,maxn,-1);
	for(int i=0;i<maxn;i++)
		gg[i]=F[i];
	for(int i=0;i<maxn;i++)F[i]=0,G[i]=0;
	for(int i=0;i<=d;i++)
		F[i]=g[1][i],
		G[i]=h[0][i];
	DFT(F,maxn,1);
	DFT(G,maxn,1);
	for(int i=0;i<maxn;i++)
		F[i]=1ll*F[i]*G[i]%mod;
	DFT(F,maxn,-1);
	for(int i=0;i<maxn;i++)
		hh[i]=F[i];
}
int calc(int D,int qaq){
	if(qaq==1)
		return (ff[D]-hh[D]+mod)%mod;
	else
		return (ff[D]-gg[D]+mod)%mod;
}
int QAQ(int x,int y){
	if(x==0)
		return min(1,y);
	return y;
}
void calcsmall(){
	matr fr,time0,time1;
	fr.init(1,8);
	fr.v[0][4]=1;
	time0.init(8,8);
	time1.init(8,8);
	for(int i=0;i<2;i++)
		for(int j=0;j<4;j++)
			for(int k=0;k<2;k++)
				for(int l=0;l<=min(1,j);l++)
					time0.v[k*4+j-l][i*4+j]=QAQ(j-l,mat[i][0][j%2][k][l]);
	for(int i=0;i<2;i++)
		for(int j=0;j<4;j++)
			for(int k=0;k<2;k++)
				for(int l=0;l<=min(1,j);l++)
					time1.v[k*4+j-l][i*4+j]=QAQ(j-l,mat[i][1][j%2][k][l]);
	time0.v[4][0]=1;
	time1.v[0][4]=1;
	if(n/2)
		fr=mul(fr,Mypow(mul(time0,time1),n/2));
	if(n%2)
		fr=mul(fr,time0);
	for(int i=0;i<8;i++)
		dpn[i/4][i%4]=fr.v[0][i];
}
int main(){
	cin>>n>>d;
	for(int i=0;i<2;i++)
		for(int j=0;j<2;j++){
			for(int k=0;k<2;k++)
				for(int l=0;l<2;l++){
					for(int times=0;times<3;times++)
						mat[k][i][j][l].push_back(0);
					for(int times=0;times<5;times++)
						dmat[k][i][j][l].push_back(0);
				}
			if(j){
				mat[0][i][j][0][0]=1;
				mat[0][i][j][1^i][1]=2;
				mat[1][i][j][1][0]=1;
				mat[1][i][j][1^i][1]=2;
			}else{
				mat[0][i][j][0][0]=1;
				mat[1][i][j][1][0]=1;
				mat[i][i][j][1^i][1]=1;
				mat[i][i][j][1^i][2]=1;
			}
		} 
	init();
	calcsmall();
	for(int i=0;i<2;i++)
		for(int j=0;j<2;j++){
			if(j){
				mat[0][i][j][1^i][1]=0;
				mat[1][i][j][1^i][1]=0;
				mat[0][i][j][i][1]=2;
				mat[1][i][j][i][1]=2;
			}else{
				mat[i][i][j][i^1][1]=0;
				mat[i][i][j][i^1][2]=0;
				mat[i^1][i][j][i][1]=1;
				mat[i^1][i][j][i][2]=1;
			}
		}
	for(int i=0;i<2;i++)
		for(int j=0;j<2;j++)
			for(int k=0;k<2;k++){
				for(int l=0;l<2;l++)
					for(int o=0;o<3;o++)
						for(int ll=0;ll<2;ll++)
							for(int oo=0;oo<3;oo++)
								dmat[i][j][k][ll][o+oo]+=mat[i][j][k][l][o]*mat[l][j^1][(k-o+2)%2][ll][oo];
				dmat[i][j][k][i][0]--;
				for(int ll=0;ll<2;ll++)
					for(int oo=1;oo<3;oo++)
						dmat[i][j][k][ll][oo]-=mat[1^i][j^1][k][ll][oo];
			}
	for(int D=4;D<=d;D++)
		for(int k=0;k<2;k++){
			for(int i=0;i<2;i++)
				for(int j=1;j<5;j++)
					if(D-j||j!=4)
						dpn[k][D]+=dpn[i][D-j]*dmat[k^1][n%2][D%2][i][j];
			dpn[k][D]-=calc(D-1,k);
			dpn[k][D]=(dpn[k][D]%mod+mod)%mod;
		}
	int ans=(dpn[0][d]+dpn[1][d])%mod;
	if(d%2)ans=1ll*ans*(mod+1)/2%mod;
	cout<<ans; 
}