埃拉托色尼筛选法 C++实现

在公元前3世纪，古希腊天文学家埃拉托色尼发现了一种找出不大于n的所有自然数中的素数的算法，即埃拉托色尼筛选法。

具体筛选步骤：

• 这种算法首先需要按顺序写出2到n中所有的数。
• 然后把第一个元素画圈，表示它是素数，然后依次对后续元素进行如下操作：

• 如果后面的元素是画圈元素的倍数，就画X，表示该数不是素数。
• 在执行完第一步后，会得到素数2，而所有是2的倍数的数将全被画掉，因为他们肯定不是素数。
• 接下来，只需要重复上述操作，把第一个既没有被圈又没有画X的元素圈起来，然后把后续的是它的倍数的数全部画X。
• 这次操作将得到素数3，而所有是3的倍数的数都被去掉。
• 以此类推，最后数组中所有的元素不是画圈就是画X。

• 所有被圈起来的元素均是素数，而所有画X的元素均是合数。

给定一个数字n, 2≤n≤2000000, 编写一个程序实现埃拉托色尼筛选法找出小于或等于n的所有素数。

最基本的实现方法：

#include <iostream>
//在外部创建全局变量解决数组过大的问题
const int m = 2000000;
using namespace std;
//给第一个赋值为2，初始化质数列表
int all_number[m] = {2, 0};

int main() {
    int n, i, j;
    cout << "input an integer: " << endl;
    cin >> n;
    //生成2~n的所有整数
    for (i = 0; i < n; ++i) all_number[i] = i + 2;
    //筛选质数
    for (i = 1; i <= n; ++i) {
        for (j = 1; j <= i; ++j) {
        	//避免出现除以前面标记为0的错误
            if (all_number[j - 1] == 0)continue;
            else {
            	//判断是否可以整除前面的数
                if (all_number[i] % all_number[j - 1] == 0) {
                        //用零替换被筛掉的数
                        all_number[i] = 0;
                        break;
                }
            }
        }
    }
    //输出剩下的质数（非零）
    for (i = 0; i <= n; ++i) {
        if (all_number[i] != 0) {
            cout << all_number[i] << " ";
        }
    }
    return 0;
}

小结：

• 用零标记
• 避免标记后再次参与运算，从而出现除以0的错误