在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;
}
答案 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;
}