素數篩選法

題目鏈接
最開始篩選素數的代碼

bool is_prime(int n )
{
    int sqrt_n=sqrt(n);
    for(int i=2;i<sqrt_n+1;i++)
    {
        if(n%i==0)
            return false;
    }
    return true;
}

然後是一般的素數篩選法

#include<cstdio>
#include <string>
#include<iostream>
#define MAX 900000
int prime[900010]={1,1};
using namespace std;
int main(){
   for(int i=2;i<MAX;i++)
   {
       if(!prime[i])
       {
           for(int j=i*2;j<=MAX;j+=i)
           {
               prime[j]=1;
           }
       }
   }
   int m,n;
   scanf("%d%d",&m,&n);
   int pi=0;
   for(int k=0;k<MAX;k++)
   {
       if(prime[k]==0)
       {
           pi++;
           if(pi>=m&&pi<n)
           {
               printf("%d",k);
               if((pi-m+1)%10==0)
                printf("\n");
               else
                printf(" ");
           }
           if(pi==n)
           {
                printf("%d",k);
           break;
           }

       }
   }
   return 0;

}

快速素數篩選法：
原理：對於每一個數i，乘上小於等於i的最小素因數的素數，就得到以i爲最大因數的合數。

#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/bjrxyz/article/details/8125913