bzoj3396[Usaco2009 Jan]Total flow 水流*

Posted YuanZiming

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj3396[Usaco2009 Jan]Total flow 水流*相关的知识,希望对你有一定的参考价值。

bzoj3396[Usaco2009 Jan]Total flow 水流

题意:

求无环图的最大流。边数≤700。

题解:

管它有没有环。注意本题的节点标号既有大写字母,也有小写字母。

代码:

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <algorithm>
 4 #include <queue>
 5 #define inc(i,j,k) for(int i=j;i<=k;i++)
 6 #define maxn 100
 7 #define INF 0x3fffffff
 8 using namespace std;
 9 
10 struct e{int t,c,n;}es[maxn*80]; int g[maxn],ess;
11 void pe(int f,int t,int c){
12     es[++ess]=(e){t,c,g[f]}; g[f]=ess; es[++ess]=(e){f,0,g[t]}; g[t]=ess;
13     es[++ess]=(e){f,c,g[t]}; g[t]=ess; es[++ess]=(e){t,0,g[f]}; g[f]=ess;
14 }
15 queue<int>q; int h[maxn];
16 bool bfs(int s,int t){
17     while(!q.empty())q.pop(); memset(h,-1,sizeof(h)); q.push(s); h[s]=0;
18     while(!q.empty()){
19         int x=q.front(); q.pop();
20         for(int i=g[x];i;i=es[i].n)if(es[i].c&&h[es[i].t]==-1){h[es[i].t]=h[x]+1; q.push(es[i].t);}
21     }
22     return h[t]!=-1;
23 }
24 int dfs(int x,int t,int f){
25     if(x==t)return f; int u=0;
26     for(int i=g[x];i;i=es[i].n)if(es[i].c&&h[es[i].t]==h[x]+1){
27         int w=dfs(es[i].t,t,min(f,es[i].c)); f-=w; u+=w;
28         es[i].c-=w; es[i^1].c+=w; if(f==0)return u;
29     }
30     if(u==0)h[x]=-1; return u;
31 }
32 int dinic(int s,int t){
33     int f=0; while(bfs(s,t))f+=dfs(s,t,INF); return f;
34 }
35 int s,t,n;
36 int main(){
37     scanf("%d",&n); ess=1; s=0; t=Z-A;
38     inc(i,1,n){char x[3],y[3]; int z; scanf("%s%s%d",x,y,&z); pe(x[0]-A,y[0]-A,z);}
39     printf("%d",dinic(s,t)); return 0;
40 }

 

20160908

以上是关于bzoj3396[Usaco2009 Jan]Total flow 水流*的主要内容,如果未能解决你的问题,请参考以下文章

bzoj 3396: [Usaco2009 Jan]Total flow 水流最大流

[BZOJ1574] [Usaco2009 Jan]地震损坏Damage(贪心 + dfs)

bzoj 1575: [Usaco2009 Jan]气象牛Baricdp

BZOJ3393 [Usaco2009 Jan]Laserphones 激光通讯 BFS

[BZOJ1575] [Usaco2009 Jan]气象牛Baric(DP)

BZOJ1576 [Usaco2009 Jan]安全路经Travel