网络流_洛谷P4001_最小割

题目：
题目描述
现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到，但抓兔子还是比较在行的，而且现在的兔子还比较笨，它们只有两个窝，现在你做为狼王，面对下面这样一个网格的地形：

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

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只狼，才能完全封锁这条道路，你需要帮助狼王安排一个伏击方案，使得在将兔子一网打尽的前提下，参与的狼的数量要最小。因为狼还要去找喜羊羊麻烦。

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

接下来分三部分

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

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

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

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

输入输出样例
输入 #1 复制
3 4
5 6 4
4 3 1
7 5 3
5 6 7 8
8 7 6 5
5 5 5
6 6 6
输出 #1 复制
14

题目思路：
构图：首先建边，然后每一条边是无向边，即2条有向边；然后跑Dinic算法。

代码注意点：
1.s 和 t 需要在 n, m初始化后再赋值
2. 最后要记得 flow = dinic(s, inf) ，别忘了给flow 赋值！！！

AC代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>

using namespace std;

const int N = 1e6 + 5, M = 12000005;
const int inf = 0x3f3f3f3f;

int n, m, s, t, maxflow;
int h[N], w[M], ne[M], e[M], idx;
int d[N];

inline void add(int a, int b, int c) {
	e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx ++;
	e[idx] = a, w[idx] = 0, ne[idx] = h[b], h[b] = idx ++;
}

inline bool bfs(void) {
	queue<int> q;
	while(q.size()) q.pop();
	memset(d, 0, sizeof d);
	
	d[s] = 1, q.push(s);
	
	while(q.size()) {
		int u = q.front(); q.pop();
		
		for(int i = h[u]; i != -1; i = ne[i]) {
			int v = e[i];
			if(w[i] && !d[v]) {
				d[v] = d[u] + 1;
				q.push(v);
			}
		}
	}
	
	if(d[t] == 0) return false;
	return true;
}

inline int dinic(int u, int flow) {
	if(u == t) return flow;
	
	for(int i = h[u]; i != -1; i = ne[i]) {
		int v = e[i];
		if(d[v] == d[u] + 1 && w[i]) {
			int k = dinic(v, min(flow, w[i]));
			if(k) {
				w[i] -= k;
				w[i ^ 1] += k;
				return k;
			} else d[v] = 0;
		}
	}
	
	return 0;
}
int main(void) {
//	freopen("in.txt", "r", stdin);

	scanf("%d%d", &n, &m);
	memset(h, -1, sizeof h);
	s = 1, t = n * m; 
	for(int i = 1; i <= n; i ++) {
		for(int j = 2; j <= m; j ++) {
			int c; scanf("%d", &c);
			int a = (i - 1) * m + j - 1;
			int b = (i - 1) * m + j;
//			printf("%d %d %d\n", a, b, c);
			add(a, b, c); add(b, a, c);
		}
	}
	
	for(int i = 2; i <= n; i ++) {
		for(int j = 1; j <= m; j ++) {
			int c; scanf("%d", &c); 
			int a = (i - 2) * m + j;
			int b = (i - 1) * m + j;
//			printf("%d %d %d\n", a, b, c);
			add(a, b, c); add(b, a, c);
		}
	}
	
	for(int i = 2; i <= n; i ++) {
		for(int j = 2; j <= m; j ++) {
			int c; scanf("%d", &c);
			int a = (i - 2) * m + j - 1;
			int b = (i - 1) * m + j;
//			printf("%d %d %d\n", a, b, c);
			add(a, b, c); add(b, a, c);
		}
	}
	
	int flow;
	while(bfs()) 
	while(flow = dinic(s, inf)) maxflow += flow;
	
	printf("%d\n", maxflow);
	
//	fclose(stdin);
	return 0;
}