平面几何最近两点（分治+sort）

A - Raid POJ - 3714

思路：

1. 本题是问station和agent之间最近的两个，不妨先想想如果是x-y平面内最近的两点怎么算呢。

对于平面最近的两点

1. 先按x坐标升序排，（如果x相同，那么y小的排到前面）
2. 之后就开始分治了**divide（l，r）【l，r】**区间内的最近点对。
3. 问题可以再缩小，考虑d1= 【l，mid】内的最近点对
   和d2=【mid+1，r】内的最近点对
4. 由于只考虑了上面左右区间最接近的点对，是不严谨的·。
   因为最近点对，可能是一个在【l，mid】，而另一个在【mid+1，r】。
5. 所以要枚举找到【l，r】内满足条件的点对，之后暴力求d3.（这里要智慧暴力）
6. 返回就是min（d1，d2，d3）；

那么解决了平面最近点对，回到这题，

其实就把station和agent标号，当两者的标号相异，就算距离，否则就返回INF

反思

1. double 的有效位是15-16，数值范围可以到1e308
2. 这次又没开longlong见了一回祖宗。qwq。

AC

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <vector>
#define For(i,x,y) for(register int i=(x); i<=(y); i++)
#define mst(x,a) memset(x,a,sizeof(x))
using namespace std;
typedef long long ll;
const int maxn=2e5+10;
const double INF=1e30;
struct point
{
    ll x,y;
    int sta;
}p[maxn];
vector<point>v;
bool cmpx(point a, point b)
{
    if(a.x!=b.x)return a.x<b.x;
    else return a.y<b.y;
}
bool cmpy(point a, point b)
{
    return a.y<b.y;
}
double dis(point a, point b)
{
    double tmp=(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
    return sqrt(tmp);
}
double divd(int l, int r)
{
    if(l==r)return INF;
    if(l+1==r)
    {
        if(p[l].sta!=p[r].sta)
            return dis(p[l],p[r]);
        else return INF;
    }
    int mid=l+(r-l)/2;
    double d1=divd(l,mid);
    double d2=divd(mid+1,r);
    double d=min(d1,d2);
    For(i,l,r)
    {
        if(fabs(p[mid].x-p[i].x)<=d)v.push_back(p[i]);
    }
    sort(v.begin(),v.end(),cmpy);
    int len=v.size();
    For(i,1,len)
    {
        For(j,i+1,len)
        {
            if(v[j].y-v[i].y<=d&&v[j].sta!=v[i].sta)
            {
                double d3=dis(v[i],v[j]);
                d=min(d,d3);
            }
        }
    }
    v.clear();
    return d;
}
int main()
{
    int t;
    scanf("%d", &t);
    while(t--)
    {
        mst(p,0);
        int n;
        scanf("%d", &n);
        For(i,1,n)
        {
            p[i].sta=1;
            scanf("%lld%lld",&p[i].x,&p[i].y);
        }
        For(i,n+1,2*n)scanf("%lld%lld",&p[i].x,&p[i].y);
        sort(p+1,p+1+2*n,cmpx);
        printf("%.3f\n",divd(1,2*n));
    }
    return 0;
}