poj2588 snakes （并查集判断连通性+几何）

最重要的写在最前面:这题要用c++编译器才能过！！！默认的G++过不了！！！坑害多少生灵，╭(╯^╰)╮
首先要发现，如果上下边界联通则一定过不去，如图1。我们把上边界叫做0，下边界叫做n+1，与他们相交的合并进一个集合，两圆相交的也合并进一个集合。看上下边界是否在同一个集合里即可判断能否通过。然后是找最大值的问题。先看起点的最大值，什么样的圆才会影响我们的起点最大值呢？至少，这个圆要与起点线（即左边界）相连，如果他还与上边界联通（注意是联通，不一定直接相连）,则他与起点线的交点以上都是不可以进入的，如图2。终点同理。
附上两张图帮助理解。(图片源自cjj490168650)

图1

图2

#include<cstdio>
#include<cstring>
#include<cmath>
#define N 1005
struct point{
    double x,y,r;
}p[N];
struct lr{
    double d;
    int i;
}l[N],r[N];
int n,fa[N],ll=0,rr=0;
inline int find(int x){
    return fa[x]==x?x:fa[x]=find(fa[x]);
}
void merge(int x,int y){
    int xx=find(x),yy=find(y);
    if(xx!=yy) fa[yy]=xx;
}
int main(){
//  freopen("a.in","r",stdin);
    scanf("%d",&n);
    for(int i=0;i<=n+1;++i) fa[i]=i;//0--上边界,n+1--下边界 
    // 输入n个圆心数据，并处理其与上下左右边界关系。
    for(int i=1;i<=n;++i){
        scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].r);
        if(p[i].y+p[i].r>=1000) merge(i,0);
        if(p[i].y-p[i].r<=0) merge(i,n+1);
        if(p[i].x-p[i].r<=0){
            l[++ll].d=p[i].y-sqrt(pow(p[i].r,2)-pow(p[i].x,2));//与起点线的下交点 
            l[ll].i=i;
        }
        if(p[i].x+p[i].r>=1000){
            r[++rr].d=p[i].y-sqrt(pow(p[i].r,2)-pow(1000-p[i].x,2));//与终点线的下交点
            r[rr].i=i;
        }
    }
    // 处理两圆相交的关系。
    for(int i=1;i<n;++i){
        for(int j=i+1;j<=n;++j){
            if(p[i].r+p[j].r>sqrt(pow(p[i].x-p[j].x,2)+pow(p[i].y-p[j].y,2))) merge(i,j);
        }
    }
    // 如果上边界与下边界联通，则一定过不去。
    if(find(0)==find(n+1)) printf("Bill will be bitten.");
    else{// 如果能过，找最高点。
        double lu=1000,ru=1000;
        for(int i=1;i<=ll;++i)//与左边界、上边界均相交,l[i].d以上一定不能过
            if(find(l[i].i)==find(0)&&l[i].d<lu) lu=l[i].d;
        for(int i=1;i<=rr;++i)
            if(find(r[i].i)==find(0)&&r[i].d<ru) ru=r[i].d;
        printf("Bill enters at (0.00, %.2lf) and leaves at (1000.00, %.2lf).",lu,ru);
    }
    return 0;
}