Min_25筛学习笔记

Preface

• 模拟赛遇到了，于是学习了一下。

• 问题一般是求积性函数$f(x)$的前缀和，如$$\sum _{i=1}^{n}F(i)$$

Solution

• 记$P$表示素数集合，记$P_i$表示第$i$个质数，其中$|P|=\sqrt{n}$。

• 则此类问题$F(i)$需要满足：

$$\{\begin{array}{l}f(p)是一个关于p的简单多项式p\in P\\ f(p^c)一般可以快速计算出otherwise\end{array}$$

• 考虑记$$S(n,j)=\sum _{i=2}^{n}F(i)[i的最小质因子大于等于P_j]$$

• Min25筛会把$F_i$拆成$\sum a_i{i}^k$的形式，对于每个$k$分别筛后最后合并到一起。

递归版本

• 那么先考虑质数的贡献，即$i$为质数，且保证括号内条件的情况下的贡献。

• 若记$min(p)$表示$i$的最小质因子，则再记$$g(n,j)=\sum _{i=1}^n[i\in P,min(p)>P_j]F(i)$$

• 由定义，不难推出$$g(n,j)=\{\begin{array}{ll}g(n,j-1)& {P_j}^2>n\\ g(n,j-1)-F(P_j)[g(\frac{n}{P_j},j-1)-{\sum }_{i=1}^{j-1}F(P_i)]& {P_j}^2\le n\end{array}$$

• 现在还有一个问题是$g(n,0)$如何求出来。

• 事实上，我们观察这整一个递推过程，实际上很类似埃氏筛，每一次都用一个$P_j$去把一些合数的贡献给筛去，最终筛到$\sqrt{n}$时，保证了还没有被质数筛去的数在$n$范围内一定不是合数。

• 那么继续考虑埃氏筛，实际上一开始我们是把所有数都看成质数的。那么$g(n,0)$的初值也就是这个${\sum }_{i=1}^n{i}^k$，例如当$k=2$时，$g(n,0)=\frac{n(n+1)(2n+1)}{6}$.

• 现在回到原问题$$\sum _{i=1}^{n}F(i)$$

• 我们要求质数的贡献，实质上就是$g(n,|P|)-g(P_{j-1},j-1)$

• 那么可以推出$$S(n,j)=g(n,|P|)-\sum _{i=1}^{j-1}f(P_i)+\sum _{k\ge j}^{|P|}\sum _e(f({P_k}^e)S(\frac{n}{{P_k}^e},k+1)+f({P_k}^{e+1}))$$

• 事实上，后面那一部分式子就是合数的贡献，枚举质数$k$，次数$e$，然后与求$g$的思想是完全类似的，注意要加上一个$f({P_k}^{e+1})$就好了.

递推版本

• 实际上递推版本比递归更慢。因为它求的是$\sqrt{n}$个$\frac{n}{i}$的前缀和。

• 具体细节与递归版本类似。我们直接设$$S^{\prime }(n,i)=\sum _{i=2}^{n}F(i)[i的最小质因子不小于p_i或i是质数]$$.

• 可以得到：$$S^{\prime }(n,i)=\{\begin{array}{ccc}{S}^{\prime }(n,i+1)+{\sum }_{e\ge 1且p_i^{e+1}\le n}F(p_i^e)[S(\lfloor \frac{n}{p_i^e}\rfloor ,i+1)-g(p_i,i)]+F(p_i^{e+1})& & p_i^2\le n\\ S^{\prime }(n,i+1)& & p_i^2\ge n\end{array}$$

• 初值$S^{\prime }(n,|P|+1)=g(n,|P|)$

Some Problems

T1

LOJ#6235. 区间素数个数

事实上就是求$1\sim n$素数，这个可以直接通过求的$g$来得到，连$s$都不用求。

注意实现的时候只需要保存$\frac{n}{i}$的$\sqrt{n}$个值即可。

#include <bits/stdc++.h>

#define F(i, a, b) for (int i = a; i <= b; i ++)

const int N = 700000;

using namespace std;

long long n, m, id1[N], id2[N], w[N], g[N], z[N], sqr;
bool bz[N];

int main() {
	scanf("%lld", &n), sqr = int(sqrt(n));

	F(i, 2, sqr) {
		if (!bz[i])
			z[++ z[0]] = i, bz[i] = 1;
		F(j, 1, z[0]) {
			if (z[j] * i > sqr) break;
			bz[z[j] * i] = 1;
			if (i % z[j] == 0) break;
		}
	}

	for (long long i = 1, j; i <= n; i = j + 1) {
		j = n / (n / i); w[++ m] = n / i;
		if (w[m] <= sqr) id1[w[m]] = m; else id2[n / w[m]] = m;
		g[m] = w[m] - 1;
	}
	F(j, 1, z[0])
		for (int i = 1; i <= m && z[j] * z[j] <= w[i]; i ++) {
			int k = (w[i] / z[j] <= sqr) ? id1[w[i] / z[j]] : id2[n / (w[i] / z[j])];
			g[i] -= g[k] - (j - 1);
		}

	printf("%lld\n", g[1]);
}

T2

#6053. 简单的函数

• 注意把$f(p)=p-1$拆成两个函数$g(p)=p,h(p)=1$，然后相减，这样可以保证积性函数的性质。然后就是板子了。

#include <bits/stdc++.h>
#define F(i, a, b) for (int i = a; i <= b; i ++)

typedef long long LL;

const int N = 3e5;
const LL ny = 5e8 + 4, Mo = 1e9 + 7;

using namespace std;

bool bz[N]; int id1[N], id2[N], sqr;
LL g[N], h[N], z[N], w[N], s[N], n, m;

LL S(LL x, int y) {
	LL k = x <= sqr ? id1[x] : id2[n / x];
	LL sum = (g[k] - h[k] - (s[y - 1] - (y - 1)) + (y == 1) * 2) % Mo;
	for (LL i = y; i <= z[0] && z[i] * z[i] <= x; i ++)
		for (LL e = 1, p = z[i]; p * z[i] <= x; e ++, p *= z[i])
			sum = (sum + (z[i] ^ e) % Mo * S(x / p, i + 1) + (z[i] ^ (e + 1))) % Mo;
	return sum;
}

int main() {
	scanf("%lld", &n), sqr = int(sqrt(n));
	F(i, 2, sqr) {
		if (!bz[i])
			z[++ z[0]] = i, s[z[0]] = (s[z[0] - 1] + i) % Mo;
		F(j, 1, z[0]) {
			if (z[j] * i > sqr) break;
			bz[z[j] * i] = 1;
			if (i % z[j] == 0) break;
		}
	}
	for (LL i = 1, j; i <= n; i = j + 1) {
		j = n / (n / i); w[++ m] = n / i;
		if (w[m] <= sqr) id1[w[m]] = m; else id2[n / w[m]] = m;
		g[m] = (((w[m] + 2) % Mo * ((w[m] - 1) % Mo)) % Mo * ny) % Mo;
		h[m] = (w[m] - 1) % Mo;
	}
	F(j, 1, z[0])
		for (LL i = 1; i <= m && z[j] * z[j] <= w[i]; i ++) {
			LL k = (w[i] / z[j] <= sqr) ? id1[w[i] / z[j]] : id2[n / (w[i] / z[j])];
			h[i] = (h[i] - (h[k] - (j - 1))) % Mo;
			g[i] = (g[i] - z[j] * (g[k] - s[j - 1])) % Mo;
		}

	printf("%lld\n", n == 1 ? 1 : ((S(n, 1) + 1) % Mo + Mo) % Mo);
}

• 两种版本的时间复杂度都是$O(\frac{n^{\frac{3}{4}}}{log(\sqrt{n})})$。经实测，递归版会稍微快一些。