HDU 3622 Bomb Game（2-SAT+二分）

Description
给n对炸弹可以放置的位置(每个位置为一个二维平面上的点),每次放置炸弹是时只能选择这一对中的其中一个点,每个炸弹爆炸的范围半径都一样,控制爆炸的半径使得所有的爆炸范围都不相交(可以相切),求解这个最大半径
Input
多组用例，每组用例第一行为一个整数n表示炸弹对数，之后n行每行四个整数表示该组两个炸弹的坐标，以文件尾结束输入
Output
对于每组用例，输出最大爆炸半径
Sample Input
2
1 1 1 -1
-1 -1 -1 1
2
1 1 -1 -1
1 -1 -1 1
Sample Output
1.41
1.00
Solution
二分爆炸半径，对于每个确定的半径，可以用2-SAT判断其可行性，对于两组炸弹pa1,pa2,pb1,pb2,若dis(pa1,pb1)<2*r,则有pa1->pb2,pb1->pa2，根据各组炸弹的相对位置构图，用tarjan算法求强联通分量并判断是否有使得布尔公式值为真的一组布尔变量赋值即可
Code

#include<iostream>  
#include<cstdio>  
#include<cstring>  
#include<algorithm>  
#include<vector>  
#include<stack>  
#include<cmath>  
#define maxn 255  
#define INF 500000005  
const double eps=1e-7;  
using namespace std;  
struct P  
{  
    int x,y;  
}pa[maxn],pb[maxn];  
vector<int>g[maxn];  
stack<int>st;  
int n,scc,index;  
int low[maxn],dfn[maxn],instack[maxn],fa[maxn];  
double dist(P a, P b)//两点间距离  
{  
    return sqrt(1.0*((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)));  
}  
void init()//初始化  
{  
    scc=index=0;  
    while(!st.empty())st.pop();  
    for(int i=0;i<maxn;i++)g[i].clear();  
    memset(dfn,0,sizeof(dfn));  
    memset(instack,0,sizeof(instack));  
    memset(low,0,sizeof(low));
}  
void tarjan(int u)//求强联通分量  
{  
    dfn[u]=low[u]=++index;  
    instack[u]=1;  
    st.push(u);  
    int v,size=g[u].size();  
    for(int i=0;i<size;i++)  
    {  
        v=g[u][i];  
        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])  
    {  
        scc++;  
        do  
        {  
            v=st.top();  
            st.pop();  
            fa[v]=scc;  
            instack[v]=0;  
        }while(v!=u);  
    }  
}  
void build(double mid)  
{  
    for(int i=1;i<=n;i++)  
        for(int j=i+1;j<=n;j++)  
        {  
            if(dist(pa[i],pa[j])+eps<2*mid) 
                g[i].push_back(j+n),g[j].push_back(i+n);  
            if(dist(pa[i],pb[j])+eps<2*mid) 
                g[i].push_back(j),g[j+n].push_back(i+n);  
            if(dist(pb[i],pa[j])+eps<2*mid) 
                g[i+n].push_back(j+n),g[j].push_back(i);  
            if(dist(pb[i],pb[j])+eps<2*mid) 
                g[i+n].push_back(j),g[j+n].push_back(i);  
        }  
}  
bool check()  
{  
    for(int i=1;i<=2*n;i++)//求强联通分量 
        if(!dfn[i]) 
            tarjan(i);  
    for(int i=1;i<=n;i++)  
        if(fa[i]==fa[i+n])//矛盾 
            return false;  
    return true;  
}  
int main()  
{  
    while(scanf("%d", &n)!=EOF)  
    {  
        for(int i=1;i<=n;i++) 
            scanf("%d%d%d%d",&pa[i].x,&pa[i].y,&pb[i].x,&pb[i].y);  
        double low=0,high=20000,ans=0;  
        while(high-low>eps)//二分爆炸范围,用ans记录答案 
        {  
            double mid=(low+high)/2;  
            init();//初始化  
            build(mid);//确定爆炸范围后构图  
            if(check()) 
            {
                ans=max(ans,mid); 
                low=mid;
            }  
            else 
                high=mid;  
        }  
        printf("%.2f\n",ans);    
    }  
    return 0;  
}