一定能让你理解的素数筛法——埃氏筛法和欧式筛法

先上代码

1. 埃氏筛法

#include<bits/stdc++.h>
using  namespace std;
//时间复杂度（nloglogn)
//适用10^6下
const int MAX = 1000000;
void FindPrime(vector<int>& prime);
int main() {
	vector<int>prime;
	FindPrime(prime);
	for (auto i : prime)
		printf("%d ", i);
	return 0;
}
void FindPrime(vector<int>& prime) {
	array<int, MAX>flag;
	int j = 0;
	fill(flag.begin(), flag.end(), 0);
	for (int i = 2; i < MAX; i++) {
		if (flag[i] == 0) {
			prime.push_back(i);
		}
		for (int k = i + i; k < MAX; k += i) {
			flag[k] = 1;
		}
	}
}

1. 欧式筛法

#include<bits/stdc++.h>
using  namespace std;
//时间复杂度（n)
//适用低于10^7
const int MAX = 10000000;
void FindPrime(vector<int>& prime);
int main() {
	vector<int>prime;

	FindPrime(prime);
	for (auto i : prime)
		printf("%d ", i);
	return 0;
}
void FindPrime(vector<int>& prime) {
	array<int, MAX>flag;
	fill(flag.begin(), flag.end(), 0);
	for (int i = 2; i < MAX; i++) {
		if (flag[i]==0) {
			prime.push_back(i);
		}
		for (int j = 0; j < prime.size() && i * prime[j] < MAX; j++) {
			flag[i * prime[j]] = 1;
			if (i % prime[j] == 0) {
				break;
			}
		}
	}
}

一些细节： 埃氏和欧式的共同点都是从已知的最小素数 2起手，并将2放入素数数组，然后开始针对2的倍数进行筛选。
其他： 假如是vs编译的话，这么大的数据量需要修改堆栈区的大小，否则一定会抛出异常。
我的修改堆栈区大小为：99999999——》100m

分析

在开始分析之前，请忘了数论的那些术语，质因数等等
知道什么是素数即可
进行手动计算过程

原因是：2分解成乘法，可有12=2×6；12=3×4
而其中最小的质数组合为2×6
也就是说欧式筛法，每次筛的是最小素数的组合，这样保证了一定不会多筛重复的数

从代码来分析

欧式筛法的脱出条件共有3条

1. j < prime.size() 

内部脱出：遍历整个素数数组，筛掉i与素数的乘法组合

2. i * prime[j] < MAX

外部进入判断：如果第一次筛选已经超过上界，则脱出

3. if (i % prime[j] == 0)

内部脱出：保证了筛掉的为最小素数的乘法组合（专业术语叫质因数）

其中1和3的组合构成素数筛选的核心判断
2其实可有可无，只为了控制边界而已

后记

如果只是想学习算法，不想做科研，大可不必深入去研究数论的那些术语，大部分算法都可通过手工模拟去理解核心原理，再去理解算法。
但是 假如有数论的基础，相对看起代码也会更快，不至于都去手工模拟。
所以如何学习，根据自己的目标量力而行即可。