筛选素数的三种方法(包含欧拉筛)

具体方法的介绍可以看下面详细代码，有注释

//选出素数方法
#include<iostream>
#include<cstring>
#include<map>
#include<vector>
#include<algorithm>
using namespace std;
const int maxn = 1010;
//方法1：试除法
bool s_prime(int a)
{
	if (a <= 1) return false;
	for (int i = 2; i*i < a; ++i) {
		if (a%i == 0) return false;
	}
	return true;
}

//方法2：埃氏筛法，时间复杂度O(nloglog2(n)),缺点，比如10000会被重复筛除很多次（2,5,10...)
bool num[maxn];
int E_sieve(int n)  //筛除0~n之间的非素数
{
	for (int i = 0; i <= n; ++i) num[i] = false;
	for (int i = 2; i*i <= n; ++i) {
		if (!num[i]) {
			for (int j = 2*i; j <= n; j = j + i) {  
				num[j] = true;
			}
			//改进,//当i为素数时，能够被i整除的数肯定不是素数,这里j=2*i可以改进为j=i*i，因为当i=5时，2*5,3*5,4*5在前面i为2,3,4时已经算过了
		}
		int k = 0;
		int prime[maxn];
		for (int i = 2; i <= n; ++i) {
			if (!num[i]) prime[k++] = i;  //筛出素数
		}
	}
}


//方法3：欧拉筛：时间复杂度O(n)
//对于每一个数，要么他是素数，要么它是素数的倍数
//所以用最小的质因子来筛，又因为主循环从2~n,所以每个数都会被筛到，而且只会被筛一次
int prime[maxn], i, j, countnum = 1;
bool num[maxn];
int oula(int n)
{
	memset(num, true, sizeof(num));
	for (i = 2; i <= n; ++i) {  //i作为本身被筛掉或者作为因子来筛选其他数
		if (num[i]) prime[countnum++] = i;
		for (int j = 1; j < countnum; ++j) {
			if (i*prime[j] > n) break;  //超出所需求解的上限，退出
			num[i*prime[j]] = false;   //排除非素数
			if (i%prime[j] == 0) break;  //这里的i就是表示倍数，6*3=2*9，因为6是2的倍数，如果6*2的时候没有退出，那么当i到9的时候，18又会被筛掉一次
		}
	}
}
int main()
{
	return 0;
}