码迷,mamicode.com
首页 > 其他好文 > 详细

Luogu2261 [CQOI2007]余数求和

时间:2020-07-20 22:36:59      阅读:69      评论:0      收藏:0      [点我收藏+]

标签:pac   stream   计算   数论   bre   long   cpp   ret   iostream   

数论分块

将原式变形,很容易得到以下式子:

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

对于\(\lfloor\frac{k}{i}\rfloor=x\),可得出\(\max i=\lfloor \frac{k}{x} \rfloor\)

计算每一块的\(l,r\),高斯求和即可

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
int n,k,r;
long long ans;
int main()
{
    scanf("%d%d",&n,&k);
    ans=(long long)n*(long long)k;
    for (int l=1;l<=n;l=r+1)
    {
        if (k/l==0)
            break;
        r=k/(k/l);
        r=min(r,n);
        ans-=(long long)(k/l)*(long long)(l+r)*(long long)(r-l+1)/2LL;
    }
    cout << ans << endl;
    return 0;
}

Luogu2261 [CQOI2007]余数求和

标签:pac   stream   计算   数论   bre   long   cpp   ret   iostream   

原文地址:https://www.cnblogs.com/GK0328/p/13347191.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!