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.

  1 #include<cstdio>
  2 #include<cstring>
  3 #include<algorithm>
  4 #include<iostream>
  5 #include<cmath>
  6 using namespace std;
  7 const int N=310;
  8 const int M=30010;
  9 const double eps=1e-4;
 10 struct Edge{
 11     int to,next;
 12 } edge[M];
 13 struct Node{
 14   double x;
 15   double y;
 16 }num1[600];
 17 int  n,m,cnt,dep,top,atype;
 18 int  dfn[N],low[N],vis[N],head[N],st[N],belong[N],in[N],out[N],sum[N];
 19 //sum[i]记录第i个连通图的点的个数,in[i],out[i],表示缩点之后点的入度和初度。
 20 void init(){
 21       cnt=dep=top=atype=0;
 22     memset(head,-1,sizeof(head));
 23     memset(dfn,0,sizeof(dfn));
 24     memset(low,0,sizeof(low));
 25     memset(vis,0,sizeof(vis));
 26     memset(belong,0,sizeof(belong));
 27     memset(in,0,sizeof(in));
 28     memset(out,0,sizeof(out));
 29     memset(sum,0,sizeof(sum));
 30 }
 31 void addedge(int u,int v){
 32     edge[cnt].to=v;
 33     edge[cnt].next=head[u];
 34     head[u]=cnt++;
 35 }
 36
 37 void Tarjan(int u){
 38     dfn[u]=low[u]=++dep;
 39     st[top++]=u;
 40     vis[u]=1;
 41     for(int i=head[u]; i!=-1; i=edge[i].next){
 42         int v=edge[i].to;
 43         if(!dfn[v]){
 44             Tarjan(v);
 45             low[u]=min(low[u],low[v]);
 46         }
 47         else if(vis[v]){
 48             low[u]=min(low[u],dfn[v]);
 49         }
 50     }
 51     int j;
 52     if(dfn[u]==low[u]){
 53         atype++;
 54         do{
 55             j=st[--top];
 56             belong[j]=atype;
 57             sum[atype]++;   //记录每个连通分量中点的个数
 58             vis[j]=0;
 59         }
 60         while(u!=j);
 61     }
 62 }
 63 bool judge(double mid){
 64     init();
 65     //int base=2*n;
 66     for(int i=1;i<=2*n;i++){
 67         for(int j=i+1;j<=2*n;j++){
 68             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));
 69             if(dis<2*mid){
 70                 if(i<=n&&j<=n){
 71                     addedge(i,j+n);
 72                     addedge(j,i+n);
 73                 }
 74                 if(i<=n&&j>n){
 75                     addedge(i,j-n);
 76                     addedge(j,i+n);
 77                 }
 78                 if(i>n&&j<=n){
 79                     addedge(i,j+n);
 80                     addedge(j,i-n);
 81                 }
 82                 if(i>n&&j>n){
 83                     addedge(i,j-n);
 84                     addedge(j,i-n);
 85                 }
 86
 87             }
 88         }
 89     }
 90     for(int i=1;i<=2*n;i++)
 91       if(!dfn[i])Tarjan(i);
 92     for(int i=1;i<=n;i++)
 93        if(belong[i]==belong[i+n])
 94           return false;
 95
 96   return true;
 97
 98 }
 99 int main(){
100     while(~scanf("%d",&n)){
101         for(int i=1;i<=n;i++){
102         scanf("%lf %lf %lf %lf",&num1[i].x,&num1[i].y,&num1[i+n].x,&num1[i+n].y);
103         }
104         double l=0,r=10000,ans=0;
105         while(abs(l-r)>eps){
106             double mid=(l+r)/2.0;
107             if(judge(mid)){
108                 ans=mid;
109                 l=mid;
110             }
111             else{
112                 r=mid;
113             }
114         }
115         printf("%.2f\n",ans);
116     }
117
118 }