首先是埃式筛,就是从2开始,把所有2的倍数全都标记,第一次遍历时有未被标记的就是素数,然后再把他的倍数给标记。
#include <bits/stdc++.h>
using namespace std;
/****
埃氏筛
****/
const int MAX=1e6+10;
int n;
int prime[MAX];
int main (){
cin>>n;
for(int i=2;i<n;i++){
if(!prime[i]){
for(int j=i*i;j<n;j+=i){
prime[j]=1;
}
}
}
for(int i=2;i<n;i++){
if(!prime[i]){
cout<<i<<" ";
}
}
return 0;
}
有一个缺点是,比如6,既要被2标记,又要被3标记,就重复了很多次。所有有一种改进方式,就是线性筛,或者又叫做欧拉筛。当然了,时间复杂度只是趋近于线性,并不是真的为o(n),只是每个数都只会被判断一次,所以效率很高。
#include <bits/stdc++.h>
using namespace std;
/****
线性筛
****/
const int MAX=1e6+10;
int n;
int prime[MAX];
int vis[MAX];
int main (){
cin>>n;
int cnt=0;
for(int i=2;i<n;i++){
if(!vis[i]){
prime[cnt++]=i;
}
for(int j=0;j<cnt&&i*prime[j]<n;j++){
vis[i*prime[j]]=1;
if(i%prime[j]==0){
break;
}
}
}
for(int i=0;i<cnt;i++){
cout<<prime[i]<<" ";
}
return 0;
}
prime数组存储的是素数。核心代码是嵌套的那个for循环。分析这个循环,如果当前判断的数字 i 是个素数,那么好说,将i与小于i的素数相乘 给筛掉,这样保证不会重复筛选。当i是个合数时,i可以分解为多个素数相乘的结果,这样,只筛 i 的最小质因子与 i 的乘积结果 ,这样筛选也是不会重复的,这样筛的数只被最小的质因子给筛了。
其它的证明可以参考
https://blog.csdn.net/dinosoft/article/details/5829550