P2261 [CQOI2007]余数求和 （整除分块）

P2261 [CQOI2007]余数求和

$$\begin{array}{rl}G(n,k)& =\sum _{i=1}^{n}kmodi\\ & =\sum _{i=1}^n(k-i\ast [\frac{k}{i}])\\ & =nk-\sum _{i=1}^{n}i\ast [\frac{k}{i}]\end{array}$$

• 以 $k=6$ 为例：

$i$	1	2	3	4	5	6	7	8	9	10	11
$[\frac{k}{i}]$	$6$	$3$	$2$	$1$	$1$	$1$	$0$	$0$	$0$	$0$	$0$

• 发现 $[\frac{k}{i}]$ 被分为了 $4$ 块，可以用整除分块来计算
• 整除分块的时间复杂度 $O(2\sqrt{k})$ , （ $[\frac{k}{i}]$ 约有 $2\sqrt{k}$ 个值 ）

#include <bits/stdc++.h>
#define int long long 
using namespace std;

signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int n,k;
    cin>>n>>k;
    int ans=n*k;
    for(int l=1,r;l<=n;l=r+1)
    {
        int t=k/l; // 即 k/i 
        if(t) r=min(k/t,n); // 确定右端点
        else r=n;
        ans -= t*(r-l+1)*(l+r)>>1; // 表示区间[l,r]的贡献，别忘了除2
    }
    cout<<ans<<endl;
    return 0;
}