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;

}