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