ZOJ 2107

来源：互联网发布：汉字笔顺演示软件编辑：程序博客网时间：2024/06/16 00:01

平面最近点对问题。。

这题以前看过一篇文章。貌似只用比较最近的6个点左右（X轴排序过的）就可以了。。

A B C G

D E F H

假设这8个点在一个矩形上。。可以知道平面最近点对。。就是X轴周围的12个点中。。 O(n)算法就能过掉了。。

#include <queue>  #include <stack>  #include <math.h>  #include <stdio.h>  #include <stdlib.h>  #include <iostream>  #include <limits.h>  #include <string.h>  #include <algorithm>  using namespace std;  const int MAX = 100010;  const double D_MAX = 1e50;  const double eps = 1e-15;  struct point{      double x,y;  };  point p[MAX];bool xy(double x,double y) // x < y   {      return x < y - eps;  }  bool dd(double x,double y)  // x == y   {      return fabs( x - y ) < eps;  }  bool cmp(point a,point b)  {      if( dd(a.x,b.x) )          return xy(a.y,b.y);      return xy(a.x,b.x);  }  double disp(point a,point b)  {      return sqrt( (a.x - b.x)*(a.x - b.x) + (a.y - b.y)*(a.y - b.y) );  }  double nearpot(int l,int r)  {  double ans = D_MAX;for(int i = l;i <= r;i++){for(int j = 1; i+j <= r && j <= 6;j++){ans = min(ans,disp(p[i],p[j+i]));}}return ans;}int main(){int n;      while( ~scanf("%d",&n) && n )      {          for(int i=0; i<n; i++)              scanf("%lf %lf",&p[i].x,&p[i].y);             sort(p,p+n,cmp);          double ans = nearpot(0,n-1);          printf("%.2lf\n",ans/2.0);      }  return 0;  }

0 0