HDU-3622 Bomb Game

题意：

放炸弹，每组炸弹有两个，必须使用其中一枚。多个炸弹爆破范围不可重叠。询问得分为-最小的爆炸半径。问得分最高

思路;

2-sat跑图，二分爆炸半径，用2-sat判别是否可以成立。

#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<iostream>
#include<math.h>
using namespace std;
const int maxn=205;
const double eps=1e-5;
struct node
{
    int v,next;
}edge[maxn*maxn],edgef[maxn*maxn];
int head[maxn*maxn];
int headf[maxn*maxn];
int tot1,tot2;
int vis1[maxn],vis2[maxn];
int belong[maxn];
int T[maxn];
int Bcnt,Tcnt;
void add(int u,int v)
{
    edge[++tot1].v=v;
    edge[tot1].next=head[u];
    head[u]=tot1;

    edgef[++tot2].v=u;
    edgef[tot2].next=headf[v];
    headf[v]=tot2;
}

void init()
{
    memset(head,-1,sizeof(head));
    memset(headf,-1,sizeof(headf));
    memset(vis1,0,sizeof(vis1));
    memset(vis2,0,sizeof(vis2));
    tot1=tot2=0;
    Bcnt=Tcnt=0;
}

void dfs1(int u) //正向图 dfs ，算结束时间
{
    vis1[u]=1;
    for(int i=head[u];i!=-1;i=edge[i].next)
    {
        int v=edge[i].v;
        if(!vis1[v])
            dfs1(v);
    }
    T[++Tcnt]=u;
}

void dfs2(int u)
{
    vis2[u]=1;
    belong[u]=Bcnt;
    for(int i=headf[u];i!=-1;i=edgef[i].next)
    {
        int v=edgef[i].v;
        if(!vis2[v])
            dfs2(v);
    }
}
struct Point{
    int x,y;
}pt[maxn];
double dis(Point a,Point b)
{
    return sqrt((double)(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int judge(int n)
{
    for(int i=0;i<2*n;i++)
    {
        if(!vis1[i])
            dfs1(i);
    }
    for(int i=Tcnt;i>=1;i--)
    {
        if(!vis2[T[i]])
        {
            dfs2(T[i]);
            Bcnt++;
        }
    }
    for(int i=0;i<=2*n-2;i+=2)
    {
        if(belong[i]==belong[i+1])
            return 0;
    }
    return 1;
}
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=0;i<n;i++)
          scanf("%d%d%d%d",&pt[2*i].x,&pt[2*i].y,&pt[2*i+1].x,&pt[2*i+1].y);
        double left,right,mid;
        left=0;
        right=40000.0;
        while(right-left>=eps)
        {
            init();
            mid=(left+right)/2;

            for(int i=0;i<2*n-2;i++)
            {
                int t;
                if(i%2==0)
                    t=i+2;
                else 
                    t=i+1;
                for(int j=t;j<2*n;j++)
                   if(dis(pt[i],pt[j])<2*mid)//因为i与j矛盾，因此在 这一步骤的判定关系->上i必定与j^1 同选， j必定与i^1同选.
                   {
                       add(i,j^1);
                       add(j,i^1);
                   }
            }
            if(judge(n))
                left=mid;
            else
                right=mid;
        }
        printf("%.2lf\n",right);
    }
    return 0;
}