bzoj1257[CQOI2007] 余数之和sum

最新推荐文章于 2020-07-09 17:11:18 发布

Euryale_

最新推荐文章于 2020-07-09 17:11:18 发布

阅读量481

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/u010336344/article/details/65631140

版权

数论专栏收录该内容

25 篇文章 0 订阅

订阅专栏

题目链接：bzoj1257
题目大意：
给出正整数n和k，计算 $J_n^k=\sum_{i=1}^{n}(k\mod i)$ 的值。

题解：
数论..
化式子变成

J k n = \sum i = 1 n (k - i \times ⌊ k i ⌋)

$J_n^k=\sum_{i=1}^{n}(k-i\times\lfloor \dfrac{k}{i}\rfloor)$
即

J k n = n \times k - \sum i = 1 n i \times ⌊ k i ⌋

$J_n^k=n\times k-\sum_{i=1}^{n}i\times\lfloor \dfrac{k}{i}\rfloor$
于是就把

⌊ki⌋ $\lfloor \dfrac{k}{i}\rfloor$ 相同的分成一块，对于同一块里的

i <script type="math/tex" id="MathJax-Element-149">i</script>一定是连续的。所以可以直接用求和公式进行计算。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;

LL mymin(LL x,LL y){return (x<y)?x:y;}
int main()
{
    //freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    LL n,k,i,l,r,mid,ret;
    scanf("%lld%lld",&n,&k);
    LL ans=n*k;
    LL lim=mymin(k,n);
    for (i=1;i<=lim;i=r+1) 
    {
        LL x=k/i;r=mymin(k/x,lim);
        ans-=x*(i+r)*(r-i+1)/2;
    }
    printf("%lld\n",ans);
    return 0;
}