质数筛法，欧拉筛与埃氏筛的比较

埃氏筛为常见的质数筛法，其原理为利用已知的质数作为因子将其他合数标记，则剩下的未标记的均为质数
其代码如下：
此算法的复杂度为O(NloglogN)

#include <bits/stdc++.h>
using namespace std;
const int N=1000;
bool b[N];
int main(){
	for(int i=2;i<=sqrt(N);i++){
		//i作为N的一个因子最大可达sqrt(N) 
		if(!b[i]){
			for(int j=i;j<=N/i;j++){
				//i为小因子，j为大因子，故j从i的值开始计数 
				b[i*j]=1;
			}
		}
	}	
	for(int i=2;i<=N;i++)
		if(!b[i])cout<<i<<" ";
	return 0;
}

该算法很明显会重复的筛选数，如12既会被2作为小因子筛出，也会被3作为小因子筛出。
故下面介绍一种以唯一最小因子为标记的筛法：欧拉筛

#include <bits/stdc++.h>
using namespace std;
int b[10000],pre[1000];
int main(){
	int n,k=0;
	cin>>n;
	for(int i=2;i<=n;i++){
		if(!b[i])pre[++k]=i;//未标记则必然为素数，素数记录数组+1
		//关键处1
		for(int j=1;j<=k;j++){
			if(pre[j]*i>n)break;//若超出n的范围则结束当前i
			b[pre[j]*i]=1;//标记合数
			if(i%pre[j]==0)break;
			//关键处2
		}
	}
	for(int i=1;i<=k;i++)
		cout<<pre[i]<<" ";
	return 0;
} 

欧拉筛的原理为：
不管为合数还是素数均可到“关键处1”
1.若i为素数的话，则一个大的素数i和任意一个小于i的素的乘积必然不会和其他的筛出的数重复，即n=p1*p2;
2.若i为合数的话，则 i=p1*p2*p3…*pn,p均为素数,则p1必然为最小的素数。
根据“关键处2”,i%pre[j]==0则退出循环，其原因为若pre[j]能整除i则说明pre[j]为i的最小因子，此时若再用i*pre[j]生成合数则pre[j]已经不是新生成的合数的最小因子了，故需要break.
因此该筛法的另一种写法为：

#include <bits/stdc++.h>
using namespace std;
const int N=1000;
int pre[N],v[N],k,n;
int main(){
	for(int i=2;i<=N;i++){
		if(v[i]==0){
			pre[++k]=i;
			v[i]=i;//当前i的最小质因子为i
			cout<<i<<" ";
		} 
		for(int j=1;j<=k;j++){
		//若pre[j]比i的最小质因子要大，则退出，因为基于pre[j]生成的合数的最小因子为v[i]这不符合预设
			if(pre[j]>v[i] || pre[j]*i>N)break;
			v[pre[j]*i]=pre[j];//新生成的合数的最小质因子为pre[j]
		}
	}

	return 0;
}

参考：
1.信息学奥赛数学一本通
2.算法进阶指南