余数求和【除法分块】

最新推荐文章于 2021-07-16 16:02:05 发布

小鱼yn

最新推荐文章于 2021-07-16 16:02:05 发布

阅读量373

点赞数 1

分类专栏：小东西文章标签：分块

本文链接：https://blog.csdn.net/qq_43803508/article/details/107228024

版权

小东西专栏收录该内容

13 篇文章 0 订阅

订阅专栏

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

题意：
给出正整数 $n$ 和 $k$ ，请计算
$G(n,k)=\sum _{i=1}^n k\bmod i$
其中 $k\bmod i$ 表示 $k$ 除以 $i$ 的余数。
例：G(10,5)=0+1+2+1+0+5+5+5+5+5=29。

因为 $k\bmod i$ 可以表示成 $k-[\frac{k}{i}]*i$ ，所以
$G(n,k)=\sum _{i=1}^n (k-[\frac{k}{i}]* i)=n*k-\sum _{i=1}^n [\frac{k}{i}]* i$
$n = k = 25$ 时， $[\frac{k}{i}]$ 的值：
在这里插入图片描述

#include<algorithm>
#include<iostream>
//#include<stdio.h>
using namespace std;

int main()
{
    ios::sync_with_stdio(false);
    LL n,k,ans;
    cin>>n>>k;
    ans=n*k;
    for(LL l=1,r;l<=n;l=r+1)//l总是倍数相同时，最小的那一个数
    {
        if(k/l==0)break;//k<i时，k%i=k
        r=min(k/(k/l),n);//通过l求r
        ans-=(k/l)*(l+r)*(r-l+1)/2;
    }
    /*
    for(int l=1,r;l<=n;l=r+1)
	{
		r=n/(n/l);
		ans+=(r-l+1)*(n/l);
	}*/
    cout<<ans<<endl;
    return 0;
}

在这里插入图片描述
题解：
因为 $x\%(k*i)$ 可以表示成 $x\%k+([\frac{x}{k}]\%i)*k$ ，所以
$ans=\sum\limits _{i=l}^r x\%(k*i)=\sum\limits _{i=l}^r (x\%k+([\frac{x}{k}]\%i)*k)$
$ans=(x\%k)*(r-l+1)+k*\sum\limits _{i=l}^r [\frac{x}{k}]\%i$
令 $n=[\frac{x}{k}]$ 通过上一题可知：
$\sum\limits _{i=l}^r [\frac{x}{k}]\%i=\sum\limits _{i=l}^r n\%i=(r-l+1)*n-\sum\limits _{i=1}^n [\frac{n}{i}]* i$

#include<algorithm>
#include<iostream>
//#include<stdio.h>
#include<vector>
#include<string>
#include<cstring>
#include<math.h>
#include<iomanip>
using namespace std;
//#define int long long
typedef long long LL;
const int manx=4e5+10;
const int manx2=1e7+10;
const int mod=1e9+7;
const double inf=1e9+7;

int main()
{
    ios::sync_with_stdio(false);
    LL t,ll,rr,x,k;;
    cin>>t;
    while(t--)
    {
        cin>>ll>>rr>>x>>k;
        LL ans=(x%k)*(rr-ll+1)%mod;
        LL n=x/k;
        LL temp=n*(rr-ll+1);
        for(LL l=1,r;l<=rr;l=r+1)
        {
            if(n/l==0)break;
            r=min(n/(n/l),rr);
            LL l2=max(l,ll);
            if(l2<=r)
                (temp-=(n/l)*(((r-l2+1)*(r+l2)/2)%mod))%=mod;
        }
        temp=(temp+mod)%mod;
        (temp*=k)%=mod;
        cout<<(ans+temp)%mod<<endl;
    }
    return 0;
}