#include #include using Fp = atcoder::modint998244353; void solve(std::istream& is, std::ostream& os) { int C, X; is >> C >> X; if (X == 0) { os << Fp(C + 1).inv().val() << '\n'; } else if (X <= C) { os << (1 - Fp(X + 1).inv()).val() << '\n'; } else { std::vector f{1}; for (int i = C + 1; i <= X; ++i) { std::vector nf(f.size() + 1); Fp r = Fp(i).inv(); for (int j = 0; j < int(f.size()); ++j) { nf[j] += f[j] * j * r; nf[j + 1] += f[j] * (i - j) * r; } f = std::move(nf); } Fp ans = 0; for (int i = 1; i < int(f.size()); ++i) { ans += f[i] * Fp(X - i + 1).inv(); } ans = 1 - ans; os << ans.val() << '\n'; } } int main() { std::ios::sync_with_stdio(false); std::cin.tie(nullptr); solve(std::cin, std::cout); }