【模板】EK求最大流

求解最大流问题的最简便方法
但是效率不是很理想

参考博客:网络流-最大流

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=505,maxe=40005;
int n,e,que[maxn],fa[maxn],id[maxn],ans;
bool vis[maxn];
struct data{
    int son[maxe],nxt[maxe],lnk[maxn],tot;
    int flw[maxe],cap[maxe];
    void add(int x,int y,int w){
        son[++tot]=y;nxt[tot]=lnk[x];lnk[x]=tot;
        flw[tot]=0;cap[tot]=w;
    }
}a;
inline int red(){
    int tot=0;char ch=getchar();
    while (ch<'0'||'9'<ch) ch=getchar();
    while ('0'<=ch&&ch<='9') tot=tot*10+ch-48,ch=getchar();
    return tot;
}
bool bfs(){
    memset(vis,0,sizeof(vis));
    int hed=0,til=1;
    que[1]=1;vis[1]=1;
    while (hed!=til)
     for (int j=a.lnk[que[++hed]];j;j=a.nxt[j])
      if (!vis[a.son[j]]&&a.cap[j]>a.flw[j]){
        que[++til]=a.son[j];vis[a.son[j]]=1;
        fa[a.son[j]]=que[hed];id[a.son[j]]=j;
        if (a.son[j]==n) return 1;
      }
    return 0;
}
int main(){
    n=red(),e=red();a.tot=1;
    for (int i=1,x,y,w;i<=e;i++)
     x=red(),y=red(),w=red(),a.add(x,y,w),a.add(y,x,0);
    while (bfs()){
        int Min=0x3f3f3f3f;
        for (int x=n;x!=1;x=fa[x])
         Min=min(Min,a.cap[id[x]]-a.flw[id[x]]);
        ans+=Min;
        for (int x=n;x!=1;x=fa[x])
         a.flw[id[x]]+=Min,a.flw[id[x]^1]-=Min;
    }
    printf("%d",ans);
    return 0;
}
以下是使用Edmonds-Karp算法实现最大流的示例C代码: ```c #include <stdio.h> #include <stdlib.h> #include <string.h> #include <limits.h> #define MAXN 1000 // 最大节点数 #define MAXM 1000000 // 最大边数 int n, m; // n为节点数,m为边数 int s, t; // 源节点和汇节点 int head[MAXN], nxt[MAXM], to[MAXM], cap[MAXM], cur[MAXN], cnt = 1; // 邻接表表示边 void add_edge(int u, int v, int c) { nxt[++cnt] = head[u]; head[u] = cnt; to[cnt] = v; cap[cnt] = c; } void init() { memset(head, 0, sizeof(head)); cnt = 1; } int bfs() { static int q[MAXN], dep[MAXN]; int l = 0, r = 0; memset(dep, -1, sizeof(dep)); q[r++] = s; dep[s] = 0; while (l < r) { int u = q[l++]; for (int i = head[u]; i; i = nxt[i]) { int v = to[i]; if (dep[v] == -1 && cap[i]) { dep[v] = dep[u] + 1; q[r++] = v; } } } return dep[t] != -1; } int dfs(int u, int flow) { if (u == t) return flow; for (int &i = cur[u]; i; i = nxt[i]) { int v = to[i]; if (cap[i] && dep[v] == dep[u] + 1) { int f = dfs(v, fmin(flow, cap[i])); if (f) { cap[i] -= f; cap[i ^ 1] += f; return f; } } } return 0; } int dinic() { int max_flow = 0; while (bfs()) { memcpy(cur, head, sizeof(cur)); int f; while (f = dfs(s, INT_MAX)) { max_flow += f; } } return max_flow; } int main() { scanf("%d %d %d %d", &n, &m, &s, &t); init(); for (int i = 1; i <= m; i++) { int u, v, c; scanf("%d %d %d", &u, &v, &c); add_edge(u, v, c); add_edge(v, u, 0); // 反向边初始容量为0 } printf("%d\n", dinic()); return 0; } ``` 代码中使用了邻接表存储图的边,使用了bfs和dfs实现了Edmonds-Karp算法。其中bfs用于搜索增广路,dfs用于在增广路上更新
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值