UVA 1146 Now or later(2-SAT问题)

2-SAT问题

有n个布尔型变量，然后同时有m个条件，每个条件都是形如“xi为真或xj为假“一类的问题。然后我们的做法就是把这些变量全都抽象成图中两个点，例如给定变量xi，我们抽象成2*i和2*i+1两个点，这两个点分别代表变量的真假两种状态。然后标记2*i点代表取xi的真状态，2*i+1代表假状态。这样的话因为有m个条件，我们就可以在不同变量抽象成的点之间连一些边。这些边代表的就是推导出的意思，即如果（2*i+1）---->（2*j）有一条边的话代表，如果xi为假的时候，xj必然为真。

然后如何判断该问题是否有解呢？我们求出这个有向图的强连通分支，然后判断某个变量所抽象成的两个点是不是同在一个强连通分支中，如果是，必然无解。因为这两个点之间有路的话，必然代表两个点可以互相推导出，也就是如果xi为真，那么根据传递性，xi同时也必然为假，这肯定是错的。

UVA 1146 Now or later

一般求最小值的最大可能，这一类最大最小的，都是用二分法来求。这一题的话我们就可以二分枚举着陆时间的最小间隔。

例如，如果第i辆飞机晚着陆时间与第j辆飞机早着陆时间中间的间隔小于现在所枚举的最小间隔t，那么我们就可以找出这么两条有向边。第一条是在2*i+1--->2*j+1，也就是说如果第i辆飞机是晚着陆的，必然第j辆飞机也得是晚着陆才满足条件。同理，我们可以找出另一条2*j---->2*i.注意如果第j辆飞机晚着陆，我们无法得知第i辆飞机到底是晚着陆还是早着陆，因为两个都可能满足条件，所以不能连接2*j+1--->2*i+1的边。

构好图之后就直接用tarjan算法求强连通分量了，如果对我的tarjan的风格不太明白的话，可以先看一下http://blog.csdn.net/u010489766/article/details/9360353，这里面我用的是同样的tarjan模版，而且注释很详细。

#include<stdio.h>
#include<string.h>
#include<vector>
#include<algorithm>
#include<stack>
#define MEM(x,y) memset(x,y,sizeof(x))
#define INF 2100
using namespace std;

int air[INF][2];	//用来存飞行器的早晚降落时间
vector<int>graph[2*INF];
stack<int>S;
int dfn[2*INF],low[2*INF],dfs_deep;
int	flag[2*INF],sub_id[2*INF],id;
int n;

void init()
{
	for(int i = 0 ; i <= 2*n; i ++)
		graph[i].clear();
	MEM(dfn,0);
	MEM(low,0);
	MEM(flag,0);
	MEM(sub_id,0);
	dfs_deep = id = 0;
}

void addedge(int x,int xval,int y,int yval)	//xval，yval代表想x,y的所取的值，是1或0
{
//x取xval时，y取yval不满足条件。
//所以x取xval时，y只能取yval^1,连一条有向边。但是y取yval^1时，x不一定取xval。
//同理，y取yval时，x只能取xval^1，连一条有向边。但是x取xval^1时，y不一定取yval.
	x = 2*x + xval;
	y = 2*y + yval;
	graph[x].push_back(y^1);
	graph[y].push_back(x^1);
}

void tarjan(int node)
{
	dfn[node] = low[node] = ++dfs_deep;
	S.push(node);
	flag[node] = 1;
	int len = graph[node].size();
	for(int i = 0 ; i < len ; i ++)
	{
		int temp = graph[node][i];
		if(!dfn[temp])
		{
			tarjan(temp);
			low[node] = min(low[node],low[temp]);
		}
		else if(flag[temp] == 1)
			low[node] = min(low[node],dfn[temp]);
	}
	if(dfn[node] == low[node])
	{
		id++;
		while(1)
		{
			int t = S.top();
			S.pop();
			flag[t] = 0;
			sub_id[t] = id;
			if(t == node)
				break;
		}
	}
}

bool check(int t)
{
	init();		//初始化
	for(int i = 0 ; i < n; i ++)
		for(int a = 0; a < 2; a ++)
			for(int j = i +1; j < n; j ++)
				for(int b = 0; b < 2; b ++)
					if(abs(air[i][a] - air[j][b]) < t)
						addedge(i,a,j,b);	//构图
	for(int i = 0; i < 2*n ; i ++)
		if(!dfn[i])
			tarjan(i);
	for(int i = 0; i < 2*n ; i += 2)
		if(sub_id[i] == sub_id[i^1])	//如果出现一个点的两个状态在同一个连通分量中，即这个点同时取两个状态，必错
			return false;
	return true;
}

int main()
{
	while(scanf("%d",&n) != EOF)
	{
		int l = 0,r = -1;
		MEM(air,0);
		for(int i = 0; i < n; i ++)	//此处只能从0开始存点，因为1^1=0,0^1=1
		{
			scanf("%d%d",&air[i][0],&air[i][1]);
			r = max(r,air[i][0]);
			r = max(r,air[i][1]);
		}
		//	二分法
		while(l < r)
		{
			int mid = l + (r-l+1)/2;
			if(check(mid))
				l = mid;
			else
				r = mid-1;
		}
		printf("%d\n",r);
	}
	return 0;
}