此题来自hdu acm网站１００７题，http://acm.hdu.edu.cn/showproblem.php?pid=1007，提交结果是ＲＥ，题目中的几个例子输入，运行结果都正确，但输入有些数据，程序就卡主没法运行了，不知道代码错哪了！
下面是代码：

#include #include #include #define M 1000 using namespace std; void sorted(double**,int); //对点的X轴系数进行排序 double get_min(double ,double ); double cal_radius(double** ,int ,int ); //计算圆环的最大半径 int main() { int N[M],n; double** pos; double radius[M]; int i=1; while(1) { cin>>N[i]; if(N[i]==0) break; if(N[i]<2||N[i]>100000) return 0; pos=new double*[N[i]]; for(n=0;n>pos[n][0]>>pos[n][1]; } sorted(pos,n); //对x系数进行排序 radius[i]=cal_radius(pos,0,n-1)/2.00; //计算所有点之间的最小距离 i++; delete []pos; //释放内存 } for(int j=1;j=0;j--) { pos[j+1][0]=pos[j][0]; pos[j+1][1]=pos[j][1]; } pos[j+1][0]=tmp[0]; pos[j+1][1]=tmp[1]; } } } double get_min(double d1,double d2) { return d1<=d2?d1:d2; } double cal_radius(double** pos,int pre,int tail) //计算所有点之间的最小距离，并返回最小距离 { int left; //left表示在中界线左边的点的个数 int n=tail-pre+1; //二维数组长度 double mid_line; //中界线 double d1,d2,dmin,tmp_d; double** tmp; tmp=new double*[n]; if(n==2) { dmin=sqrt(pow(pos[pre][0]-pos[tail][0],2)+pow(pos[pre][1]-pos[tail][1],2)); return dmin; } if(n==3) { double t1,t2,t3; t1=sqrt(pow(pos[pre][0]-pos[pre+1][0],2)+pow(pos[pre][1]-pos[pre+1][1],2)); t2=sqrt(pow(pos[pre][0]-pos[pre+2][0],2)+pow(pos[pre][1]-pos[pre+2][1],2)); t3=sqrt(pow(pos[pre+1][0]-pos[pre+2][0],2)+pow(pos[pre+1][1]-pos[pre+2][1],2)); dmin=get_min(get_min(t1,t2),t3); return dmin; } mid_line=pos[n/2+pre][0]; if(n%2==0) { d1=cal_radius(pos,pre,n/2+pre-1); d2=cal_radius(pos,n/2+pre,tail); } else { d1=cal_radius(pos,pre,n/2+pre); d2=cal_radius(pos,n/2+pre+1,tail); } dmin=get_min(d1,d2); for(int i=n/2+pre,j=0;pos[i][0]>(mid_line-dmin);i--,j++) { tmp[j]=new double[2]; tmp[j][0]=pos[i][0]; tmp[j][1]=pos[i][1]; } left=j; left--; for(int k=n/2+pre+1;pos[k][0]<(mid_line+dmin);k++,j++) { tmp[j]=new double[2]; tmp[j][0]=pos[k][0]; tmp[j][1]=pos[k][1]; } j--; for(int h=0;h=left;j--) { if(tmp[j][1]-tmp[h][1]>-dmin&&tmp[j][1]-tmp[h][1]