素数筛法

平常的素数求法太慢了，所以需要加一下速。

第一种，筛法的最基本做法，但是效率不是很高：

bool prime[MAXN+5];
int ans[MAXN+5],cnt;
void getPrimeSlower(int n){
    memset(prime,1,sizeof(prime));
    cnt=0, prime[0]=prime[1]=0;
    for(int i=2;i<=n;i++)
    {
        if(prime[i])
        {
            ans[cnt++]=i;
            for(int j=i*2;j<=n;j+=i)
            {
                prime[j]=0;
            }
        }
    }
}

第二种，第一种方法做了一点优化，只需要判断n^(1/2)个数即可，时间复杂度大概为O(NlogN)：

bool prime[MAXN+5];
int ans[MAXN+5],cnt;
void getPrime(int n){
    memset(prime,1,sizeof(prime));
    cnt=0, prime[0]=prime[1]=0;
    for(int i=2;i<=n;i++)
    {
        if(prime[i])
        {
            if(n/i<i) break;
            for(int j=i*i;j<=n;j+=i)
            {
                prime[j]=0;
            }
        }
    }
    for(int i=2;i<=n;i++)
        if(prime[i])    
            ans[cnt++]=i;
}

第三种，比较快的一种，思路是用已经筛出的素数进行筛选，时间复杂度O(N)：

bool prime[MAXN+5];
int ans[MAXN+5],cnt;
void getPrimeFaster(int n){
    memset(prime,1,sizeof(prime));
    cnt=0, prime[0]=prime[1]=0;
    for(int i=2;i<=n;i++)
    {
        if(prime[i])
            ans[cnt++]=i;
        for(int j=0;(j<cnt)&&(i*ans[j]<=n);j++)
        {
            prime[i*ans[j]]=0;
            if(i%ans[j]==0)
                break;
        }
    }
}

其实还有一种很恶心的筛法，此处不列举了。