Joseph's Problem UVA - 1363 取模求和

题目链接：https://vjudge.net/problem/UVA-1363
题意：计算${\sum }_{i=1}^n$（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;
}