算法学习笔记（23）：杜教筛

合集 - 算法(31)

1.算法学习笔记（1）：CDQ分治2023-11-172.算法学习笔记（2）：分块2023-11-303.算法学习笔记（3）：莫队算法2024-01-014.算法学习笔记（4）：FHQ平衡树（无旋平衡树）2024-01-015.算法学习笔记（5）：AC自动机2024-01-016.算法学习笔记（6）：优秀trick和性质合集2024-03-097.算法学习笔记（7）：数论2024-03-098.算法学习笔记（9）：第k大问题合集2024-03-169.算法学习笔记（10）：各种序的美好性质2024-03-1610.算法学习笔记（11）：历史版本和线段树2024-03-3111.算法学习笔记（12）：左偏树2024-04-0812.算法学习笔记（13）：同余最短路2024-04-1213.算法学习笔记（14）：区间最值操作和历史最值问题2024-04-2814.算法学习笔记（15）： Splay树2024-05-0615.算法学习笔记（16）：Link Cut Tree2024-05-0716.算法学习笔记（17）：Slope trick2024-05-1017.算法学习笔记（18）：珂朵莉树2024-05-1318.算法学习笔记（20）：网络流2024-05-1719.算法学习笔记（21）：数论分块2024-06-0620.算法学习笔记（22）：莫比乌斯反演2024-06-06

21.算法学习笔记（23）：杜教筛2024-06-07

22.算法学习笔记（24）：筛法2024-06-0823.暑假集训学习笔记（1）：lxl DS Day 12024-07-0224.暑假集训学习笔记（2）：lxl DS Day 22024-07-0225.暑假集训学习笔记（3）：lxl DS Day 32024-07-0426.多项式笔记2024-07-0827.生成函数笔记2024-07-0828.插头DP2024-07-2529.DP选讲做题记录 by 付乙淼2024-07-2630.拓展摩尔投票2024-09-1831.图论知识总结2024-09-25

收起

杜教筛

参考来源： OI-Wiki， 网上博客
线性筛可以在线性时间求积性函数前缀和， 而杜教筛可以用低于线性时间求解积性函数前缀和。

我们考虑 $S(n)$ 就是积性函数的前缀和， 所以我们尝试构造关于 $S(n)$ 关于 $S(\lfloor \frac{n}{i}\rfloor )$ 的递推式。

对于任意一个数论函数 g，必满足：

$\begin{array}{rl}\sum _{i=1}^n(f\ast g)(i)& =\sum _{i=1}^n\sum _{d\mid i}g(d)f\left(\frac{i}{d}\right)\\ & =\sum _{i=1}^{n}g(i)S\left(\lfloor \frac{n}{i}\rfloor \right)\end{array}$

将 $S(n)$ 提出来就可以得到：

$\begin{array}{rl}g(1)S(n)& =\sum _{i=1}^{n}g(i)S\left(\lfloor \frac{n}{i}\rfloor \right)-\sum _{i=2}^{n}g(i)S\left(\lfloor \frac{n}{i}\rfloor \right)\\ & =\sum _{i=1}^n(f\ast g)(i)-\sum _{i=2}^{n}g(i)S\left(\lfloor \frac{n}{i}\rfloor \right)\end{array}$

再将 $g(1)$ 移过去即可得到 $S(n)$ 的递推式， 考虑我们可以自己构造 $g(n)$ 使得计算变快。

当我们构造出这样的 $g(n)$ 时：

1. 可以快速计算 $\sum _{i=1}^n(f\ast g)(i)$；
2. 可以快速计算 g 的前缀和，以用数论分块求解 $\sum _{i=2}^{n}g(i)S\left(\lfloor {\displaystyle \frac{n}{i}}\rfloor \right)$。

则我们可以在较短时间内求得 $g(1)S(n)$。

例题

P4213 【模板】杜教筛

求 $\mu$ 和 $\varphi$ 的前缀和， 考虑一个性质， $\mu \ast 1=\epsilon$， $\varphi \ast 1=id$

考虑 $\sum _{i=1}^n(f\ast g)(i)$ 也就是 $\sum id$ 和 $\sum \epsilon$， 可以快速计算。

考虑 $g(i)$ 也就是 $1$， 我们可以快速计算， 所以可以直接数论分块了。

点击查看代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned int uint;
const int N = 2e6 + 10;
unordered_map<int, ll> sphi, smu; 
ll phi[N], mu[N];
int pri[N], tot, n, T;
bool p[N];
void init() {
	phi[1] = 1; mu[1] = 1;
	for(int i = 2; i <= 2e6; i++) {
		if(!p[i]) pri[++tot] = i, phi[i] = i - 1, mu[i] = -1;
		for(int j = 1; j <= tot && i * pri[j] <= 2e6; j++) {
			p[i * pri[j]] = 1;
			if(i % pri[j] == 0) {
				phi[i * pri[j]] = phi[i] * pri[j];
				mu[i * pri[j]] = 0;
				break;
			} 
			else {
				phi[i * pri[j]] = phi[i] * (pri[j] - 1);
				mu[i * pri[j]] = -mu[i];
			}
		}
	}
	for(int i = 2; i <= 2e6; i++) phi[i] += phi[i - 1], mu[i] += mu[i - 1];
}
ll calcphi(int n) {
	if(n <= 2e6) return phi[n];
	if(sphi[n]) return sphi[n];
	ll res = (ll)n * ((ll)n + 1) / 2;
	for(uint l = 2, r; l <= n; l = r + 1) {
		r = n / (n / l);
		res -= 1ll * (r - l + 1) * calcphi(n / l);
	}	
	return sphi[n] = res;
}
ll calcmu(int n) {
	if(n <= 2e6) return mu[n];
	if(smu[n]) return smu[n];
	ll res = 1;
	for(uint l = 2, r; l <= n; l = r + 1) {
		r = n / (n / l);
		res -= 1ll * (r - l + 1) * calcmu(n / l);
	} 
	return smu[n] = res;
}
int main() {
	scanf("%d", &T); init();
	while(T--) {
		scanf("%d", &n);
		printf("%lld %lld\n", calcphi(n), calcmu(n));
	}
	return 0;
}