acm 1007 Quoit Design

这道题用到分治，以前只在数据结构上学到过，今天终于给碰上了。
一开始我用暴力，结果当然是超时。看了别人的代码，是用分治做的，我还是写不来，接着我就边理解边敲。里面还是有一些小细节要注意的。
1.注意循环的范围
2.注意优化
3.多用库函数

#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
double t;

struct point{
    double x,y;
}p[100001];

double power(double x){
    return x*x;
}

bool fun(point p1,point p2){
    if(p1.x<p2.x)
        return true;
    else if(p1.x==p2.x&&p1.y<p2.y)
        return true;
    else
        return false;
}

double min(double a,double b){
    return a<b?a:b;
}

double dis(point p1,point p2){
    return power(p1.x-p2.x)+power(p1.y-p2.y);
}

void mind(int l,int r){

    if(l+5>=r){

        for(int i=l;i<=r;i++){

            for(int j=i+1;j<=r;j++){
            //  printf("%lf\n",dis(p[i],p[j]));
                t=min(t,dis(p[i],p[j]));
            }

        }
        return; 
    }
    int mid=(l+r)/2;
    mind(l,mid);
    mind(mid+1,r);
    for(int i=l;i<=mid;i++){

        for(int j=mid+1;j<=r;j++){

            double d=dis(p[i],p[j]);
            if(d>t)
                break;
            t=d;
        }

    }

    return;
}

int main(){
    int n,i;

    while(scanf("%d",&n)!=EOF){ 
        t=100000000;
        if(n==0)
            break;
        for( i=0;i<n;i++)
            scanf("%lf%lf",&p[i].x,&p[i].y);
        sort(p,p+n,fun);

        mind(0,n-1);
        t=sqrt(t);
        printf("%.2lf\n",t/2);  
    }
    return 0;
}