простые числа — Получение SIGABRT в c ++ в следующей программе

Я представил эту проблему в спой, и это показывает runtime error(SIGABRT), Он работает правильно на моей машине и на Ideone.com, но показывает ошибку там. какой-либо причине? Я пишу свой код здесь:
я пытаюсь вычислить простые числа порядка 10 ^ 8 и обработать некоторые операции над ними.
вот ссылка на проблему: http://www.spoj.com/problems/CPRIME

#include<iostream>
#include<vector>
#include<cmath>
#include<cstdio>

using namespace std;
int n=100000001;
int d = 10001;//3163;
vector<int>num(100000001,0);

int main (){
//atkin's seive for generating prime numbers of order 10^8
for (int x = 1; x <= d; x++){
for (int y = 1; y <= d; y++){
long long sqx=x*x,sqy=y*y;
long long m =(3*sqx)-sqy;
//m=m-sqx;
if(m<n){
if(x>y&&m%12==11)num[m]=num[m]^1;
m=m+2*sqy;
if(m<n){
if(m%12==7)num[m]=num[m]^1;
m=m+sqx;
if(m<n&&(m%12==1||m%12==5))num[m]=num[m]^1;
}
}

}
}
for (int i=5;i<d;i++){
if(num[i]){
for(int j=i*i;j<n;j+=i*i)num[j]=0;
}
}   //sieve finished here
//main code start from here and for loop is for counting number of primes
//less than or equal to that number.
int add=2;
num[2]=1;num[3]=2;
for(int a=4;a<n;a++){
if(num[a]){
//if(a==2||a==3||a==5)cout<<num[a]<<" "<<a<<" ";
add++;
num[a]=add;
}
else num[a]=add;
}
//cout<<num[2]<<" "<<num[3]<<endl;
while(1){
int t;
scanf("%d",&t);
if(t==0)break;
double ans=(num[t]-t*1.0/(log(t)))*100.0/num[t];
if(ans<0)ans*=-1.0;
printf("%.1lf\n",ans);
}
return 0;}

1

Решение

100000001*4 = 400000004 bytes
400000.004 kilo bytes
400.000004 mega bytes

Вам нужна половина концерта минимум! Поскольку вектор добавляет некоторые накладные расходы тоже.

ОС пойдет bannanas пытаясь поместить это в кучу.


Как примечание: вычисление простых чисел за O (n ^ 2) является наиболее наивным и медленным подходом. Рассмотреть возможность использования Сито Аткинса.


Убедитесь, что вы не делите на ноль. 100.0/num[t], Обычно деление на ноль генерирует SIGFPE но аппаратное обеспечение также может генерировать SIGABRT на такое разделение.

0

Другие решения

Вы должны выделить num как массив из кучи вместо стека
int* num=new int[size]; это будет работать нормально.

0

По вопросам рекламы [email protected]