poj 2749 Building roads 2-SAT

同样我们二分距离，我们可以这样考虑，对每个仓库a用a表示到第一个结合点连边，非a到第二个转接点连边，接下来我们对于相互恨的和相互喜欢的可以直接加边，对于约束条件，我们暴力每两个仓库，判断这两个仓库分别连在两个转接点时是否可行，用2-SAT作为判断条件。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<iostream>
#include<cmath>
#define maxn 10105
using namespace std;
int low[maxn],dnf[maxn],que[maxn],tear,head[maxn],c[maxn],mark[maxn];
vector<int>g[maxn];
struct pi{
    int a;
    int b;
}pp[maxn],pp1[maxn],pp2[maxn];
void init(int n){
    for(int i=1;i<=n;i++) g[i].clear();
}
int cal,tot;
void tarjin(int p){
    dnf[p]=low[p]=cal++;
    int i,k;
    k=(int)g[p].size();
    que[tear++]=p;
    for(i=0;i<k;i++){
        int v=g[p][i];
        if(!dnf[v]){
            tarjin(v);
            low[p]=min(low[p],low[v]);//如果儿子就比较low
        }
        else if(!mark[v]){
            low[p]=min(low[p],dnf[v]);//与栈里面的点进行比较，如果mark[i]有值就代表已经出栈了。
        }
    }
    if(low[p]==dnf[p]){
        while(tear>0){
            mark[que[tear-1]]=tot;
            if(que[tear-1]==p){
                tear--;
                break;
            }
            tear--;
        }
        tot++;
    }
}
int solve(int n){
    int i;
    tear=0;
    cal=1;
    tot=1;
    memset(dnf,0,sizeof(dnf));
    memset(c,0,sizeof(c));
    memset(mark,0,sizeof(mark));
    for(i=1;i<=2*n;i++){
        tear=0;
        if(!dnf[i]) tarjin(i);
    }
    for(i=1;i<=n;i++){
        if(mark[i]==mark[n+i]) return 0;
    }
    return 1;
}
void add(int a,int b){
    g[a].push_back(b);
}
int a[maxn],b[maxn];
int main()
{
    int i,j,n,m,p,t;
    int x1,y1,x2,y2;
    while(cin>>n>>m>>p){
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        t=abs(x1-x2)+abs(y1-y2);
        for(i=1;i<=n;i++){
            scanf("%d%d",&pp[i].a,&pp[i].b);
            a[i]=abs(pp[i].a-x1)+abs(pp[i].b-y1);
            b[i]=abs(pp[i].a-x2)+abs(pp[i].b-y2);
        }
        for(i=1;i<=m;i++){
            scanf("%d%d",&pp1[i].a,&pp1[i].b);
        }
        for(i=1;i<=p;i++) scanf("%d%d",&pp2[i].a,&pp2[i].b);
        int le,ri,mid;
        le=0;
        ri=1000000000;
        while(le<=ri){
            mid=(le+ri)/2;
            init(2*n);
            for(i=1;i<=m;i++){
                add(pp1[i].a,n+pp1[i].b);
                add(n+pp1[i].a,pp1[i].b);
                add(pp1[i].b,n+pp1[i].a);
                add(n+pp1[i].b,pp1[i].a);
            }
            for(i=1;i<=p;i++){
                add(pp2[i].a,pp2[i].b);
                add(n+pp2[i].b,n+pp2[i].a);
                add(n+pp2[i].a,n+pp2[i].b);
                add(pp2[i].b,pp2[i].a);
            }
            for(i=1;i<=n;i++){
                for(j=1;j<i;j++){
                    if((a[i]+a[j])>mid){
                        add(i,n+j);
                        add(j,n+i);
                    }
                    if(b[i]+b[j]>mid){
                        add(n+i,j);
                        add(n+j,i);
                    }
                    if(a[i]+t+b[j]>mid){
                        add(i,j);
                        add(n+j,n+i);
                    }
                    if(b[i]+a[j]+t>mid){
                        add(n+i,n+j);
                        add(j,i);
                    }
                }
            }
            if(solve(n))  ri=mid-1;
            else le=mid+1;
        }
        if(le>1000000000) printf("-1\n");
        else printf("%d\n",le);
    }
}