Powerful Number筛

Powerful Number筛

要求积性函数$f(n)$的前缀和$F(n)$。

Powerful Number(PN)

定义：一个数$n={\prod }_{i=1}^m{p}_i^{c_i}$是PN当且仅当$\forall 1\le i\le m,c_i>1$。

性质1：所有PN都可以写成$a^2{b}^3$的形式

证明1：当$c_i$为偶数时$p_i^{c_i}$归入$a$，否则$p_i^3$归入$b$，$p_i^{c_i-3}$归入$a$

性质2：$n$以内的PN只有$O(n)$个

证明2：根据证明1，一组$(a,b)$可以唯一的还原出一个PN，并且每个PN都可以唯一对应一组$(a,b)$，所以只需枚举$a$然后算$b$个个数的和即可证明：
$${\int }_1^{\sqrt{n}}\sqrt[3]{\frac{n}{x^2}}\mathrm{d}x=O(\sqrt{n})$$
找PN可以先把$\sqrt{n}$内的质数筛出来然后爆搜，由于PN只有$O(\sqrt{n})$个，所以复杂度就是$O(\sqrt{n})$。

PN筛

构造一个易于求前缀和的积性函数$g$，且满足$\forall p\in P,g(p)=f(p)$。

假定有函数$h$满足$f=g\ast h$，这里的$\ast$是狄利克雷卷积，显然$h$也是积性函数。

根据$h$的假设，$\forall p\in P,f(p)=g(1)h(p)+g(p)h(1)=g(p)+h(p)$，再根据$g$的构造，可以得知$\forall p\in P,h(p)=0$。由于$h$是积性函数，所以有$\forall n\notin \mathbf{PN},h(n)=0$，那么$h$只有$O(\sqrt{n})$个有效值。

前面都是假设存在这样的$h$，现在考虑怎么计算$h$。

显然只要能算$h(p^c)$就行了，根据假设$f(p^c)={\sum }_{i=0}^{c}g(p^i)h(p^{c-i})$，移项得$h(p^c)=f(p^c)-{\sum }_{i=1}^{c}g(p^i)h(p^{c-i})$。当然，也有可能$h(p^c)$存在只用$p,c$表达的式子，可以先表后猜再证。综上，$h$是一定存在的，假设成立。

那么现在唯一剩下的就是怎么算$F$了，根据$h$的假设，有如下：
$$\begin{array}{rl}F(n)& =\sum _{i=1}^{n}f(i)\\ & =\sum _{i=1}^n\sum _{d|i}g(d)h(\frac{i}{d})\\ & =\sum _{d=1}^{n}h(d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }g(i)\\ & =\sum _{1\le d\le n,d\in \mathbf{PN}}h(d)G(\lfloor \frac{n}{d}\rfloor )\end{array}$$

复杂度

在$\sqrt{n}$内大约有$O(\frac{\sqrt{n}}{\log \sqrt{n}})$个质数，对于一个质数$p$，指数最多是$O(\log \sqrt{n})$，计算$h(p^c)$时复杂度为$O(c)$，所以前面处理$h$的复杂度为$O(\frac{\sqrt{n}}{\log \sqrt{n}}(\log \sqrt{n}{)}^2)=O(\sqrt{n}\log \sqrt{n})$。

对于$g$的计算，可能$O(1)$算，如果比较麻烦，那么套杜教筛$O(n^{\frac{2}{3}})$也是可以的。

总结

PN筛一般步骤如下：

1. 构造$g$和快速计算$G$的方法
2. 找PN并计算$h$
3. 套式子计算$F$

计算$G$时可能要用杜教筛，复杂度似乎比直接杜教筛更劣，但是考虑到对$f$杜教筛有时候可能构造不出来合适的$g,h$所以PN筛还是很有用的。