P4001 [ICPC-Beijing 2006]狼抓兔子（ISAP模板题）

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

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

(x,y)\rightleftharpoons(x+1,y)(x,y)⇌(x+1,y)

(x,y)\rightleftharpoons(x,y+1)(x,y)⇌(x,y+1)

(x,y)\rightleftharpoons(x+1,y+1)(x,y)⇌(x+1,y+1)

道路上的权值表示这条路上最多能够通过的兔子数，道路是无向的。左上角和右下角为兔子的两个窝，开始时所有的兔子都聚集在左上角 (1,1)(1,1) 的窝里，现在它们要跑到右下角 (N,M)(N,M) 的窝中去，狼王开始伏击这些兔子。当然为了保险起见，如果一条道路上最多通过的兔子数为 KK，狼王需要安排同样数量的 KK 只狼，才能完全封锁这条道路，你需要帮助狼王安排一个伏击方案，使得在将兔子一网打尽的前提下，参与的狼的数量要最小。因为狼还要去找喜羊羊麻烦。

输入格式
第一行两个整数 N,MN,M，表示网格的大小。

接下来分三部分。

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

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

第三部分共 N-1N−1 行，每行 M-1M−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
说明/提示
数据规模与约定
对于全部的测试点，保证 3 \leq N,M \leq 10003≤N,M≤1000，所有道路的权值均为不超过 10^610
6
的正整数。

思路：很明显的最大流问题

#include <bits/stdc++.h>
#define inf 0x7fffffff
//#define ll long long
#define int long long
//#define double long double
#define eps 1e-8
//#define mod 1e9+7
using namespace std;
const int mod=1e9+7;
const int M=1e5;
const int N=1e7+5;//空间最大限制 4e8
int n,m,s,t;
struct node
{
	int ver,edge,next;
}e[N];
int tot=1,head[N];
int d[N],gap[N],now[N];
void add(int x,int y,int z)
{
	e[++tot].ver=y;
	e[tot].edge=z;
	e[tot].next=head[x];
	head[x]=tot;
}
void addedge(int x,int y,int z)
{
	add(x,y,z);add(y,x,z);
}
int get(int x,int y)
{
	return (x-1)*m+y;
}
void bfs()
{
	d[t]=gap[1]=1;
	queue < int > q;
	q.push(t);
	while(q.size())
	{	
		int x=q.front();q.pop();
		for(int i=head[x];i;i=e[i].next)
		{
			int y=e[i].ver;
			int z=e[i].edge;
			if(d[y])  continue;
			d[y]=d[x]+1;
			q.push(y); 
			gap[d[y]]++;
		}
	}
}
int dfs(int x,int flow)
{
	if(x==t)  return flow;
	int i,rest=0,k;
	for(int i=head[x];i;i=e[i].next)
	{
		int y=e[i].ver;
		int z=e[i].edge;
		if(!z)  continue;
		if(d[x]!=d[y]+1)  continue;
		k=dfs(y,min(z,flow-rest));
		if(k<=0)  continue;
		e[i].edge-=k;
		e[i^1].edge+=k;
		rest+=k;
		if(flow==rest)  return flow;  
	}
	if(--gap[d[x]]==0)  d[s]=n*m+1;
	++gap[++d[x]];
	return rest;
}
int isap()
{
	int maxflow=0;
	bfs();
	while(d[s]<n*m)  maxflow+=dfs(s,1e18);
	return maxflow;
}
signed main()
{
//	ios::sync_with_stdio(false);
	cin>>n>>m;
	s=1,t=get(n,m);
	for(int i=1;i<=n;i++)  for(int j=1;j<m;j++)
	{
		int x;
		scanf("%lld",&x);
		addedge(get(i,j),get(i,j+1),x);
	}
	for(int i=1;i<n;i++)  for(int j=1;j<=m;j++)
	{
		int x;
		scanf("%lld",&x);
		addedge(get(i,j),get(i+1,j),x);
	}
	for(int i=1;i<n;i++)  for(int j=1;j<m;j++)
	{
		int x;
		scanf("%lld",&x);
		addedge(get(i,j),get(i+1,j+1),x);
	}
	cout<<isap()<<endl;
    return 0;
}