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

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

图1

图2

#include<cstdio>#include<cstring>#include<cmath>#define N 1005struct 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;}