POJ 2749 Building roads 2-sat+二分答案

把爱恨和最大距离视为限制条件,可以知道,最大距离和限制条件多少具有单调性

所以可以二分最大距离,加边+check

 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<vector>
 #include<stack>
 #define N 5010
 #define INF 4000000
 using namespace std;
 int dis1[N],a,b,n,dis2[N],sx1,sy1,sx2,sy2,head[N],dfn[N],low[N],ecnt,x[N],y[N],tmp,ok,belong[N],hate[N][],love[N][],indx,inst[N],cnt;
 stack <int> st;
 int read()
 {
     int ret=,neg=;
     char j=getchar();
     for (;j>'' || j<'';j=getchar())
     if (j=='-') neg=-;
     for (;j>='' && j<='';j=getchar())
     ret=ret*+j-'';
     return ret*neg;
 }
 struct edge
 {
     int nxt,v;
 }e[N*N];
 int ABS(int x)
 {
     return x>?x:-x;
 }
 void add(int u,int v)
 {
     e[++ecnt].v=v;
     e[ecnt].nxt=head[u];
     head[u]=ecnt;
 }
 void init()
 {
     memset(head,,sizeof(head));
     memset(dfn,,sizeof(dfn));
     ecnt=;
     indx=;
     cnt=;
 }
 void buildG(int lim)
 {
     for (int i=;i<=a;i++)
     {
     int u=hate[i][],v=hate[i][];
     add(u,v+n),add(v,u+n),add(u+n,v),add(v+n,u);
     }
     for (int i=;i<=b;i++)
     {
     int u=love[i][],v=love[i][];
     add(u,v),add(v+n,u+n),add(u+n,v+n),add(v,u);
     }
     for (int i=;i<=n;i++)
     for (int j=i+;j<=n;j++)
     {
         if (dis1[i]+dis2[j]+tmp>lim)
         add(i,j),add(j+n,i+n);
         if (dis2[i]+dis1[j]+tmp>lim)
         add(i+n,j+n),add(j,i);
         if (dis1[i]+dis1[j]>lim)
         add(i,j+n),add(j,i+n);
         if (dis2[i]+dis2[j]>lim)
         add(i+n,j),add(j+n,i);
     }
 }
 void tar(int u)
 {
     dfn[u]=low[u]=++indx;
     inst[u]=;
     st.push(u);
     for (int i=head[u];i;i=e[i].nxt)
     {
     int v=e[i].v;
     if (!dfn[v])
     {
         tar(v);
         low[u]=min(low[v],low[u]);
     }
     else if (inst[v])
         low[u]=min(dfn[v],low[u]);
     }
     if (low[u]==dfn[u])
     {
     int t;
     ++cnt;
     while ()
     {
         t=st.top();
         inst[t]=;
         st.pop();
         belong[t]=cnt;
         if (t==u) break;
     }
     }
 }
 int check(int lim)
 {
     init();
     buildG(lim);
     for (int i=;i<=*n;i++)
     if (!dfn[i]) tar(i);
     for (int i=;i<=n;i++)
     if (belong[i]==belong[i+n]) return ;
     return ;
 }
 int main()
 {
     n=read(),a=read(),b=read();
     sx1=read(),sy1=read(),sx2=read(),sy2=read();
     tmp=ABS(sx1-sx2)+ABS(sy1-sy2);
     for (int i=;i<=n;i++)
     {
     x[i]=read(),y[i]=read();
     dis1[i]=ABS(x[i]-sx1)+ABS(y[i]-sy1);
     dis2[i]=ABS(x[i]-sx2)+ABS(y[i]-sy2);
     }
     for (int i=;i<=a;i++)
     hate[i][]=read(),hate[i][]=read();
     for (int i=;i<=b;i++)
     love[i][]=read(),love[i][]=read();
     int l=,r=INF,mid;
     while (l<r)
     {
     mid=(l+r)>>;
     if (check(mid)==) r=mid;
     else l=mid+;
     }
     printf("%d\n",l==INF?-:l);
     return ;
 }