埃氏筛法全名:埃拉托斯特尼筛法
以下是如何编程实现该算法:
算法的具体实现说明如下:
(1)定义一个数组a,每一个元素的下标对应一个整数,它的值表示这个整数是否为素数,值1表示是素数,值0表示不是素数;例如:a[12] = 0,表示在数组中存的数为12,a[12]的值表示12是不是素数,显然12不是素数,故a[12]中存的是0;
(2)下面输入你要求的上边界n,即你要从2筛到哪个数;
(3)在还没有使用埃氏筛法之前,我们显然要先默认每一个a[i]是素数,即令每一个a[i] = 1;这样我们才能操作后面的算法,筛去不是素数的数;
(4)选定第一个整数2,从它对应的数组元素a[2]的下一个元素a[3]开始判断数组元素的下标值是否能被2整除,若可以,则把这个数组元素的值改为0,即a[i] = 0,这个数 i 不是素数;如果不可以则不改变a[i]的值,即还是1;这样一直判断到n, 即元素a[n]为止;
(5)再选定下一个比2大,且a[i]值为1的数组元素,它的下标 i 就是比2大的第一个素数,从它对应的数组元素a[i]的下一个元素a[i + 1]开始判断数组元素的下标值是否能被 i 整除,若可以,则把这个数组元素的值改为0,这个数不是素数;如果不可以,则不改变这个数组元素的值,即还是1;这样一直判断到n, 即元素a[n]为止;
(6)再继续向后:选定,判断,筛去三个重复步骤,直到本次选定的数是大于n的算数平方根即可结束;
经过以上步骤:就可以得到一个a[]数组,数组中的每一个元素值为0或1,其中0表示这个元素的下标值不是素数,而1表示这个元素的下标值是素数;
以下给出模板代码:
#include <iostream>
#include <math.h>
#define N 10000
using namespace std;
int a[N];
void Sieve(int n)
{
int i;
for(i = 0; i <= n; i++)
{
a[i] = 1;
}
a[0] = 0;
a[1] = 0;
for(i = 2; i <= sqrt(n); i++)
{
if(a[i] == 1)
{
for(int j = i*i; j <= n; j = j + i)
{
a[j] = 0;
}
}
}
}
int main()
{
int n;
cin >> n;
Sieve(n);
for(int i = 2; i <= n; i++)
{
if(a[i] == 1)
{
cout << i << endl;
}
}
return 0;
}