就是要计算sum{k%i|1<=i<=n},对于i>k的情况很容易处理。所以只要可以快速计算出sum{k%i|1<=i<=k}即可,额……式子太难打了,具体的看代码吧,反正可以通过公式把复杂度降低到sqrt(k)。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
typedef long long LL;
LL n,k,m;
int main()
{
while(scanf("%lld%lld",&n,&k)!=EOF)
{
LL ans;
if(n>k) ans=(n-k)*k;
else ans=0;
m=sqrt(1.0*k);
for(int i=1;i<=min(n,k/(m+1));i++) ans+=k%i;
for(int i=1;i<=m;i++)
{
LL up=k/i,down=k/(i+1)+1;
if(down>n) continue;
up=min(up,n);
ans+=(up-down+1)*k-(up-down+1)*(up+down)*i/2;
}
printf("%lld\n",ans);
}
return 0;
}