杜教筛入门

以下主要的话都用无序列表表示。

诶，是不是应该先讲背景

有什么好讲的？

问一个积性函数的前缀和，项数到1e10。

前置知识

线性筛，积性函数。

正文

钦定你已经可以再$O(\sqrt{n})$的复杂度内求出：

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

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

对于第一个，先枚举小于$\sqrt{n}$的i，得出这段的值；又因为$i$在一段区间内$\lfloor \frac{n}{i}\rfloor$都是$i_0\le \sqrt{n}$，因此可以求前式和后式。

• 形式化地讲，若${\sum }_{i=1}^{n}f(i)\ast g(\lfloor \frac{n}{i}\rfloor )$中，若能方便的求得$f(i)$的前缀和以及$g(x)$，就能方便的求得原式。

然后对于积性函数$f(i)$，我们想求$S_n={\sum }_{i=1}^{n}f(i)$。

那么，我们找一个积性函数$g$，令$T_n={\sum }_{d\mid n}g(d)\ast f(\frac{n}{d})$（就是狄利克雷卷积）。则有：
$$\begin{array}{rl}\sum _{i=1}^n{T}_i& =\sum _{i=1}^n\sum _{d\mid i}g(d)\ast f(\frac{i}{d})\\ & =\sum _{d=1}^{n}g(d)\ast \sum _{d\mid i,i\le n}f(\frac{i}{d})\\ & =\sum _{d=1}^{n}g(d)\ast S_{\lfloor \frac{n}{d}\rfloor }\end{array}$$

然后，钦定$g(1)=1$，那么就有
$$\begin{array}{rl}{S}_n& =\sum _{d=1}^{n}g(d)\ast S_{\lfloor \frac{n}{d}\rfloor }-\sum _{d=2}^{n}g(d)\ast S_{\lfloor \frac{n}{d}\rfloor }\\ & =\sum _{i=1}^n{T}_i-\sum _{d=2}^{n}g(d)\ast S_{\lfloor \frac{n}{d}\rfloor }\end{array}$$

如果不管怎么求$S_{\lfloor \frac{n}{d}\rfloor }$的话，发现满足以下两个条件就可以求$S_n$了。

• 可以求g的前缀和
• 可以求T的前缀和

考虑怎么求上式后面的$S_{\lfloor \frac{n}{d}\rfloor }$，由于f是积性函数，必定可以用线性筛筛出前n项，于是可以用线性筛筛了f求出$S_{1..\sqrt{n}}$，至于大于$\sqrt{n}$的下标，可以记搜：因为$\lfloor \frac{\lfloor \frac{n}{d}\rfloor }{e}\rfloor =\lfloor \frac{n}{d\ast e}\rfloor$，于是可以记录n/d的d，当$d>\sqrt{n}$时直接返回结果就行。

• 于是，求前缀和成功转化成用人类智慧求一个g。

模板题

对于$\varphi$，g为1（常值函数），T为i（自然数序列），用到的结论是$x={\sum }_{d\mid x}ph{i}_d$。

证明有这样几个方向：

1. 证明$\forall i\mid n,1\le j<i,(i,j)=1$，都有$i_1\ast j_1̸=i_2\ast j_2$。显然是要被叉翻的（其实是我在用3时不知道干嘛了）。
2. 直接通过$\varphi$的计算公式和积性通过一波推理得到一些式子，化简得到n。（巨佬做法）
3. 证明$\forall i\mid n,1\le j<i,(i,j)=1$，都有$n/i_1\ast j_1̸=n/i_2\ast j_2$。显然，因为互质，所以两个都是最简分数，于是当i,j不同时，分数不可能相等，于是证明了任意一个$x\in [1,n]$都只能由一个(i,j)对转移而来，即一一对应，证毕。

以上2的证明：
$$\begin{array}{rl}\sum _{d\mid n}\varphi (d)& =\sum \varphi (\prod _{i=1}^m{p}_i^j)\\ & =\sum \prod _{i=1}^m\varphi (p_i^j)\\ & =\prod _{i=1}^m\sum _{j=0}^{a_i}\varphi (p_i^j)\\ & =\prod _{i=1}^m(\sum _{j=1}^{a_i}(p_i^j-p_i^{j-1})+1)\\ & =\prod _{i=1}^m{p}_i^{a_i}=n\end{array}$$

然后写出来就可以了。

【代码代填坑】