如何使用Sieve打印素数直到10 ^ 8?

Fahim Hassan:

我最近学习了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;
}
cigien:

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

而是使用a vector将数据存储在堆上,如下所示:

vector<bool> isPrime(MAX+1);

这是一个演示

另外,请注意,由于要在索引处进行索引MAX,因此会有一个错误,因此向量应为size MAX+1

此外,还应避免使用using namespace std;以及像一样的typedef ll,它们会使代码更难阅读。

本文收集自互联网,转载请注明来源。

如有侵权,请联系 [email protected] 删除。

编辑于
0

我来说两句

0 条评论
登录 后参与评论

相关文章