AK F.*ing leetcode 流浪计划之欧拉筛选素数

欢迎关注更多精彩
关注我，学习常用算法与数据结构，一题多解，降维打击。

素数筛选

素数筛选的基本过程就是从小到大遍历数字，如果之前没有被标记过，就是素数。然后拿当前遍历数字的倍数去标记成合数。效率的核心就是如何拿已经遍历过的数字去筛选掉合数。

简单倍数法

bool isPrime[100000010] = { 0 };
//isPrime[i] == 1表示：i是素数
int Prime[6000010], cnt = 0;
//Prime存质数

void simplePrime(int n) {
	cnt = 0;
	memset(isPrime, 1, sizeof(isPrime));
	isPrime[1] = 0;
	for (lld i = 2; i <= n; ++i) {
		if (isPrime[i]) Prime[cnt++] = i;
		for (lld j = 2; j*i<=n; ++j) {
			isPrime[i*j] = 0;
		}
	}
}

利用素数优化步长

上述算法中里面一层循环会对同一个合数筛选多次。
比如n=20;

i	i*j
2	4,6,8,10,12,14,16,18,20
3	6,9,12,15,18
4	8,12,16,20
5	10,15,20
6	12,18
7	14
8	16
9	18
10	20

i超过10以后就没有筛选到合数了。
上面表格中有很多重复的合数。
可以减少j的枚举量来优化。

我们知道一个合数x肯定可以拆解成一个素数m乘以另一个数p.
$x=m\ast p(其m是素数，m<=\sqrt{n},p>1)$
所以，j可以只枚举已经得到的素数。

void simplePrime(int n) {
	cnt = 0;
	memset(isPrime, 1, sizeof(isPrime));
	isPrime[1] = 0;
	for (lld i = 2; i <= n; ++i) {
		if (isPrime[i]) Prime[cnt++] = i;
		// 只枚举已知的素数
		for (int j = 0; j < cnt && 1LL * i * Prime[j] <= n;++j) {
			isPrime[1LL * i * Prime[j]] = 0;
		}
	}
}

欧拉筛思路

再次对上述算法进行模拟
n=20;

i	i*Prime[j]
2	4
3	6,9
4	8,12
5	10,15
6	12,18
7	14
8	16
9	18
10	20

这次改进后，筛选合数的次数比直接倍数法减少了很多，但是还有重复筛选，一旦数量多了以后也会影响时间复杂度。

上述算法会导致重复的原因如下：
$设x=p_0^{s_0}\cdot p_1^{s_1}...p_n^{s_n},p_i为质数，那么x就有可能被每个质数筛选到$

欧拉筛选法希望每一个合数只被最小的质数筛选到，也就是上述过程中最小的p0筛选。
$把x拆成x=p_0\cdot p_0^{s_0-1}\cdot p_1^{s_1}...p_n^{s_n}$

对应到上述代码中就是Prime[j]不能在i作为因子中出现过。

bool isPrime[100000010] = { 0 };
//isPrime[i] == 1表示：i是素数
int Prime[6000010], cnt = 0;
//Prime存质数

void initPrime(int n) {
	cnt = 0;
	memset(isPrime, 1, sizeof(isPrime));
	isPrime[1] = 0;
	for (int i = 2; i <= n; ++i) {
		if (isPrime[i]) Prime[cnt++] = i;
		for (int j = 0; j < cnt && i * Prime[j] <= n; ++j) {
			isPrime[i * Prime[j]] = 0;
			if (i % Prime[j] == 0)break;
		}
	}
}

正确性证明：

完备性：每个合数都会被筛选到
假设x=p*h, p 是x最小质因子，h>1。

当i=h时，Prime[j]在到达p之前肯定不会退出循环，保证了完备性。

唯一性：每个合数只会被筛选一次
反证法：
假设x=ph, p 是x最小质因子，h>1，p1>p。
x是不可能被p1h1筛选到，因为h1中必然有p因子，所以Prime[j]会到达p时就会退出。

故x只会被p*h筛选掉，只会筛选一次。根据这个结论可以得出欧拉筛的最复杂度是O(n)的。

求解积性函数应用

积性函数指对于所有互质的整数a和b有性质f(ab)=f(a)f(b)的数论函数。

欧拉函数

定义与公式

$欧拉函数\phi (x)定义为所有小于等x正整数中与x最大公约数为1的个数$
欧拉函数公式

$\phi (x)=\{\begin{array}{lc}1& x=1\\ x-1& x是素数\\ x\cdot \prod _{p_i|x}\frac{p_i-1}{p_i}& x>1,p_i是x的质因子\end{array}$

欧拉筛过程中求解欧拉函数

在程序运行过程中，我们可以知道每个合数被拆解的情况。

$令P(x)=\prod _{p_i|x}\frac{p_i-1}{p_i}$

$x=p\ast h,p是质数,\phi (h)=h\ast P(h)是已经算出来了,也可以求得P(h)=\phi (h)/h$

$如果p是h质因子，那么P(x)=P(h),\phi (x)=\phi (h)/h\ast x=\phi (h)\ast p$

$如果p不是h质因子，那么P(x)=P(h)\ast \frac{p-1}{p},\phi (x)=\phi (h)/h\ast \frac{p-1}{p}\ast (p\ast h)=\phi (h)\ast (p-1)$

bool isPrime[100000010] = { 0 };
//isPrime[i] == 1表示：i是素数
int Prime[6000010], cnt = 0;
//Prime存质数
lld phi[N] = { 0 };

void initPrime(int n) {
	cnt = 0;
	memset(isPrime, 1, sizeof(isPrime));
	isPrime[1] = 0;
	phi[1] = 1;
	for (int i = 2; i <= n; ++i) {
		if (isPrime[i]) {
			Prime[cnt++] = i;
			phi[i] = i - 1;
		}
		for (int j = 0; j < cnt && i * Prime[j] <= n; ++j) {
			isPrime[i * Prime[j]] = 0;
			if (i % Prime[j] == 0) {
				phi[i * Prime[j]] = phi[i] * (Prime[j]);
				break;
			}
			else {
				phi[i * Prime[j]] = phi[i] * (Prime[j]-1);
			}
		}
	}
}

---

本人码农，希望通过自己的分享，让大家更容易学懂计算机知识。创作不易，帮忙点击公众号的链接。