HDU1007 平面最近点对

只判六个点的抽屉原理好像不太懂，，打了个类似暴力的东西但是听说复杂度可证？反正也过去了。。

先将区间分成左右两个子区间，递归的计算这两个子区间内的点对

的最近距离，记这两个最小值中较小的一个为res。

然后计算跨越这两个子区间分界线的点对的最近距离，显然，枚举的

时候不需要枚举坐标之差的绝对值大于res的点对，同理，y坐标之差的绝

对值大于res的点对也不需要进行枚举。

因此，我们在递归的同时用归并排序的方法将区间[l,r]内的点按y坐标

从小到大排序，将x坐标不在[mid−res,mid+ees]范围内的点剔除(mid为x坐

标分界线)，枚举剩下的y坐标差的绝对值不大于res的点对，更新答案。

以上为st的讲课题解

这道题简单代码粗暴写得很快

//QWsin
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const double INF=1e10;
const int maxn=100000+10;

struct point{
double x,y;int t;
inline void input(){scanf("%lf%lf",&x,&y);}
inline bool operator < (const point &rhs)const{
return y<rhs.y||(y==rhs.y&&x<rhs.x);
}
}a[maxn],t[maxn];

inline int cmp(point a,point b){return a.x<b.x;}

inline double dis(point a,point b)
{
double dx=a.x-b.x,dy=a.y-b.y;
return sqrt(dx*dx+dy*dy);
}

int n;

double solve(int l,int r)
{
if(l==r) return INF;
if(l+1==r) return dis(a[l],a[r]);

int mid=(l+r)>>1;
double d=min(solve(l,mid),solve(mid+1,r));

int k=0;
for(int i=l;i<=r;++i)
if(fabs(a[i].x-a[mid].x) <= d) t[++k]=a[i];

sort(t+1,t+k+1);

for(int i=1;i<=k;++i)
for(int j=i+1;j<=k&&t[j].y-t[i].y <= d;++j)
d=min(d,dis(t[i],t[j]));

return d;
}

inline void solve()
{
for(int i=1;i<=n;++i) a[i].input();
sort(a+1,a+n+1,cmp);
printf("%.2f\n",solve(1,n)/2);
}

int main()
{
while(scanf("%d",&n)==1&&n) solve();
return 0;
}