如何使用 Sieve 打印 10^8 之前的质数？

Question

我最近学习了Sieve算法并开始使用它来学习如何在问题中使用该算法。我已经正确地编写了代码，因为我在其中找不到任何错误，但它在没有显示任何输出的情况下关闭了。找不到问题所在。帮助将不胜感激。

#include <iostream>
#include <vector>
//#define MAX 10000
typedef long long int ll;
using namespace std;
vector <ll> primes;
void sieve(){
    ll MAX = 100000000;
    bool isPrime [MAX];
    for(ll i = 0;i < MAX; ++i)isPrime[i] = true;
    //isPrime[0] = isPrime[1] = false;
    
    for(ll i=3; i*i <= MAX; i += 2){
        if(isPrime[i]){
            for(ll j = i*i; j <= MAX; j += i){
                isPrime[j] = false;
            }
        }
    }
    primes.push_back(2);
    for(ll i = 3; i <= MAX; i += 2){
        if(isPrime[i]){
            primes.push_back(i);
        }
    }
    for(ll i = 0; i <= 10; ++i){
        cout<<primes[i]<<endl;
    }
}
int main(){
    sieve();
    return 0;
}

Answer 1

您正在创建一个大小为 10^8 的静态数组，它存储在堆栈中。这对于堆栈来说太大了，可能会导致堆栈溢出。

相反，使用 vector 将数据存储在堆上，如下所示：

vector<bool> isPrime(MAX+1);

这里是 demo。

此外，请注意您有一个错误，因为您是在索引 MAX 处建立索引，所以向量的大小应为 MAX+1.

此外，您应该避免 using namespace std;，以及像 ll 这样的类型定义，它们会使代码更难阅读。

如何使用 Sieve 打印 10^8 之前的质数？

How can I print prime numbers till 10^8 using Sieve?

c++

algorithm

numbers

sieve-of-eratosthenes