#include using namespace std; using ll = long long; const ll MOD = 998244353; const int N = 1000100; ll fac[N], invfac[N]; ll modpow(ll a, ll b, ll MOD){ ll ret = 1; while (b > 0){ if(b & 1){ ret *= a; ret %= MOD; } a *= a; a %= MOD; b >>= 1; } return ret; } void init(){ fac[0] = 1; for(int i = 1; i < N; i++) fac[i] = fac[i - 1] * i % MOD; invfac[N - 1] = modpow(fac[N - 1], MOD - 2, MOD); for(int i = N - 2; i >= 0; i--) invfac[i] = invfac[i + 1] * (i + 1) % MOD; } ll nCk(int n, int k){ if(k < 0 || k > n) return 0; else return (fac[n] * invfac[k] % MOD) * invfac[n - k] % MOD; } void solve(){ init(); ll n, m; cin >> n >> m; ll ans = 1; int mi = m / n; int ma = mi + 1; int mac = m % n; int mic = n - m % n; while(mic--){ ans *= nCk(m, mi); ans %= MOD; m -= mi; } while(mac--){ ans *= nCk(m, ma); ans %= MOD; m -= ma; } cout << ans << endl; } int main(){ solve(); }