hdu 3549 Flow Problem（最大流入门）

Problem Description

Network flow is a well-known difficult problem for ACMers. Given a graph, your task is to find out the maximum flow for the weighted directed graph.

Input

The first line of input contains an integer T, denoting the number of test cases.
For each test case, the first line contains two integers N and M, denoting the number of vertexes and edges in the graph. (2 <= N <= 15, 0 <= M <= 1000)
Next M lines, each line contains three integers X, Y and C, there is an edge from X to Y and the capacity of it is C. (1 <= X, Y <= N, 1 <= C <= 1000)

Output

For each test cases, you should output the maximum flow from source 1 to sink N.

Sample Input

2

3 2

1 2 1

2 3 1

3 3

1 2 1

2 3 1

1 3 1

Sample Output

Case 1: 1

Case 2: 2

解题思路：Dicic实现，即每个阶段先进行1次bfs给图分层，然后在该图上进行1次或多次寻找增广路，如果当前层次图中已找不到增广路，就重新给图分层，然后继续找增广路，只要t的level小于0就说明当前网络已达到最大流。时间复杂度大概为O(|E||V|²)。

AC代码一(312ms)：

#include<bits/stdc++.h>
using namespace std;
const int INF=0x3f3f3f3f;
const int maxn=1005;//int cnt=1;
struct edge{ int to, cap; size_t rev;
    edge(int _to, int _cap, size_t _rev) :to(_to), cap(_cap), rev(_rev){}//构造函数，初始化结构体变量
};//指向节点to，边容量是cap，rev是记录为当前邻接点to反向边的编号
vector<edge> G[maxn];//邻接表，G[i][j]表示节点i连接的第j条边包含的所有信息
int t,n,m,x,y,c,level[maxn];//level数组在bfs时为分层图所用
void add_edge(int from,int to,int cap){//向图中增加一条从s到t容量为cap的边
    G[from].push_back(edge( to, cap, G[to].size() ));
    G[to].push_back(edge( from, 0, G[from].size() - 1));//关键：反向建边，初始流量为0
}
//bfs给图分层次
void bfs(int s){
    memset(level,-1,sizeof(level));//刚开始每个节点的层次置为-1
    queue<int> que;//队列实现bfs
    level[s]=0;//源点s为第0层
    que.push(s);
    while(!que.empty()){//给图分层
        int v=que.front();que.pop();
        for(size_t i=0;i<G[v].size();++i){//遍历节点v与之相连的每条边
            edge &e=G[v][i];//取出与节点v相连的第i条边
            if(e.cap>0&&level[e.to]<0){//如果边残余流量大于0，且节点e.to还未分层
                level[e.to]=level[v]+1;//节点e.to的层次为指向它的节点v所在层次数加1
                que.push(e.to);
            }
        }
    }
}
//dfs寻找增广路，寻找当前图中s-->t的一条增广路
int dfs(int v,int t,int f){//v->t（t为汇点），当前增广路径上的最小剩余流量为f
    if(v==t)return f;//到达汇点t
    for(size_t i=0;i<G[v].size();++i){//遍历节点v与之相连的每条边
        edge &e=G[v][i];//取出与节点v相连的第i条边
        if(e.cap>0 && level[v]<level[e.to]){//如果该边残流量大于0，且邻接点e.to是v的下一级，就增广下去
            //cout<<"当前遍历到的点v："<<v<<",邻接点to："<<e.to<<endl;
            int d=dfs(e.to,t,min(f,e.cap));//维护当前增广路上最小的剩余流量f
            if(d>0){//若f>0，说明找到了一条增广路
                //cout<<v<<"-->"<<e.to<<",边容量为:"<<e.cap<<endl;
                e.cap-=d;//正向边流量减去f
                G[e.to][e.rev].cap+=d;//反向边流量加上f
                return d;//沿着增广路回溯到源点s，不会在途中去深搜其他点，返回当前增广路上的最小剩余流量
            }
        }
    }
    return 0;//否则说明没有增广路，返回0
}
//Dinic算法实现最大流，每个阶段执行完一次bfs分层之后，只需查找当前层次图中是否还增广路即可
int max_flow(int s,int t){
    int flow=0;
    while(1){
        bfs(s);//每个阶段先bfs将图分层标记
        if(level[t]<0)return flow;
        //如果分层之后，汇点t的层次小于0，即未被分层，说明再无增广路，则直接返回当前最大流量
        int f=dfs(s,t,INF);
        //cout<<"第"<<cnt++<<"次的最小剩余容量为:"<<f<<endl;
        while(f>0){//在该层次图中找到增广路
            flow+=f;//累加到最大流中
            f=dfs(s,t,INF);//继续找该层次图中是否还有增广路，直到f为0，
            //cout<<"第"<<cnt++<<"次的最小剩余容量为:"<<f<<endl;
        }
    }
}
int main(){
    while(~scanf("%d",&t)){
        for(int cas=1;cas<=t;++cas){
            scanf("%d%d",&n,&m);//cnt=1;
            for(int i=0;i<=n;++i)G[i].clear();
            while(m--){
                scanf("%d%d%d",&x,&y,&c);
                add_edge(x,y,c);
            }
            printf("Case %d: %d
",cas,max_flow(1,n));
        }
    }
    return 0;
}

AC代码二(93ms)：当前弧优化Dinic算法。

#include<bits/stdc++.h>
using namespace std;
const int INF=0x3f3f3f3f;
const int maxn=1005;
struct edge{ int to,cap;size_t rev;
    edge(int _to, int _cap, size_t _rev):to(_to),cap(_cap),rev(_rev){}
};
int T,n,m,x,y,c,level[maxn];queue<int> que;vector<edge> G[maxn];size_t curfir[maxn];
void add_edge(int from,int to,int cap){
    G[from].push_back(edge(to,cap,G[to].size()));
    G[to].push_back(edge(from,0,G[from].size()-1));
}
bool bfs(int s,int t){
    memset(level,-1,sizeof(level));
    while(!que.empty())que.pop();
    level[s]=0;
    que.push(s);
    while(!que.empty()){
        int v=que.front();que.pop();
        for(size_t i=0;i<G[v].size();++i){
            edge &e=G[v][i];
            if(e.cap>0&&level[e.to]<0){
                level[e.to]=level[v]+1;
                que.push(e.to);
            }
        }
    }
    return level[t]<0?false:true;
}
int dfs(int v,int t,int f){
    if(v==t)return f;
    for(size_t &i=curfir[v];i<G[v].size();++i){
        edge &e=G[v][i];
        if(e.cap>0&&(level[v]+1==level[e.to])){
            int d=dfs(e.to,t,min(f,e.cap));
            if(d>0){
                e.cap-=d;
                G[e.to][e.rev].cap+=d;
                return d;
            }
        }
    }
    return 0;
}
int max_flow(int s,int t){
    int f,flow=0;
    while(bfs(s,t)){
        memset(curfir,0,sizeof(curfir));
        while((f=dfs(s,t,INF))>0)flow+=f;
    }
    return flow;
}
int main(){
    while(~scanf("%d",&T)){
        for(int cas=1;cas<=T;++cas){
            scanf("%d%d",&n,&m);
            for(int i=0;i<=n;++i)G[i].clear();
            while(m--){
                scanf("%d%d%d",&x,&y,&c);
                add_edge(x,y,c);
            }
            printf("Case %d: %d
",cas,max_flow(1,n));
        }
    }
    return 0;
}