题目描述:
有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;
}