HDU 3046 Pleasant sheep and big big wolf | 最小割

刷了一周的dp，现在又回来到图论上来。之前只是学了最大流，而现在学的最小割等同于最大流。因此两类题目的求法是基本不变的，需要的是思维上的改变。

可以构成二分图的图，它的最小点权覆盖，最大点权独立集都是和最小割，最大流紧密联系，详情可以看我转载的一篇文章《网络流之--最小点权覆盖和最大点权独立集》。

————————————————————————————————————————————————————————————————

题意：

方格上有羊和狼，问要最少要建多少个fence，才能使狼无法干扰到羊。

————————————————————————————————————————————————————————————————

方法：

格子与其周围的鸽子（有公共边）连一条容量为1的边，有狼的格子与源点s相连，有羊的格子与汇点t相连，求出最大流即是题目要求的答案。

一个grid最多能容纳一个animal，于是边的容量为1，总感觉有点不太对劲。一般来说，一只狼走了那个格子，但先来的狼可以走开，然后让另一只狼进来，那容量为1就有问题。但回到题目的本意，是让羊不受到狼的干扰。那多只狼能到达的格子与一只狼能到达没什么区别。

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <queue>
#include <vector>
#define debug cout<<"debug"<<endl;
using namespace std;
const int MAXN = 40050;
const int MAXEDGE = 500005;
const int INF = 0x3F3F3F3F;
int n, m;
struct Edge
{
        int to, cap, flow, next;
};
Edge edge[MAXEDGE];
struct Dinic
{
        int pp, s, t;
        int d[MAXN], cur[MAXN];
        bool vis[MAXN];
        int head[MAXN];
        Dinic(int ts, int tt)
        {
                s = ts, t = tt;
                pp = 0;
                memset(head, -1, sizeof(head));
        }
        void AddEdge(int u, int v, int c)
        {
                edge[pp] = (Edge){v, c, 0, head[u]};
                head[u] = pp++;

                edge[pp] = (Edge){u, 0, 0, head[v]};
                head[v] = pp++;
        }

        bool bfs()
        {
                memset(vis, false, sizeof(vis));
                queue <int> q;
                q.push(s);
                d[s] = 0;
                vis[s] = true;
                while(!q.empty())
                {
                        int u = q.front();
                        q.pop();
                        int next = head[u];
                        while(next != -1)
                        {
                                Edge &e = edge[next];
                                if(!vis[e.to] && e.cap > e.flow)
                                {
                                        vis[e.to] = true;
                                        q.push(e.to);
                                        d[e.to] = d[u] + 1;
                                }
                                next = e.next;
                        }
                }
                return vis[t];
        }

        int dfs(int u,int a)
        {
                if(u == t || a == 0)    return a;
                int flow = 0, f;
                int &next = cur[u];
                while(next != -1)
                {
                        Edge &e = edge[next];
                        if(d[e.to] == d[u] + 1 && (f = dfs(e.to, min(a, e.cap - e.flow))) > 0)
                        {
                                e.flow += f;
                                edge[next^1].flow -= f;
                                flow += f;
                                a -= f;
                                if(a == 0)      break;
                        }
                        next = e.next;
                }
                return flow;
        }
        int getMaxFlow()
        {
                int res = 0;
                while(bfs())
                {
                        for(int i = 0;i < n*m + 2; i++)
                        {
                                cur[i] = head[i];
                        }
                        res += dfs(s, INF);
                }
                return res;
        }
};
int main()
{
        int cas = 0;
        while(scanf("%d%d", &n, &m) != EOF)
        {
                int tmp, s = 0, t = n * m + 1;
                Dinic dinic(s, t);
                for(int i = 0;i < n; i++)
                {
                        for(int j = 1;j <= m; j++)
                        {
                                scanf("%d",&tmp);
                                if(i != 0)      dinic.AddEdge( i*m + j, (i-1)*m + j, 1);
                                if(i != n-1)    dinic.AddEdge( i*m + j, (i+1)*m + j, 1);
                                if(j != 1)      dinic.AddEdge( i*m + j, i * m + j - 1, 1);
                                if(j != m)      dinic.AddEdge( i*m + j, i * m + j + 1, 1);
                                if(tmp == 1)    dinic.AddEdge( i*m + j, t, INF);
                                if(tmp == 2)    dinic.AddEdge( s, i*m + j, INF);
                        }
                }
                printf("Case %d:\n%d\n", ++cas, dinic.getMaxFlow());
        }
        return 0;
}