POJ 2749 || HDU 1815 Building roads 2-sat

来源：互联网发布：网络在线测速编辑：程序博客网时间：2024/05/01 05:01

题目：

http://poj.org/problem?id=2749

题意：

有两个中转站并给出坐标，有n个牛舍并给出n个坐标，要把n个牛舍连接到两个中转站上使任意两个牛舍想通。每个牛舍里有一头牛，然后有A对牛之间是敌对关系，着意味着他们不能连接在同一个中转站上，有B对牛之间是朋友关系，意味着他们必须连在同一个中转站上，问满足上述条件的基础上，任意两个牛舍的最远距离最小是多少

思路：

最小化最大距离，二分枚举距离，然后2-sat判定是否可行。对于敌对关系的牛，意味着两只牛不能连接到同一个中转站，也不能同时不连接到同一个中转站（此时说明两只牛同时连接到另外一个中转站），即连边(i,~j)(j,~i)(~i,j)(~j,i)。对于朋友关系的牛，必须在同一个中转站，即连边(i,j)(j,i)(~i,~j)(~j,~i)。然后因为要满足任意两个牛舍之间的距离小于枚举值，所以对任意两个牛舍，如果距离大于枚举值，那么不能共存，即(i,~j)(j,~i)。另外注意枚举距离时的右界，虽然点坐标的范围是[-1000000, 1000000]，最大距离貌似是根号2乘以1000000，但是我们计算过程最大距离是这种情况：i-s1,s1-s2,s2-j，此时最大距离应该是乘以3的。。。

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#include <queue>#include <cmath>using namespace std;const int N = 1010, INF = 0x3f3f3f3f;const double eps = 1e-8;struct edge{    int to, next;} g[N*N*2];int cnt, head[N];int dfn[N], low[N], scc[N], st[N], top, num, idx;int x[N], y[N], a[N], b[N], c[N], d[N], dis[N*2];bool vis[N];int n, m, k, len;void add_edge(int v, int u){    g[cnt].to = u, g[cnt].next = head[v], head[v] = cnt++;}void init(){    memset(head, -1, sizeof head);    memset(dfn, -1, sizeof dfn);    memset(vis, 0, sizeof vis);    top = num = idx = cnt = 0;}void tarjan(int v){    dfn[v] = low[v] = ++idx;    vis[v] = true, st[top++] = v;    int u;    for(int i = head[v]; i != -1; i = g[i].next)    {        u = g[i].to;        if(dfn[u] == -1)        {            tarjan(u);            low[v] = min(low[v], low[u]);        }        else if(vis[u]) low[v] = min(low[v], dfn[u]);    }    if(dfn[v] == low[v])    {        num++;        do        {            u = st[--top], vis[u] = false, scc[u] = num;        }        while(u != v);    }}bool solve(int mid){    init();    for(int i = 1; i <= m; i++)    {        add_edge(a[i], b[i]+n), add_edge(b[i], a[i]+n);        add_edge(a[i]+n, b[i]), add_edge(b[i]+n, a[i]);    }    for(int i = 1; i <= k; i++)    {        add_edge(c[i], d[i]), add_edge(d[i], c[i]);        add_edge(c[i]+n, d[i]+n), add_edge(d[i]+n, c[i]+n);    }    for(int i = 1; i <= n; i++)        for(int j = i + 1; j <= n; j++)         {            if(dis[i] + dis[j] > mid) add_edge(i, j+n), add_edge(j, i+n);            if(dis[i+n] + dis[j+n] > mid) add_edge(i+n, j), add_edge(j+n, i);            if(dis[i+n] + dis[j] + len > mid) add_edge(i+n, j+n), add_edge(j, i);            if(dis[i] + dis[j+n] + len > mid) add_edge(i, j), add_edge(j+n, i+n);        }    for(int i = 1; i <= 2*n; i++)        if(dfn[i] == -1) tarjan(i);    for(int i = 1; i <= n; i++)        if(scc[i] == scc[i+n]) return false;    return true;}int main(){    while(~ scanf("%d%d%d", &n, &m, &k))    {        int sx1, sy1, sx2, sy2;        scanf("%d%d%d%d", &sx1, &sy1, &sx2, &sy2);        for(int i = 1; i <= n; i++) scanf("%d%d", &x[i], &y[i]);        for(int i = 1; i <= m; i++) scanf("%d%d", &a[i], &b[i]);        for(int i = 1; i <= k; i++) scanf("%d%d", &c[i], &d[i]);        for(int i = 1; i <= n; i++)            dis[i] = abs(x[i] - sx1) + abs(y[i] - sy1), dis[i+n] = abs(x[i] - sx2) + abs(y[i] - sy2);        len = abs(sx1 - sx2) + abs(sy1 - sy2);        int l = 0, r = 8000000, res = -1;        while(l <= r)        {            int mid = (l + r) >> 1;            if(solve(mid)) r = mid - 1, res = mid;            else l = mid + 1;        }        printf("%d\n", res);    }    return 0;}

0 0