埃拉托色尼的经典素数计算方法

看到埃拉托色尼的经典素数计算方法：

第一步，列出如下这样以2开头的序列：

2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

第二步，标出序列中的第一个素数，主序列变成：(其实什么都没变，我们知道第一个素数是2)

2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

第三步，将剩下的序列中，第二项开始每隔一项划掉(2的倍数，用红色标出)，主序列变成：

2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

第四步，如果现在剩下的序列中，最大数小于第一个素数的平方，那么剩下的序列全部是素数，否则返回第二步。本例中，因为25大于2的平方，我们返回第二步：

         剩下的序列中第一个素数是3 (大概前面2已经被认定为素数了，不记在其中), 接下来划出3的倍数。

2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

         然后再比较3的平方跟25的大小。然后再确定5，再划出5的倍数。一直做下去。

在本例子中，5的平方等于25了，所以跳出循环，否则继续划出7的倍数，11的倍数，一直做下去。

 

综上所述，就是在一个序列中，一次划出2,3,5,7……，等素数的倍数，直到整个序列全部都是素数，什么时候停止，即某个素数的平方等于最后一个数。

 

C程序实现，在这里利用一个数组来表示，该位置1表示素数，0表示合数。

 

#define N 10000

main()

{

         int i ,j,a[N];

         for(i=2;i<N;i++) a[i]=1;

         for(i=2;i<N;i++)

            if(a[i])

               for(j=i;i*j<n;j++) a[j*i]=0;

         for(i=2;i<N;j++)

            if(a[i]) printf("M",i);

 

}