uva 1363 公式求和(优化)

点击打开链接

//计算 累加k%i (i:1~n)
//n,k<=1e9直接算超时,k固定除数增加,下取整,商相同的有很多 
//若k/i=p &&k/(i+1)=p
//则k%(i+1)=k-(i+1)*p=k-ip-p=k%i-p
//即商的整数部分相同的i~j的余数呈等差数列 

#include <iostream>
#include <cstring>
using namespace std;
const int N=1e3+20;
typedef long long ll;
const ll mod=1e9+7;
int main()
{
	ll n,k; 
	//计算 累加k%i (i:1~n)
	//n,k<=1e9直接算超时,k固定除数增加,下取整,商相同的有很多 
	//若k/i=p &&k/(i+1)=p
	//则k%(i+1)=k-(i+1)*p=k-ip-p=k%i-p
	//即商的整数部分相同的i~j的余数呈等差数列 
	while(cin>>n>>k)
	{
		int i=1;
		ll ans=0;
		while(i<=n)
		{
			ll p=k/i;
			ll j;
			if(p==0)//k<i k%i=k
			{
				ans+=k*(n-i+1);
				break;
			}
			//i从小到大,k/i为p的最大i是多少？ k/i=p.0 i=k/p   
			j=min(k/p,n);
			ll x=k%i;
			ll y=k%j;
		//	ans+=(x+y)*(j-i+1)/2;
			ans+=(x+x-p*(j-i))*(j-i+1)/2;
			i=j+1;			
		} 
		printf("%lld\n",ans); 
	}
	return 0;
}