C ++:分段错误(核心转储)问题

时间:2018-02-17 13:34:22

标签: c++ segmentation-fault

在C ++中实现Sieve of eratosthenes

当我运行我的C ++程序时,我得到了

  

“分段错误(核心转储)”

编译没有任何错误。

在这个程序中,我试图打印两个数字a和b之间的所有素数。

#include <iostream>
#include <string.h>
#define MAX 1000000

using namespace std;

// Print all primes s.t. a <= prime <= b
int main()
{
    int t; // no of test cases
    cin>>t;

    bool prime[MAX + 1];    // a[i] = true for i = prime
    long int count_primes_lte_me[MAX + 1]; // a[i] = Count ( primes ) <= i
    long int counter_of_visited_primes;

    prime[0] = prime[1] = false;
    for(int i = 2 ; i <= MAX ; i++)
    {
        if(prime[i] == true)
            count_primes_lte_me[i] = ++counter_of_visited_primes;

        for(int j = i*i ; j <= MAX ; j += i)  // sieve of eratosthenes
            prime[j] = false;
    }

    long int a , b;
    while(t--)
    {
        cin>>a>>b;
        cout<<count_primes_lte_me[b] - count_primes_lte_me[a - 1]<<endl;
    }

    return 0;
}

3 个答案:

答案 0 :(得分:4)

尝试使用MAX的较小值。在堆栈上分配如此庞大的数组会产生(初始)问题。然后用动态分配替换:

bool *prime = new bool[MAX+1];
if (prime==nullptr) // error

...

delete [] prime;

您也可以使用静态分配(将变量定义为全局)。

最好可能是使用一些合适的容器,例如bitset

小心 i*i可能会溢出,因此可能会出现其他问题......

答案 1 :(得分:1)

i * i会在i <= MAX的内循环中溢出。由于负j,这可能会导致段错误。使用更大的整数类型:long long j = (long long)i * i

为了正确起见,您需要将prime数组初始化为true,将计数数组初始化为0值。

答案 2 :(得分:1)

你最多只有一百万......我试图减少它并且没有错误但仍然无效......试试这个

/************************************
***Array names as pointers***********
************************************/
#include <iostream>
#include <iomanip>
using namespace std;

int main()
{
    const int MAX = 10000;
    long prime[MAX] = { 2, 3, 5 };
    long trial = 5;
    int count = 3;
    int found = 0;

    do
    {
        trial += 2;
        found = 0;
        for (int i = 0; i < count; i++)
        {
            found = (trial % *(prime + i)/* prime[i] */) == 0;

            if (found)
                break;
        }
        if (found == 0)
            *(prime + count++)/* prime[count++] */ = trial;
    } while (count < MAX);

    for (int i = 0; i < MAX; i++)
    {
        if (i % 5 == 0)
            cout << endl;
        cout << setw(10) << *(prime + i)/* prime[i] */;
    }
    cout << endl;
    return 0;
}