Acwing 2171-EK求最大流

//EK算法
#include<bits/stdc++.h>

using namespace std;

const int N=1010,M=20010,INF=1e8;

int n,m ,S,T;
int h[N],e[M],f[M],ne[M],idx;
int d[N],pre[N];
bool st[N];

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(st,false,sizeof st);
    q.push(S),st[S]=true,d[S]=INF;

    while(!q.empty())
    {
        int t=q.front();
        for(int i=h[t];~i;i=ne[i])
        {
            int ver=e[i];
            if(!st[ver]&&f[i])
            {
                st[ver]=true;
                d[ver]=min(d[t],f[i]);
                pre[ver]=i;//前驱边
                if(ver==T) return true;
                q.push(ver);
            }
        }
        q.pop();
    }
    return false;
}

int EK()
{
    int r=0;
    while(bfs())
    {
        r+=d[T];
        for(int i=T;i!=S;i=e[pre[i]^1])
        {
            f[pre[i]]-=d[T],f[pre[i]^1]+=d[T];
        }
    }
    return r;
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie();cout.tie();

    cin>>n>>m>>S>>T;
    memset(h,-1,sizeof h);
    while(m--)
    {
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
    }

    cout<<EK()<<endl;
    return 0;
}
/*
 最大流最小割定理
(1)流f是最大流
(2)残留网络中不存在增广路径
(3)存在某一个割(S,T),使得 |f|=c(S,T)
 */

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 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用于在增广路上更新流。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值