poj3714 Raid 题解报告

题目传送门

【题目大意】

有$N$个核电站和$N$个特工，均已知坐标，求特工与核电站之间的最短距离。

【思路分析】

平面最近点对板子题，注意特判一下两个点之间的关系，要求的是特工与核电站之间的距离，一共有$2N$个点。

平面最近点对相关知识走链接

【代码实现】

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define rg register
#define go(i,a,b) for(rg int i=a;i<=b;i++)
#define db double
using namespace std;
const int N=200002;
const db INF=1e9;
int T,n;
struct point{
    db x,y;
    bool mark;
}p[N];
db D(int a,int b){
    db X=p[a].x-p[b].x;
    db Y=p[a].y-p[b].y;
    return sqrt(X*X+Y*Y);
}
bool cmp(point &a,point &b){
    if(a.x!=b.x) return a.x<b.x;
    return a.y<b.y;
}
bool cm(int &a,int &b){
    return p[a].y<p[b].y;
}
int mid_id[N];
db work(int l,int r){
    if(l==r) return INF;
    if(l+1==r&&p[l].mark==p[r].mark) return INF;
    if(l+1==r) return D(l,r);
    int mid=(l+r)/2;db d=INF;
    db d1=work(l,mid),d2=work(mid+1,r);
    d=min(d1,d2);
    int num=0;
    go(i,l,r)
        if((db)abs(p[i].x-p[mid].x)<=d) mid_id[++num]=i;
    sort(mid_id+1,mid_id+1+num,cm);
    go(i,1,num-1)
        go(j,i+1,num){
        if(p[mid_id[j]].y-p[mid_id[i]].y>=d) break;
        if(p[mid_id[j]].mark==p[mid_id[i]].mark) continue;
        db d3=D(mid_id[i],mid_id[j]);
        d=min(d,d3);
    }
    return d;
}
int main(){
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        go(i,1,n) scanf("%lf%lf",&p[i].x,&p[i].y),p[i].mark=0;
        go(i,1,n) scanf("%lf%lf",&p[i+n].x,&p[i+n].y),p[i+n].mark=1;
        sort(p+1,p+1+n*2,cmp);
        printf("%.3lf
",work(1,n*2));
    }
    return 0;
}