HDU 1007 Quoit Design_二分法求最近点对

题目：~click here~

具体推导过程见王晓东《计算机算法设计与分析》第四版29页，一大堆图，公式，我就不抄上来了。

#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
int n;
struct node
{
  double x;
  double y;
}p[100005];
int a[100005];
double cmpx(node a,node b)
{
  return a.x<b.x;
}
double cmpy(int a,int b)
{
  return p[a].y<p[b].y;
}
double min(double a,double b)
{
  return a<b?a:b;
}
double dis(node a,node b)
{
  return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double find(int l,int r)
{
     if(r==l+1)
      return dis(p[l],p[r]);
     if(l+2==r)
      return min(dis(p[l],p[r]),min(dis(p[l],p[l+1]),dis(p[l+1],p[r])));
     int mid=(l+r)>>1;
     double ans=min(find(l,mid),find(mid+1,r));
     int i,j,cnt=0;
     for(i=l;i<=r;i++)
     {
       if( p[i].x>=p[mid].x-ans && p[i].x<=p[mid].x+ans)
          a[cnt++]=i;
     }
     sort(a,a+cnt,cmpy);
     for(i=0;i<cnt;i++)
     {
       for(j=i+1;j<cnt;j++)
       {
         if(p[a[j]].y-p[a[i]].y>=ans) break;
         ans=min(ans,dis(p[a[i]],p[a[j]]));
       }
     } 
     return ans;
}
int main()
{
    int i;
   
    while(scanf("%d",&n)!=EOF)
    {
      if(!n) break;
      for(i=0;i<n;i++)
       scanf("%lf %lf",&p[i].x,&p[i].y);
      sort(p,p+n,cmpx);
      printf("%.2lf%\n",find(0,n-1)/2);
    }
    return 0;
}