图论模板——最大流问题

最新推荐文章于 2024-07-23 22:06:37 发布
最新推荐文章于 2024-07-23 22:06:37 发布
阅读量288 收藏 1
点赞数
分类专栏： 图论板子 文章标签： 图论 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_51771856/article/details/125371624
版权

Dinic算法的时间复杂度为O（n^2m）。实际运用远远小于这个上界。
特别的，Dinic算法求解二分图最大匹配的时间复杂度为O（msqrt（n））
最大流问题模板

#include <bits/stdc++.h>
using namespace std;
const int N = 10010,M = 200010, INF = 1e8;
int n,m,S,T;
int h[N],ne[M],e[M],f[M],idx;
int d[N],q[N],cur[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()
{
    int hh = 0, tt = 0;
    memset(d, -1, sizeof d);
    q[0] = S, d[S] = 0, cur[S] = h[S];
    while (hh <= tt)
    {
        int t = q[hh ++ ];
        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[ ++ tt] = ver;
            }
        }
    }
    return false;
}
int find(int u,int limit)
{
    if (u==T) return limit;
    int flow = 0;
    for (int i = cur[u];~i&&flow<limit;i = ne[i])
    {
        cur[u] = i;//当前弧优化
        int ver = e[i];
        if (d[ver] == d[u]+1&&f[i])
        {
            int t = find(ver,min(f[i],limit-flow));
            if (!t) d[ver] = -1;
            f[i] -= t,f[i^1] += t,flow+=t;
        }
    }
    return flow;
}
int dinic()
{
    int r = 0,flow;
    while (bfs()) while (flow = find(S,INF)) r += flow;
    return r;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&S,&T);
    memset(h,-1,sizeof h);
    while (m--)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
    }
    printf("%d\n",dinic());
    return 0;
}

应用一：二分图匹配问题
飞行员配对方案问题 https://www.acwing.com/problem/content/2177/

int main()
{
    scanf("%d%d",&m,&n);
    S = 0,T = n+1;
    memset(h,-1,sizeof h);
    for (int i = 1;i<=m;i++) add(S,i,1);
    for (int i = m+1;i<=n;i++) add(i,T,1);
    int a,b;
    while (scanf("%d%d",&a,&b),a!=-1) add(a,b,1);
    printf("%d\n",dinic());
    for (int i=0;i<idx;i+=2)
        if (e[i]>m&&e[i]<=n&&!f[i]) 
            printf("%d %d\n",e[i^1],e[i]);
    return 0;
}

应用二：无源汇的最大流问题
无源汇上下界可行流 https://www.acwing.com/problem/content/2190/)
注意流量守恒

int main()
{
    memset(h,-1,sizeof h);
    scanf("%d%d",&n,&m);
    S = 0,T = n+1,tot = 0;
    for (int i=0;i<m;i++)
    {
        int a,b,c,d;
        scanf("%d%d%d%d",&a,&b,&c,&d);
        add(a,b,c,d);
        A[a] -= c,A[b] += c;//出去c，进来c
    }
    for (int i=1;i<=n;i++)
    {
        if (A[i]>0) add(S,i,0,A[i]),tot += A[i];//如果进来大于0
        else if (A[i]<0) add(i,T,0,-A[i]);//如果出去
    }
    if (dinic()!=tot) printf("NO\n");//如果最低的流量都不行
    else 
    {
        printf("YES\n");
        for (int i=0;i<m*2;i+=2)
            printf("%d\n",f[i^1]+l[i]);//输出方案
    }
    return 0;
}

应用三：有源汇的最大流问题
有源汇上下界最大流 https://www.acwing.com/problem/content/2191/
先链接一个t到s的边，看是否能满足所有边的最小下限
如果满足再做一次dinic()

