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; 
}