poj2749 Building roads 2-sat

链接：http://poj.org/problem?id=2749

题意：有n个奶牛，每个奶牛有坐标。有两个集散地，给出坐标，每个奶牛要链接到其中一个集散地上，让最远的两个奶牛的距离（曼哈顿距离）最近。而且有的奶牛间互相讨厌，不能连在一个集散地上，有的互相喜欢，必须连在一个上。

思路：很明显的2-sat，不过还是跑偏了。。首先给出的限制关系还有集散地数量为2都符合2-sat，不过2-sat中没有距离，只是判定。不过只要二分最大值，判断是否可行就可以了。

//开始还以为要让距离和最小，完全读错题了，后来发现了= =

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <stack>
using namespace std;
#define mem(a,b) memset((a),(b),sizeof((a)))
#define For(i,a,b) for(int (i)=(a);(i) < (b);(i)++)
#define Ror(i,a,b) for(int (i)=(a);(i) > (b);(i)--)
#define mp make_pair
#define pb push_back
#define inf 0x3f3f3f3f

void RI (int& x){
x = 0;
char c = getchar ();
while (c == ' '||c == '\n') c = getchar ();
bool flag = 1;
if (c == '-'){
flag = 0;
c = getchar ();
}
while (c >= '0' && c <= '9'){
x = x * 10 + c - '0';
c = getchar ();
}
if (!flag) x = -x;
}
void RII (int& x, int& y){RI (x), RI (y);}
void RIII (int& x, int& y, int& z){RI (x), RI (y), RI (z);}

const int maxn = 1010;
int n,m;
struct Side{
int v,next;
}side[maxn*maxn];
int top,node[maxn];
int t,dfn[maxn],low[maxn],sum;
stack<int>s;
void add_side(int u,int v){
side[top] = (Side){v,node[u]};
node[u] = top++;
}
void init(){
memset(node,-1,sizeof(node));
memset(dfn,0,sizeof(dfn));
top=sum=t=0;
}
void dfs(int u){
dfn[u]=low[u]=++t;
s.push(u);
for(int i=node[u];i!=-1;i=side[i].next){
int v=side[i].v;
if(!dfn[v])dfs(v);
if(dfn[v]!=-1)low[u]=min(low[u],low[v]);
}
if(low[u]==dfn[u]){
int v;
do{
v=s.top();s.pop();
dfn[v]=-1;
low[v]=sum;
}while(v!=u);
sum++;
}
}
bool solve(){
for(int i=0;i<n;i++)//判断是否有可行解
if(low[2*i]==low[2*i+1])return false;
return true;
}

struct PP{
int x,y;
}pp[2],cow[maxn],hate[maxn],like[maxn];
int get_dis(PP a,PP b){
return abs(a.x - b.x) + abs(a.y - b.y);
}
int dis[maxn][2],ddd;
int a,b;

void build(int mid){
For(i,0,a){
int a = hate[i].x*2,b = hate[i].y*2;
int fa = a+1,fb = b+1;
add_side(a,fb);
add_side(fa,b);
add_side(b,fa);
add_side(fb,a);
}
For(i,0,b){
int a = like[i].x*2,b = like[i].y*2;
int fa = a+1,fb = b+1;
add_side(a,b);
add_side(fa,fb);
add_side(b,a);
add_side(fb,fa);
}
For(i,0,n){
For(j,i+1,n){
For(a,0,2){
For(b,0,2){
if(a == b){
if(dis[i][a] + dis[j][b] > mid){
add_side(i*2+a,j*2+(!b));
add_side(j*2+b,i*2+(!a));
}
}
else {
if(dis[i][a] + dis[j][b] + ddd > mid){
add_side(i*2+a,j*2+(!b));
add_side(j*2+b,i*2+(!a));
}
}
}
}
}
}
}
int get_ans(){
int l = 0,r = 6000010,mid;
while(l < r){
init();
mid = (l + r)>>1;
build(mid);
For(i,0,2*n)
if(!dfn[i])dfs(i);
if(solve())r = mid;
else l = mid + 1;
}
if(r == 6000010)return -1;
return r;
}
int main()
{
//freopen("test.txt","r",stdin);
while(~scanf("%d%d%d",&n,&a,&b)){
RII(pp[0].x,pp[0].y);
RII(pp[1].x,pp[1].y);
ddd = get_dis(pp[0],pp[1]);
For(i,0,n){
RII(cow[i].x,cow[i].y);
}
For(i,0,n){
For(j,0,2){
dis[i][j] = get_dis(pp[j],cow[i]);
}
}
For(i,0,a){
RII(hate[i].x,hate[i].y);
hate[i].x --;
hate[i].y --;
}
For(i,0,b){
RII(like[i].x,like[i].y);
like[i].x --;
like[i].y --;
}
int ans = get_ans();
cout<<ans<<endl;
}
return 0;
}