1

この問題を spoj に提出しましたが、これは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;


}
4

2 に答える 2

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

最低でもハーフ ギグが必要です。ベクトルはオーバーヘッドも追加するためです。

OSは、それをヒープに収めようと大騒ぎします。

補足として: O(n^2) 時間で素数を計算することは、最も素朴で遅いアプローチです。アトキンスふるいの使用を検討してください 。

ゼロで割っていないことを確認してください。100.0/num[t]. 通常、ゼロによる除算が生成されますSIGFPEが、ハードウェアもSIGABRTそのような除算で生成できます。

于 2013-08-15T09:39:32.337 に答える
0

スタックの代わりにヒープから配列として num を割り当てる必要があります int* num=new int[size];。これは正常に機能します。

于 2013-10-29T13:46:12.960 に答える