poj2749 Building roads（二分答案+2-sat）

求任意两点间最大距离最小。我们二分答案，每次用2-sat去判定是否可行。
三种限制：
1、i,j不能连同一个。
2、i,j必须连同一个
3、i连k1,j连k2就超了的话，这俩互斥。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <stack>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 520
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m1,m2,id[N][2],tot=-1,w[N][2],X[2],Y[2],h[N<<1],num=0,h1[N<<1],num1;
int dfn[N<<1],low[N<<1],dfnum,scc,bel[N<<1],W;
bool inq[N<<1];
struct edge{
    int to,next;
}data[1000010];
inline void add(int x,int y){
    data[++num].to=y;data[num].next=h[x];h[x]=num;
}stack<int>qq;
inline void tarjan(int x){
    low[x]=dfn[x]=++dfnum;qq.push(x);inq[x]=1;
    for(int i=h[x];i;i=data[i].next){
        int y=data[i].to;
        if(!dfn[y]) tarjan(y),low[x]=min(low[x],low[y]);
        else if(inq[y]) low[x]=min(low[x],dfn[y]);
    }if(low[x]==dfn[x]){
        ++scc;while(1){
            int y=qq.top();qq.pop();inq[y]=0;
            bel[y]=scc;if(y==x) break;
        }
    }
}
inline bool jud(int mid){
    memcpy(h,h1,sizeof(h));num=num1;
    memset(dfn,0,sizeof(dfn));dfnum=0;scc=0;
    for(int i=1;i<=n;++i)
        for(int j=i+1;j<=n;++j)
            for(int k1=0;k1<2;++k1)
                for(int k2=0;k2<2;++k2){
                    int res=w[i][k1]+w[j][k2]+W*(k1^k2);
                    if(res>mid)
                        add(id[i][k1],id[j][k2^1]),add(id[j][k2],id[i][k1^1]);
                }   
    for(int i=0;i<=tot;++i) if(!dfn[i]) tarjan(i);
    for(int i=1;i<=n;++i) if(bel[id[i][0]]==bel[id[i][1]]) return 0;return 1;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m1=read();m2=read();
    X[0]=read();Y[0]=read();X[1]=read();Y[1]=read();W=abs(X[0]-X[1])+abs(Y[0]-Y[1]);
    for(int i=1;i<=n;++i){
        id[i][0]=++tot;id[i][1]=++tot;
        int x=read(),y=read();w[i][0]=abs(X[0]-x)+abs(Y[0]-y);
        w[i][1]=abs(X[1]-x)+abs(Y[1]-y);
    }while(m1--){
        int x=read(),y=read();
        add(id[x][0],id[y][1]);add(id[x][1],id[y][0]);
        add(id[y][1],id[x][0]);add(id[y][0],id[x][1]);
    }while(m2--){
        int x=read(),y=read();
        add(id[x][0],id[y][0]);add(id[x][1],id[y][1]);
        add(id[y][0],id[x][0]);add(id[y][1],id[x][1]);
    }for(int i=0;i<=tot;++i) if(!dfn[i]) tarjan(i);
    for(int i=1;i<=n;++i) if(bel[id[i][0]]==bel[id[i][1]]){puts("-1");return 0;}
    int l=0,r=12e6;memcpy(h1,h,sizeof(h1));num1=num;
    while(l<=r){
        int mid=l+r>>1;
        if(jud(mid)) r=mid-1;
        else l=mid+1;
    }printf("%d\n",r+1);
    return 0;
}