P3376 【模板】网络最大流当前弧优化

最新推荐文章于 2021-09-13 09:16:44 发布

箱@子

最新推荐文章于 2021-09-13 09:16:44 发布

阅读量129

点赞数

CC 4.0 BY-SA版权

分类专栏：网络流文章标签：网络流

本文链接：https://blog.csdn.net/qq_33914105/article/details/100976311

网络流专栏收录该内容

8 篇文章

订阅专栏

本文深入探讨了Dinic算法的原理与应用，这是一种高效的网络流算法，用于解决最大流问题。通过详细的代码解析，展示了如何在C++中实现Dinic算法，并提供了完整的示例代码。适用于对网络流算法感兴趣的读者。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<bits/stdc++.h>?
#include<iostream>
#include<cstring>
#include<stack>
#include<stdio.h>
const int maxn = 1e4 + 10;
int head[maxn];
int dft[maxn];
int cur[maxn];
int tot;
int n, m, S, T;
using namespace std;
struct node{
    int to;
    int next;
    int w;
    node() {}
    node(int a, int b, int c) : to(a), next(b), w(c) {}
}edge[maxn * 2];
 
void edgeadd(int a, int b, int c){
    edge[tot] = node(b, head[a], c);
    head[a] = tot++;
    edge[tot] = node(a, head[b], 0);
    head[b] = tot++;
}
 
int bfs(int s, int e){
    memset(dft, 0, sizeof(dft));
    queue<int> q;
    q.push(s);
    dft[s] = 1;
    while(!q.empty()){
        int nw = q.front();
        q.pop();
        for(int i = head[nw]; i != -1; i = edge[i].next){
            int w = edge[i].w;
            int v = edge[i].to;
            if(w <= 0 || dft[v]) continue;
            dft[v] = dft[nw] + 1;
            q.push(v);
        }
    }
    return dft[e]; 
} 
 
int dfs(int s, int e, int f){
    if(s == e) return f;
    int ans = 0;
    for(int i = cur[s]; i != -1; i = edge[i].next){
        int v = edge[i].to;
        int w = edge[i].w;
        cur[s] = i;
        if(dft[v] != dft[s] + 1 || w <= 0)
            continue;
        int we = dfs(v, e, min(w, f - ans));//由于多次增广， 不能为f
        edge[i].w -= we;
        edge[i ^ 1].w += we;
        ans += we;
    }
    return ans;
}
void init(){
    memset(head, -1, sizeof(head));
    tot = 0;
}
 
int dinic(int s, int e){
    int cnt = 0;
    while(bfs(s, e)){
        for(int i = 1; i <= n; i++){
            cur[i] = head[i];
        } 
        cnt += dfs(s, e, 0x3f3f3f3f);
    }
    return cnt;
}
 
int main(){
    ios::sync_with_stdio(false);
    cin >> n >> m >> S >> T;
    init();
    for(int i = 1; i <= m; i++){
        int a, b, c;
        cin >> a >> b >> c;
        edgeadd(a, b , c);
    }
    cout << dinic(S, T) << endl;
    return 0;
}