[hdu3549]Flow Problem(最大流模板题)

解题关键：使用的挑战程序设计竞赛上的模板，第一道网络流题目，效率比较低，且用不习惯的vector来建图。

看到网上其他人说此题有重边，需要注意下，此问题只在邻接矩阵建图时会出问题，邻接表不会存在的，也体现了邻接表的优越性？

edge结构体的第三个变量为from的下标。

模板一：

#include<bits/stdc++.h>
#define MAX_V 17
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
int n,m,s,t;
struct edge{
    int to,cap,rev;
};
vector<edge>G[MAX_V];
bool used[MAX_V];
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,,G[from].size()-});
}
int dfs(int v,int t,int f){
    if(v==t)return f;
    used[v]=true;
    for(int i=;i<G[v].size();i++){
        edge &e=G[v][i];
        if(!used[e.to]&&e.cap>){
            int d=dfs(e.to,t,min(f,e.cap));
            if(d>){
                e.cap-=d;
                G[e.to][e.rev].cap+=d;
                return d;
            }
        }
    }
    return ;
}
 
int max_flow(int s,int t){
    int flow=;
    while(){
        memset(used,,sizeof used);
        int f=dfs(s,t,inf);
        if(f==) return flow;
        flow+=f;
    }
    return flow;
}
 
int main(){
    int T,u,v,f;
    scanf("%d",&T);
    for(int ca=;ca<=T;ca++){
        memset(G,,sizeof G);
        scanf("%d%d",&n,&m);
        for(int i=;i<m;i++){
            scanf("%d%d%d",&u,&v,&f);
            add_edge(u,v,f);
        }
        s=,t=n;
        int ans=max_flow(s,t);
        printf("Case %d: %d\n",ca,ans);
    }
    return ;
}

模板二：dinic，187ms，比第一个快，在层次图上进行增广，且进行了当前弧优化。

#include<bits/stdc++.h>
#define inf 0x3f3f3f3f
#define MAX_V 17
using namespace std;
typedef long long ll;
struct edge{int to,cap,rev;};//终点，容量，反向边
vector<edge>G[MAX_V];
int level[MAX_V],iter[MAX_V];
int n,m,s,t;
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,,G[from].size()-});
}
void bfs(int s){
    memset(level,-,sizeof level);
    queue<int>que;
    level[s]=;
    que.push(s);
    while(!que.empty()){
        int v=que.front();que.pop();
        for(int i=;i<G[v].size();i++){
            edge &e=G[v][i];
            if(e.cap>&&level[e.to]<){
                level[e.to]=level[v]+;
                que.push(e.to);
            }
        }
    }
}
 
int dfs(int v,int t,int f){
    if(v==t) return f;
    for(int &i=iter[v];i<G[v].size();i++){
        edge &e=G[v][i];
        if(e.cap>&&level[v]<level[e.to]){
            int d=dfs(e.to,t,min(f,e.cap));
            if(d>){
                e.cap-=d;
                G[e.to][e.rev].cap+=d;
                return d;
            }
        }
    }
    return ;
}
 
int max_flow(int s,int t){
    int flow=;
    while(){
        bfs(s);
        if(level[t]<) return flow;
        memset(iter,,sizeof iter);
        int f;
        while((f=dfs(s,t,inf))>){
            flow+=f;
        }
    }
    return flow;
}
 
int main(){
    int T,u,v,f;
    scanf("%d",&T);
    for(int ca=;ca<=T;ca++){
        memset(G,,sizeof G);
        memset(iter,,sizeof iter);
        scanf("%d%d",&n,&m);
        for(int i=;i<m;i++){
            scanf("%d%d%d",&u,&v,&f);
            add_edge(u,v,f);
        }
        s=,t=n;
        int ans=max_flow(s,t);
        printf("Case %d: %d\n",ca,ans);
    }
    return ;
}