数论 —— 整除分块，常见经典例题。

什么是整除分块？

在我们求 $\lfloor \frac{n}{i}\rfloor$ 时会发现，当 $i$ 不同时 $\lfloor \frac{n}{i}\rfloor$ 却可能相同。

i	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15
$\frac{15}{i}$	15	7	4	3	3	2	2	1	1	1	1	1	1	1	1

如果我们要求 ${\sum }_{i=1}^{15}\frac{15}{i}$，如果直接枚举 i，那么时间复杂度是 O(n) 的。

而如上图所示，存在相邻的若干项的值是相同的，我们是否可以每次把相同值的区间算出来，用 区间长度 * 区间中的值 便是这段区间中所有值之和了。
这样就不需要遍历区间中重复值的位置，而 $\lfloor \frac{n}{i}\rfloor$ 中不同值有 $2\sqrt{n}$ 种，故复杂度降为 $O(\sqrt{n})$。

证明：
对于 $i\le \lfloor \sqrt{n}\rfloor ，\lfloor \frac{n}{i}\rfloor$ 最多有 $\lfloor \sqrt{n}\rfloor$ 种取值；
对于 $i>\lfloor \sqrt{n}\rfloor$ ，有 $\lfloor \frac{n}{i}\rfloor \le \lfloor \sqrt{n}\rfloor$ ，最多也只有 $\lfloor \sqrt{n}\rfloor$ 种取值；
所以 $\lfloor \frac{n}{i}\rfloor$ 最多只有 $2\sqrt{n}$ 种取值。
（准确来说是 $\sqrt{4n+1}-1$ 种取值）

那么我们就可以从前往后枚举，每到一个位置，找到其所在相同值区间的最右端位置，两位置相减+1便是区间长度，然后跳过这个区间。

那么如何得到相同值区间最右端位置呢？
设一段区间 $[l,r]$ 中所有元素的值都为 $k=\lfloor \frac{n}{l}\rfloor$。
现在要找到最大的 $i$ 使得 $k=\lfloor \frac{n}{i}\rfloor$ 满足。
也就是，找到满足 $i\le \frac{n}{k}$ 的最大的 $i$，那么这个 $i$ 为 $\lfloor \frac{n}{k}\rfloor$，即 $r=\lfloor \frac{n}{k}\rfloor$。
将 $k=\lfloor \frac{n}{l}\rfloor$ 代入上式即得：$r=\lfloor \frac{n}{\lfloor \frac{n}{l}\rfloor }\rfloor$。

所以我们只要在位置 $l$，就能得到和 $\lfloor \frac{n}{l}\rfloor$ 值相同的区间最右端位置：$$r=\lfloor \frac{n}{\lfloor \frac{n}{l}\rfloor }\rfloor$$

那么对于求 $$\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor$$ 就可以这样实现，复杂度$O(\sqrt{n})$：

int ans = 0;
for(int l = 1, r; l <= n; l = r + 1) //l是枚举的i，r是相同值的最大i
{
	r = n/(n/l);
	ans += (r-l+1) * (n/l); //注意n/l的括号一定要加上，先整除再相乘
}
cout << ans;

---

下面来看几道例题：

例1：约数研究

Link
定义 $f(N)$ 为 $N$ 的约数的个数。
求 ${\sum }_{i=1}^{n}f(i)$。

$1\le n\le 2\times 10^{14}$

思路
如果挨个求每个数的约数的话，复杂度是 $O(n\sqrt{n})$，会超时。
转化思路，考虑 1~n 每个约数分别出现了多少次？
对于一个数 x 来说，肯定作为其所有倍数的约数，那 1~n 中 x 的倍数有 $\lfloor \frac{n}{x}\rfloor$ 个，那这个约数就出现 $\lfloor \frac{n}{x}\rfloor$ 次。遍历所有数累加其出现的次数。
所以题目就转化为求：
$$\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor$$
就是模板题了。

Code

#include<bits/stdc++.h>
using namespace std;

#define Ios ios::sync_with_stdio(false),cin.tie(0)
#define int long long

const int N = 200010, mod = 1e9+7;
int T, n, m;
int a[N];

signed main(){
	Ios;
	cin >> n;
	
	int ans = 0;
	for(int l = 1, r; l <= n; l = r + 1)
	{
		r = n/(n/l);
		ans += (r-l+1) * (n/l);
	}
	cout << ans;
	
	return 0;
}

---

例2：P2424 约数和

Link
定义 $f(N)$ 为 $N$ 的约数之和。
给定 $X,Y$，求 $f(X)+f(X+1)+\cdots +f(Y)$。

$1\le n\le 2\times 10^{14}$

思路
我们可以仿照上一题的思路求 ${\sum }_{i=1}^{n}f(i)$，然后用 ${\sum }_{i=1}^{Y}f(i)$ 减去 ${\sum }_{i=1}^{X-1}f(i)$ 即可。

