数论--整除分块

1.概念

在介绍整除分块之前，我们先来看一道算数题：已知正整数n，求
$$\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor$$

我们写一个表格看一看1-20的整除是什么样子的

表中同样的值会连续出现，而相同的值所划分的区间积是整出分块。整除的性质使得从1到n的数组表可根据数值划分为不同的分块，且分块数远远小于n。利用这种性质，我们如果能推导出每个分块具体的左右端点位置在哪，这个问题就可以快速求解出来了。

2.整除分块公式推导

搬出例题：$$\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor$$

假设我们已知某一个分块的左端点$l$，要求解出该分块的右端点$r$。设该分块的数值为$k$，对于该分块中的每个数$i$，有$k=\lfloor \frac{n}{l}\rfloor =\lfloor \frac{n}{i}\rfloor$即$i\ast k\le n$，也就是说我们找到可得使$i\ast k\le n$成立的最大的$i$的值即是我们所求的右端点$r$，因此我们可以得到下列式子:
$$k=\lfloor \frac{n}{l}\rfloor$$

$$r=max(i),i\ast k\le n$$

推导可得：
$$r=\lfloor \frac{n}{k}\rfloor =\lfloor \frac{n}{\lfloor \frac{n}{l}\rfloor }\rfloor$$

转换成代码就是：

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

还没有结束，我们再看这样一题：
$$p=\sum _{i=1}^n\lfloor \frac{n}{ax+b}\rfloor$$
第一种方法：

第二种方法：

例题1

$$\sum _{i=1}^n\lfloor \frac{n}{i^2}\rfloor$$

例题2

$$\sum _{i=1}^n\lceil \frac{n}{i}\rceil$$

我们就不推导向上取整了，这里只需要一个小转化，将向上取整转化为向下取整。
我们考虑没有整除的时候是不是就有
$$\lceil \frac{n}{i}\rceil =\lfloor \frac{n}{i}\rfloor +1$$如果整除的时候就相等了
那么我们只要不加1，我们加上$\frac{i-1}{i}$就可以了

最后放一道模板题：

余数求和

$PS$：需要注意$n/k$等于零时的情况
ACcode：

#include<bits/stdc++.h>
#define LL long long
using namespace std;
LL n,t,k,ans;

void solve()
{
	cin>>n>>k;
	ans = n*k;
	for(int l=1,r;l<=n;l=r+1)
	{
		if( k/l==0 ) r = n;
		else r = min( k/(k/l)  , n );
		ans -= (k/l)*(l+r)*(r-l+1)/2;
	}
	cout<<ans;
}

int main()
{
	solve();
	return 0;
}

参考大佬的博客