质数超过 600 万

Question

我在 Hackerrank 中解决了一个问题，问题是找到一个范围内的素数计数。由于使用常规方法面临超时，我使用了埃拉托色尼筛。除了两个隐藏的测试用例外，大多数测试用例都有效。我在 GDB 编译器中运行代码，发现代码只支持高达 600 万的值。我该怎么办？代码如下：

#include<cstring>
#include<cstdio>
#include <iostream>
#include <algorithm>
using namespace std;


void SieveOfEratosthenes(unsigned long long int a,unsigned long long int b) 
{ 
    unsigned long long int count=0; 
    bool prime[b+1]; 
    memset(prime, true, sizeof(prime)); 
  
    for (unsigned long long int p=2; p*p<=b; p++) 
    { 
        // If prime[p] is not changed, then it is a prime 
        if (prime[p] == true) 
        { 
            for (unsigned long long int i=p*p; i<=b; i += p) 
                prime[i] = false; 
        } 
    } 
  
    for (unsigned long long int p=a; p<b; p++) 
       if (prime[p] &&p!=1) 
           count++;
    cout<<count;
          
} 
  
int main() 
{ 
    unsigned long long int a,b;
    cin>>a>>b;
    SieveOfEratosthenes(a,b); 
    return 0; 
} 

Answer 1

看起来像经典的堆栈溢出。 bool prime[b+1]; 分配在堆栈上，您已达到限制。

如果它在 Linux 上运行，则允许的最大堆栈大小通常总计约为 8MB 或更小，因此很有可能您刚刚超过了这个值。

将其移出堆栈，或执行位打包而不是完整的bool，它应该会再次正常工作。

Answer 2

您正在函数中创建一个 bool 数组，该数组将存储在堆栈中。在 Windows 上，堆栈的典型最大大小为 1MB，而在典型的现代 Linux 上为 8MB。您正在创建一个包含 600 万条记录的数组，将近 6MB。

要解决这个问题，您可以创建一个vector 而不是数组，这将是堆中的stored。