转载自:腾讯云 - 杜教筛入门
ps:其实我就是看看里面的公式,原文公式不显示了。
Orz OO0OOO00O0OOO0O00OOO0OO - 杜教筛的发明者
前置知识
狄利克雷卷积
杜教筛
套路
杜教筛是用来求一类积性函数的前缀和
它通过各种转化,最终利用数论分块的思想来降低复杂度
假设我们现在要求\(S(n) = \sum_{i = 1}^n f(i)\),\(f(i)\)为积性函数,\(n \leqslant 10^{12}\)
直接求肯定是不好求的,不过现在假设有另一个积性函数\(g\)
我们来求它们狄利克雷卷积的前缀和
\[(g * f) = \sum_{i = 1}^n \sum_{d \mid i} g(d) f(\frac{i}{d})\]
\[= \sum_{d = 1}^n g(d) \sum_{d|i} f(\frac{i}{d})\]
\[= \sum_{d = 1}^n g(d) \sum_{i = 1}^{\frac{n}{d}}f(i)\]
\[= \sum_{d = 1}^n g(d) S(\frac{n}{d})\]
然后就化不动了,不过我们发现我们化出了\(\frac{n}{d}\)!excited
但是\(S(n)\)怎么求呢?
容斥一下
\[g(1)S(n) = \sum_{d = 1}^n g(d)S(\frac{n}{d}) - \sum_{d = 2}^n g(d)S(\frac{n}{d})\]
前半部分是狄利克雷卷积的前缀和的形式
后半部分可以数论分块。这样看起来就好搞多了
现在我们的问题是,如何选择\(g\)才能使得上面这个式子好算
这个就要因情况而定了
下面煮几个典型栗子
\(\mu\)函数
定理:\(\sum_{d \mid n} \mu(d) = [n = 1]\)
那么我们如果选择\(g = e\),\(e\)为原函数,\(e = [n = 1]\)
\(g\)与\(\mu\)的卷积的前缀和肯定为\(1\)
上面的式子变为
\(S(n) = 1 - \sum_{d = 2}^n S(\frac{n}{i})\)
后半部分直接数论分块就好
\(\varphi\)函数
定理:\(\sum_{d \mid n}\varphi(d) = n\)
我们的\(g\)还是选\(e\)做卷积
我们要求得式子变为
\[S(n) = \frac{n * (n + 1)}{2} - \sum_{d = 2}^n S(\frac{n}{i})\]
前半部分\(O(1)\)算,后半部分数论分块