杜教筛详解

前置技能　狄利克雷卷积　线性筛函数

杜教筛

设要求积性函数 $f$ 的前缀和， 设 $S(n)=\sum _{i=1}^{n}f(i)$。
设另有一函数 $g$，则它们的狄利克雷卷积的前缀和为
$$\begin{array}{rl}\sum _{i=1}^n(f\ast g)& =\sum _{i=1}^n\sum _{d|i}f(d)g(\frac{i}{d})\\ & =\sum _{d=1}^n\sum _{i|d}g(d)f(\frac{i}{d})\\ & =\sum _{d=1}^{n}g(d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }f(i)\\ & =\sum _{d=1}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )\end{array}$$

接下来考虑 $g(1)S(n)$
$$\begin{array}{rl}g(1)S(n)& =\sum _{i=1}^{1}g(i)S(\lfloor \frac{n}{i}\rfloor )\\ & =\sum _{i=1}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )-\sum _{i=2}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )\\ & =\sum _{i=1}^n(f\ast g)(i)-\sum _{i=2}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )\end{array}$$

于是，我们得到了需要的核心结论
$$g(1)S(n)=\sum _{i=1}^n(f\ast g)(i)-\sum _{i=2}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )$$

这让我们开始考虑是否有某些十分容易得到前缀和的积性函数，如果有这样的积性函数的话，那右边那个式子可以使用除法分块解决。令求出 $S(n)$ 的时间复杂度为 $T(n)$，若左边的式子可以 $\Theta (1)$ 求解，考虑到除法分块的时间复杂度为$\Theta (\sqrt{n})$，有：
$$T(n)=\sum _{i=1}^{\sqrt{n}}\Theta (\sqrt{i})+\Theta \left(\sqrt{\frac{n}{i}}\right)=\Theta \left(n^{\frac{3}{4}}\right)$$

如果我们先预处理出前 $m$ 个前缀和，则时间复杂度为
$$T(n)=\sum _{i=1}^{\lfloor \frac{n}{m}\rfloor }\Theta \left(\sqrt{\frac{n}{i}}\right)=\Theta \left(\frac{n}{\sqrt{m}}\right)$$

当 $m=n^{\frac{2}{3}}$ 时，$T(n)=\Theta (n^{\frac{2}{3}})$ ，时间复杂度达到最低。
可惜上述时间复杂度都必须建立在一个基础条件下：记忆化 ，如果不这样做，时间复杂度无法保证，可是程序不允许开如此大的数组，如果使用 $map$，则程序的时间复杂度会变为多一个 $lo{g}_{2}n$。
其实并不需要使用 $map$，考虑到每次递归求解时，我们需要的并不是 $S(n)$ 的值，而是 $S(\lfloor \frac{n}{i}\rfloor )$ 的值，因此可以仅仅用 $S[i]$ 保存 $S(\lfloor \frac{n}{i}\rfloor )$ 。
代码如下：

long long S(int n) {
	if(n<maxm) return mSphi[n];//线性筛
	int x=maxn/n;//求出S(n)应该放在的位置
	if(vis[x]) return Sphi[x];//记忆化
	vis[x]=true;
	long long &ans=Sphi[x]=calc_left();//左边的答案
	for (int i=2,j;i<=n;i=j+1) {
		j=n/(n/i);
		ans-=(j-i+1)*S(n/i)*g(i,j);//g(i,j)表示 $\sum g(i)$
	} return ans;
}

应用

$\mu$ 的前缀和

求 $\begin{array}{r}\sum _{i=1}^n\mu (i)\end{array}$，$1⩽n⩽2^{31}$
首先，线筛肯定炸。
考虑到 $\mu \ast 1=ϵ$，令 $f=\mu ,g=1$，则 $f\ast g=ϵ$。
代入到上面的结论中。
$$\begin{array}{rl}\mu (1)S(n)& =\sum _{i=1}^nϵ(i)-\sum _{i=2}^{n}1(i)S(\lfloor \frac{n}{i}\rfloor )\\ S(n)& =1-\sum _{i=2}^{n}S(\lfloor \frac{n}{i}\rfloor )\end{array}$$
然后就很愉快了。

const int maxn=2147483647;
const int maxm=2000000;
long long mSmu[maxm+10],Smu[1300];
bool vis[1300];
long long S(int n) {
	if(n<maxm) return mSmu[n];
	int x=maxn/n;
	if(vis[x]) return Smu[x];
	vis[x]=true;
	long long &ans=Smu[x]=1;
	for (int l=2,r;l<=n;l=r+1) {
		r=n/(n/l);
		ans-=(r-l+1)*S(n/l);
	} return ans;
}

