洛谷 P2261 [CQOI2007]余数求和 解题报告

P2261 [CQOI2007]余数求和

题意:

\(G(n,k)=\sum_{i=1}^n k \ mod \ i\)

数据范围:

\(1 \le n,k \le 10^9\)


\(G(n,k)\)
\(=\sum_{i=1}^n k-i*\lfloor \frac{k}{i} \rfloor\)
\(=n*k-\sum_{i=1}^n i*\lfloor \frac{k}{i} \rfloor\)

显然,\(\lfloor \frac{k}{i} \rfloor\)的分布可能会有重复。

根据除法分块(别在意它只是一个名字),这些值不重复的个数大约是\(\sqrt k\)

我们只需要统计每一块的值即可,注意到值在区间上的出现是单调递增的。

如果我们得到某一块最开始的下标\(l\)(可以从上一块的\(r\)得到),如何推得这一块的\(r\)呢?

其实很简单,\(\frac{k}{l}\)的余数是最大的,而\(\frac{k}{r}\)的余数显然得为0

\(t=\frac{k}{l}\),则\(r=\frac{k}{t}\)

加上不能越界的判断,完整的即为
\(if \ t==0\)
\(r=n\)
\(else\)
\(r=min(n,\frac{k}{t})\)

对每一块直接统计即可


Code:

#include <cstdio>
#define ll long long
ll min(ll x,ll y){return x<y?x:y;}
ll k,ans,n,l,r,t;
int main()
{
    scanf("%lld%lld",&n,&k);
    ans=n*k;
    l=1;
    while(r!=n)
    {
        ll t=k/l;
        if(!t) r=n;
        else r=min(n,k/t);
        ans-=(r+1-l)*(l+r)*t/2;
        l=r+1;
    }
    printf("%lld\n",ans);
    return 0;
}


2018.7.23

posted @ 2018-07-23 21:35  露迭月  阅读(164)  评论(0编辑  收藏  举报