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

题目链接：点击进入

题目

题意

n组，每组2个炸弹的坐标，每一组选且只选一个炸弹，每个炸弹的爆炸半径相同，求出最大的爆炸范围半径。（不能存在有炸弹重复炸到的地方）

思路

二分求最大半径，用2-sat检查是否存在合理的答案，如果两个点（不是同一对的）的距离小于2*r（我用的是两边平方），那么说明两圆相交，有冲突，那么就连边到它们各自取反的点，剩下的就是常规的2-sat

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<stack>
#include<map>
#include<algorithm>
using namespace std;
const int maxn=1e3+10;
const double eps=1e-7;
stack<int>s;
int t,n,m,cnt,tot,flag,ant;
int dfn[maxn],low[maxn],belong[maxn],head[maxn];
int x[maxn],y[maxn];
bool vis[maxn];
double dis[310][310];
struct node
{
	int to;
	int next;
}edge[100010];
void init()
{
	memset(head,-1,sizeof(head));
	memset(vis,0,sizeof(vis));
    memset(dfn,0,sizeof(dfn));
    //memset(low,0,sizeof(low));
   // memset(belong,0,sizeof(belong));
    cnt=tot=ant=0;
    while(!s.empty())
        s.pop();
}
void add(int u,int v)
{
	edge[ant].to=v;
	edge[ant].next=head[u];
	head[u]=ant++;
}
void tarjan(int pos)
{
    dfn[pos]=low[pos]=++cnt;
    s.push(pos);
    vis[pos]=1;
    for(int i=head[pos];i!=-1;i=edge[i].next)
    {
        int v=edge[i].to;
        if(!dfn[v])
        {
            tarjan(v);
            low[pos]=min(low[pos],low[v]);
        }
        else if(vis[v]==1)
            low[pos]=min(low[pos],dfn[v]);
    }
    if(low[pos]==dfn[pos])
    {
    	int v;tot++;
    	do{
    		v=s.top();
    		s.pop();
    		vis[v]=0;
    		belong[v]=tot;
		}while(pos!=v);
    }
}
double dist(int i,int j)
{
    return (x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]);
}
bool check(double r)
{
	init();
	double tmp=r*r*4;
	for(int i=0;i<n<<1;i++)
	{
		for(int j=0;j<i;j++)
		{
			if(dist(i,j)<tmp)
				add(i,j^1),add(j,i^1);
		}
	}
	for(int i=0;i<n<<1;i++)
        if(dfn[i]==0)
            tarjan(i);
    for(int i=0;i<n;i++)
        if(belong[i<<1]==belong[i<<1|1])
            return 0;
    return 1;
}
int main()
{
//	ios::sync_with_stdio(false);
	while(~scanf("%d",&n))
	{
		double ans=0.0;
		for(int i=0;i<n;i++)
			scanf("%d%d%d%d",&x[i<<1],&y[i<<1],&x[i<<1|1],&y[i<<1|1]);
		double l=0,r=1e9,mid;
        while(r-l>eps)
        {
            mid=(l+r)/2;
            if(check(mid))
                l=mid,ans=mid;
            else
                r=mid;
        }
        printf("%.2lf\n",ans);
	}
    return 0;
}