其他分享
首页 > 其他分享> > 进阶数论artalter级服务第三弹 杜教筛

进阶数论artalter级服务第三弹 杜教筛

作者:互联网

1.杜教筛

杜教筛是用来在低于线性的时间复杂度\((O(n^\frac{2}{3} )?)\)内求出积性函数的前缀和的算法

根据杜教筛的定义,我们设

\[S(n)=\sum_{i=1}^nf(i) \]

\[g是一个积性函数 \]

\[h(n)=f \times g \]

\[H(n)=\sum_{i=1}^nh(i) \]

那么有

\[\begin{aligned} H(n)&=\sum_{i=1}^nh(i)\\ &=\sum_{i=1}^n(f\times g)(i)\\ &=\sum_{i=1}^n\sum_{d|i}{g(d)f(\frac{i}{d})}\\ &=\sum_{i=1}^n\sum_{d|i}{g(d)f(\frac{i}{d})}\\ &=\sum_{d=1}^ng(d)\sum_{i=1}^nf(\frac{i}{d})[d|i] \\ &=\sum_{d=1}^ng(d)\sum_{i=1}^{\lfloor\frac{i}{d}\rfloor}f(i) \\ &=\sum_{d=1}^ng(d)S(\lfloor\frac{i}{d}\rfloor)\\ \end{aligned} \]

变化一下可得:

\[\begin{aligned} g(1)S(n)&=H(n)-\sum_{d=2}^ng(d)S(\lfloor\frac{i}{d}\rfloor) \end{aligned} \]

这就是杜教筛最根本的式子。
如果我们能找到一个g,使得g的前缀和\(g\times h\)的前缀和都很好求,那我们就能用整数分块快速求出S了

应用

\[\sum_{i=1}^n\mu(i) \]

我们知道

\[\mu \times I =\epsilon \]

\(\epsilon\)的前缀和永远是1,\(I\)的前缀和就是n

那我们就可以推出\(\mu\)的前缀和式子了

标签:artalter,frac,前缀,sum,杜教,ng,aligned,进阶
来源: https://www.cnblogs.com/artalter/p/16491424.html