const int maxn=2147483647;
const int maxm=2000000;
long long mSmu[maxm+10];
bool vis[maxm];map<int,int> Smu;
long long S(int n) {
	if(n<maxm) return mSmu[n];
	if(Smu.find(n)!=Smu.end()) return Smu[n];
	int &ans=(Smu[n]=1);
	for (int l=2,r;l<=n;l=r+1) {
		r=n/(n/l);
		ans-=(r-l+1)*S(n/l);
	} return ans;
}

$\phi$ 的前缀和

求 $\begin{array}{r}\sum _{i=1}^n\phi (i)\end{array}$，$1⩽n⩽2^{31}$
考虑到 $\phi \ast 1=id$，则选择 $f=\phi ,g=1$，则 $f\ast g=id$，代入。
$$\begin{array}{rl}g(1)S(n)& =\sum _{i=1}^n(f\ast g)(i)-\sum _{i=2}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )\\ S(n)& =\sum _{i=1}^{n}i-\sum _{i=2}^{n}S(\lfloor \frac{n}{i}\rfloor )\end{array}$$
前面可以根据小学奥数得到，后面除法分块即可，代码如下：

const int maxn=2147483647;
const int maxm=2000000;
long long mSphi[maxm+10],Sphi[1300];
bool vis[1300];
long long Sum_phi(int n) {
	if(n<maxm) return mSphi[n];
	int x=maxn/n;
	if(vis[x]) return Sphi[x];
	vis[x]=true;
	long long &ans=Sphi[x]=(long long)(1+n)*n/2;
	for (int l=2,r;l<=n;l=r+1) {
		r=n/(n/l);
		ans-=(r-l+1)*S(n/l);
	} return ans;
}

const int maxn=2147483647;
const int maxm=2000000;
long long mSphi[maxm+10];
bool vis[maxm];map<int,long long> phi;
long long Sum_phi(int n) {
	if(n<maxm) return mSphi[n];
	if(Sphi.count(n)!=Sphi.end()) return Sphi[n];
	long long &ans=Sphi[n]=(long long)(1+n)*n/2;
	for (int l=2,r;l<=n;l=r+1) {
		r=n/(n/l);
		ans-=(r-l+1)*S(n/l);
	} return ans;
}

如果题目同时要求 $\mu$ 和 $\phi$ 呢？用杜教筛求出 $\mu$ 之后其实不需要再用杜教筛求 $\phi$。
引理　$\begin{array}{r}\sum _{i=1}^n\phi (i)=\frac{1}{2}(1+\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor \lfloor \frac{n}{i}\rfloor \mu (i))\end{array}$
证明　可以发现，$\begin{array}{r}\sum _{i=1}^n\phi (i)\end{array}$ 其实就是 $1\sim n$ 中互质的数的对数。去除重复计算和只出现了一次的 $(1,1)$，可以得到下面这样的式子：
$$\sum _{i=1}^n\phi (i)=\frac{1}{2}(1+\sum _{i=1}^n\sum _{j=1}^n[(i,j)=1])$$

根据莫比乌斯反演，得：
$$\sum _{i=1}^n\sum _{j=1}^n[(i,j)=1]=\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor \lfloor \frac{n}{i}\rfloor \mu (i)$$

回代即可。
因此我们得到了一种更好的计算 $\phi$ 的方法——整除分块。
注意不需要记忆化，因为本来就没有递归计算。

long long Sum_phi(int n) {
    if(n<maxm) return mSphi[n];
    long long phi=0;
    for(long long l=1,r;l<=n;l=r+1){
    	r=n/(n/l);
    	phi+=(long long)(n/l)*(n/l)*(Sum_mu(r)-Sum_mu(l-1));
	} return (phi+1)/2;
}

上面代码的时间复杂度似乎是 $\Theta (\sqrt{n}\times n^{\frac{2}{3}})=\Theta (n^{\frac{1}{2}}\times n^{\frac{2}{3}})=\Theta (n^{\frac{7}{6}})>\Theta (n)$ 的，根本无法通过测试。但是考虑到每次杜教筛是记忆化的。因而相当于只考虑了一次杜教筛，即 $\Theta (n^{\frac{1}{2}}+n^{\frac{2}{3}})<\Theta (n)$。这也奠定了杜教筛在莫比乌斯反演中的作用，保证下面题表中所有的算法的时间复杂度的正确。

题表

题目	题解
Luogu P3768 简单的数学题	欧拉反演，杜教筛 $\phi \ast i{d}^2$