【整除分块】洛谷P2261 余数求和题解

搂鱼114514

已于 2024-03-10 10:36:45 修改

阅读量243

点赞数 13

文章标签：算法

于 2024-03-09 13:27:19 首次发布

本文链接：https://blog.csdn.net/m0_60506105/article/details/136581097

版权

题意

给定n和k，求

$\sum_{i=1}^{n}k\: mod\; i$

思路

这一题求的是余数，观察余数是怎么来的：

$a\div b=c\cdot \cdot \cdot \cdot \cdot \cdot d$

那么

$d=a-bc$

则原始可以化为

$\sum_{i=1}^{n}k-i*\left \lfloor \frac{k}i{} \right \rfloor$

把k提出来：

$nk-\sum_{i=1}^{n}i*\left \lfloor \frac{k}{i} \right \rfloor$

则可以考虑对其中的 $\left \lfloor \frac{k}{i} \right \rfloor$ 进行整除分块，计算同一分块对答案的贡献：对于 $\left \lfloor \frac{k}{i} \right \rfloor$ 相同的 $i$ ，将其归到一块统一处理，从而达到减小时间复杂度的目的。

拿 $n=10\: k=5$ 举例如图为 $i$ 与 $\left \lfloor \frac{k}{i} \right \rfloor$ 的关系：

众所周知，在分块中，是通过双指针l和r实现分块的截取。

对左右指针l和r的计算：对于有序数组 $\left \{ 1,2,\cdot \cdot \cdot \cdot \cdot \cdot n \right \}$ ，让 $l$ 为分块的左边界，那么分块的单个值为 $p=\left \lfloor \frac{n}{l} \right \rfloor$ ，那么右边界就是值为 $k$ 的最大下标，即寻找满足 $i\leq \frac{n}{p}$ 的 $i_{max}$ ，即：

$r=i_{max}=\left \lfloor \frac{n}{p} \right \rfloor=\left \lfloor \frac{n}{\left \lfloor \frac{n}{l} \right \rfloor} \right \rfloor$

下一次截取令 $l'=r+1$ 即可。

得知l和r怎么计算后，就要将 $[l,r]$ 中的 $i*\left \lfloor \frac{k}{i} \right \rfloor$ 算出来。观察双指针运行中的原理如图：

对于 $i=[3,5]$ 时的 $i*\left \lfloor \frac{k}{i} \right \rfloor=(3+4+5)\times 1=12$

$3\: 4\: 5$ ？那不就是等差数列求和吗！

那么对于每个分块，其对答案的贡献为 $(\frac{n}{l})(l+r)\times (r-l+1)\div 2$ ！

还要注意，因为分块的单个值出现了等于0的情况，所以计算右端点时要特判，否则会报错！

具体细节见代码。

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const ll z=1e6+9;
void fastread()
{
	ios::sync_with_stdio(0);
}
ll n,k,ans; 
int main()
{
	fastread();
	cin>>n>>k;
	ans=n*k;
	for(ll l=1,r;l<=n;l=r+1)
	{
		if(k<l)r=n;//k/l已经全为0 
		else r=min(n,k/(k/l));
		ans-=(k/l)*(l+r)*(r-l+1)/2;
		//当前块*(首项+末项)*项数/2 
	}
	cout<<ans;
	return 0;
}

搂鱼114514

关注

13
点赞
踩
5

收藏

觉得还不错? 一键收藏
0
评论
【整除分块】洛谷P2261 余数求和题解

还要注意，因为分块的单个值出现了等于0的情况，所以计算右端点时要特判，否则会报错！，将其归到一块统一处理，从而达到减小时间复杂度的目的。众所周知，在分块中，是通过双指针l和r实现分块的截取。对左右指针l和r的计算：对于有序数组。为分块的左边界，那么分块的单个值为。那么对于每个分块，其对答案的贡献为。得知l和r怎么计算后，就要将。的最大下标，即寻找满足。，那么右边界就是值为。
复制链接

扫一扫