hdu 1007 /ZOJ 2107 Quoit Design （二分，最近点距离）

题目地址：hdu：点击打开链接

    ZOJ：点击打开链接

该代码在ZOJ能AC，但hdu就不行，不过可以作为一个经典案例分析：

#include<stdio.h>
#include<math.h>
#include <stdlib.h>

/*结构体尽量写成全局变量*/
struct point{
    double x,y;
}node[100005];

int tmp[100005];

double min(double a,double b){
  if(a<=b)
    return a;
  if(a>b)
    return b;
}
/*结构体的快排，注意返回值是int，如果直接两double相减会出错*/
int cmpx( const void *a , const void *b ){
    struct point *c = (struct point *)a;
    struct point *d = (struct point *)b;
    if(c->x != d->x) return c->x < d->x;
    else return d->y < c->y;
}
/*传的数组作为下标时可以引入中间变量，只要控制传回的-1，0，1正确即可*/
int cmp( const void *a ,const void *b){
    double  temp = ((struct point*)a)->y-((struct point*)b)->y;
      if(temp>0)
          return 1;
      else if(temp==0)
          return 0;
      else
          return -1;
}

double dist(struct point a,struct point b){
    double dis=(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
    return sqrt(dis);
}
double min_dist(int left,int right){
  int i,j;
  double d=1e100;
  /*分的只有一个点时肯定不能算进去，算作无穷大*/
  if(left==right)
    return d;
  /*两个点时传回距离*/
  else if(left+1==right)
    return dist(node[left],node[right]);
  int mid=(left+right)/2;
  d=min(min_dist(left,mid),min_dist(mid+1,right));

  /*搞完上面的才能搞下面的，也就是说多次递归后直到计算完d才进行下面的
  而计算中间矩形区域最小值是不需要递归的*/
  int pos=0;
  /*一个O(n)复杂度的循环真的不算什么，O(n)的循环尽管用就行*/
  for(i=left;i<=right;i++){
    if(fabs(node[i].x-node[mid].x)<=d)   /*把以mid为中线，左右两边距离为d的矩形里的点都搞进来*/
      tmp[pos++]=i;                      /*记录有几个点和其位置*/
  }
  qsort(tmp,pos,sizeof(tmp[0]),cmp);     /*把记录的这几个数按y排了*/
  for(i=0;i<pos;i++)
  for(j=i+1;j<pos&&node[tmp[j]].y-node[tmp[i]].y<=d;j++){  /*把y轴上的距离判断一下，主要是若单看y上的投影都比d大，距离肯定更大了*/
    double dd=dist(node[tmp[i]],node[tmp[j]]);            /*其实每个点只要和对面判断6次就行了。符合条件比较下距离*/
    d=min(d,dd);
  }
  return d;
}

int main(){
  int i,j,n,a,b;
  while(scanf("%d",&n)&&n){
    for(i=0;i<n;i++)
      scanf("%lf%lf",&node[i].x,&node[i].y);
    /*超时的原因可能是快排最坏时间复杂度太高*/
    qsort(node,n,sizeof(node[0]),cmpx);
    double ans=min_dist(0,n-1)/2.0;
    printf("%.2lf\n",ans);

  }
    return 0;
}