2-SAT

1．构图
2．求图的极大强连通子图
3．把每个子图收缩成单个节点，根据原图关系构造一个有向无环图
4．判断是否有解，无解则输出（退出）
5．对新图进行拓扑排序
6．自底向上进行选择、删除

HDU 3622

题意:给n对炸弹可以放置的位置(每个位置为一个二维平面上的点),每次放置炸弹是时只能选择这一对中的其中一个点,

每个炸弹爆炸的范围半径都一样,控制爆炸的半径使得所有的爆炸范围都不相交(可以相切),求解这个最大半径.

Solution：

首先二分最大半径值,然后2-sat构图判断其可行性,对于每两队位置(u,uu)和(v,vv),

如果u和v之间的距离小于2*id,也就是说位置u和位置v处不能同时防止炸弹(两范围相交)

所以连边(u,vv)和(v,uu),求解强连通分量判断可行性.

#include<iostream>
#include<cstdio>
#include<cmath>
#include<stack>
#include<cstring>
using namespace std;
const double Eps=1e-5;
int n,m,Cnt,Scc[210],Dfs_t;
int Pre[210],Low[210],Scc_Cnt;
stack<int> S;

struct Node
{
	int to;
	Node *nxt;
}G[210],Poor[50010];

struct Pos
{
	int x,y;
}p[210];

inline void Push(int u,int v)
{
	 Node *t=&Poor[++Cnt];
	 t->to=v;
	 t->nxt=G[u].nxt;
	 G[u].nxt=t;
}

inline double Dis(int i,int j)
{
	return sqrt((double)(p[i].x-p[j].x)*
	 (p[i].x-p[j].x)+(p[i].y-p[j].y)*(p[i].y-p[j].y));
}

void Init()
{
	Cnt=Dfs_t=Scc_Cnt=0;
	for(int i=0;i<2*n;++i)G[i].nxt=NULL;
	memset(Pre,0,sizeof(Pre));
	memset(Low,0,sizeof(Low));
	memset(Scc,0,sizeof(Scc));
}

void Dfs(int u)
{
	Pre[u]=Low[u]=++Dfs_t;
	S.push(u);
	for(Node *i=G[u].nxt;i;i=i->nxt)
	{
		int v=i->to;
		if(!Pre[v])
		{
			Dfs(v);
			Low[u]=min(Low[u],Low[v]);
		}
		else
		if(!Scc[v])
		Low[u]=min(Low[u],Pre[v]);
	}
	if(Low[u]==Pre[u])
	{
		Scc_Cnt++;
		while(1)
		{
			int x=S.top();
			S.pop();
			Scc[x]=Scc_Cnt;
			if(x==u)break;
		}
	}
}

void Tarjan(double M)
{
	for(int i=0;i<2*n;++i)
	{
		int t;
		if(i&1)t=i+1;
		else t=i+2;
		for(int j=t;j<2*n;++j)
		if(Dis(i,j)<2*M)
		{
			Push(i,j^1);
			Push(j,i^1);
		}
	}
	for(int i=0;i<2*n;++i)
	if(!Pre[i])Dfs(i);
}

bool Check(double M)
{
	Init();
	Tarjan(M);
	for(int i=0;i<2*n;i+=2)
	if(Scc[i]==Scc[i^1])return false;
	return true;
}

int main()
{
	while(~scanf("%d",&n))
	{
		for(int i=0;i<2*n;i+=2)
		{
			scanf("%d %d",&p[i].x,&p[i].y);
			scanf("%d %d",&p[i^1].x,&p[i^1].y);
		}
		double L=0,R=40000.0,Res=0;
		while(fabs(R-L)>Eps)
		{
			double M=(L+R)/2.0;
			if(Check(M))
			{
				Res=M;
				L=M+Eps;
			}
			else R=M-Eps;
		}
		printf("%.2f\n",Res);
	}
	return 0;
}