题意:
给出平面上的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;
}