#include #include #include using namespace std; typedef long long ll; template struct static_modint { using mint = static_modint; public: int val; static_modint() : val(0) {} static_modint(long long v) { if (abs(v) >= mod()) { v %= mod(); } if (v < 0) { v += mod(); } val = v; } mint &operator++() { val++; if (val == mod()) { val = 0; } return *this; } mint &operator--() { if (val == 0) { val = mod(); } val--; return *this; } mint &operator+=(const mint &x) { val += x.val; if (val >= mod()) { val -= mod(); } return *this; } mint &operator-=(const mint &x) { val -= x.val; if (val < 0) { val += mod(); } return *this; } mint &operator*=(const mint &x) { val = (int)((long long)val * x.val % mod()); return *this; } mint &operator/=(const mint &x) { *this *= x.inv(); return *this; } mint operator-() { return mint() - *this; } mint pow(long long n) const { mint x = 1, r = *this; while (n) { if (n & 1) { x *= r; } r *= r; n >>= 1; } return x; } mint inv() const { return pow(mod() - 2); } friend mint operator+(const mint &x, const mint &y) { return mint(x) += y; } friend mint operator-(const mint &x, const mint &y) { return mint(x) -= y; } friend mint operator*(const mint &x, const mint &y) { return mint(x) *= y; } friend mint operator/(const mint &x, const mint &y) { return mint(x) /= y; } friend bool operator==(const mint &x, const mint &y) { return x.val == y.val; } friend bool operator!=(const mint &x, const mint &y) { return x.val != y.val; } friend std::ostream &operator<<(std::ostream &os, const mint &x) { return os << x.val; } friend std::istream &operator>>(std::istream &is, mint &x) { long long v; is >> v; x = mint(v); return is; } private: static constexpr int mod() { return MOD; } }; template struct combination { public: combination(int maxn) : maxn(maxn) { fact.resize(maxn + 1); ifact.resize(maxn + 1); fact[0] = 1; for (int i = 1; i <= maxn; i++) { fact[i] = fact[i - 1] * i; } ifact[maxn] = fact[maxn].inv(); for (int i = maxn - 1; i >= 0; i--) { ifact[i] = ifact[i + 1] * (i + 1); } } M operator()(int n, int k) { assert(n <= maxn); if (k > n || k < 0) { return 0; } return fact[n] * ifact[k] * ifact[n - k]; } private: int maxn; std::vector fact; std::vector ifact; }; using mint = static_modint<998244353>; int main() { int n, m; cin >> n >> m; combination com(m); int s = m; mint ans = 1; for (int i = 0; i < n; i++) { if (i < m % n) { ans *= com(s, m / n + 1); s -= m / n + 1; } else { ans *= com(s, m / n); s -= m / n; } } cout << ans << endl; }