进阶数论第三弹 杜教筛

1.杜教筛

杜教筛是用来在低于线性的时间复杂度$(O(n^{\frac{2}{3}})?)$内求出积性函数的前缀和的算法

根据杜教筛的定义，我们设

$$S(n)=\sum _{i=1}^{n}f(i)$$
$$g是一个积性函数$$
$$h(n)=f\times g$$
$$H(n)=\sum _{i=1}^{n}h(i)$$
那么有
$$\begin{array}{rl}H(n)& =\sum _{i=1}^{n}h(i)\\ & =\sum _{i=1}^n(f\times g)(i)\\ & =\sum _{i=1}^n\sum _{d|i}g(d)f(\frac{i}{d})\\ & =\sum _{i=1}^n\sum _{d|i}g(d)f(\frac{i}{d})\\ & =\sum _{d=1}^{n}g(d)\sum _{i=1}^{n}f(\frac{i}{d})[d|i]\\ & =\sum _{d=1}^{n}g(d)\sum _{i=1}^{\lfloor \frac{i}{d}\rfloor }f(i)\\ & =\sum _{d=1}^{n}g(d)S(\lfloor \frac{i}{d}\rfloor )\end{array}$$
变化一下可得：
$$\begin{array}{rl}g(1)S(n)& =H(n)-\sum _{d=2}^{n}g(d)S(\lfloor \frac{i}{d}\rfloor )\end{array}$$
这就是杜教筛最根本的式子。
如果我们能找到一个g，使得g的前缀和$g\times h$的前缀和都很好求，那我们就能用整数分块快速求出S了

应用

$$\sum _{i=1}^n\mu (i)$$

我们知道
$$\mu \times I=ϵ$$

$ϵ$的前缀和永远是1，$I$的前缀和就是n

那我们就可以推出$\mu$的前缀和式子了