题目链接:https://vjudge.net/problem/UVA-1363
题意:计算
∑
i
=
1
n
\sum_{i=1}^n
∑i=1n(k%i)。
思路:紫书p338。假设k / i整数部分为p,那么k % i= k - i * p。如果k / (i+1)的整数部分也是p, 那么k % ( i + 1 ) = k - ( i + 1 ) * p = k - i * p-p=k%i-p。k % i ,k % ( i + 1 )是一个以-p为公差的等差序列。假设j为等差序列最后一项,因为k/j>=p,则j<=k / p,如此可计算项数。
#include<cstdio>
#include<iostream>
#include<algorithm>
#include <cstring>
#include <cmath>
#include <cstdio>
using namespace std;
int main()
{
long long n,k;
long long ans;
while(~scanf("%lld%lld",&n,&k))
{
ans=0;
for(int i=1;i<=n;i++)
{
long long p=k/i;//计算整数部分
if(p==0)//提前结束
{
ans=ans+(n-i+1)*k;
break;
}
long long x=k%i;//首项
long long y=min(n,k/p);//项数
long long z=x-p*(y-i);//末项
ans=ans+(x+z)*(y-i+1)/2;
i=y;
}
printf("%lld\n",ans);
}
return 0;
}