算法竞赛宝典 分治算法 最接近点对问题(二维二分)

问题描述

给定平面上n个点，找其中一对点，使得在n个点组成的所有点对中，该点对间的距离最小。

二维的情形:

选取一垂直线l：x=m来作为分割直线。其中m为s中各点x坐标的中位数。 S1={p∈s∣x(p)≤m}   S2={p∈s∣x(p)>m}      ， S=s1∪s2

   递归地在S1和S2上解最接近点对问题，我们分别得到S1和S2中的最小距离δl和δ2。现设δ＝min(δl，δ2)。若S的最接近点对(p，q)之间的距离d(p，q)< δ，则p和q必分属于Sl和S2。

设p∈S1，q∈S2那么，p和q距直线l的距离均小于δ。因此，我们若用P1和P2分别表示直线l左边和右边的宽为δ的2个垂直长条,则p∈P1,q∈P2,   如下图所示

考虑P1中任意一点p，它若与P2中的点q 构成最接近点对的候选者，则必有d(p，g)< δ。这样的点一定落在一个δ×2δ的矩形R中，如图所示

由δ的意义可知,P2中任何两个S中的点的距离都不小于δ.由此矩形R中最多只有6 个S中的点.如图

证明:将矩形R的长为2d的边3等分，将它的长为d的边2等分，由此导出6个(d/2)×(2d/3)的矩形。若矩形R中有多于6个S中的点，则由鸽舍原理易知至少有一个(d/2)×(2d/3)的小矩形中有2个以上S中的点。设u，v是位于同一小矩形中的2个点，则

distance(u,v)<d。这与d的意义相矛盾。分治法合并时，检验6×n/2=3n对侯选者 ，在O(n)时间内完成合并。

#include<iostream>
#include<bits/stdc++.h>
#define min(a,b) ((a)>(b)?(b):(a))
#define ABS(a) (((a)>0)?(a):(-a))
const int N=100000+10;
using namespace std;
const int INF=0x3f3f3f3f;
int n;
struct node
{
	double x;
	double y;
}s[N];
bool cmp(node a,node b)
{
	if(a.x==b.x)
		return a.y<b.y;
	else
		return a.x<b.x;
}
double dis(node a,node b)
{
	double xx=(b.x-a.x);
	double yy=(b.y-a.y);
	return sqrt(xx*xx+yy*yy);
}
bool cmpy(int a,int b)
{
	return s[a].y<s[b].y;
}
double binary_search(int left,int right)
{
	if(left+1==right)
		return dis(s[left],s[right]);
	if(left>=right)
		return INF;
	double dl,dr,d;
	int mid;
	int tmp[1000];
	mid=(right-left)/2+left;
	dl=binary_search(left,mid);
	dr=binary_search(mid+1,right);
	d=min(dl,dr);
	int k=0;
	for(int i=left; i<=right; i++)
	{
		if(ABS(s[i].x-s[mid].x)<=d)
			tmp[k++]=i;
	}
	sort(tmp,tmp+k,cmpy);
	for(int i=0; i<k; i++)
	{
		for(int j=i+1; j<k&&(s[j].y-s[i].y)<=d; j++)
		{
			d=min(d,dis(s[tmp[i]],s[tmp[j]]));        //合并的时候考虑的需要多一点
		}        //不一定最接近的点在两边中的一点与mid上的点的距离，可能是两边之间的某一个距离
	}                //这里的d用来线性划分出可能出现的最近点对，通过扫描，可以得到存在的情况，进而划分，求出最小
	return d;
}
int main()
{
	while(cin>>n)
	{
		for(int i=1; i<=n; i++)
		{
			scanf("%lf%lf",&s[i].x,&s[i].y);
		}
		sort(s+1,s+n+1,cmp);
		printf("%.2lf\n",binary_search(1,n)/2.0);
	}
	return 0;
}