#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
const int mod=998244353;
int n,m;
struct Mat {
	int a[3][3];
	void init() {memset(a,0,sizeof a);}
};
Mat T={
	0,1,1,
	1,0,0,
	0,1,0
};
Mat Mul(Mat A,Mat B) {
	Mat C; C.init();
	for(int i=0; i<3; i++) for(int k=0; k<3; k++) for(int j=0; j<3; j++)
		C.a[i][j]=(C.a[i][j]+1ll*A.a[i][k]*B.a[k][j])%mod;
	return C;
}
int dp[31][2];
int main() {
	scanf("%d%d",&n,&m);
	int k=n-2,res=0;
	Mat C; C.init(); for(int i=0; i<3; i++) C.a[i][i]=1;
	for(; k; k>>=1,T=Mul(T,T)) if(k&1) C=Mul(C,T);
	res=(1ll*C.a[0][0]+C.a[0][2]+C.a[1][0]+C.a[1][1]+C.a[2][1])%mod;
	if(m&1) dp[0][0]=res;
	else dp[0][0]=1,dp[0][1]=1;
	for(int i=1; i<=30; i++) {
		if(m&(1<<i)) dp[i][0]=(1ll*res*dp[i-1][0]%mod+dp[i-1][1])%mod,dp[i][1]=dp[i-1][1];
		else dp[i][0]=dp[i-1][0],dp[i][1]=(1ll*res*dp[i-1][1]+dp[i-1][0])%mod;
	}
	printf("%d\n",dp[30][0]);
	return 0;
}