BZOJ1001 洛谷4001 [BJOI2006]狼抓兔子 题解

题目

这个题目有多种解法，这个题也是一个比较经典的题了，正是因为他的多样的做法，这个题主要难在建图和优化，因为这是一个网格图，所以spfa肯定过不去，所以用最短路解法的话，只能用dij，而网络流也是要加上几个必不可少的优化，具体在代码中解释：

代码：

// luogu-judger-enable-o2
#include <bits/stdc++.h>
#define M 6000100
#define N 200100
#define inf 1 << 30
using namespace std;
queue <int> q;
struct cym {
    int to, next, len;
}e[M * ];
int lin[M * ], cnt = ;
int n, m, s, t, deep[M];
inline void add(int u, int v, int w)
{
    e[++cnt].to = v;
    e[cnt].len = w;
    e[cnt].next = lin[u];
    lin[u] = cnt;
}
inline bool bfs()
{
    int i, now;
    memset(deep,  ,sizeof(deep));
    q.push(s);
    deep[s] = ;
    while(!q.empty())
    {
        int now = q.front();
        q.pop();
        for(int i = lin[now]; i; i = e[i].next)
        {
            if(e[i].len && !deep[e[i].to])
            {
                deep[e[i].to] = deep[now] + ;
                q.push(e[i].to);
            }
        }
    }
    if(deep[t] > )
        return ;
    return ;
}
int dfs(int x, int maxflow)//maxflow是当前最多可以流多少流量
{
    if(x == t || maxflow == )
        return maxflow;
    int w, used = ;
    for(int i = lin[x]; i; i = e[i].next)
    {
        if(deep[e[i].to] == deep[x] +  && e[i].len)
        {
            w = dfs(e[i].to, min(maxflow - used, e[i].len));//w=这条边权的最小值或当前最多可以流的流量减去当前增广的流量
            if(!w)//如果当前不能走了，就把这条路封掉，封掉的方法就是把层数设为0
            {
                deep[e[i].to] = ;
                continue;
            }
            used += w;//当前可以流的流量
            e[i].len -= w;
            e[i ^ ].len += w;
            if(used == maxflow)//相当于当前用光了 ，没有可以增广的流量了
                return used;
        }
    }
    if(!used)//used 是现在可以增广的流量
        deep[x] = ;
    return used;
}
int dinic()
{
    int ans = ;
    while(bfs())
    {
        ans += dfs(s, inf);
    }
    return ans;
}
inline int read() {
    char ch = getchar(); int x = , f = ;
    while(ch < '' || ch > '') {
        if(ch == '-') f = -;
        ch = getchar();
    } while('' <= ch && ch <= '') {
        x = x *  + ch - '';
        ch = getchar();
    } return x * f;
}
int main()
{
    n = read(), m = read();
    s = , t = n * m;
    for (int i = ; i <= n; i++)
        for (int j = ; j < m; j++)
        {
            int w;
            w = read();
            add ( (i - ) * m + j, (i - ) * m + j + , w);
            add ( (i - ) * m + j + , (i - ) * m + j, w);
        }
    for (int i = ; i < n; i++)
        for (int j = ; j <= m; j++)
        {
            int w;
            w = read();
            add ( (i - ) * m + j, (i - ) * m + j + m, w);
            add ( (i - ) * m + j + m, (i - ) * m + j, w);
        }
    for (int i = ; i < n; i++)
        for (int j = ; j < m; j++)
        {
            int w;
            w = read();
            add ( (i - ) * m + j, (i - ) * m + j + m + , w);
            add ( (i - ) * m + j + m + , (i - ) * m + j, w);
        }
    /*
    for(int i = 1; i <= cnt; i++)
    {
        printf("%d %d %d\n", e[i].to, e[i].len, e[i].next);
    }*/
    printf("%d", dinic());
    return ;
}