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[100000+10];
struct point SCL[100000+10];
double find(int l,int r){
    int i,j,k,sn,mid=(l+r)/2;
    double res=1000000,lres,lx,rx,rres,tmp;
    if(r-l+1<=3){
        for(i=l;i<=r;++i){
            for(j=i+1;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+1,r);
    res=min(lres,rres);
    lx=P[mid].x-res;    rx=P[mid].x+res;
    sn=0;
    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
",P[i].x,P[i].y);
        }
    }
    sort(SCL,SCL+sn,cmpy);
    for(i=0;i<sn;++i){
        for(j=i+1;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=0;i<n;++i){
            scanf("%lf%lf",&P[i].x,&P[i].y);
        }
        sort(P,P+n,cmp);
        printf("%.2lf
",find(0,n-1)/2);
    }
    return 0;
}

　　然后是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[100000+10],R[100000+10];
struct point SCL[100000+10];
double find(int n, struct point P[]){
    int i,j,k,ln=0,rn=0,sn,mid=n/2;
    double res=1000000,lx,rx,lres,rres,tmp;
//    for(i=0;i<n;++i) printf("%lf %lf 
",P[i].x,P[i].y);
//    printf("
");
    if(n<=3){
        for(i=0;i<n;++i){
            for(j=i+1;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=0;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=0;
    for(i=0;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
",P[i].x,P[i].y);
        }
    }
    sort(SCL,SCL+sn,cmpy);
    for(i=0;i<sn;++i){
        for(j=i+1;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[100000+10];
    while(~scanf("%d",&n)&&n){
        for(i=0;i<n;++i){
            scanf("%lf%lf",&INIT[i].x,&INIT[i].y);
        }
        sort(INIT,INIT+n,cmp);
        printf("%.2lf
",find(n,INIT)/2);
    }
    return 0;
}