学步园

一般筛法求素数代码

思路：素数的倍数不是素数。

#include<iostream>
using namespace std;
const int M = 2000000;
int a[M] = {1,1};
int prime[M]={0};
int k=0;
void init()
{
    for(int i=2;i<=M;i++)
    {
        if(!a[i])
        {
            prime[k++] = i;
            for(int j=i+i;j<=M;j+=i)
                    a[j]= 1;
        }
    }
}
int main()
{
    int n;
    while(cin>>n)
    {
        for(int i=0;i<n;i++)
        cout<<prime[i]<<" ";
        cout<<endl;
    }
}

快速线性筛法

原理：1.一个合数是由n个素数的乘积所组成的

2.素数的倍数不是素数。

建议无论手模拟一个两个过程 就明白谁原理和优势了。

#include<iostream>
using namespace std;

const int Max=2000000;
long long prime[Max] = {0};
int k=0;
int a[Max]={1,1};

void init()
{
    for(long long i=2;i<Max;i++)
    {
        if(!a[i])//
        prime[k++] = i;
        for(long long j=0;j<k&&i*prime[j]<Max;j++)
        {
            a[i*prime[j]] = 1;
            if(i%prime[j]==0)//
                break;
        }
    }
}

int main()
{
    int n;
    init();
    while(cin>>n)
    {
        for(int i=0;i<n;i++)
       cout<<prime[i]<<" ";
        cout<<endl;
    }
}

推荐：http://blog.csdn.net/dinosoft/article/details/5829550