埃氏筛法
/*用“埃氏筛法”求2~1/*00以内的素数。2~100以内的数,先去掉2的倍数,再去掉3的倍数,
再去掉4的倍数,……依此类推,
最后剩下的就是素数。*/
一传统:
#include<stdio.h>
#include<iostream>
using namespace std;
#define N 10000000
int prime[N];
bool flag[N];
int f(int n)
{
for(int i=0;i<=n;i++){
if(i%2==0)
flag[i]=false;
else
flag[i]=true;
}
int cnt=0;
flag[2]=true;
for(int i=2;i<=n;i++){
if(flag[i]==true){
prime[cnt++]=i;
for(int j=i*2;j<=n;j+=i)
flag[j]=false;
}
}
return cnt;
}
void print(int count)
{
for(int i=0;i<count;i++)
cout<<prime[i]<<"\t";
cout<<endl;
}
int main()
{
int n;
while(scanf("%d",&n)!=EOF){
int count=f(n);
print(count);
}
}
二改进:
#include<iostream>
#include<cstdio>
using namespace std;
const int SIZE = 1e7;
int prime[SIZE]; // 第i个素数
bool is_prime[SIZE]; //true表示i是素数
int slove(int n)
{
int p = 0;
for(int i = 0; i <= n; i++)
is_prime[i] = true; //初始化
is_prime[0] = is_prime[1] = false; //0,1不是素数
for(int i = 2; i <= n; i++){
if(is_prime[i]) //这里比较巧妙, 我只是意会{
prime[p++]=i; //计算素数的个数,也记录下了素数
for(int k=0,j=i*prime[k];k<p && j<=n;k++,j=i*prime[k]) {
is_prime[j] = false;
} }
return p;
}
int main()
{
int n;
while(cin >> n){
int res = slove(n);
for(int i= 0; i < res; i++)
cout<<prime[i]<<"\t";
}}