和上一题的思路一样，考虑每个约数的贡献。
对于每个约数 x，其出现了 $\lfloor \frac{n}{x}\rfloor$ 次，那么这个约数之和就是 $\lfloor \frac{n}{x}\rfloor \ast x$。
所以这题就转化为求：
$$\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor \ast i$$
对于一段区间 [l, r] 来说，其 $\lfloor \frac{n}{i}\rfloor$ 值相同，但是 i 值不同。
但由于 i 在同一区间，可以利用等差数列算出这一区间的 i 之和 sum，用 sum * $\lfloor \frac{n}{i}\rfloor$ 便是这段区间的贡献 ${\sum }_{i=l}^r\lfloor \frac{n}{i}\rfloor \ast i$。

Code

#include<bits/stdc++.h>
using namespace std;

#define Ios ios::sync_with_stdio(false),cin.tie(0)
#define int long long

const int N = 200010, mod = 1e9+7;
int T, n, m;
int a[N];

int sum_1(int n)
{
	return n * (n+1) / 2;
}

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

signed main(){
	Ios;
	cin >> n >> m;
	
	cout << ans(m) - ans(n-1);
	
	return 0;
}

---

例3：几番烟雾，只有花难护

Link
题目转化为，求：
$$\sum _{i=1}^n\lceil \frac{n}{i}\rceil \ast i^2$$
答案对 $998244353$ 取模。

思路
上取整转化为下取整：
$$\lceil \frac{n}{i}\rceil =\lfloor \frac{n-1}{i}\rfloor +1$$
注意这个公式只能在保证 n 不为 0 的情况下使用，n = 0 时出现错误为 1。
当 n 可能为 0 时可以用公式 $\lceil \frac{n}{i}\rceil =\lfloor \frac{n+i-1}{i}\rfloor$。

于是原式转化为：
$$\sum _{i=1}^n\lfloor \frac{n-1}{i}\rfloor \ast i^2+\sum _{i=1}^n{i}^2$$
后半部分可以直接求出，前半部分和上一题类似，这题需要处理出区间的 $i^2$ 之和，乘以区间相同值 $\lfloor \frac{n}{i}\rfloor$ 即可。

需要知道：
$$1^2+2^2+3^2+\dots +{\mathrm{n}}^2=\frac{\mathrm{n}(\mathrm{n}+1)(2\mathrm{n}+1)}{6}$$

• 除法要用逆元。

Code

#include<bits/stdc++.h>
using namespace std;

#define Ios ios::sync_with_stdio(false),cin.tie(0)
#define int long long

const int N = 200010, mod = 998244353;
int T, n, m;
int a[N];

int qmi(int x, int y)
{
	int ans = 1;
	while(y)
	{
		if(y & 1) ans = ans * x % mod;
		x = x * x % mod;
		y >>= 1;
	}
	return ans;
}

int sum_2(int n)
{
	return n * (n + 1) % mod * (2*n + 1) % mod * qmi(6, mod-2) % mod;
}

signed main(){
	Ios;
	
	cin >> T;
	while(T--)
	{
		cin >> n;
		
		int ans = sum_2(n);
		
		int x = n-1;
		for(int l=1, r; l <= x; l = r+1)
		{
			r = x/(x/l);
			int t = (sum_2(r) - sum_2(l-1) + mod) % mod;
			ans = (ans + t * (x/l) % mod) % mod;
		}
		cout << ans << endl;
	}
	
	return 0;
}

此外还有
$$1^3+2^3+\dots +{\mathrm{n}}^3={\left[\frac{\mathrm{n}(\mathrm{n}+1)}{2}\right]}^2=(1+2+\dots +\mathrm{n}{)}^2$$

---

例4：P2261 [CQOI2007] 余数求和

Link
给出正整数 $n$ 和 $k$ ，请计算 $$G(n,k)=\sum _{i=1}^{n}k\mathrm{mod}i$$

思路
$k\mathrm{mod}i$ 不好处理，将其转化，等价于 $k-i\ast \lfloor \frac{k}{i}\rfloor$。
原式转化为：
$$\sum _{i=1}^{n}k-i\ast \lfloor \frac{k}{i}\rfloor$$
即：
$$\sum _{i=1}^{n}k-\sum _{i=1}^{n}i\ast \lfloor \frac{k}{i}\rfloor$$
前半部分可以直接算出为 n*k，后半部分和例2一样，利用前缀和相减 * 区间值。

Code

#include<bits/stdc++.h>
using namespace std;

#define Ios ios::sync_with_stdio(false),cin.tie(0)
#define int long long

const int N = 200010, mod = 1e9+7;
int T, n, m;
int a[N];
int k;

int sum(int l, int r)
{
	return (l+r)*(r-l+1)/2;
}

