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;
}