hdoj 1007 Quoit Design(分治法)

题意:给你n个点，让你求出n个点里距离最近的两点距离。

思路:分治法先把点按照坐标x进行排序，划分为左右两部分，

最近点对有三种情况，都在左半部分，都在右半部分，两点分别在左边和右边

递归求出都在左边和都在右边的情况，选一个最小值m1。

然后求出两边都有的情况。分别在两边的情况可以对点进行进一步的筛选，将点的下标

保存在t[]，对筛选出的点按照y坐标进行排序，方便求出最小值。

代码如下:

#include<stdio.h>#include<iostream>#include<algorithm>#include<math.h>using namespace std;struct point{double x,y;};point p[100005];int t[100005]; double min(double a,double b){return (a>b)?b:a;}double dis(point a,point b){return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}bool cmp(const point &a,const point &b){return a.x<b.x;}bool cmp1(int a,int b){return p[a].y<p[b].y;}double spit(int a,int b){int i,j,c;double m1;if(a+1==b)return dis(p[a],p[b]);if(a+2==b)return min(dis(p[a],p[a+1]),min(dis(p[a],p[b]),dis(p[b],p[a+1])));    c=(a+b)>>1;      m1=min(spit(a,c),spit(c+1,b));  int k=0;for(i=a;i<=b;i++) {if(p[i].x>=p[c].x-m1&&p[i].x<=p[c].x+m1)        t[k++]=i;} sort(t,t+k,cmp1);    for(i=0;i<k;i++)          for(j=i+1;j<k;j++){                                   if(p[t[j]].y-p[t[i]].y>m1) break;              m1=min(m1,dis(p[t[i]],p[t[j]]));        }      return m1; }int main(){      int n;      int i;      double r;      while(scanf("%d",&n)){          if(n==0) break;          for(i=0;i<n;i++) scanf("%lf %lf",&p[i].x,&p[i].y);          sort(p,p+n,cmp);          r=spit(0,n-1);          r=r/2;          printf("%.2lf\n",r);      }      return 0;  }