HDU 1532 Drainage Ditches(基础最大流)

最近开始学习图论,到了网络流部分 。

该部分比较难,难度在于有点看不懂, 不过这没有关系,慢慢看,先刷几道水题,将模板中需要注意的地方考虑好   。   之后再刷紫书提高 。 逐步加深对网络流的理解 。 

该题非常简单,套上模板就过了 。

我第一次写最大流模板,理解的还不深刻,下面的注释是我对代码的理解。

细节参见代码:

#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;
}


你可能感兴趣的:(图论,最大流)