埃拉托色尼筛法用于求一定范围内的素数。

第一种方法：

原理是在数组元素中存放数值，若判定为非素数则将其值修改为0，输出时会跳过值为0的元素。

 
  
  1. //埃拉托色尼筛法  
  2.   
  3. #include <iostream> 
  4. #include <iomanip>  
  5. #define NUM 100 
  6. #define W 2 
  7. using namespace std; 
  8.   
  9. int main(int argc,char * argv[]) 
  10. { 
  11.     int array[NUM/2]; 
  12.     int sign,flag; 
  13.   
  14.     //给数组赋值 
  15.     for(sign=0,flag=3;sign<NUM/2-1;sign++) 
  16.     { 
  17.         array[sign]=flag; 
  18.         flag+=2; 
  19.     }  
  20.   
  21.     //对数组进行筛选 
  22.     for(sign=0;sign<NUM/2-1;sign++) 
  23.     { 
  24.         if(array[sign]) 
  25.         { 
  26.             for(flag=sign+1;flag<49;flag++) 
  27.             { 
  28.                 if(!(array[flag]%array[sign])) 
  29.                 { 
  30.                     array[flag]=0; 
  31.                 } 
  32.             } 
  33.         } 
  34.     }  
  35.   
  36.     //输出素数组 
  37.     cout<<"100以内的素数有：\n 2\t"; 
  38.     flag=1; 
  39.     for(sign=0;sign<NUM/2-1;sign++) 
  40.     { 
  41.         if(array[sign]) 
  42.         { 
  43.             cout<<setw(W)<<array[sign]<<'\t'; 
  44.             flag+=1; 
  45.             if(!(flag%5)) 
  46.             { 
  47.                 cout<<endl; 
  48.             } 
  49.         } 
  50.     }  
  51.     cout<<endl; 
  52.   
  53.     return 0; 
  54. } 
 
  
 
 
  
代码中会反复出现 NUM/2-1 的原因：
 if sign<=50 最后判断的数字为 3+2*50=103 超出了上限100
 if sign<=49 最后判断的数字为 3+2*49=101 超出了上限100
 if sign<49（即sign<=48） 最后判断的数字为 3+2*48=99 在范围之内
 至于为什么是 NUM/2-1 ，在原理上想不通怎么解释，只是在运行过程中观察结果得出的结论。
 
  
另一种方法：
 
  
用数组的下标存放数值，数组中各元素的值用0和1表示是不是素数。最终输出值为1的数组元素的下标。
 原理是，发现一个数，然后在后面的所有数中删除它的倍数（赋值为0）。
 
  
 
  
  1. //埃拉托色尼筛法  
  2.   
  3. #include <iostream> 
  4. #include <iomanip> 
  5. #define N 100 
  6. #define W 2 
  7. using namespace std; 
  8.   
  9. int main(int argc,char *argv[]) 
  10. { 
  11.     int array[N]={0},sign,flag,i=0; 
  12.   
  13.     cout<<N<<"以内的素数有："<<endl;  
  14.     for(sign=2;sign<N;sign++) 
  15.     { 
  16.         if(array[sign]==0) 
  17.         { 
  18.             cout<<setw(W)<<sign<<" "; 
  19.             i++; 
  20.             //每输出5个值后换行  
  21.             if(i%5==0) 
  22.             { 
  23.                 cout<<endl; 
  24.             } 
  25.         } 
  26.         for(flag=1;sign*flag<=N;flag++) 
  27.         { 
  28.             array[sign*flag]=1; 
  29.         } 
  30.     } 
  31.   
  32.     return 0; 
  33. }