BZOJ1257余数之和sum

1257: [CQOI2007]余数之和sum
Time Limit: 5 Sec Memory Limit: 162 MB
Submit: 2968 Solved: 1370
[Submit][Status][Discuss]
Description
给出正整数n和k，计算j(n, k)=k mod 1 + k mod 2 + k mod 3 + … + k mod n的值，其中k mod i表示k除以i的余数。例如j(5, 3)=3 mod 1 + 3 mod 2 + 3 mod 3 + 3 mod 4 + 3 mod 5=0+1+0+3+3=7
Input
输入仅一行，包含两个整数n, k。
Output
输出仅一行，即j(n, k)。
Sample Input
5 3
Sample Output
7
HINT
50%的数据满足：1<=n, k<=1000 100%的数据满足：1<=n ,k<=10^9
k%i=k-trunc(k/i)*i，当n>k时，直接把长度为(n-k)的区间长度乘以k加入答案，然后就处理i≤k的时候就可以了，然后就用j=n/(n/i)直接跳到（k/i都相等）区间末尾，然后需要对这个区间的等差数列求和加入答案，输出即可。。
注意点：输出%lld，加入答案的时候强制类型转换为long long。。
精髓点：①j=n/（n/i）
②已知区间始末和trunc(k/i)时，求k%i的等差数列的和的公式：
∑=(j-i+1)k-(j-i+1)(i+j)/2*(k/i);
附上本蒟蒻的代码：

#include<cstdio>
using namespace std;
int n,k,i,j;
long long ans=0;
int main()
{
    scanf("%d %d",&n,&k);
    if (n>k)
      {
        ans+=(long long)(n-k)*k;
        n=k;
      }
    j=0;
    for (i=1;i<=n;i=j+1)
      {
        j=k/(k/i);
        if (j>=n)
          j=n;
        ans+=(long long)(j-i+1)*k-(long long)(j-i+1)*(i+j)/2*(k/i);
      }
    printf("%lld",ans);
    return 0;
}