poj3714 最近点对

最近点对，采用分治方法。过程：
1对原数组依据x左标从小到大排序。
2二分数组，左边求出最小值，右边求出最小值，我们求最小的。
3找出对于左右两边的可能小于当前最小值的最近点对，更新最小值。
这题目需要区分一下点，让我们求的是闪兵到任意一个核电站的最短距离，加一个标志就可以了。
代码：

#include <iostream>
#include <stdio.h>
#include <iomanip>
#include <algorithm>
#include <cmath>
using namespace std;
const int N=200000+1;
const double INF=1e100;
struct Point{
	double x,y;
	bool flag;
};
Point m[N+1];
int tmp[N+1];
double dis(Point p1,Point p2)
{
	return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
double minL(double l1,double l2)
{
	return l1<=l2?l1:l2;
}
bool cmpy(int a,int b)
{
	return m[a].y<m[b].y;
}
bool cmpx(Point a,Point b)
{
	return a.x<b.x;
}
double getMinLen(Point *s,int left,int right)
{
	double rs=INF;
	if(left == right)
		return rs;
	if(left+1 == right)
	{
		if(s[left].flag == s[right].flag)
			return rs;
		return dis(s[left],s[right]);
	}
	int mid=(left+right)>>1;
	rs=getMinLen(s,left,mid);
	rs=minL(rs,getMinLen(s,mid+1,right));//rs现在是两边中最点值了
	int i,j,num=0;
	for(i=left;i<=right;++i)
	{
		if(fabs(s[i].x-s[mid].x)<=rs)//过滤掉不可能是最近点对的点
			tmp[num++]=i;
	}
	sort(tmp,tmp+num,cmpy);//按y坐标排序，减少比较次数
	double d=INF;
	for(i=0;i<num;++i)
	{
		for(j=i+1;j<num;++j)
		{
			if(fabs(s[tmp[i]].y-s[tmp[j]].y)>=rs)//i，j两点的y坐标已经大于rs,不需要在遍历了
				break;
			if(s[tmp[i]].flag!=s[tmp[j]].flag && (d=dis(s[tmp[i]],s[tmp[j]]))<rs)
				rs=d;
		}
	}
//	cout<<left<<" "<<right<<" "<<rs<<endl;
	return rs;
}
int main()
{
	int t,n,i;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d",&n);
		for(i=0;i<n;++i)
		{
			scanf("%lf%lf",&m[i].x,&m[i].y);
			m[i].flag=0;
		}
		for(i=0;i<n;++i)
		{
			scanf("%lf%lf",&m[i+n].x,&m[i+n].y);
			m[i+n].flag=1;
		}
		n<<=1;
		sort(m,m+n,cmpx);
		double ans=getMinLen(m,0,n-1);
		cout<<setiosflags(ios::fixed)<<setprecision(3)<<ans<<endl;
	}
	return 0;
}