min_25筛

min_25筛

设有积性函数$f(x)$，满足$f(P_i)$是多项式，且能够快速计算出$f(P_i^k)$。求其前缀和。

一个函数

令$minp(x)$表示$x$的最小质因子。令$P$为质数集，其中第$1$个元素为$2$。设$g(n,j)={\sum }_{i=1}^{n}f(i)[i\in P\vee minp(i)>P_j]$。

考虑从$g(n,j-1)$推到$g(n,j)$。这个过程相当于删去那些最小质因数为$P_j$的合数。也就是说，这些数除以$P_j$后的最小质因数还是大于$P_{j-1}$的。那要删去的数的答案就是$f(P_j)\cdot g(\frac{n}{P_j},j-1)$吗？发现这时还多删了除以$P_j$后比$P_j$小的质数的答案，因此把这部分补回去。则转移就是$$g(n,j)=g(n,j-1)-f(P_j)\cdot [g(\frac{n}{P_j},j-1)-g(P_j-1,j-1)]$$
这个是没有特判各种情况时的转移。首先，若$P_j^2>n$，则显然这次什么也删不了，即$g(n,j)=g(n,j-1)$。另外还有边界，即$g(P_j,j)={\sum }_{i=1}^{j}f(P_j)$。

如果我们要计算$g(n,m)$，则我们只需要对$\sqrt{n}$内的质数计算即可。该部分的复杂度为$O(\frac{n^{\frac{3}{4}}}{\log n})$。

另一个函数

设$S(n,j)={\sum }_{i=1}^{n}f(i)[minp(i)\ge P_j]$。则$f$的前缀和就是$S(n,1)+f(1)$。

先预处理要求范围内的$g$。考虑如何转移$S$，将答案为质数的答案和合数的答案，其中质数直接计算，合数递归计算。质数的答案可以由$g$算出，是$g(n,j)-{\sum }_{i=1}^{j-1}f(P_i)$。对于合数，枚举其最小质因子及其出现次数。转移为$$S(n,j)=g(n,j)-\sum _{i=1}^{j-1}f(P_i)+\sum _{k=j}^{P_k^2\le n}\sum _{l=1}^{P_k^{l+1}\le n}f(P_k^l)\cdot S(\frac{n}{P_k^l},k+1)+f(P_k^{l+1})$$
总复杂度为$O(\frac{n^{\frac{3}{4}}}{\log n})$。