#pragma GCC optimize("trapv")
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int mod = 998244353;
int inv(int x) {
    return x < 2 ? x : mod - (mod / x) * 1ll * inv(mod % x) % mod;
}
vector<ll> get_pent(int n) {
    vector<ll> penta(n + 1);
    penta[0] = 1;
    penta[1] = -1;
    for(ll x = 0; x <= n; x++) {
        ll a = (3*x*x+7*x+4) / 2;
        ll b = a + 2*x + 3;
        // cout << a << " " << b << endl;
        if(a > n && b > n) break;
        if(a <= n) penta[a] = x % 2 ? 1 : -1;
        if(b <= n) penta[b] = (x+1) % 2 ? 1 : -1;
    }
    return penta;
}
int main() {
    ios_base::sync_with_stdio(0);
    cin.tie(0);
    int n, m;
    cin >> n >> m;
    // vector<ll> dp(m); dp[0] = 1;
    // for(int k = n; k >= 1; k--) {
    //     int X = inv(k) % mod;
    //     vector<ll> ndp(m);
    //     for(int i = m; i--;) {
    //         for(int t = 0; t < k && i - t >= 0; t++)
    //             ndp[i] = (ndp[i] + dp[i - t] * X) % mod;
    //     }
    //     dp = ndp;
    // }
    // ll fact = 1, v = 0;
    // for(int i = 0; i < m; i++) v = (v + dp[i]) % mod;
    // cout << v << endl;
    // vector<ll> dp(n + 1);dp[0] = 1;
    // for(int i = 1; i <= n; i++) {
    //     for(int j = n; j - i >= 0; j--)
    //         dp[j] = (dp[j] + mod - dp[j - i]) % mod;
    // }
    auto dp = get_pent(n + 1);

    vector<ll> fact(2*n + 1, 1), ifact(2*n + 1, 1);
    for(int i = 2; i <= 2*n; i++) {
        ifact[i] = mod - (mod / i) * ifact[mod % i] % mod;
    }
    for(int i = 2; i <= 2*n; i++) {
        fact[i] = fact[i-1] * 1ll * i % mod;
        ifact[i] = ifact[i - 1] * 1ll * ifact[i] % mod;
    }
    auto C = [&](int n, int k) -> ll{
        if(k < 0 || k > n) return 0;
        ll res = fact[n] * ifact[k] % mod;
        return res * ifact[n - k] % mod;
    };

    int T = n + 1;
    // for(int T = n + 1; T--;) {
    //     for(int i = 1; i <= n; i++) {
    //         // assert(i < dp.size());
    //         (dp[i] += dp[i - 1]) %=mod;
    //     }
    // }
    ll v = 0;m--;
    for(int i = 0; i <= m; i++) {
        // cout << C(i+T-1, i) << " " << dp[m - i] << endl;
        v = (v + C(i+T-1, i) * dp[m - i]) % mod;
    }
    v = v * ifact[n] % mod;
    v = (mod + v) % mod;
    cout << v << endl;
}