int main()
{
    memset(h,-1,sizeof h);
    int s,t,tot = 0;
    scanf("%d%d%d%d",&n,&m,&s,&t);
    S = 0,T = n+1;
    while (m--)
    {
        int a,b,c,d;
        scanf("%d%d%d%d",&a,&b,&c,&d);
        add(a,b,d-c);
        A[a] -= c,A[b] += c;
    }
    for (int i=1;i<=n;i++)
        if (A[i]>0) add(S,i,A[i]),tot += A[i];
        else if (A[i]<0) add(i,T,-A[i]);
    add(t,s,INF);
    if (dinic()<tot) puts("No Solution");
    else 
    {
        int res = f[idx-1];
        f[idx-1] = f[idx-2] = 0;
        S = s,T = t;
        printf("%d",res+dinic());
    }
}

应用四：有源汇的最小流问题
有源汇上下界最小流 https://www.acwing.com/problem/content/2192/
加t到s的边，然后dinic，然后再反向dinic

int main()
{
    memset(h,-1,sizeof h);
    int s,t,tot = 0;
    scanf("%d%d%d%d",&n,&m,&s,&t);
    S = 0,T = n+1;
    while (m--)
    {
        int a,b,c,d;
        scanf("%d%d%d%d",&a,&b,&c,&d);
        add(a,b,d-c);
        A[a] -= c,A[b] += c;
    }
    for (int i=1;i<=n;i++)
        if (A[i]>0) add(S,i,A[i]),tot += A[i];
        else if (A[i]<0) add(i,T,-A[i]);
    add(t,s,INF);
    if (dinic()<tot) puts("No Solution");
    else 
    {
        int res = f[idx-1];
        f[idx-1] = f[idx-2] = 0;
        S = t,T = s;
        printf("%d",res-dinic());
    }
}

应用五：多源汇的最大流问题
多源汇最大流 https://www.acwing.com/problem/content/2236/

int main()
{
    memset(h,-1,sizeof h);
    int ss,tt;
    scanf("%d%d%d%d",&n,&m,&ss,&tt);
    S = 0,T = n+1;
    while (ss--)
    {
        int x;scanf("%d",&x);
        add(S,x,INF);
    }
    while (tt--)
    {
        int x;scanf("%d",&x);
        add(x,T,INF);
    }
    while (m--)
    {
        int a,b,c;scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
    }
    printf("%d\n",dinic());
}

应用六：关键边
伊基的故事 I - 道路重建 https://www.acwing.com/problem/content/2238/
dinic后从S和T遍历

void dfs(int u, bool st[], int t)
{
    st[u] = true;
    for (int i = h[u]; ~i; i = ne[i])
    {
        int j = i ^ t, ver = e[i];
        if (f[j] && !st[ver])
            dfs(ver, st, t);
    }
}
int main()
{
    memset(h,-1,sizeof h);
    scanf("%d%d",&n,&m);
    for (int i = 0;i<m;i++)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
    }
    S = 0,T = n-1;
    dinic();
    dfs(S, vis_s, 0);
    dfs(T, vis_t, 1);

    int res = 0;
    for (int i = 0; i < m * 2; i += 2)
        if (!f[i] && vis_s[e[i ^ 1]] && vis_t[e[i]])
            res ++ ;

    printf("%d\n", res);
    return 0;
}

应用七：有边大小使用限制的最大流判定
秘密挤奶机 https://www.acwing.com/problem/content/2279/

bool check(int x)
{
    for (int i=0;i<idx;i++) 
        if (w[i]>x) f[i] = 0;
        else f[i] = 1;
    return dinic()>=k;
}
int main()
{
    memset(h,-1,sizeof h);
    scanf("%d%d%d",&n,&m,&k);
    S = 1,T = n;
    for (int i = 0;i<m;i++)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
    }
    int l = 1,r = 1e6;
    while (l<r)
    {
        int mid = l+r>>1;
        if (check(mid)) r = mid;
        else l = mid+1;
    }
    printf("%d\n",r);
    return 0;
}

应用八：拆点
餐饮 https://www.acwing.com/problem/content/2242/

