【ZJOI 2009】狼和羊的故事（最小割）

最新推荐文章于 2021-03-30 13:15:57 发布

Galaxy Coder

最新推荐文章于 2021-03-30 13:15:57 发布

阅读量348

点赞数 2

CC 4.0 BY-SA版权

分类专栏：最小割文章标签： BZOJ

本文链接：https://blog.youkuaiyun.com/weixin_42068627/article/details/80789106

最小割专栏收录该内容

2 篇文章

订阅专栏

本文介绍了ZJOI2009竞赛中一道名为“狼和羊的故事”的题目解决方法。该题通过建立最大流模型求解最优方案，特别讨论了非狼非羊领地的处理方式，即与周围所有节点相连，实现自由切割。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接

【ZJOI 2009】狼和羊的故事

题解

如果只有 $1$ 和 $2$ 的话，那么就将源点连羊，狼连汇点（容量 $\inf$ ），狼羊之间连一条容量为 $1$ 的边。

但是，有的领地既不是羊的，也不是狼的。
你可能会想把 $0$ 归为 $1$ （或 $2$ ），但是仔细想想可以发现还有更优的解。

最终的解决方案： $0$ 点与四周的所有节点连一条边，这样就可以在 $0$ 的四周自由切割，没有附加限制。

代码

#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
const int maxl = 105;
const int maxn = 10005;
const int maxm = 60005;
const int inf = 1000000000;
int n, m, src, snk, c[maxl][maxl], dis[maxn];
int tot, ter[maxm], len[maxm], nxt[maxm], lnk[maxn];
int id(int i, int j) {
    return (i - 1) * m + j;
}
void add(int u, int v, int w) {
    ter[tot] = v;
    len[tot] = w;
    nxt[tot] = lnk[u];
    lnk[u] = tot++;
}
void adde(int u, int v, int w) {
    add(u, v, w);
    add(v, u, w);
}
int min(int x, int y) {
    return x < y ? x : y;
}
bool bfs(int s, int t) {
    queue<int> que;
    que.push(s);
    memset(dis, -1, sizeof(dis));
    dis[s] = 0;
    for (int u, v, w; !que.empty(); ) {
        u = que.front();
        que.pop();
        for (int i = lnk[u]; ~i; i = nxt[i]) {
            v = ter[i], w = len[i];
            if (w && dis[v] == -1) {
                dis[v] = dis[u] + 1;
                que.push(v);
            }
        }
    }
    return ~dis[t];
}
int find(int u, int lft) {
    if (u == snk) {
        return lft;
    }
    int tmp, res = 0;
    for (int v, w, i = lnk[u]; ~i && res < lft; i = nxt[i]) {
        v = ter[i], w = len[i];
        if (w && dis[u] + 1 == dis[v]) {
            tmp = find(v, min(w, lft - res));
            if (tmp) {
                len[i] -= tmp;
                len[i ^ 1] += tmp;
                res += tmp;
            }
        }
    }
    if (res < lft) {
        dis[u] = -1;
    }
    return res;
}
int dinic() {
    int tmp, res = 0;
    while (bfs(src, snk)) {
        tmp = find(src, inf);
        while (tmp) {
            res += tmp;
            tmp = find(src, inf);
        }
    }
    return res;
}
int main() {
    memset(lnk, -1, sizeof(lnk));
    scanf("%d %d", &n, &m);
    src = 0, snk = n * m + 1;
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= m; j++) {
            scanf("%d", c[i] + j);
        }
    }
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= m; j++) {
            if (c[i][j] == 1) {
                adde(src, id(i, j), inf);
            } else if (c[i][j] == 2) {
                adde(id(i, j), snk, inf);
            }
        }
    }
    for (int i = 1; i < n; i++) {
        for (int j = 1; j <= m; j++) {
            if (c[i][j] == c[i + 1][j] && c[i][j]) {
                continue;
            }
            adde(id(i, j), id(i + 1, j), 1);
        }
    }
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j < m; j++) {
            if (c[i][j] == c[i][j + 1] && c[i][j]) {
                continue;
            }
            adde(id(i, j), id(i, j + 1), 1);
        }
    }
    printf("%d\n", dinic());
    return 0;
}