#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;

constexpr long long mod = 998244353;

long long pow_mod(long long x, long long y, long long mod) {
    long long result = 1;
    while (y > 0) {
        if (y & 1) {
            result *= x;
            result %= mod;
        }
        x *= x;
        x %= mod;
        y >>= 1;
    }
    return result;
}

template<typename T>
void print(const T t, bool ln = true) {
    cout << t;
    if (ln) {
        cout << '\n';
    } else {
        cout << ' ';
    }
}

template<typename T>
void vprint(const vector<T> &v) {
    for (int i = 0; i < (int)v.size(); i++) {
        print(v[i], false);
    }
    cout << '\n';
}

template<typename T>
void vvprint(const vector<vector<T>> &vv) {
    for (int i = 0; i < (int)vv.size(); i++) {
        vprint(vv[i]);
    }
}

void solve() {
    long long n, m;
    cin >> n >> m;
    long long mx = max(n, m);
    vector<long long> fact(mx + 1, 1), fact_inv(mx + 1, 1);

    for (int i = 2; i <= mx; i++) {
        fact[i] = fact[i - 1] * i;
        fact[i] %= mod;
    }

    fact_inv[mx] = pow_mod(fact[mx], mod - 2, mod);
    for (int i = mx - 1; i > 1; i--) {
        fact_inv[i] = fact_inv[i + 1] * (i + 1);
        fact_inv[i] %= mod;
    }

    auto comb = [&](long long n, long long k) {
        if (n < k) {
            return 1LL;
        } else {
            return fact[n] * fact_inv[k] % mod * fact_inv[n - k] % mod;
        }
    };

    long long ans = 1;
    while (n > 0 && m > 0) {
        long long k = (m / n) + ((m % n) > 0);
        ans *= comb(m, k);
        ans %= mod;
        m -= k;
        n--;
    }

    print(ans);
}

int main() {
    int q = 1;
    // cin >> q;
    while (q--) {
        solve();
    }
}