BZOJ 1001 [BeiJing2006]狼抓兔子
题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1001
题意: ...
很容易想到求的是一个最小割=最大流。
之前一直用的刘汝佳的模板STL过题，很久没用过数组模拟了。
再次熟悉一下写法，first数组是索引数组，标记的结点的最后一条边，利用next数组找到上一条边。
1 #include <bits/stdc++.h> 2 3 using namespace std; 4 5 #define N 1100000 6 int n,m,first[N],next[N*6],v[N*6],w[N*6],tot,jy,ans,vis[N],S,T,xx,pos[N]; 7 8 void Add(int x,int y,int z) 9 { 10 w[tot] = z; 11 v[tot] = y; 12 next[tot] = first[x]; 13 first[x] = tot++; 14 } 15 16 void add(int x,int y,int z) 17 { 18 Add(x,y,z); 19 Add(y,x,z); 20 } 21 22 bool bfs() 23 { 24 memset(vis,-1,sizeof(vis)),vis[S]=0; 25 queue<int>q; 26 q.push(S); 27 while(!q.empty()) 28 { 29 int t=q.front(); 30 q.pop(); 31 for(int i=first[t]; ~i; i=next[i]) 32 if(w[i]&&vis[v[i]]==-1) 33 vis[v[i]]=vis[t]+1,q.push(v[i]); 34 } 35 return vis[T]!=-1; 36 } 37 38 int dfs(int x,int y) 39 { 40 if(x==T)return y; 41 int r=0; 42 for(int i=first[x]; ~i&&y>r; i=next[i]) 43 if(w[i]&&vis[v[i]]==vis[x]+1) 44 { 45 int t=dfs(v[i],min(y-r,w[i])); 46 w[i]-=t,w[i^1]+=t,r+=t; 47 } 48 if(!r)vis[x]=-1; 49 return r; 50 } 51 52 int main() 53 { 54 memset(first,-1,sizeof(first)); 55 scanf("%d%d",&n,&m); 56 S=m+1,T=n*m+m; 57 for(int i=1; i<=n; i++) 58 for(int j=1; j<m; j++) 59 scanf("%d",&xx),add(i*m+j,i*m+j+1,xx); 60 for(int i=1; i<n; i++) 61 for(int j=1; j<=m; j++) 62 scanf("%d",&xx),add(i*m+j,(i+1)*m+j,xx); 63 for(int i=1; i<n; i++) 64 for(int j=1; j<m; j++) 65 scanf("%d",&xx),add(i*m+j,(i+1)*m+j+1,xx); 66 67 while(bfs()) 68 while(jy=dfs(S,0x3fffffff)) 69 ans+=jy; 70 printf("%d\n",ans); 71 }
