筛素数

ACM比赛中经常能遇到关于素数的题目，可能需要我们预处理求出100W，1000W以内的所有素数。笔者在这里介绍共三种筛素数的方法，代码如下：

#include <cstdio>  
#include <cstring>  
#include <algorithm>  
using namespace std;  

const int maxn=100000000;  
int prime[maxn];  
int primeNum;  
bool h[maxn];  

// 普通版  
void calPrime()  
{  
    for(int i=2;i<maxn;i++) if(!h[i])  
    {  
        prime[primeNum++]=i;  
        for(int j=i+i;j<maxn;j+=i)  
            h[j]=true;  
    }  
    printf("%d\n", primeNum);  
}  

// 优化版  
void calPrime2()  
{  
    for(int i=2;i<maxn;i++) if(!h[i])  
    {  
        prime[primeNum++]=i;  
        if(maxn>(long long)i*i) for(int j=i*i;j<maxn;j+=i)  
            h[j]=true;  
    }  
    printf("%d\n", primeNum);  
}  

// 欧拉筛法  
void calPrime3()  
{  
    for(int i=2;i<maxn;i++)  
    {  
        if(!h[i]) prime[primeNum++]=i;  
        for(int j=0;j<primeNum;j++)  
        {  
            if(i*prime[j]>maxn) break;  
            h[i*prime[j]]=true;  
            if(i%prime[j]==0) break;  
        }  
    }  
    printf("%d\n", primeNum);  
}  

int main()  
{  
    calPrime();  
//    calPrime2();  
//    calPrime3();  
}  

求1亿以内的所有素数，普通方法大概4.4S，优化版大概2.9S，欧拉筛法大概1.5S，优化程度还是很高的。
当然，只求100W以内的所有素数三种方法基本无差别。