最近开始学习图论,到了网络流部分 。
该部分比较难,难度在于有点看不懂, 不过这没有关系,慢慢看,先刷几道水题,将模板中需要注意的地方考虑好 。 之后再刷紫书提高 。 逐步加深对网络流的理解 。
该题非常简单,套上模板就过了 。
我第一次写最大流模板,理解的还不深刻,下面的注释是我对代码的理解。
细节参见代码:
#include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #include<vector> #include<queue> using namespace std; const int maxn = 205; const int INF = 1000000000; int n,nn,b,c,bb,m,a[maxn],p[maxn]; struct Edge{ int from,to,cap,flow; Edge(int u,int v,int c,int f):from(u),to(v),cap(c),flow(f) {} }; vector<Edge> edges; vector<int> g[maxn]; void init(int n) { for(int i=0;i<n;i++) g[i].clear(); edges.clear(); } void addedge(int from,int to,int cap) { //增加边并将每个节点对应的边保存在g中 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; //最大流初始化为0 for(;;) { //核心算法,需要注意,我们一开始加进来的边的流量都是0,通过求最小残量逐步增广,更新最大流 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; //终点残量为0,不能再增广,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; } int main() { while(~scanf("%d%d",&n,&nn)) { init(nn); for(int i=0;i<n;i++) { scanf("%d%d%d",&bb,&b,&c); addedge(bb,b,c); } printf("%d\n",maxflow(1,nn)); } return 0; }