bzoj 1001 狼抓兔子 网络流

现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到，但抓兔子还是比较在行的，

而且现在的兔子还比较笨，它们只有两个窝，现在你做为狼王，面对下面这样一个网格的地形：

 

左上角点为(1,1),右下角点为(N,M)(上图中N=4,M=5).有以下三种类型的道路 

1:(x,y)<==>(x+1,y) 

2:(x,y)<==>(x,y+1) 

3:(x,y)<==>(x+1,y+1) 

道路上的权值表示这条路上最多能够通过的兔子数，道路是无向的. 左上角和右下角为兔子的两个窝，

开始时所有的兔子都聚集在左上角(1,1)的窝里，现在它们要跑到右下解(N,M)的窝中去，狼王开始伏击

这些兔子.当然为了保险起见，如果一条道路上最多通过的兔子数为K，狼王需要安排同样数量的K只狼，

才能完全封锁这条道路，你需要帮助狼王安排一个伏击方案，使得在将兔子一网打尽的前提下，参与的

狼的数量要最小。因为狼还要去找喜羊羊麻烦.

Input

第一行为N,M.表示网格的大小，N,M均小于等于1000.

接下来分三部分

第一部分共N行，每行M-1个数，表示横向道路的权值. 

第二部分共N-1行，每行M个数，表示纵向道路的权值. 

第三部分共N-1行，每行M-1个数，表示斜向道路的权值. 

输入文件保证不超过10M

Output

输出一个整数，表示参与伏击的狼的最小数量.

Sample Input 3 45 6 44 3 17 5 35 6 7 88 7 6 55 5 56 6 6

Sample Output 14

题意&&题解：题目的意思很明确了，就是求从起点到终点能走多少头羊，每一条边都有一个权值代表了最多能通过多少头羊。

裸的网络流，直接跑就行了，建边需要稍微想一想，也很容易。弱鸡的第一发网络流。。。菜的一批。还有一种做法是建对偶图，然后跑最短路，表示对偶图不会建2333，日后再补（应该是不会了）吧。。。

AC代码：

#include <bits/stdc++.h>
using namespace std;
const int INF = 1e9;
const int maxn = 4000000,maxnm = 5e5;

struct Edge{
    int v,f,nxt;
    Edge(){}
    Edge(int _v,int _f,int _nxt):v(_v),f(_f),nxt(_nxt){}
};
int n,scr,sink;
int head[maxn + 10];
int nume;
Edge e[maxn * 2 + 10];

void addedge(int u,int v,int c){
    e[++nume].v = v;
    e[nume].f = c;
    e[nume].nxt = head[u];
    head[u] = nume;

    e[++nume].v = u;
    e[nume].f = c;
    e[nume].nxt = head[v];
    head[v] = nume;
}

void init(){
    memset(head,0,sizeof(head));
    nume = 1;
}

queue <int> que;
bool vis[maxn + 10];
int dis[maxn + 10];

bool bfs(){
    memset(dis,0,sizeof(dis));
    memset(vis,false,sizeof(vis));
    while(!que.empty()) que.pop();
    vis[scr] = true;
    que.push(scr);
    while(!que.empty()){
        int u = que.front();
        que.pop();
        for(int i=head[u];i;i=e[i].nxt){
            if(e[i].f && !vis[e[i].v]){
                que.push(e[i].v);
                dis[e[i].v] = dis[u] + 1;
                vis[e[i].v] = true;
            }
        }
    }
    return dis[sink] != 0;
}

int dfs(int u,int delta){
    if(u == sink || !delta){
        return delta;
    }else{
        int ret = 0;
        for(int i=head[u];delta && i;i=e[i].nxt){
            if(e[i].f && dis[e[i].v] == dis[u] + 1){
                int dd = dfs(e[i].v,min(e[i].f,delta));
                e[i].f -= dd;
                e[i ^ 1].f += dd;
                delta -= dd;
                ret += dd;
                if(!delta) return ret;
            }
        }
        if(!ret) dis[u] = -1;
        return ret;
    }
}

int maxflow(){
    int ret = 0;
    while(bfs()){
        ret += dfs(scr,INF);
    }
    return ret;
}
int m;
int main(){
    scanf("%d%d",&n,&m);
    int cnt = 0;
    init();
    for(int i=1;i<=n;i++){
        cnt++;
        for(int j=1;j<m;j++){
            int v;
            scanf("%d",&v);
            addedge(cnt,cnt+1,v);
            cnt += 1;
        }
    }
    cnt = 0;
    for(int i=1;i<n;i++){
        cnt = cnt == 0 ? 1 : cnt;
        for(int j=1;j<=m;j++){
            int v;
            scanf("%d",&v);
            addedge(cnt,cnt+m,v);
            cnt += 1;
        }
    }
    cnt = 0;
    for(int i=1;i<n;i++){
        cnt++;
        for(int j=1;j<m;j++){
            int v;
            scanf("%d",&v);
            addedge(cnt,cnt+m+1,v);
            cnt += 1;
        }
    }
    scr = 1;
    sink = n * m;
    int ans = maxflow();
    printf("%d\n",ans);
    return 0;
}