杜教筛学习笔记

看着身边的大佬们这么快就学了各种筛，很是紧张啊……
接下来强行学习一下

简介

用途：$O(n^{\frac{2}{3}})$求积性函数的前缀和。

本质就是推式子。
设有某个函数$f(i)$，我们要求${\sum }_{i=1}^{n}f(i)$
根据具体情况建出辅助函数$g(i)$。求狄利克雷卷积$h=g\ast f$
记$S(n)={\sum }_{i=1}^{n}f(i)$
开始推式子：
$$\begin{gathered} \sum _{i=1}^{n}h(i)=\sum _{i=1}^n\sum _{d|i}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{gathered}$$
将右边式子第一项拆出来：
$$=g(1)S(n)+\sum _{d=2}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )$$
移项得$g(1)S(n)={\sum }_{i=1}^{n}h(i)-{\sum }_{d=2}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )$
这个就是杜教筛的套路式了。
使用杜教筛的时候，注意${\sum }_{i=1}^{n}h(i)$和$g(i)$要能够快速地求出来。
用整除分块来搞，时间复杂度据说是$O(n^{\frac{2}{3}})$

于是，学会这个套路之后，问题主要就是如何选取$g$（以及$h$）。

---

实现

洛谷P4213 【模板】杜教筛（Sum）
注意实现的时候，最好把$n$比较小的$S(n)$都预处理出来。
否则哈希容易爆炸。

using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <unordered_map>
#define ll long long
#define M 1000000
int n;
int p[M+10],np;
bool inp[M+10];
int phi[M+10],mu[M+10];
ll sum1[M+10],sum2[M+10];
unordered_map<int,ll> ans1,ans2;
ll get1(int n){
	if (n<=M)
		return sum1[n];
	auto p=ans1.find(n);
	if (p!=ans1.end())
		return p->second;
	ll res=(ll)(n+1)*n>>1;
	for (unsigned i=2,j=2;i<=n;i=j+1){
		j=n/(n/i);
		res-=(j-i+1)*get1(n/i);
	}
	return ans1[n]=res;
}
ll get2(int n){
	if (n<=M)
		return sum2[n];
	auto p=ans2.find(n);
	if (p!=ans2.end())
		return p->second;
	ll res=1;
	for (unsigned i=2,j=2;i<=n;i=j+1){
		j=n/(n/i);
		res-=(j-i+1)*get2(n/i);
	}
	return ans2[n]=res;
}
int main(){
	phi[1]=mu[1]=1;
	for (int i=2;i<=M;++i){
		if (!inp[i]){
			p[++np]=i;
			phi[i]=i-1;
			mu[i]=-1;
		}
		for (int j=1;j<=np && i*p[j]<=M;++j){
			inp[i*p[j]]=1;
			if (i%p[j]==0){
				phi[i*p[j]]=phi[i]*p[j];
				mu[i*p[j]]=0;
				break;
			}
			phi[i*p[j]]=phi[i]*(p[j]-1);
			mu[i*p[j]]=-mu[i];
		}
	}
	for (int i=1;i<=M;++i){
		sum1[i]=sum1[i-1]+phi[i];
		sum2[i]=sum2[i-1]+mu[i];
	}
	int T;
	scanf("%d",&T);
	while (T--){
		scanf("%d",&n);
		printf("%lld %lld\n",get1(n),get2(n));
	}
	return 0;
}

---

参考资料
杜教筛