P2261 [CQOI2007]余数求和

原题传送门

题目大意:给出正整数 n n n, k k k,请计算 G ( n , k ) = ∑ i = 1 n   k   m o d   i G(n,k)=\sum_{i=1}^n\ k\ mod\ i G(n,k)=i=1n k mod i
思路:一道整除分块的入门题


整除分块的基本形式为: ∑ i = 1 n ⌊ n i ⌋ \sum_{i=1}^n\lfloor \frac{n}{i}\rfloor i=1nin
对于任意一个 i ( i ≤ n ) i(i\le n) i(in),我们要找到一个最大的 j ( i ≤ j ≤ n ) j(i\le j\le n) j(ijn),使得 ⌊ n i ⌋ = ⌊ n j ⌋ \lfloor\frac{n}{i}\rfloor=\lfloor\frac{n}{j}\rfloor in=jn,那么此时 j = ⌊ n ⌊ n i ⌋ ⌋ j=\lfloor\frac{n}{\lfloor\frac{n}{i}\rfloor}\rfloor j=inn

显然 i ≤ j ≤ n i\le j\le n ijn(证明从略),设 k = ⌊ n i ⌋ k=\lfloor\frac{n}{i}\rfloor k=in,当 ⌊ n j ⌋ = k \lfloor\frac{n}{j}\rfloor=k jn=k j m a x = ⌊ n k ⌋ j_{max}=\lfloor\frac{n}{k}\rfloor jmax=kn(证明从略)

那么对于每次以 [ i , j ] [i,j] [i,j]为一块,分块求和即可

Code:

for(ll l=1,r;l<=n;l=r+1){
        r=n/(n/l);
        ans+=(r-l+1)*(n/l);
}

对于本题来说, ∑ i = 1 n   k   m o d   i = ∑ i = 1 n   k − i ∗ ⌊ k i ⌋ = n ∗ k − ∑ i = 1 n   i ∗ ⌊ k i ⌋ \sum_{i=1}^n\ k\ mod\ i=\sum_{i=1}^n\ k-i*\lfloor\frac{k}{i}\rfloor=n*k-\sum_{i=1}^n\ i*\lfloor\frac{k}{i}\rfloor i=1n k mod i=i=1n kiik=nki=1n iik

对于 ∑ i = 1 n   i ∗ ⌊ k i ⌋ \sum_{i=1}^n\ i*\lfloor\frac{k}{i}\rfloor i=1n iik,我们把每块 [ l , r ] [l,r] [l,r]对应的 ∑ i = l r   i ∗ ⌊ k i ⌋ \sum_{i=l}^r\ i*\lfloor\frac{k}{i}\rfloor i=lr iik中,令 T = ∑ i = l r ⌊ k i ⌋ T=\sum_{i=l}^r\lfloor\frac{k}{i}\rfloor T=i=lrik ∑ i = l r   i ∗ ⌊ k i ⌋ = T ∗ ∑ i = l r   i = T ∗ l + r 2 \sum_{i=l}^r\ i*\lfloor\frac{k}{i}\rfloor=T*\sum_{i=l}^r\ i=T*\frac{l+r}{2} i=lr iik=Ti=lr i=T2l+r

复杂度 O ( n ) O(\sqrt n) O(n )

Code:

#include <bits/stdc++.h>
using namespace std;
const int N=2e3+10,mod=1e9+7;
typedef long long ll;
ll n,k;
int main(){
	cin>>n>>k;
	ll ans=n*k;
	for(ll l=1,r;l<=n;l=r+1){
		if(k/l!=0) r=min(n,k/(k/l));
		else r=n;
		ans-=(l+r)*(k/l)*(r-l+1)/2;
	}
	cout<<ans<<endl;
	return 0;
}
 
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

学不会数据库

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值