BZOJ 1001狼抓兔子（最大流）

Description


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


 


左上角点为(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 4


5 6 4


4 3 1


7 5 3


5 6 7 8


8 7 6 5


5 5 5


6 6 6
Sample Output


14


我用的是网络流大暴力，这道题可以转为对偶图跑最短路。
下面贴一份最大流过了的代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
using namespace std;

#define maxn (1000000+1)
int N,M,inf=1000000,S,T,ans=1000000,col=0;
int dep[maxn],head[maxn];

struct node{
	int u,v,next,c;
}e[6000001];//数组要开够,我第一次就是数组开小了导致RE 

int k=0;
void adde(int u,int v,int c){
	e[k].u=u; e[k].v=v; e[k].c=c; e[k].next=head[u]; head[u]=k++;
	e[k].u=v; e[k].v=u; e[k].c=c; e[k].next=head[v]; head[v]=k++;
}

void input(){
	memset(head,-1,sizeof(head));
	scanf("%d%d",&N,&M);
	int x;
	if(N==1||M==1)col=1;
	for(int i=1;i<=N;i++)
		for(int j=1;j<M;j++){//横向边 
			scanf("%d",&x);
			ans=min(ans,x);
			adde((i-1)*M+j,(i-1)*M+j+1,x);
		}
	for(int i=1;i<N;i++)
		for(int j=1;j<=M;j++){
			scanf("%d",&x);
			ans=min(ans,x);
			adde((i-1)*M+j,i*M+j,x);//纵向边 
		}
	for(int i=1;i<N;i++)
		for(int j=1;j<M;j++){
			scanf("%d",&x);
			ans=min(ans,x);
			adde((i-1)*M+j,i*M+j+1,x);//对角线 
		}
	S=1; T=N*M;
}

bool bfs(){
	queue<int>q;
	memset(dep,-1,sizeof(dep));
	dep[S]=0; q.push(S);
	while(!q.empty()){
		int u=q.front();
		q.pop();
		
		for(int i=head[u];i!=-1;i=e[i].next){
			int v=e[i].v;
			if(dep[v]==-1&&e[i].c>0){
				dep[v]=dep[u]+1;
				q.push(v);
			}
		}
	}
	return dep[T]!=-1;
}

int dfs(int a,int b){
	if(a==T)return b;
	int ret=0;
	for(int i=head[a];i!=-1;i=e[i].next){
		int v=e[i].v;
		if(dep[v]==dep[a]+1&&e[i].c>0){
			int x=dfs(v,min(b-ret,e[i].c));
			e[i].c-=x;
			e[i^1].c+=x;
			ret+=x;
			if(ret==b)return ret;
		}
	}
	if(ret==0)dep[a]=-1;//一个小优化
	return ret;
}

int dinic(){
	int ret=0;
	while(bfs()){
		ret+=dfs(S,inf);
	}
	return ret;
}

void solve(){
	if(col==0)ans=dinic();
	printf("%d",ans);
}

int main(){
	input();
	solve();
	return 0;
}