埃氏筛法

埃氏筛法全名：埃拉托斯特尼筛法

以下是如何编程实现该算法：
算法的具体实现说明如下：

（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;
}