最小点对

最小点对

给定平面上N个点，找出其中的一对点的距离，使得在这N个点的所有点的对中，该距离为所有点对中最小的。（2<=n<50000）

输入格式：

输入第一行一个数N，点的个数，以下N行，每行二个数，表示一个点的横坐标和纵坐标

 

输出格式：

一个数，最小点对的距离。

 

样例输入：

3
1 1
1 2
2 2

样例输出：

1.00

 

思路:

这道题就是 分治算法，先按照x排序 ，log（n）枚举其中的点 作为中间点，找在这个中间点两边距离最近的点 记为s，再在左边和右边找距离这个中间点为s的点，用一个数组记录（刚好形成了一个2d * d 的一个矩形），再以y排序，寻找距离最小点。

 

#include<bits/stdc++.h>
using namespace std;
const double inf=1e10;
struct node
{
	double x,y;
} a[1000010],b[1000010];
bool cmpx(node x,node y)
{
	return x.x<y.x;
}
bool cmpy(node x,node y)
{
	return x.y<y.y;
}
double work(node x,node y)
{
	return sqrt((x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y));
}
double sol(int x,int y)
{
	int mid,i,j,cnt=0;
	double s;
	if(x==y) return inf;
	mid=(x+y)/2;
	s=min(sol(x,mid),sol(mid+1,y));
	for(i=mid; i>=x&&a[mid].x-a[i].x<s; i--) b[++cnt]=a[i];
	for(i=mid+1; i<y&&a[i].x-a[mid].x<s; i++) b[++cnt]=a[i];
	sort(b+1,b+1+cnt,cmpy);
	for(i=1; i<=cnt; i++)
	{
		for(j=i+1; j<=cnt&&b[j].y-b[i].y<s; j++) s=min(s,work(b[i],b[j]));
	}
	return s;
}
int main()
{
	freopen("a.txt","r",stdin);
	int i,n,l,r;
	scanf("%d",&n);
	for(i=1; i<=n; i++) scanf("%lf%lf",&a[i].x,&a[i].y);
	sort(a+1,a+1+n,cmpx);
	double x=sol(1,n+1);
	printf("%.4lf",x);
	return 0;
}