hdu3622：Bomb Game

给n<=100对点，从每对点里面挑一个并以这些挑出的点为圆心画圆，并且这些圆不能互相覆盖，找出一种方案使得这些圆半径中最小的那个最大。

“最小值最大”就是二分答案啦！考虑现在每个点都画出半径x的圆，如何选点呢？

可以发现选了一个点P之后与其距离相差2x内的点Q都不能被选，也就是“与P在同一对的另一个点P'    或    与Q在同一对的另一个点Q'”

有了或语句描述选点冲突，接下来就是2-SAT构图啦！

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<cstdlib>
 #include<cmath>
 //#include<iostream>
 using namespace std;
 
 int n;
 #define maxn 111*2
 struct Edge{int to,next;};
 #define maxe maxn*maxn
 struct Graph
 {
     Edge edge[maxe];int le;
     int first[maxn],vis[maxn];
     void clear()
     {
         le=;
         memset(first,,sizeof(first));
     }
     void insert(int x,int y)
     {
         edge[le].to=y;
         edge[le].next=first[x];
         first[x]=le++;
     }
     void add_clause(int x,int y)
     {
         insert(x^,y);
         insert(y^,x);
     }
     int sta[maxn],top;
     bool dfs(int x)
     {
         if (vis[x^]) return ;
         if (vis[x]) return ;
         vis[x]=;
         sta[++top]=x;
         for (int i=first[x];i;i=edge[i].next)
             if (!dfs(edge[i].to)) return ;
         return ;
     }
     bool twosat()
     {
         memset(vis,,sizeof(vis));
         for (int i=;i<=n;i++)
             if (!vis[i*] && !vis[i*+])
             {
                 top=;
                 if (!dfs(i*))
                 {
                     for (;top;top--) vis[sta[top]]=;
                     if (!dfs(i*+)) return ;
                 }
             }
         return ;
     }
 }G;
 struct Point
 {
     int x,y;
 }a[maxn];
 double dis[maxn][maxn];
 #define eps 1e-5
 double sqr(double x) {return x*x;}
 double caldis(int p,int q)
 {
     return sqrt(sqr(a[p].x-a[q].x)+sqr(a[p].y-a[q].y));
 }
 bool judge(double x)
 {
     for (int i=;i<=n;i++)
         for (int j=i+;j<=n;j++)
         {
             if (dis[i*][j*]<x*) G.add_clause(i*+,j*+);
             if (dis[i*+][j*]<x*) G.add_clause(i*,j*+);
             if (dis[i*][j*+]<x*) G.add_clause(i*+,j*);
             if (dis[i*+][j*+]<x*) G.add_clause(i*,j*);
         }
     return G.twosat();
 }
 int main()
 {
     while (~scanf("%d",&n))
     {
         for (int i=;i<=n;i++)
             scanf("%d%d%d%d",&a[i*].x,&a[i*].y,&a[i*+].x,&a[i*+].y);
         double L=0.0,R=0.0;
         for (int i=;i<=n*+;i++)
             for (int j=;j<=n*+;j++)
                 dis[i][j]=caldis(i,j),R=max(R,dis[i][j]);
         while (R-L>=eps)
         {
             G.clear();
             double mid=(L+R+eps)/;
             if (judge(mid)) L=mid;
             else R=mid-eps;
         }
         printf("%.2lf\n",L);
     }
     return ;
 }