给出正整数n和k,计算G(n, k)=k mod 1 + k mod 2 + k mod 3 + … + k mod n的值,其中k mod i表示k除以i的余数。例如G(10, 5)=5 mod 1 + 5 mod 2 + 5 mod 3 + 5 mod 4 + 5 mod 5 …… + 5 mod 10=0+1+2+1+0+5+5+5+5+5=29
设G(n,k) =
∑
i
=
1
n
k
m
o
d
i
\sum_{ i = 1}^{n} {k mod i}
∑i=1nkmodi =
(那个符号实在打不出来。。)
(去外面写了一个截个图)
因为后面那个东西(k/i什么的)的值呈块状分布(很明显吧),就可以通过整除分块在O(sqrt(k))的时间复杂度内求出
for (long long l = 1 , r, s ; l <= n ; l = r + 1 )
{
r = (s = k / l) ? min( k / s, n ) : n ;
ans -= s * (r - l + 1) * (l + r) >> 1;
}
printf("%lld\n", ans);