P2261 [CQOI2007]余数求和 【整除分块】

一、题面

　　P2261 [CQOI2007]余数求和

二、分析

　　参考文章：click here

　　对于整除分块，最重要的是弄清楚怎样求的分得的每个块的范围。

　　假设$ n = 10 ，k = 5 $　　

$$   i : 1   2 3 4 5 6 7 8 9 10  \  lfloor frac{k}{i} floor :  5 2 1 1 1 0 0 0 0 0   $$ 

　　我们推导出假设$ L = i $，那么，对应的 $ lfloor frac{k}{i}  floor $ 相等的最右边界为 $  R =  lfloor frac{k}{  lfloor frac{k}{i}  floor }  floor $.（具体证明可以看参考文章。）

　　需要注意的细节是

　　1　$R$可能超过$n$,所以要限制一下。

　　2　一定要用$long long$。

三、AC代码

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;

int main()
{
    //freopen("input.txt", "r", stdin);
    ll n, k;
    while(scanf("%lld%lld", &n, &k) != EOF)
    {
        ll ans = n * k;
        ll L, R;
        for(L = 1; L <= n; L = R + 1)
        {
            ll res = k/L;
            if(res)
            {
                // 必须加min,因为k/res可能超过n，例如 k = 10, n = 6
                R = min(k/res, n);
            }
            else
                R = n;
            ans -= res * (R - L + 1) * (R + L) / 2;
        }
        printf("%lld
", ans);
    }
    return 0;
}