【模板】 计几平面最近圆对

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3124

思路：二分答案，每个圆扩大len，找到最大使得所有圆都不想交的len即可。用到了扫描线，x轴排序，y轴用set维护。

参考了博客：https://blog.csdn.net/morejarphone/article/details/50114867

也不知道为什么删除的时候哦还要查看一边有没有相交。。。。。

#include<bits/stdc++.h>
using namespace std;
const int N = 5e4+9;
#define eps 1e-8
#define pii pair<double,int> 
struct Circle{
    double x,y,r;
}p[N],tp[N];
pii L[N*3];
int n;
bool inter(Circle a,Circle b){
    double dis = sqrt( (a.x - b.x)*(a.x - b.x) + (a.y - b.y)*(a.y - b.y) );
    return dis - a.r - b.r <= eps;
}
bool judge(double len){
    int cnt =0;
    for(int i = 1;i<=n;++i){
        tp[i] = p[i];
        tp[i].r += len;
        L[++cnt] = make_pair( tp[i].x - tp[i].r,i);
        L[++cnt] = make_pair( tp[i].x + tp[i].r,i+n);
    }
    sort(L+1,L+1+cnt);
    set< pii > st;
    for(int i = 1;i<=cnt;++i){
        int id = L[i].second;
        if(id <= n){
            set< pii >::iterator it = st.lower_bound( make_pair(tp[id].y,id) );
            if( it != st.end() ){
                int tid = it->second;
                if( inter(tp[id],tp[tid]) ) return 0;
            }
            if( it != st.begin() ){
                --it; int tid = it->second;
                if( inter(tp[id],tp[tid]) ) return 0;
            }
            st.insert( make_pair(tp[id].y,id) );
        }
        else st.erase( make_pair(tp[id-n].y,id-n) );
    }
    return 1;
}
int main(){
    int T; scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        for(int i = 1; i<=n;++i) scanf("%lf %lf %lf",&p[i].x,&p[i].y,&p[i].r);
        double ans,l=0,r=1000000;
        while(r-l >= eps){
            double m = (l+r)/2;
            if(judge(m)){
                l = m;
                ans = m;
            }
            else r = m;
        }
        printf("%lf\n",ans*2);
    }
    return 0;
}