signed main(){
	cin >> n >> k;
	
	int ans = n * k;
	for(int l = 1, r; l <= min(n, k); l = r + 1) //注意这里只需要循环到min(n,k)。否则当l>k时，k/l=0, k/0无意义
	{
		r = min(n, k/(k/l)); //这里的右端点要和n取min
		ans -= sum(l, r) * (k/l);
	}
	cout << ans;
	
	return 0;
}

当枚举的除数 i 的范围和原数 x 不相同的话，要注意两个地方：

• 枚举的除数 i 最大应该到 min(n, x)。因为当 $i>x$ 时，x/i = 0，x/(x/i) 便会无意义。
• r 的范围最大不超过 n，要和 n 取 min。当 n 和 x 相同时，n/(n/l) 一定不会超过 n，没问题。但是当 n 和 x 不同时，x/(x/l) 却可能大于 n，但是除数只需要枚举到 n，就加多了。

---

例5：求和公式

Link
给定 $n,m$ 求：
$$\sum _{i=1}^n\sum _{j=1}^m{i}^{2}j(n\%i)(m\%j)$$

$10\le n,m\le 10^9$，答案取模 $1e9+7$。

思路
首先观察出 i 和 j 没有联系，将其分开：
$$\sum _{i=1}^n{i}^2(n\%i)\ast \sum _{j=1}^{m}j(m\%j)$$
将取模操作转化：
$$\sum _{i=1}^n{i}^2\left(n-i\ast \lfloor \frac{n}{i}\rfloor \right)\ast \sum _{j=1}^{m}j\left(m-j\ast \lfloor \frac{m}{j}\rfloor \right)$$
即：
$$\left[n\ast \sum _{i=1}^n{i}^2-\sum _{i=1}^n{i}^3\ast \lfloor \frac{n}{i}\rfloor \right]\ast \left[m\ast \sum _{j=1}^{m}j-\sum _{j=1}^m{j}^2\ast \lfloor \frac{m}{j}\rfloor \right]$$
而对于左右两个部分的前半部分都可以直接求出，后半部分和上例一样，利用前缀和相减 * 区间相同值。

Code

//https://ac.nowcoder.com/acm/contest/2891/C
#include<bits/stdc++.h>
using namespace std;

#define Ios ios::sync_with_stdio(false),cin.tie(0)
#define int long long

/*
	尽量保证运算过程中的每一个小块在参与下一步运算时都是正数。
	可以把每个小块单独拉出来写。 
*/

const int N = 200010, mod = 1e9+7;
int T, n, m;
int a[N];

int qmi(int x, int y)
{
	int ans = 1;
	while(y)
	{
		if(y & 1) ans = ans * x % mod;
		x = x * x % mod;
		y >>= 1;
	}
	return ans;
}

int sum_1(int n)
{
	return n * (n+1) % mod * qmi(2, mod - 2) % mod;
}

int sum_2(int n)
{
	return n * (n+1) % mod * (2*n + 1) % mod * qmi(6, mod - 2) % mod;
}

int sum_3(int n)
{
	int x = n * (n+1) % mod * qmi(2, mod - 2) % mod;
	return x * x % mod;
}

signed main(){
	Ios;
	cin >> n >> m;
	
	int ans1 = n * sum_2(n) % mod;
	for(int l = 1, r; l <= n; l = r + 1)
	{
		r = min(n, n/(n/l));
		int t = (sum_3(r) - sum_3(l-1) + mod) % mod;
		ans1 = (ans1 - t * (n/l) % mod + mod) % mod; //两数相乘之后一定先取模再进行其他操作 
	}
	
	int ans2 = m * sum_1(m) % mod;
	for(int l = 1, r; l <= m; l = r + 1)
	{
		r = min(m, m/(m/l));
		int t = (sum_2(r) - sum_2(l-1) + mod) % mod;
		ans2 = (ans2 - t * (m/l) % mod + mod) % mod;
	}
	
	cout << ans1 * ans2 % mod;
	
	return 0;
}

---

多维整除分块（N 维数论分块）

求形如下式的和式的值： $$\sum _{i=1}^n\lfloor \frac{a_1}{i}\rfloor \lfloor \frac{a_2}{i}\rfloor \cdots \lfloor \frac{a_m}{i}\rfloor$$

