网络流-最大流 DINIC算法及例题luogu P3376

DINIC是建立在增广路算法上的一种优化算法。
它的优点有两个：
多路增广 ：
在找到一条增广路后，继续利用它剩余的流量找增广路，以达到在一次DFS中找到多条增广路，提高算法效率。
当前弧优化 ：
标记当前路是否被访问过，若被访问过则跳过，这样可以省去原本EK算法中将多余处理的数据。
DINIC算法可分为两部分，一部分是BFS广搜寻找分层图（判断是否还有到达终点的路径即增广路存在的最基本条件），另一部分是DFS深搜查找增广路的值。

BFS部分代码

bool bfs() {
    queue<int> Q;
    memset(dep,0,sizeof(dep)); //用dep数组来存储点的层次
    dep[s]=1; //起点在第一层
    Q.push(s);
    while(!Q.empty()) {
        int u=Q.front();Q.pop();
        for(int i=head[u];i!=-1;i=edge[i].next) { //链式星向前的遍历
            int v=edge[i].v,f=edge[i].w;
            if(!dep[v] && f) { //如果当前点的层次还是0（即未被访问过）并且剩余流量大于0则符合存在增广路的基本条件
                dep[v]=dep[u]+1;
                Q.push(v);
            }
        }
    }
    return dep[t]!=0;
}

DFS部分代码

int dfs(int u,int flow) {
    if(u==t) return flow; //当前点为终点，那么此时的最小流量就等于这条增广路的值，将其返回
    for(int i=head[u];i!=-1;i=edge[i].next) {
        int v=edge[i].v;
        int f=edge[i].w;
        if(dep[v]==dep[u]+1 && f) {
            int df=dfs(v,min(f,flow)); //用DFS算法对下一个点进行增广
            if(df) { //若增广成功，更新路上的剩余流量，正向路要减，反向路要加。具体原因详见增广路算法
                edge[i].w-=df;
                edge[i^1].w+=df;
                return df;
            }
        }
    }
    return 0;
}

在得出以上的BFS部分和DFS部分后，不难发现我们只需里用BFS去寻找分层图，如果存在能到达终点分层图，那么就有可能存在增广路，然后再用DFS去求得增广路的值。一直重复利用 BFS->DFS->BFS->DFS->… 的步骤直到BFS找不到分层图即不满足增广路存在的基本条件为止。
那么只需再写出一个dinic函数去调用这两个函数即可，dinic算法也就显现了出来。在学习过最基本的增广路算法后再这样分步去理解dinic就简单易懂了。

DINIC算法调用函数

int dinic() {
    int ans=0;
    while(bfs()) {
        int delta;
        do {
            delta=dfs(s,0x3f3f3f3f);
            ans+=delta;
        }while(delta);
    }
    return ans;
}

最后来看一下最大流的例题 luogu P3376 : https://www.luogu.com.cn/problem/P3376

题目描述

如题，给出一个网络图，以及其源点和汇点，求出其网络最大流。

输入格式

第一行包含四个正整数 n,m,s,t，分别表示点的个数、有向边的个数、源点序号、汇点序号。 接下来M行每行包含三个正整数 u_i, v_i,
w_i ，表示第 i 条有向边从 u_i 出发，到达 v_i，边权为 w_i（即该边最大流量为 w_i ）。

输出格式

一行，包含一个正整数，即为该网络的最大流。

这道题是最大流的模板例题，套入上述代码即可，在下面我放出完整的dinic算法代码，希望帮助到各位理解dinic的解题。

#include<bits/stdc++.h>
using namespace std;
int cnt=0,n,m;
int head[10005];
int dep[10005];
int s,t;
struct Edge {
    int v,w,next;
}edge[200005];
void add(int u,int v,int w) {
    edge[cnt].v=v;
    edge[cnt].w=w;
    edge[cnt].next=head[u];
    head[u]=cnt++;
}
bool bfs() {
    queue<int> Q;
    memset(dep,0,sizeof(dep));
    dep[s]=1;
    Q.push(s);
    while(!Q.empty()) {
        int u=Q.front();Q.pop();
        for(int i=head[u];i!=-1;i=edge[i].next) {
            int v=edge[i].v,f=edge[i].w;
            if(!dep[v] && f) {
                dep[v]=dep[u]+1;
                Q.push(v);
            }
        }
    }
    return dep[t]!=0;
}
int dfs(int u,int flow) {
    if(u==t) return flow;
    for(int i=head[u];i!=-1;i=edge[i].next) {
        int v=edge[i].v;
        int f=edge[i].w;
        if(dep[v]==dep[u]+1 && f) {
            int df=dfs(v,min(f,flow));
            if(df) {
                edge[i].w-=df;
                edge[i^1].w+=df;
                return df;
            }
        }
    }
    return 0;
}
void dinic() {
    int ans=0;
    while(bfs()) {
        int delta;
        do {
            delta=dfs(s,0x3f3f3f3f);
            ans+=delta;
        }while(delta);
    }
    cout<<ans;
    return ;
}
int main() {
    memset(head,-1,sizeof(head));
    cin>>n>>m>>s>>t;
    for(int i=0,u,v,w;i<m;i++) {
        cin>>u>>v>>w;
        add(u,v,w);
        add(v,u,0); //添加反向路，反向路的初始剩余流量是0
    }
    dinic();
    return 0;
}