Dinic算法——最大流

最新推荐文章于 2024-07-24 22:13:06 发布

green qwq

最新推荐文章于 2024-07-24 22:13:06 发布

阅读量286

点赞数 4

分类专栏：洛谷算法题合集文章标签：算法深度优先

本文链接：https://blog.csdn.net/m0_60738889/article/details/138233831

版权

洛谷算法题合集专栏收录该内容

142 篇文章 0 订阅

订阅专栏

网络最大流

题目描述

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

输入格式

第一行包含四个正整数 $n, m, s, t$ ，分别表示点的个数、有向边的个数、源点序号、汇点序号。

接下来 $m$ 行每行包含三个正整数 $u_i,v_i,w_i$ ，表示第 $i$ 条有向边从 $u_i$ 出发，到达 $v_i$ ，边权为 $w_i$ （即该边最大流量为 $w_i$ ）。

输出格式

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

样例 #1

样例输入 #1

样例输出 #1

提示

样例输入输出 1 解释

题目中存在 $3$ 条路径：

$4\to 2\to 3$ ，该路线可通过 $20$ 的流量。
$4\to 3$ ，可通过 $20$ 的流量。
$4\to 2\to 1\to 3$ ，可通过 $10$ 的流量（边 $4\to 2$ 之前已经耗费了 $20$ 的流量）。

故流量总计 $20 + 20 + 10 = 50$ 。输出 $50$ 。

数据规模与约定

对于 $30\%$ 的数据，保证 $n\leq10$ ， $m\leq25$ 。
对于 $100\%$ 的数据，保证 $\leq n\leq200$ ， $\leq m\leq 5000$ ， $\leq w\lt 2^{31}$ 。

模板

//dinic算法

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
// #define int long long
using namespace std;

const int N = 10210,M = 120010*2,INT = 1e9;

int e[M],ne[M],h[N],f[M],idx;
bool st[N];
int d[N],cur[N];//当前弧优化
//这个是拿来干什么的呢？用于find函数，就是碰到后面的流量都大于limits的时候，就不搜索了。
int n,m,S,T;

void add(int a,int b,int c){
    e[idx]=b,f[idx]=c,ne[idx]=h[a],h[a]=idx++;
    e[idx]=a,f[idx]=0,ne[idx]=h[b],h[b]=idx++;
}

bool bfs(){
    queue<int>q;
    memset(d,-1,sizeof d);//这个-1可以用来st数组
    q.push(S);
    d[S]=0;
    cur[S]=h[S];
    // 将源点加入队列，源点深度为0，初始化源点当前弧为表头
    
    while(q.size()){
        int t=q.front();
        q.pop();
        
        for(int i=h[t];~i;i=ne[i]){
            int ver=e[i];
            if(d[ver]==-1&&f[i]){
                // 存在增广路才能从这里分层
                d[ver]=d[t]+1;// 更新深度
                cur[ver]=h[ver];// 所有点当前弧最开始都是表头
                if(ver==T)return true;
                q.push(ver);
            }
        }
    }
    return false;
}

int find(int u,int lim){// DFS，u是当前点，lim是u之前的路径能流过的最大值
    if(u==T)return lim;// 如果已经流到汇点了就返回这个最大值
     // 当前往下流的流量
     
    int flow=0;
     
    for(int i=cur[u];~i&&flow<lim;i=ne[i]){
         // 注意应从当前弧开始遍历，并且每次要更新。没有剩余流量也应该退出
        cur[u]=i;//这样子就是为了你前面肯定不能走了，才往这里走，因此把cur更新成当前值
        int ver=e[i];
         
        if(d[ver]==d[u]+1&&f[i]){
            //只有在某一层的下边才能搜索（为了防止出现环）
             
            int t=find(ver,min(f[i],lim-flow));
            // 如果没有流量那么说明后面没有增广路了，这个点不会用到，将深度设为-1
            if(!t)d[ver]=-1;
            // 当前边减流量，反向边加流量，实际流量加流量
            f[i]-=t,f[i^1]+=t,flow+=t;
        }
    }
    return flow;
}

int dinic(){
    int r=0,flow=0;
    while(bfs())while((flow=find(S,INT)))r+=flow;
    // 只要存在增广路就一直DFS，直到DFS出的流量为0
    
    return r;
}

signed main(){
    //cin>>n>>m>>S>>T;
    
    scanf("%d%d%d%d",&n,&m,&S,&T);
    
    memset(h,-1,sizeof h);
    
    for(int i=1;i<=m;i++){
        int a,b,c;
        // cin>>a>>b>>c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
    }
    
    printf("%d",dinic());
    
    return 0;
}