HDU 3622 Bomb Game (二分+2-SAT)

题目地址：HDU 3622

先二分半径。然后小于该半径的不能选，对这些不能选的点对进行加边。然后推断可行性就可以。

代码例如以下：

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <ctype.h>
#include <queue>
#include <map>
#include <set>
#include <algorithm>

using namespace std;
#define LL __int64
const int INF=0x3f3f3f3f;
const double eqs=1e-3;
int head[410], cnt, top, ans, index;
int dfn[410], low[410], instack[410], stak[410], belong[410];
struct node
{
    int u, v, next;
}edge[1000000];
struct Point
{
    int x, y;
}dian[10000];
void add(int u, int v)
{
    edge[cnt].v=v;
    edge[cnt].next=head[u];
    head[u]=cnt++;
}
void tarjan(int u)
{
    dfn[u]=low[u]=++index;
    instack[u]=1;
    stak[++top]=u;
    for(int i=head[u];i!=-1;i=edge[i].next)
    {
        int v=edge[i].v;
        if(!dfn[v])
        {
            tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if(instack[v])
        {
            low[u]=min(low[u],dfn[v]);
        }
    }
    if(dfn[u]==low[u])
    {
        ans++;
        while(1)
        {
            int v=stak[top--];
            instack[v]=0;
            belong[v]=ans;
            if(u==v) break;
        }
    }
}
void init()
{
    memset(head,-1,sizeof(head));
    cnt=0;
    memset(dfn,0,sizeof(dfn));
    memset(instack,0,sizeof(instack));
    top=ans=index=0;
}
double dist(Point x, Point y)
{
    return sqrt((x.x-y.x)*(x.x-y.x)*1.0+(x.y-y.y)*(x.y-y.y));
}
int solve(double mid, int n)
{
    int i, j;
    init();
    for(i=0;i<n<<1;i++)
    {
        for(j=0;j<i;j++)
        {
            if(dist(dian[i],dian[j])<mid)
            {
                add(i,j^1);
                add(j,i^1);
                //add(j^1,i);
                //add(i^1,j);
                //printf("%d %d %d %d\n",i<<1,j<<1|1,j<<1,i<<1|1);
            }
        }
    }
    for(i=0;i<n<<1;i++)
    {
        if(!dfn[i])
            tarjan(i);
    }
    for(i=0;i<n;i++)
    {
        if(belong[i<<1]==belong[i<<1|1])
            return 0;
    }
    return 1;
}
int main()
{
    int n, i, j;
    while(scanf("%d",&n)!=EOF)
    {
        for(i=0;i<n;i++)
        {
            scanf("%d%d%d%d",&dian[i<<1].x,&dian[i<<1].y,&dian[i<<1|1].x,&dian[i<<1|1].y);
        }
        double r=40000.0, l=0, mid, ans;
        while(r-l>eqs)
        {
            mid=(l+r)/2;
            //printf("%.2lf\n",mid);
            if(solve(mid, n))
            {
                ans=mid;
                l=mid;
            }
            else
                r=mid;
        }
        printf("%.2lf\n",ans/2.0);
    }
    return 0;
}