几种筛选求素数的实现

1. #include <iostream>  
2. using namespace std;  
3. bool flag[10000];  
4. int prime[10000];  
5. int main()  
6. {  
7.     memset(flag,0,sizeof(flag));  
8.     int t=0;  
9.     for(int i=2;i<=10000;i++)  
10.     {  
11.         if(!flag[i]){  
12.             prime[t++]=i;  
13.             for(int j=2;i*j<=10000;j++)  
14.                 flag[i*j]=1;  
15.         }  
16.     }  
17.     for(int i=0;i<=100;i++)  
18.         cout<<prime[i]<<"\t";  
19.     cout<<endl;     
20. }
21. 另一种：
    1. 我推荐这个算法！ 易于理解。 只算奇数部分，时空效率都还不错！  
    2. half=SIZE/2;   
    3. int sn = (int) sqrt(SIZE);   
    4. for (i = 0; i < half; i++)   
    5.    p[i] = true;// 初始化全部奇数为素数。p[0]对应3，即p[i]对应2*i+3   
    6. for (i = 0; i < sn; i++) {      
    7. if(p[i])//如果 i+i+3 是素数  
    8. {       
    9.     for(k=i+i+3, j=k*i+k+i; j < half; j+=k)   
    10.     // 筛法起点是 p[i]所对应素数的平方 k^2                                          
    11.     // k^2在 p 中的位置是 k*i+k+i  
    12.     //    下标 i         k*i+k+i  
    13.     //对应数值 k=i+i+3   k^2           
    14.        p[j]=false;   
    15. }   
    16. }   
    17. //素数都存放在 p 数组中，p[i]=true代表 i+i+2 是素数。  
    18. //举例，3是素数，按3*3,3*5,3*7...的次序筛选，因为只保存奇数，所以不用删3*4，3*6....