#include #include using namespace std; using lint = long long; using mint = atcoder::modint998244353; int main() { lint n, m; cin >> n >> m; lint sq = sqrtl(m); mint ans = 0; for (lint i = 1; i <= min(sq, n); i++) { ans += mint(m % i); } auto calc = [](lint fst, lint lst, lint q) -> mint { lint cnt = (fst - lst) / q + 1; lint res = (fst + lst) * cnt / 2; return mint(res); }; for (lint q = m / (sq + 1); q >= 1 && n > sq; q--) { if (q == m / n) { lint fst = m % (m / (q + 1) + 1); lint lst = m % n; ans += calc(fst, lst, q); break; } else { lint fst = m % (m / (q + 1) + 1); lint lst = m % q; ans += calc(fst, lst, q); } } ans += mint(max(0LL, n - m)) * mint(m); cout << ans.val() << endl; }