埃氏筛选素数及其变形



埃氏筛法

/*用“埃氏筛法”求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";
    }}