/*
用筛选法求100以内的素数
所谓“筛选法”指的是“埃拉托色尼(Eratosthenes)筛法”。他是古希腊的著名数学家。他采取的方法是,在一张纸上写上1到100全部整数,然后逐个判断它们是否是素数,找出一个非素数,就把它挖掉,最后剩下的就是素数。
具体做法如下:
<1> 先将1挖掉(因为1不是素数)。
<2> 用2去除它后面的各个数,把能被2整除的数挖掉,即把2的倍数挖掉。
<3> 用3去除它后面的各数,把3的倍数挖掉。
<4> 分别用4、5…各数作为除数去除这些数以后的各数。这个过程一直进行到在除数后面的数已全被挖掉为止。例如找1~50的素数,要一直进行到除数为47为止(事实上,可以简化,如果需要找1~n范围内素数表,只需进行到除数为n^2(根号n),取其整数即可。例如对1~50,只需进行到将50^2作为除数即可。)
如上算法可表示为:
<1> 挖去1;
<2> 用刚才被挖去的数的下一个数p去除p后面各数,把p的倍数挖掉;
<3> 检查p是否小于n^2的整数部分(如果n=1000, 则检查p<31?),如果是,则返回(2)继续执行,否则就结束;
<4> 纸上剩下的数就是素数。
定义1个数组a, a[1]~a[100] 分别代表1~100这100个数。
如果检查出数组a的某一元素的值不是素数,就将其值设为0,最后剩下不为0的就是素数。
*/
#include <iostream>
#include <iomanip>
using namespace std;
#include <math.h>
int main()
{int i,j,n,a[101]; //定义a数组包含101个元素
for (i=1;i<=100;i++) //a[0]不用,只用a[1]~a[100]
a[i]=i; //使a[1]~a[100]得值分别为1~100
a[1]=0; //先“挖掉”a[1]
for (i=2;i<sqrt(100);i++) //i表示被除数下标,从2~sqrt(n) n=100
for (j=i+1;j<=100;j++) //j表示除数下标,从i后面的i+1到100
{if(a[i]!=0 && a[j]!=0)
if (a[j]%a[i]==0) //如果整除,表示a[j]不是素数,则挖掉它
a[j]=0; }
cout<<endl;
for (i=1,n=0;i<=100;i++)
{if (a[i]!=0) //没有被挖掉的数,即值不为0的数,为素数
{cout<<setw(5)<<a[i]<<" ";
n++;} //累计素数个数
if(n==10) //输出10个数后换行
{cout<<endl;
n=0;}
}
cout<<endl;
return 0;
}