最大流Edmonds-Karp模板
Posted 033000-
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了最大流Edmonds-Karp模板相关的知识,希望对你有一定的参考价值。
const int maxn=1e4+5; struct Edge{ int from,to,cap,flow; Edge(int u,int v,int c,int f):from(u),to(v),cap(c),flow(f){} }; struct EdmodsKarp{ int n,m; vector<Edge> edges; vector<int>G[maxn];//邻接表,G[i][j]表示节点i的第j条边在edges数组中的需要 int a[maxn];//当起点到i的可改进量 int p[maxn];//最短路树上p的入弧编号 void init(int n){ for(int i=0;i<n;i++) G[i].clear(); edges.clear(); } void AddEdge(int from,int to,int cap){ edges.push_back(Edge(from,to,cap,0)); edges.push_back(Edge(to,from,0,0));//反向弧 m=edges.size(); G[from].push_back(m-2); G[to].push_back(m-1); } int Maxflow(int s,int t){ int flow=0; for(;;){ //开始寻找增广路 memset(a,0,sizeof(a)); queue<int> Q; Q.push(s); a[s]=inf; while(!Q.empty()){ int x=Q.front();Q.pop(); for(int i=0;i<G[x].size();i++){ Edge &e=edges[G[x][i]]; if(!a[e.to]&&e.cap>e.flow){ p[e.to]=G[x][i]; a[e.to]=min(a[x],e.cap-e.flow); Q.push(e.to); } } if(a[t]) break;//找到了一条增广路 } if(!a[t]) break;//不存在增广路,当前流就是最大流 for(int u=t;u!=s;u=edges[p[u]].from){ edges[p[u]].flow+=a[t]; edges[p[u]^1].flow-=a[t]; } flow+=a[t]; } return flow; } };
以上是关于最大流Edmonds-Karp模板的主要内容,如果未能解决你的问题,请参考以下文章
Poj 1273 Drainage Ditches(最大流 Edmonds-Karp )