poj-2749 Building roads

题意：

给出平面上的n个点和两个特殊点s1,s2，定义距离为曼哈顿距离；

现想将这些点一些连到s1上，另外的连在s2上，并且连接s1,s2将整个图连成一颗树；

有一些规则，使x,y不能连到一个点上，或者x,y必须连到一个点上；

求所有可行方案中树的最小直径，如果没有合法方案输出-1；

题解：

对于这样一个最小化最大距离的问题我们很容易想到二分答案；

进行二分之后，问题转化成了是否存在一组解使树的直径不超过mid；

我们进行这样的转化之后，枚举所有的点对判断它们能否分别连到s1,s2上共四种情况；

这样也将二分的限制转化到和题目本身限制一样的形式啦；

那么这样的限制用2-SAT算法来解决，时间复杂度是O(n^2logn)的；

代码：

#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 1100
#define pr pair<int,int>
using namespace std;
pr p,s[2];
int dis[N][2],D;
int next[N*N],to[N*N],head[N],thead[N],ce,tce;
int n,dfn[N],low[N],belong[N],tim,tot;
int st[N],top;
bool ins[N];
void add(int x,int y)
{
	to[++ce]=y;
	next[ce]=head[x];
	head[x]=ce;
}
void dfs(int x)
{
	dfn[x]=low[x]=++tim;
	st[++top]=x;
	ins[x]=1;
	for(int i=head[x];i;i=next[i])
	{
		if(!dfn[to[i]])
			dfs(to[i]),low[x]=min(low[x],low[to[i]]);
		else if(ins[to[i]])
			low[x]=min(low[x],dfn[to[i]]);
	}
	if(dfn[x]==low[x])
	{
		tot++;
		int k;
		do
		{
			k=st[top--];
			ins[k]=0;
			belong[k]=tot;
		}while(k!=x);
	}
}
int calc(pr x,pr y)
{
	return abs(x.first-y.first)+abs(x.second-y.second);
}
bool check(int mid)
{
	memcpy(head,thead,sizeof(head));
	memset(dfn,0,sizeof(dfn));
	ce=tce;
	tot=top=tim=0;
	int i,j;
	for(i=1;i<=n;i++)
	{
		for(j=i+1;j<=n;j++)
		{
			if(dis[i][0]+dis[j][0]>mid)
			{
				add(i<<1,j<<1|1);
				add(j<<1,i<<1|1);
			}
			if(dis[i][1]+dis[j][1]>mid)
			{
				add(i<<1|1,j<<1);
				add(j<<1|1,i<<1);
			}
			if(dis[i][0]+dis[j][1]+D>mid)
			{
				add(i<<1,j<<1);
				add(j<<1|1,i<<1|1);
			}
			if(dis[i][1]+dis[j][0]+D>mid)
			{
				add(i<<1|1,j<<1|1);
				add(j<<1,i<<1);
			}
		}
	}
	for(i=1;i<=n;i++)
	{
		if(!dfn[i<<1])
			dfs(i<<1);
		if(!dfn[i<<1|1])
			dfs(i<<1|1);
		if(belong[i<<1]==belong[i<<1|1])
			return 0;
	}
	return 1;
}
int main()
{
	int A,B,i,j,k,x,y,l,r,mid;
	scanf("%d%d%d",&n,&A,&B);
	scanf("%d%d%d%d",&s[0].first,&s[0].second,&s[1].first,&s[1].second);
	D=calc(s[0],s[1]);
	for(i=1;i<=n;i++)
	{
		scanf("%d%d",&p.first,&p.second);
		dis[i][0]=calc(p,s[0]);
		dis[i][1]=calc(p,s[1]);
	}
	for(i=1;i<=A;i++)
	{
		scanf("%d%d",&x,&y);
		add(x<<1,y<<1|1);
		add(x<<1|1,y<<1);
		add(y<<1,x<<1|1);
		add(y<<1|1,x<<1);
	}
	for(i=1;i<=B;i++)
	{
		scanf("%d%d",&x,&y);
		add(x<<1,y<<1);
		add(x<<1|1,y<<1|1);
		add(y<<1,x<<1);
		add(y<<1|1,x<<1|1);
	}
	memcpy(thead,head,sizeof(head));
	tce=ce;
	l=0,r=4000000;
	while(l<=r)
	{
		mid=l+r>>1;
		if(check(mid))
			r=mid-1;
		else
			l=mid+1;
	}
	if(l>4000000)
		puts("-1");
	else
		printf("%d\n",l);
	return 0;
}