每段区间的右端点 $r={\min }_{i=1}^m\left\{[\frac{a_i}{\lfloor \frac{a_i}{l}\rfloor }\rfloor \right\}$ 即可，即取每一区间的右端点中最小的那个作为整体的右端点。

较常用的是 2 维数论分块：求 $$\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{i}\rfloor$$ 区间右端点 $r=\min \left(\lfloor \frac{n}{\lfloor \frac{n}{l}\rfloor }\rfloor ,\lfloor \frac{m}{\lfloor \frac{m}{l}\rfloor }\rfloor \right)$，即 r = min(n, n/(n/l), m/(m/l))。

时间复杂度仍为 $O(\sqrt{n})$。

---

例题：P2260 [清华集训2012] 模积和

Link
题意
求：
$$\sum _{i=1}^n\sum _{j=1}^m(n\mathrm{mod}i)\times (m\mathrm{mod}j),i\ne j$$

$\mathrm{mod}19940417$ 的值。

思路
乍一看和上一题的做法一样，将 i 和 j 分开分别处理就行了，但是要注意到一个条件：$i\ne j$。
可以用一个容斥，先不考虑相等的情况，分开 i 和 j 分别算出答案相乘之后再减去 i 和 j 相同情况的答案：
$$\sum _{i=1}^n(n\mathrm{mod}i)\ast \sum _{j=1}^m\left(m\mathrm{mod}j\right)-\sum _{i=1}^{min(n,m)}\left(n\mathrm{mod}i\right)\left(m\mathrm{mod}i\right)$$

前两项直接求解，最后一项化简：
$$\begin{gathered} \sum _{i=1}^n(n\mathrm{mod}i)(m\mathrm{mod}i) \\ =\sum _{i=1}^n\left(n-\lfloor \frac{n}{i}\rfloor i\right)\left(m-\lfloor \frac{m}{i}\rfloor i\right) \\ =\sum _{i=1}^n\left(nm-\lfloor \frac{n}{i}\rfloor i\cdot m-\lfloor \frac{m}{i}\rfloor i\cdot n+\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{i}\rfloor i^2\right) \\ =n^{2}m-m\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor i-n\sum _{i=1}^n\lfloor \frac{m}{i}\rfloor i+\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{i}\rfloor i^2 \end{gathered}$$

第 2,3 项都是模板，第 4 项就是一个 2 维数论分块。

Code

#include<bits/stdc++.h>
using namespace std;

#define Ios ios::sync_with_stdio(false),cin.tie(0)
#define int long long

const int N = 200010, mod = 19940417;
int T, n, m;
int a[N];
int inv2 = 9970209, inv6 = 3323403;

int x, y;
void exgcd(int a, int b)
{
	if (!b)
	{
		x = 1, y = 0;
		return;
	}
	exgcd(b, a % b);
	int tmp = x;
	x = y;
	y = tmp - a / b * y;
}

int inv(int a)
{
	exgcd(a, mod);
	x = (x % mod + mod) % mod;
	return x;
}

int sum_1(int n){
	return n * (n+1) % mod * inv2 % mod;
}

int sum_2(int n){
	return n * (n+1) % mod * (2*n+1) % mod * inv6 % mod;
}

int ans_1(int n, int x)
{
	int ans = 0;
	for(int l=1, r; l <= min(n, x); l = r+1)
	{
		r = min(n, x/(x/l));
		int t = (sum_1(r) - sum_1(l-1) + mod) % mod;
		ans = (ans + t * (x/l) % mod + mod) % mod;
	}
	return ans;
}

int ans_2(int n, int x, int y)
{
	int ans = 0;
	for(int l = 1, r; l <= n; l = r+1)
	{
		r = min(n, min(x/(x/l), y/(y/l)));
		int t = (sum_2(r) - sum_2(l-1) + mod) % mod;
		ans = (ans + t * (x/l) % mod * (y/l) % mod) % mod;
	}
	return ans;
}

signed main(){
	Ios;
	
	cin >> n >> m;
	
	int res1 = (n * n % mod - ans_1(n, n) + mod) % mod;
	int res2 = (m * m % mod - ans_1(m, m) + mod) % mod;
	int mina = min(n, m);
	int res3 = (n * m % mod * mina % mod - n * ans_1(mina, m) % mod - m * ans_1(mina, n) % mod
		+ ans_2(mina, n, m) + mod) % mod;
		
	cout << (res1 * res2 % mod - res3 + mod) % mod;
	
	return 0;
}

注意在该题中，模数 19940417 不是素数，则不能用快速幂求逆元，需要用扩展欧几里得。
逆元只求了 2 和 6 两个数，为了降低复杂度，可以预处理出来。

---

总结：

总之，整除分块就是一个可以以 $O(\sqrt{n})$ 的复杂度求解形如 ${\sum }_{i=1}^n\lfloor \frac{a_1}{i}\rfloor \lfloor \frac{a_2}{i}\rfloor \cdots \lfloor \frac{a_m}{i}\rfloor$ 的问题的算法。

哪里有问题的话可以留言评论~

---

参考文章：
ACM竞赛中常见的整数公式与整数分块
OI Wiki - 数论分块
整数分块 —— 完美的根号n算法
【数论】整除分块（数论分块）
C -求和公式(整除分块)
【数学】数论分块（整除分块）
1³⁺²3+3^3+…+n3通项公式