1、首先定义大小为n+1的bool型数组prime[n+1]
2、将下标为奇数的赋值为true,偶数赋值为false
3、.然后:
for( i=3; i<=sqrt(n); i+=2 )
{ if(prime[i]) //将素数的所有倍数赋值为false
for( j=i+i; j<=n; j+=i ) prime[j]=false;
}
一个简单的筛素数的过程:n=30。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
第 1 步过后2 4 ... 28 30这15个单元被标成false,其余为true。
第 2 步开始:
i=3; 由于prime[3]=true, 把prime[6], [9], [12], [15], [18], [21], [24], [27], [30]标为false.
i=4; 由于prime[4]=false,不在继续筛法步骤。
i=5; 由于prime[5]=true, 把prime[10],[15],[20],[25],[30]标为false.
i=6>sqrt(30)算法结束。
4、输出:
先将prime[2]赋值为true
输出值为true的下标
把prime[]值为true的下标输出来:
for(i=2; i<=30; i++)
if(prime[i]) printf("%d ",i);
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
using namespace std;
int main(){
int n,i,j;
cin>>n;
bool prime[n+1];
for(i=0; i<=n; ++i){
if(i%2==0)
prime[i]=false;
else
prime[i]=true;
}
for(i=3; i<sqrt(n); i+=2){
if(prime[i])
for(j=i+i; j<=n; j+=i) prime[j]=false;
}
prime[2]=true;
for(i=2; i<=n; ++i)
if(prime[i])
cout<<i<<' ';
cout<<endl;
return 0;
}
输入:50
输出:2 3 5 7 11 13 17 19 23 29 31 37 41 43 47