HDU2620：数论+分块

HDU2620

题解

Σ1<=i<=n(k mod i)
k mod i = k - (k/i)*i
Σ1<=i<=n(k mod i)
 = Σ1<=i<=n[k - (k/i)*i]
 = nk - Σ1<=i<=n (k/i)*i
 存在一段区间的k/i相同，只要求出数量。

设k / i = k / j = d，i为区间的左端，j为区间的最右端。

如果已知区间的最多端i，d = k / i，可以证明j = k / d。

代码

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll n,k;
int main(){
	while(~scanf("%lld%lld",&n,&k)){
		ll ans = n * k;
		if(n > k)	n = k;  //n > k，大于n的部分模值是k本身。
		for(ll i=1;i<=n;){
			ll d = k / i;    //如果i > k,d = 0。下面就不合法了
			ll j = k / d;
			if(j > n)	j = n;
			ans -= (i + j) * (j - i + 1) / 2 * d;
			i = j + 1;
		}
		printf("%lld\n",ans);
	}
	return 0;
}