筛素数

一：埃拉托斯尼斯筛法.

#define maxn 1000000
bool isPrime[maxn+10]={true};
int prime[maxn],k=0;
void selPrime()
{
       isPrime[0]=isPrime[1]=false;
       for(int i=2;i<=maxn;i++)
       {
               if(isPrime[i])
                {
                       prime[k++]=i;
                       for(int j=i*i;j<=maxn;j+=i)
                       isPrime[j]=false;
                 }
        }
}

二：6N(+-)1法.

#define maxn 1000
int prime[maxn+10],k=0;
bool isPrime(int k)
{
    if(k%2==0)return false;
    for(int i=3;i*i<=k;i+=2)
    if(!(k%i))return false;
    return true;
}
void doPrime()
{
    prime[k++]=2;
    prime[k++]=3;
    for(int i=3;i<=maxn;i+=3)
    for(int j=0;j<2;j++)
    if(isPrime(2*(i+j)-1))
    prime[k++]=2*(i+j)-1;
}