#include<stdio.h>
#include<algorithm>
typedef unsigned int uint;
typedef unsigned long long ull;
constexpr uint mod{998244353};
constexpr uint plus(const uint &x,const uint &y) {
    if(x+y>=mod) {
        return x+y-mod;
    }
    return x+y;
}
template<int N>
struct vector {
    uint a[N];
    constexpr vector():a{} {}
};
template<uint N>
struct matrix {
    uint a[N][N];
    constexpr matrix():a{} {}
    constexpr matrix<N> operator*(const matrix<N> &x) const {
        matrix<N> y;
        for(uint i=0;i!=N;i++) {
            for(uint j=0;j!=N;j++) {
                for(uint k=0;k!=N;k++) {
                    y.a[i][k]=(y.a[i][k]+ull(this->a[i][j])*x.a[j][k])%mod;
                }
            }
        }
        return y;
    }
    constexpr vector<N> operator*(const vector<N> &x) const {
        vector<N> y;
        for(uint i=0;i!=N;i++) {
            for(uint j=0;j!=N;j++) {
                y.a[i]=(y.a[i]+ull(this->a[i][j])*x.a[j])%mod;
            }
        }
        return y;
    }
};
int main() {
    int n,m;
    scanf("%d%d",&n,&m);
    uint coef[3]{1,0,1};
    matrix<4> a,b;
    for(int i=0;i!=4;i++) {
        for(int j=0;j!=2;j++) {
            if((i&j&1)+(i>>1)==1) {
                a.a[i][i>>1|j<<1]=1;
            }
        }
        b.a[i][i]=1;
    }
    for(int i=n-2;i;i>>=1) {
        if(i&1) {
            b=a*b;
        }
        a=a*a;
    }
    for(int i=0;i!=4;i++) {
        for(int j=0;j!=4;j++) {
            if((i&j&1)+(j>>1)==1&&(j>>1&i>>1)+(i&1)==1) {
                coef[1]=plus(coef[1],b.a[i][j]);
            }
        }
    }
    uint dp[2]{};
    dp[0]=1;
    for(int i=0;1<<i<=m;i++) {
        uint dp2[2]{};
        for(int j=0;j!=2;j++) {
            for(int k=0;k<=2;k++) {
                if((j+k&1)==(m>>i&1)) {
                    dp2[j+k>>1]=(dp2[j+k>>1]+ull(coef[k])*dp[j])%mod;
                }
            }
        }
        std::copy(dp2,dp2+2,dp);
    }
    printf("%u\n",dp[0]);
    return 0;
}