筛选质数

首先，先简单介绍一下质数。
如果一个数能被除1和其本身以外的数整除，那它就是合数，否则就是质数。
1 既不是质数也不是合数。
接下来，介绍三种筛选质数的方法。

一、普通筛选法

最简单的筛选方法，直接通过定义然后枚举写出代码

//普通
bool isPrime(int x) {
	if (x <= 1) return false;
	int sqr = (int)sqrt(x * 1.0);
	for (int i = 2; i <= sqr; i++)
		if (x % i == 0) return false;
	return true;
}
void FindPrime(vector<int>& prime, int maxn) {
	bool isP[1005];
	for (int i = 2; i <= maxn; i++) {
		if (isPrime(i)) {
			prime.push_back(i);
			/*isP[i] = true;*/
		}
	}
}

二、埃式筛选法

也很简单，利用任何一个数的倍数是合数的定理得出代码,复杂度O(nloglogn)

//埃拉托斯特尼筛
void FindPrime2(vector<int>& prime, int maxn) {
	bool isP[1005];
	memset(isP, true, sizeof(isP)); //先把所有数都定义为质数
	for (int i = 2; i <= maxn; i++) {
		if (isP) {
			prime.push_back(i);
			for (int j = i * 2; j <= maxn; j++)
				isP[j] = false;
		}
	}
}

三、欧拉筛选法

上面的线性筛还是有一些细节需要优化的，比如一个数可能被筛了很多遍，如6既被2筛了一遍，又被3筛了一遍。
所以我们可以用 “最小质因数 × 最大因数（非自己） = 这个合数” 的途径删掉。
即对于每一个数（无论质数合数）x，筛掉所有小于x最小质因子的质数乘以x的数。
代码中的if (i % prime[j] == 0) break;即实现了这个核心思想。
具体证明不再这论述，感兴趣的小朋友可自行wiki或百度百科。
由于每个数只被筛一次，时间复杂度为 O(n)，所以是真正的线性筛。

//欧拉筛
void FindPrime3(vector<int>& prime, int maxn) {
	bool isP[1005];
	memset(isP, true, sizeof(isP)); //先把所有数都定义为质数
	for (int i = 2; i <= maxn; i++) {
		if (isP[i]) {
			prime.push_back(i);
			for (int j = 0; j < prime.size() && i * prime[j] <= maxn; j++) {
				isP[i * prime[j]] = false;
				if (i % prime[j] == 0) break; //保证了一个数只被筛一次
			}
		}
	}
}