素数筛总结

原始筛：

        就是根据素数的定义：素数的因子只有1和它本身。

        所以就只需要从2开始遍历到根号n，看定义是否成立。

                若符合定义，则是素数，反之，就不是素数。代码如下：

void isprime1(int n) {
	for(int i=2; i<=sqrt(n); i++)
		if(n%i==0) return false;
	return true;
}

---

埃氏筛：

        埃氏筛又称为普通筛，其原理就是从2开始打表，若2为质数，则2的倍数必然不是质数，因此便让表中的4 6 8 10。。。。全部删除，最后剩下的就是所要的素数。其运行速度要远高于普通筛。

bool is_prime[MAX+5];//标记是否为素数
int prime[MAX/2];//存放素数
int cnt=0;//已存放的素数个数
void isprime2() {
	memset(is_prime,true,sizeof is_prime);
	for(int i=2; i<=MAX; i++) {
		if(is_prime[i]) {
			prime[cnt++]=i;
			for(int j=i+i; j<=MAX; j+=i) //二次筛选法:i是素数，则i的倍数必不为素数，筛掉
				is_prime[j]=false;
		}
	}
}

欧拉筛：

        欧拉筛的原理就是利用：最小质因数 × 最大因数（非自己） = 这个合数。

        只需利用最小质因子即可。原理和埃氏筛较为类似，但并不是单纯的一个数一个数的找倍数，其效率又远远高于埃氏筛。我认为其中也利用打表的思路。

#include <cstdio>
#include <cstring>

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

void GetPrime(int n)//筛到n
{
	memset(isPrime, 1, sizeof(isPrime));
	//以“每个数都是素数”为初始状态，逐个删去
	isPrime[1] = 0;//1不是素数
	
	for(int i = 2; i <= n; i++)
	{
		if(isPrime[i])//没筛掉 
			Prime[++cnt] = i; //i成为下一个素数
			
		for(int j = 1; j <= cnt && i*Prime[j] <= n/*不超上限*/; j++)//注意j要<=cnt, 至少要把自身的平方给筛掉把 
		{
        	//从Prime[1]，即最小质数2开始，逐个枚举已知的质数，并期望Prime[j]是(i*Prime[j])的最小质因数
            //当然，i肯定比Prime[j]大，因为Prime[j]是在i之前得出的
			isPrime[i*Prime[j]] = 0;
            
			if(i % Prime[j] == 0)//i中也含有Prime[j]这个因子
				break; //重要步骤。见原理
                //若满足i%Prime[j]==0条件，就可以保证Prime[j]是i的因子，而欧拉筛是利用最小质因子，因而就不需要i了，所以break即可。
		}
	}
}

int main()
{
	int n, q;
	scanf("%d %d", &n, &q);
	GetPrime(n);
	while (q--)
	{
		int k;
		scanf("%d", &k);
		printf("%d\n", Prime[k]);
	}
	return 0;
}

最近看了一篇大神的题解，似乎发现了更高效简洁的筛法：

void init(int n){
	ntp[1]=1;
	for(int i=2;i<=n;++i){
		if(!ntp[i]){
			for(int j=i+i;j<=n;j+=i)ntp[j]=1;	
		}	
	}	
}

ntp数组的值若为1则说明该数不是素数，否则是