Bomb Game

hdu3622:http://acm.hdu.edu.cn/showproblem.php?pid=3622

题意：你有n次，每次你可以在平面上放置一个点，并且每一次都会有两个位置可以选，每一次只能选择其中一个。然后在自己位置上以该点为圆心画圆，这n个圆不能相交，问你最后最小的圆的半径的最大值是多少。

题解：二分+2-sat。对于2-sat的建图始终很迷糊。这里的一次有两个点，并且每一次能选择其中一个，这就相当于把n个点弄成了2*n点。如果i--j距离太近 ，则建边i-->~j,j—>~i.

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<iostream>
 #include<cmath>
 using namespace std;
 const int N=;
 const int M=;
 const double eps=1e-;
 struct Edge{
     int to,next;
 } edge[M];
 struct Node{
   double x;
   double y;
 }num1[];
 int  n,m,cnt,dep,top,atype;
 int  dfn[N],low[N],vis[N],head[N],st[N],belong[N],in[N],out[N],sum[N];
 //sum[i]记录第i个连通图的点的个数,in[i],out[i],表示缩点之后点的入度和初度。
 void init(){
       cnt=dep=top=atype=;
     memset(head,-,sizeof(head));
     memset(dfn,,sizeof(dfn));
     memset(low,,sizeof(low));
     memset(vis,,sizeof(vis));
     memset(belong,,sizeof(belong));
     memset(in,,sizeof(in));
     memset(out,,sizeof(out));
     memset(sum,,sizeof(sum));
 }
 void addedge(int u,int v){
     edge[cnt].to=v;
     edge[cnt].next=head[u];
     head[u]=cnt++;
 }

 void Tarjan(int u){
     dfn[u]=low[u]=++dep;
     st[top++]=u;
     vis[u]=;
     for(int i=head[u]; i!=-; i=edge[i].next){
         int v=edge[i].to;
         if(!dfn[v]){
             Tarjan(v);
             low[u]=min(low[u],low[v]);
         }
         else if(vis[v]){
             low[u]=min(low[u],dfn[v]);
         }
     }
     int j;
     if(dfn[u]==low[u]){
         atype++;
         do{
             j=st[--top];
             belong[j]=atype;
             sum[atype]++;   //记录每个连通分量中点的个数
             vis[j]=;
         }
         while(u!=j);
     }
 }
 bool judge(double mid){
     init();
     //int base=2*n;
     for(int i=;i<=*n;i++){
         for(int j=i+;j<=*n;j++){
             double dis=sqrt((num1[i].x-num1[j].x)*(num1[i].x-num1[j].x)+(num1[i].y-num1[j].y)*(num1[i].y-num1[j].y));
             if(dis<*mid){
                 if(i<=n&&j<=n){
                     addedge(i,j+n);
                     addedge(j,i+n);
                 }
                 if(i<=n&&j>n){
                     addedge(i,j-n);
                     addedge(j,i+n);
                 }
                 if(i>n&&j<=n){
                     addedge(i,j+n);
                     addedge(j,i-n);
                 }
                 if(i>n&&j>n){
                     addedge(i,j-n);
                     addedge(j,i-n);
                 }

             }
         }
     }
     for(int i=;i<=*n;i++)
       if(!dfn[i])Tarjan(i);
     for(int i=;i<=n;i++)
        if(belong[i]==belong[i+n])
           return false;

   return true;

 }
 int main(){
     while(~scanf("%d",&n)){
         for(int i=;i<=n;i++){
         scanf("%lf %lf %lf %lf",&num1[i].x,&num1[i].y,&num1[i+n].x,&num1[i+n].y);
         }
         double l=,r=,ans=;
         while(abs(l-r)>eps){
             double mid=(l+r)/2.0;
             if(judge(mid)){
                 ans=mid;
                 l=mid;
             }
             else{
                 r=mid;
             }
         }
         printf("%.2f\n",ans);
     }

 }