int main()
{
    memset(h,-1,sizeof h);
    scanf("%d%d%d",&n,&m,&k);
    S = 0,T = 2*n+m+k+1;
    for (int i = 2*n+1;i<=2*n+m;i++) add(S,i,1);
    for (int i = 1;i<=n;i++) add(i,i+n,1);
    for (int i = 2*n+m+1;i<=2*n+m+k;i++) add(i,T,1);
    for (int i = 1;i<=n;i++)
    {
        int a,b,x;
        scanf("%d%d",&a,&b);
        for (int j = 0;j<a;j++)
        {
            scanf("%d",&x);
            add(2*n+x,i,1);
        }
        for (int j = 0;j<b;j++)
        {
            scanf("%d",&x);
            add(n+i,2*n+m+x,1);
        }
    }
    printf("%d\n",dinic());
    return 0;
}
卐邾翊
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
常用代码模板3——搜索与图论.html
06-13
常用代码模板3——搜索与图论.html
ACM模板——清华大学
05-30
清华大学ACM模板，汇集了很多算法的标准模式，而且有详细的算法讲解，图论、数论、网络流等等，从浅到深，彻底解剖，是ACM进阶者不二的选择~！
图论模板——最大流及费用流模板
weixin_30851409的博客
08-09 67
图论模板——最大流及费用流模板 最大流——SAP 时间复杂度：O(v^2*e) const int MAXN=1010;//点数的最大值 const int MAXM=1010;//边数的最大值 const int INF=0x3f3f3f3f; struct Node { int from,to,next; int cap; }edge[MAXM]; int tol; int...
图论模板——费用流（无法处理负环）
weixin_51771856的博客
06-23 209
费用流
图论——网络流问题
victor_hq的博客
05-05 647
网络流最大流算法模板题EK算法dinic算法二分图匹配飞行员的配对问题圆桌问题上下界可行流无源汇上下界可行流有源汇上下界最大流例题——东方文花帖有源汇上下界最小流多源汇最大流多源汇最大流（最大流，超级源汇点）关键边伊基的故事 I - 道路重建最大流判定秘密挤奶机（最大流，最大流判定，二分）星际转移问题（最大流，最大流判定，分层图）拆点（不相交路径，只能用k次）餐饮（最大流，拆点）最长递增子序列问题（最大流，拆点，DP）企鹅游行（最大流+拆点+枚举汇点）建图实战猪（最大流，分层图，困难题）最小割模板Dinic
【图论】网络流问题——最大流入门（Dinic算法）
nefu_ljw的博客
04-04 1200
参考文章： 1. 博客园：Dinic算法（研究总结，网络流） 2. 洛谷博客：网络最大流-从入门开始，详细讲到实用易懂的Dinic算法 前置知识：DFS、BFS、建图（链式前向星） 本文主要是用 Dinic算法 解决最大流问题。 洛谷 P3376 【模板】网络最大流 最大流问题，Dinic算法，最坏时间复杂度为O(V2*E)，然而一般情况下是达不到这么大的，要不这题早超时了。 算法流程 1、根据残量网络计算层次图。 2、在层次图中使用DFS进行增广直到不存在增广路。 3、重复以上步骤直到无法增广。
图论--网络流最大流问题
日拱一卒,功不唐捐
10-28 7903
问题表述：给定一幅图（n个结点，m条边），每一条边有一个容量，现在需要将一些物品从结点s（称为源点）运送到结点t（称为汇点），可以从其他结点中转，求最大的运送量。 在介绍最大流问题的解决方法之前，先介绍几个概念. 网络：网络是一个有向带权图，包含一个源点和一个汇点，没有反向平行边。 网络流：网络流即网上的流，是定义在网络边集E上的一个非负函数flow={flow（u，v）}, flow（u,...
图论进阶——网络流
y13438985462的博客
09-16 235
网络流基本概念 基本概念 1.1 流网络，不考虑反向边 1.2 可行流，不考虑反向边 1.2.1 两个条件：容量限制、流量守恒 1.2.2 可行流的流量指从源点流出的流量 −-− 流入源点的流量 1.2.3 最大流是指最大可行流 1.3 残留网络，考虑反向边，残留网络的可行流f′f'f′ +++ 原图的可行流fff = 原题的另一个可行流 (1) ∣f′+f∣=∣f′∣+∣f∣|f' + f| = |f'| + |f|∣f′+f∣=∣f′∣+∣f∣ (2) ∣f′∣|f'|∣f′∣ 可能是负数 1.4 增
图论 —— 最大团问题
热门推荐
Alex_McAvoy的博客
01-23 2万+
【问题描述】 当 G′是图 G的子图，且 G′是关于 V′ 的完全图时，子图 G' 为图 G 的团；当 G' 是团，且不是其他团的子集时，G' 为图 G 的极大团；当 G' 是极大团时，且点数最多，G' 为图 G 最大团 当 G′ 中所有点不相邻，最大点集最大的图 G′ 为图 G 的最大独立集，且最大独立集数=补图的最大团 当用个数最少的团覆盖图 G 所有的点时，称为最小团覆盖...
代码模板——图论
shuijiaozbc的博客
08-04 147
图论
图论建模——分层图详解
m0_74315382的博客
06-01 1118
分层图，顾名思义。是将原图按不同状态分为若干与原图连接方式相同的图层，图层之间以特定方式连接的一类建图方式。如果画成立体图，大概是这样的：假设原图位于0{0}0层，总共有kkk层图。那么对于任意i∈1ki\in[1,k]i∈1k，层iii内的节点之间的连接方式与0{0}0层是完全相同的（与原图连接方式相同）；但是层与层之间的连接方式不一定相同，具体取决于题意假设不考虑节点所在层，仅考虑它在原图上对应的编号。假设有一条无向边SV(S,V)SV。
最大流_最大流模板_
09-30
最大流问题是一个经典的图论问题，它在计算机科学和网络流理论中占有重要地位，尤其在运筹学、算法设计和优化等领域应用广泛。在信息学奥林匹克竞赛中，选手们经常需要解决这类问题来求解网络中的最大传输能力。本文...
建议收藏算法基础课模板大全
08-12
搜索与图论 —— 代码模板链接 常用代码模板3——搜索与图论 DFS与BFS 树与图的遍历：拓扑排序 最短路 最小生成树 二分图：染色法、匈牙利算法 数学知识 —— 代码模板链接 常用代码模板4——数学知识 质数 约数 ...
ACM.rar_图论
09-23
6. **网络流**：如Ford-Fulkerson方法和Edmonds-Karp算法，用于在图中找到最大流量，常用于资源分配问题。 7. **二分图**：图中任意两条边不相交的顶点集，二分图匹配问题可通过匈牙利算法解决，广泛应用于匹配问题...
力扣刷题-图论-岛屿类问题-集合实现（c++实现）
07-23 622
【代码】力扣刷题-图论-岛屿类问题-集合实现（c++实现）
算法-KMP字符串匹配
最新发布
xx_xb的博客
07-23 230
next[i] 就是使子串 s[0…i] 有最长相等前后缀的前缀的最后一位的下标。总体来说解next数组和模板串匹配的过程很相似，触类旁通。
算法题目整合5
qq_62571233的博客
07-23 879
思路：利用left和right记录连续相同的字符串索引，计算长度len，如果是0,1,2,4这种情况直接输出-1并return。如果是奇数直接输出，如果是偶数的话，可以拆成3和len-3两个部分。的最小值是101100，求其第4554小即在从后数的四个问号位置分别加上‘3’，‘5’，‘5’，‘4’即可，即551153；表示当前卡牌正面（或反面）加上上一个取余结果再对3取余的值，举个例子，比如说。），由于不同的取余方案搭配上正反面的组合，共有6中结果，所以在。，表示当前卡牌正面为1，当遍历的。
【Hot100】LeetCode—152. 乘积最大子数组
weixin_44382896的博客
07-23 277
【代码】【Hot100】LeetCode—152. 乘积最大子数组。
吉林大学ACM竞赛代码库：图论与网络流算法概览
有上下界限制的最小（最大）流问题，以及Dinic和HLPP算法分别用于求解最大流，它们的时间复杂度分别为O(V^2*E)和O(V^3)；最小费用流问题考虑了边上的成本，有多种高效实现，包括两种不同的O(V*E*F)复杂度的算法；...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值