#include<bits/stdc++.h>
using namespace std;
const int p=998244353;
int n,m;
long long k;
struct matrix{int a[3][3];}res,X;
int pls(int x,int y){return (x+=y)>=p?x-p:x;}
matrix operator*(const matrix &x,const matrix &y){
	matrix z;
	for(int i=0;i<3;i++)
		for(int j=0;j<3;j++){
			z.a[i][j]=0;
			for(int k=0;k<3;k++)
				z.a[i][j]=(z.a[i][j]+1ll*x.a[i][k]*y.a[k][j])%p;
		}
	return z;
}
map<int,int>mp;
int dfs(int x){
	if(x==0)return 1;
	if(mp.count(x))return mp[x];
	return mp[x]=(x&1?1ll*k*dfs(x/2)%p:pls(dfs(x/2),dfs(x/2-1)));
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	cin>>n>>m;
	res.a[0][0]=res.a[1][1]=res.a[2][2]=1;
	X.a[0][1]=X.a[1][0]=X.a[1][2]=X.a[2][0]=1;
	while(n){if(n&1)res=res*X;X=X*X,n>>=1;}
	for(int i=0;i<3;i++)k=pls(k,res.a[i][i]);
	return cout<<dfs(m),0;
}