HDU1007 TLE代码和AC代码对比

这题卡了一天，上午开始看算法导论，然后实现了，一开始是wa，后来TLE，由于我开始的实现方式比较笨，而且在递归调用的时候很是混乱，用了好多数组。导致我的代码不断的出问题。具体是算法导论33-4.

后来改动了一点，也是看到别人代码改的，我把两个代码都写上，大家一看便知！

在实现分治的时候，我是把左右两个分组都复制了，然后传递的是一个数组，好慢啊。。。改动就是把参数传递改成传的是地址了，这样就不用复制结构体数组啦，不知道节省了多少力气，(⊙o⊙)…

首先是AC代码

 #include <iostream>
 #include <cstdio>
 #include <algorithm>
 #include <cmath>
 #include <string.h>
 #define  eps 0.000000001
 #define  S(x) (x)*(x)
 #define  aabs(x) (x)>0?(x):-(x)
 using namespace std;
 struct point{
     double x,y;
 };
 bool cmp(point a,point b){
     if( a.x==b.x) return a.y>b.y;
     return a.x<b.x;
 }
 bool cmpy(point a,point b){
     return (a.y>b.y);
 }
 struct point P[+];
 struct point SCL[+];
 double find(int l,int r){
     int i,j,k,sn,mid=(l+r)/;
     double res=,lres,lx,rx,rres,tmp;
     if(r-l+<=){
         for(i=l;i<=r;++i){
             for(j=i+;j<=r;++j){
                 tmp=sqrt( S(P[i].x-P[j].x) + S(P[i].y-P[j].y) );
                 if(res-tmp>eps) res=tmp;
             }
         }
         return res;
     }
     lres=find(l,mid);    rres=find(mid+,r);
     res=min(lres,rres);
     lx=P[mid].x-res;    rx=P[mid].x+res;
     sn=;
     for(i=l;i<=r;++i){
         if(P[i].x-lx>eps && eps<rx-P[i].x){
             SCL[sn].x=P[i].x;    SCL[sn].y=P[i].y;    sn++;
 //            printf("~~~~~~%lf,%lf\n",P[i].x,P[i].y);
         }
     }
     sort(SCL,SCL+sn,cmpy);
     for(i=;i<sn;++i){
         for(j=i+;j<sn;++j){
             if(SCL[i].y-SCL[j].y>res) break;
             if(res- sqrt( S(SCL[j].x-SCL[i].x) + S(SCL[j].y-SCL[i].y) ) > eps)
                 res = sqrt( S(SCL[j].x-SCL[i].x) + S(SCL[j].y-SCL[i].y) );
         }
     }
     return res;
 }

 int main(){
     int n;
     int i,j;
     while(~scanf("%d",&n)&&n){
         for(i=;i<n;++i){
             scanf("%lf%lf",&P[i].x,&P[i].y);
         }
         sort(P,P+n,cmp);
         printf("%.2lf\n",find(,n-)/);
     }
     return ;
 }

　　然后是TLE代码！

！！！！！！！！！！！！

 #include <iostream>
 #include <cstdio>
 #include <algorithm>
 #include <cmath>
 #include <string.h>
 #define  eps 0.000000001
 #define  S(x) (x)*(x)
 #define  aabs(x) (x)>0?(x):-(x)
 using namespace std;
 struct point{
     double x,y;
 };
 bool cmp(point a,point b){
     if( aabs(a.x-b.x) < eps) return (b.y-a.y)>eps;
     return (b.x-a.x)>eps;
 }
 bool cmpy(point a,point b){
     return (a.y>b.y);
 }
 struct point L[+],R[+];
 struct point SCL[+];
 double find(int n, struct point P[]){
     int i,j,k,ln=,rn=,sn,mid=n/;
     double res=,lx,rx,lres,rres,tmp;
 //    for(i=0;i<n;++i) printf("%lf %lf \n",P[i].x,P[i].y);
 //    printf("\n");
     if(n<=){
         for(i=;i<n;++i){
             for(j=i+;j<n;++j){
                 tmp=sqrt( S(P[i].x-P[j].x) + S(P[i].y-P[j].y) );
                 if(res-tmp>eps) res=tmp;
             }
         }
         return res;
     }
     for(i=;i<n;++i){                //init the L , R
         if(i<mid){
             L[ln].x=P[i].x;    L[ln].y=P[i].y;    ln++;
         }else{
             R[rn].x=P[i].x; R[rn].y=P[i].y; rn++;
         }
     }
     lres=find(ln,L);    rres=find(rn,R);
     res=min(lres,rres);
     lx=P[mid].x-res;    rx=P[mid].x+res;
     sn=;
     for(i=;i<n;++i){
         if(P[i].x-lx>eps && eps<rx-P[i].x){
             SCL[sn].x=P[i].x;    SCL[sn].y=P[i].y;    sn++;
 //            printf("~~~~~~%lf,%lf\n",P[i].x,P[i].y);
         }
     }
     sort(SCL,SCL+sn,cmpy);
     for(i=;i<sn;++i){
         for(j=i+;j<sn;++j){
             if(SCL[i].y-SCL[j].y>res) break;
             if(res- sqrt( S(SCL[j].x-SCL[i].x) + S(SCL[j].y-SCL[i].y) ) > eps)
                 res = sqrt( S(SCL[j].x-SCL[i].x) + S(SCL[j].y-SCL[i].y) );
         }
     }
     return res;
 }

 int main(){
     int n;
     int i,j;
     point INIT[+];
     while(~scanf("%d",&n)&&n){
         for(i=;i<n;++i){
             scanf("%lf%lf",&INIT[i].x,&INIT[i].y);
         }
         sort(INIT,INIT+n,cmp);
         printf("%.2lf\n",find(n,INIT)/);
     }
     return ;
 }