#include #include using namespace std; using namespace atcoder; using ll=long long; ll n,m; modint998244353 ans=0,j_sum; //解法 Σ[j=1→m]j×(n/j)×(n/j+3)/2 n/jの値で場合分け void clac(ll x) { // j_sum=n/jがxとなるjの和 j_sum=(min(m+1,n/(x+1)+1)+min(m,n/x))*(min(m,n/x)-min(m,n/(x+1)))/2; ans+=j_sum*x*(x+3)/2; } int main() { cin>>n>>m; for(ll i=1;i*i<=n;i++) { // n/j=iの場合 clac(i); // n/j=n/iの場合 if(n/i!=i)clac(n/i); } cout<