1. 程式人生 > >BZOJ1001: [BeiJing2006]狼抓兔子

BZOJ1001: [BeiJing2006]狼抓兔子

div tar name family 設有 amp 兔子 ast mil

【傳送門:BZOJ1001


簡要題意:

有一個n*m大小的矩陣,假設有一個點(x,y),那麽這個點與(x+1,y)、(x,y+1)、(x+1,y+1)三條邊都連有一條有流量單向邊,且方向是以(x,y)為起點的,求出從(1,1)為起點,並且起點一開始有無限的可流量,以(n,m)為終點的最大流量


題解:

這就是一道典型的網絡流,可以說是很裸,所以只要建好邊,直接網絡流就可以了


參考程序:

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include
<cmath> using namespace std; struct node { int x,y,c,next,other; }a[6100000];int len,last[1100000]; void ins(int x,int y,int c) { int k1=++len,k2=++len; a[k1].x=x;a[k1].y=y;a[k1].c=c; a[k1].next=last[x];last[x]=k1; a[k2].x=y;a[k2].y=x;a[k2].c=c; a[k2].next=last[y];last[y]=k2; a[k1].other
=k2; a[k2].other=k1; } int h[1100000],list[1100000],head,tail; int st,ed; bool bfs() { memset(h,0,sizeof(h)); h[st]=1;head=1;tail=2; memset(list,0,sizeof(list)); list[1]=st; while(head!=tail) { int x=list[head]; for(int k=last[x];k;k=a[k].next) {
int y=a[k].y; if(h[y]==0&&a[k].c>0) { h[y]=h[x]+1; list[tail++]=y; } } head++; } if(h[ed]>0) return true; return false; } int dfs(int x,int flow) { if(x==ed) return flow; int t=0,minf; for(int k=last[x];k;k=a[k].next) { int y=a[k].y; if(h[y]==(h[x]+1)&&flow>t&&a[k].c>0) { minf=dfs(y,min(a[k].c,flow-t)); t+=minf; a[k].c-=minf;a[a[k].other].c+=minf; } } if(t==0) h[x]=0; return t; } int main() { int n,m; scanf("%d%d",&n,&m); st=1;ed=n*m; len=0;memset(last,0,sizeof(last)); for(int i=1;i<=n;i++) { for(int j=1;j<m;j++) { int c; scanf("%d",&c); ins((i-1)*m+j,(i-1)*m+j+1,c); } } for(int i=1;i<n;i++) { for(int j=1;j<=m;j++) { int c; scanf("%d",&c); ins((i-1)*m+j,i*m+j,c); } } for(int i=1;i<n;i++) { for(int j=1;j<m;j++) { int c; scanf("%d",&c); ins((i-1)*m+j,i*m+j+1,c); } } int ans=0; while(bfs()) { ans+=dfs(st,999999999); } printf("%d\n",ans); return 0; }

BZOJ1001: [BeiJing2006]狼抓兔子