POJ 2800

#include<iostream>
#include<cmath>
using namespace std;

int main()
{
long long N,K;
while(cin>>N>>K)
{
long long S,T,i,s,e,ans=0;
S=sqrt((double)K);
T=K/S;
for(i=1,ans=0;i<=N&&i<=T;i++)
ans+=K%i;
if(N>K)
ans+=(N-K)*K;
for(i=S;i>1;i--)
{
s=K/i;
e=K/(i-1);
if(s>N)
break;
if(e>N)
e=N;
ans=ans+(e-s)*(K%e+K%(s+1))/2;
}
cout<<ans<<endl;
}
return 0;
}