POJ--2749[Building roads] 2-SAT第五题，2—SAT+二分判定

题目描述：

有n个农场，每个农场有坐标x，y。

有两个集合点s1和s2（也有坐标），每个农场必须连接其中的一个（有且仅有一个）。

然后有A个条件，每个条件a，b表示a农场不能和b农场连接在一个集合点。

然后再有B个条件，每个条件a，b表示a农场必须和b农场连接在一个集合点。

问你，在各种合法的连接情况中，任何两个农场间的距离的最大值的最小值是多少。

 

思路分析：

每个农场i分成两个点，i和i + n，前面表示连接左侧集合点，后面的表示连接右侧集合点。

对于条件A中的ab，连接

<a, b + n> <a + n, b> <b, a + n> <b + n, a>

这样就保证了ab不再一个集合点。同理，B条件也很好写。

然后就是用二分枚举答案（距离的最大值）。

对于每一次枚举的答案key，枚举任意两个农场a,b, 如果他俩的4种距离（a->s1->s1>b, a->s2->s2->b, a->s1->s2->b, a->s2->s1->b)有大于key的，就加入判定条件，不能这样连接。

这个key能够成功，二分右边界等于mid - 1，否则，左边界等于mid + 1.

 

CODE:

/*2-SAT第五题*/
/*AC代码：375ms*/
#include <iostream>
#include <cmath>
#define MAXN 1005
#define INF 0x7fffffff
#define max(a,b) (a>b?a:b)
#define min(a,b) (a<b?a:b)
struct edge
{
	int from,to,next;
}E[2000000];
int head[MAXN],ecnt;
int Low[MAXN],DFN[MAXN],Stack[MAXN],Belong[MAXN];
bool Instack[MAXN];
int Index,scc,top,N,A,B,len,Min,Max;
struct point
{
	int x,y;
};
struct point P[MAXN],S[3],A1[MAXN],B1[MAXN];
void Insert(int from,int to)
{
	E[ecnt].from=from;
	E[ecnt].to=to;
	E[ecnt].next=head[from];
	head[from]=ecnt++;
}
void Tarjan(int u)
{
	int i,v;
	Low[u]=DFN[u]=++Index;
	Stack[++top]=u;
	Instack[u]=true;
	for(i=head[u];i!=-1;i=E[i].next)
	{
		v=E[i].to;
		if(!DFN[v])
		{
			Tarjan(v);
			if(Low[u]>Low[v])
				Low[u]=Low[v];
		}
		else if(Instack[v]&&Low[u]>DFN[v])
			Low[u]=DFN[v];
	}
	if(Low[u]==DFN[u])
	{
		scc++;
		do{
			v=Stack[top--];
			Instack[v]=false;
			Belong[v]=scc;
		}while(u!=v);
	}
	return;
}
int get_dis(int i,int j)
{
	return abs(P[i].x-S[j].x)+abs(P[i].y-S[j].y);
}
void Init()
{
	int i;
	Max=-1;Min=INF;
	scanf("%d%d%d%d",&S[1].x,&S[1].y,&S[2].x,&S[2].y);
	len=abs(S[1].x-S[2].x)+abs(S[1].y-S[2].y);
	for(i=1;i<=N;i++)
	{
		scanf("%d%d",&P[i].x,&P[i].y);
		Max=max(Max,get_dis(i,1));Max=max(Max,get_dis(i,2));
		Min=min(Min,get_dis(i,1));Min=min(Min,get_dis(i,2));
	}
	for(i=1;i<=A;i++)
		scanf("%d%d",&A1[i].x,&A1[i].y);
	for(i=1;i<=B;i++)
		scanf("%d%d",&B1[i].x,&B1[i].y);
}
bool Judge(int key)
{
	int i,j,x,y;
	memset(head,-1,sizeof(head));ecnt=0;
	for(i=1;i<=A;i++)
	{
		x=A1[i].x;y=A1[i].y;
		Insert(x,y+N);
		Insert(y+N,x);
		Insert(y,x+N);
		Insert(x+N,y);
	}
	for(i=1;i<=B;i++)
	{
		x=B1[i].x;y=B1[i].y;
		Insert(x,y);
		Insert(y,x);
		Insert(x+N,y+N);
		Insert(y+N,x+N);
	}
	for(i=1;i<=N;i++)//枚举任意两个点(寻找矛盾)
	{
		for(j=i+1;j<=N;j++)
		{
			if(get_dis(i,1)+get_dis(j,1)>key)//(1,1)
			{
				Insert(i,j+N);
				Insert(j,i+N);
			}
			if(get_dis(i,2)+get_dis(j,2)>key)//(2,2)
			{
				Insert(i+N,j);
				Insert(j+N,i);
			}
			if(get_dis(i,1)+get_dis(j,2)+len>key)//(1,2)
			{
				Insert(i,j);
				Insert(j+N,i+N);
			}
			if(get_dis(i,2)+get_dis(j,1)+len>key)//(2,1)
			{
				Insert(i+N,j+N);
				Insert(j,i);
			}
		}
	}
	//Trajan+判定
	memset(Instack,false,sizeof(Instack));
	memset(DFN,0,sizeof(DFN));
	memset(Low,0,sizeof(Low));
	Index=scc=top=0;
	for(i=1;i<=2*N;i++)
	{if(!DFN[i]) Tarjan(i);}
	for(i=1;i<=N;i++)
		if(Belong[i]==Belong[i+N])
			return false;
		return true;
}
int main()
{
	int l,r,mid;
	while(scanf("%d%d%d",&N,&A,&B)!=EOF)
	{
		Init();
		l=2*Min;r=2*Max+len;
		int ans=INF;
		while(r>=l)//不能写成while(r>l)
		{
			mid=(l+r)>>1;
			if(Judge(mid))
			{
				r=mid-1;
				if(mid<ans)
					ans=mid;
			}
			else
				l=mid+1;
		}
		if(ans==INF)
			printf("-1\n");
		else
			printf("%d\n",ans);
	}
	return 0;
}