[洛谷]最大流 模板

最新推荐文章于 2022-03-19 16:31:58 发布
最新推荐文章于 2022-03-19 16:31:58 发布
阅读量338 收藏
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ureaster/article/details/78428085
版权

 P3376 【模板】网络最大流


开始我打的是假的网络流,深搜宽搜的临时变量存的都是边

明显变比点要多得多


这是EK

#include<stdio.h>
#include<string.h>
#define INF 1000000009
#define N 10010
int n,m,s,t,x,y,z,pre[N*20],way[N*20];
int head[N],next[N*20],to[N*20],flow[N*20],size;
int Q[N*20],front,rear;
int min(int a,int b) {
    return a<b ? a : b;
}
void add(int a,int b,int c) {
    size++;
    to[size]=b; flow[size]=c;
    next[size]=head[a];
    head[a]=size;
}
int bfs() {
    memset(pre,0,sizeof pre);
    front=0; rear=0;
    Q[rear++]=1;
    while(front<rear) {
        int d=Q[front++];
        for(int i=head[to[d]];i;i=next[i]) {
            if(!pre[i]&&flow[i]) {
                pre[i]=d;
                way[i]=min(flow[i],way[d]);
                if(to[i]==t) return i;
                Q[rear++]=i;
            }
        }
    }
    return 0;
}
int maxflow() {
    int ans=0;
    for(int delt,e;e=bfs();ans+=delt) {
        for(delt=way[e];e!=1;e=pre[e]) {
            flow[e]-=delt; flow[e^1]+=delt;
        }
    }
    return ans;
}
int main() {
    scanf("%d%d%d%d",&n,&m,&s,&t);
    add(0,s,INF); way[1]=INF;
    for(int i=1;i<=m;i++) {
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z); add(y,x,0);
    }
    printf("%d\n",maxflow());
}


..假的dinic 

#include<stdio.h>
#include<string.h>
#define INF 1000000009
#define N 10010
int n,m,s,t,x,y,z,level[N*20];
int head[N],next[N*20],to[N*20],flow[N*20],size;
int Q[N*20],front,rear;
int min(int a,int b) {
    return a<b ? a : b;
}
void add(int a,int b,int c) {
    size++;
    to[size]=b; flow[size]=c;
    next[size]=head[a];
    head[a]=size;
}
int bfs() {
    memset(level,0,sizeof level);
    front=0; rear=0;
    level[1]=1;
    Q[rear++]=1;
    int link=0;
    while(front<rear) {
        int d=Q[front++];
        for(int i=head[to[d]];i;i=next[i]) {
            if(!level[i]&&flow[i]) {
                level[i]=level[d]+1;
                Q[rear++]=i;
                if(to[i]==t) link=1;
            }
        }
    }
    return link;
}
int dfs(int e,int minf) {
    if(to[e]==t||!minf) return minf;
    int f=0,delt;
    for(int i=head[to[e]];i;i=next[i]) {
        if(level[i]==level[e]+1&&(delt=dfs(i,min(minf,flow[i])))) {
            flow[i]-=delt; flow[i^1]+=delt;
            f+=delt; minf-=delt;
            if(!minf) break;
        }
    }
    return f;
}
int maxflow() {
    int ans=0;
    while(bfs()) {
        ans+=dfs(1,INF);
    }
    return ans;
}
int main() {
    scanf("%d%d%d%d",&n,&m,&s,&t);
    add(0,s,INF);
    for(int i=1;i<=m;i++) {
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z); add(y,x,0);
    }
    printf("%d\n",maxflow());
}


..真的dinic 

#include<stdio.h>
#include<string.h>
#define INF 1000000009
#define N 10010
int n,m,s,t,x,y,z,level[N];
int head[N],next[N*20],to[N*20],flow[N*20],size=1;
int Q[N*20],front,rear;
int min(int a,int b) {
    return a<b ? a : b;
}
void add(int a,int b,int c) {
    size++;
    to[size]=b; flow[size]=c;
    next[size]=head[a];
    head[a]=size;
}
int bfs() {
    memset(level,0,sizeof level);
    front=0; rear=0;
    level[s]=1;
    Q[rear++]=s;
    while(front<rear) {
        int d=Q[front++];
        for(int i=head[d];i;i=next[i]) {
            if(!level[to[i]]&&flow[i]) {
                level[to[i]]=level[d]+1;
                if(to[i]==t) return 1;
                Q[rear++]=to[i];
            }
        }
    }
    return 0;
}
int dfs(int u,int minf) {
    if(u==t) return minf;
    int f=0,delt;
    for(int i=head[u];i;i=next[i]) {
        if(level[to[i]]==level[u]+1&&(delt=dfs(to[i],min(minf,flow[i])))) {
            flow[i]-=delt; flow[i^1]+=delt;
            f+=delt; minf-=delt;
            if(!minf) break;
        }
    }
    return f;
}
int maxflow() {
    int ans=0;
    while(bfs()) {
        ans+=dfs(s,INF);
    }
    return ans;
}
int main() {
    scanf("%d%d%d%d",&n,&m,&s,&t);
    for(int i=1;i<=m;i++) {
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z); add(y,x,0);
    }
    printf("%d\n",maxflow());
}


然后用dinic A了一些水题


 P2756 飞行员配对方案问题

#include<stdio.h>
#include<string.h>
#define N 105
#define INF 1000000007
int m,n,x,y,level[N];
int head[N],next[N*N*2],to[N*N*2],flow[N*N*2],size=1;
int Q[N*N*2],front,rear;
int min(int a,int b) {
    return a<b ? a : b;
}
void add(int a,int b,int c) {
    size++;
    to[size]=b; flow[size]=c;
    next[size]=head[a];
    head[a]=size;
}
int bfs() {
    memset(level,0,sizeof level);
    front=0; rear=0;
    level[0]=1;
    Q[rear++]=0;
    while(front<rear) {
        int d=Q[front++];
        for(int i=head[d];i;i=next[i]) {
            if(!level[to[i]]&&flow[i]) {
                level[to[i]]=level[d]+1;
                if(to[i]==n+1) return 1;
                Q[rear++]=to[i];
            }
        }
    }
    return 0;
}
int dfs(int u,int maxf) {
    if(u==n+1||!maxf) return maxf;
    int f=0,delt;
    for(int i=head[u];i;i=next[i]) {
        if(level[to[i]]==level[u]+1&&(delt=dfs(to[i],min(maxf,flow[i])))) {
            flow[i]-=delt; flow[i^1]+=delt;
            f+=delt; maxf-=delt;
            if(!maxf) break;
        }
    }
    return f;
}
int maxflow() {
    int ans=0;
    while(bfs()) {
        ans+=dfs(0,INF);
    }
    return ans;
}
int main() {
    scanf("%d%d",&m,&n);
    for(int i=1;i<=m;i++) {
        add(0,i,1); add(i,0,0);
    }
    for(int i=m+1;i<=n;i++) {
        add(i,n+1,1); add(n+1,i,0);
    }
    while(scanf("%d%d",&x,&y)&&x!=-1) {
        add(x,y,1); add(y,x,0);
    }
    printf("%d\n",maxflow());
    for(int i=1;i<=m;i++) {
        for(int j=head[i];j;j=next[j]) {
            if(!flow[j]&&to[j]>m&&to[j]<=n) {
                printf("%d %d\n",i,to[j]);
                break;
            }
        }
    }
}

 P1343 地震逃生

#include<stdio.h>
#include<string.h>
#define N 205
#define INF 1000000007
int n,m,k,x,y,z,level[N];
int head[N],next[N*20],to[N*20],flow[N*20],size=1;
int Q[N*20],front,rear;
int min(int a,int b) {
    return a<b ? a : b;
}
void add(int a,int b,int c) {
    size++;
    to[size]=b; flow[size]=c;
    next[size]=head[a];
    head[a]=size;
}
int bfs() {
    memset(level,0,sizeof level);
    front=0; rear=0;
    level[1]=1;
    Q[rear++]=1;
    while(front<rear) {
        int d=Q[front++];
        for(int i=head[d];i;i=next[i]) {
            if(!level[to[i]]&&flow[i]) {
                level[to[i]]=level[d]+1;
                if(to[i]==n) return 1;
                Q[rear++]=to[i];
            }
        }
    }
    return 0;
}
int dfs(int u,int minf) {
    if(u==n||!minf) return minf;
    int f=0,delt;
    for(int i=head[u];i;i=next[i]) {
        if(level[to[i]]==level[u]+1&&(delt=dfs(to[i],min(minf,flow[i])))) {
            flow[i]-=delt; flow[i^1]+=delt;
            f+=delt; minf-=delt;
            if(!minf) break;
        }
    }
    return f;
}
int maxflow() {
    int ans=0;
    while(bfs()) {
        ans+=dfs(1,INF);
    }
    return ans;
}
int main() {
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=m;i++) {
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z); add(y,x,0);
    }
    int once=maxflow();
    if(!once) puts("Orz Ni Jinan Saint Cow!");
    else printf("%d %d\n",once,(k-1)/once+1);
}

 P3386 【模板】二分图匹配

嗯这用的是匈牙利算法.. 
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 1005
using namespace std;
int n,m,e,x,y;
int g[N*2][N*2],vis[N*2],link[N*2];
int dfs(int u) {
    for(int i=n+1;i<=n+m;i++) {
        if(!vis[i]&&g[u][i]) {
            vis[i]=1;
            if(!link[i]||dfs(link[i])) {
                link[i]=u; return 1;
            }
        }
    }
    return 0;
}
int hungary() {
    int ret=0;
    for(int i=1;i<=n;i++) {
        memset(vis,0,sizeof vis);
        ret+=dfs(i);
    }
    return ret;
}
void input() {
    scanf("%d%d%d",&n,&m,&e);
    for(int i=1;i<=e;i++) {
        scanf("%d%d",&x,&y);
        if(x<=n&&y<=m) {
            g[x][n+y]=1;
        }
    }
}
int main() {
    input();
    printf("%d\n",hungary());
}




ureaster
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
洛谷 [模板]KMP字符串匹配
01-03
又水了一篇 题目链接：https://www.luogu.com.cn/problem/P3375 #include using namespace std; int Next[1000010]={0}; void getNEXT(string str) { Next[0] = -1; int i=-1, j=0;... Next[j] =
最大流模板(Maximum Flow)
τ's blog
06-02 977
最大流模板(Maximum Flow) 一篇写得通俗易懂介绍最大流的文章:最大流模板【EdmondsKarp算法，简称EK算法，O(m^2n)】 残留网络 残存容量 c~f~(u, v) = c(u, v) - f(u, v); 算法导论中对残存网络的一些解释 一条边所能允许的反向流量最多将其正向流量抵消. 残存网络中的这些反向边允许算法将发送出来的流量发送回去. 而将流量从...
网络流最大流入门（洛谷P3376）
forezxl的博客
09-28 908
所谓网络或容量网络指的是一个连通的赋权有向图 D= （V、E、C） ， 其中V 是该图的顶点集，E是有向边(即弧)集，C是弧上的容量。此外顶点集中包括一个起点和一个终点。
[洛谷P3376题解]网络流（最大流）的实现算法讲解与代码
qq_16204149的博客
07-20 591
[洛谷P3376题解]网络流（最大流）的实现算法讲解与代码 定义 对于给定的一个网络，有向图中每个的边权表示可以通过的最大流量。假设出发点S水流无限大，求水流到终点T后的最大流量。 起点我们一般称为源点，终点一般称为汇点 内容前置 1.增广路 ​ 在一个网络从源点S到汇点T的一条各边剩余流量都大于0(还能让水流通过,没有堵住)的一条路。 2.分层 ​ 预处理出源点到每个点的距离(每次寻找增广路都要,因为以前原本能走的路可能因为水灌满了,导致不能走了).作用是保证只往更远的地方放水,避免兜圈子或者是没事就走回
最大流模板
You and me
08-21 294
最近学习了最大流虽然学的一般都是会简单的使用，一下是我用过的几个模板 模板一sap： #include &lt;bits/stdc++.h&gt; /* 用maze来存图 */ typedef long long ll; using namespace std; const int MAXN=1100; int maze[MAXN][MAXN]; int gap[MAXN],dis[MAXN...
高斯列主消元
左岸的博客
04-20 372
#include<stdio.h> #define N 10010 double A[N][N]; double X[N]; int n,m; int GSLX(int m,int n)//列主消元函数 { int p; int i,j,k; for(i=1;i<n;i++) { double Max=A[i][i]; p=i;
2020洛谷初赛模拟.pdf
10-10
2020洛谷初赛模拟, 简介 针对 2019 年 CSP 非专业组第一轮（相当于原 NOIP 初赛）的题目改革，我们命制了一套第一轮的模拟试卷。 本次模拟测试形式和 CSP 第一轮 完 全 一 致，难度介于 J 组（基础组）和 S 组...
洛谷B3623枚举排列.cpp
最新发布
12-16
洛谷真题源代码【附注释】
最小费用最大流（洛谷-P3381）.rar
09-16
最小费用最大流（洛谷-P3381）.rar
洛谷P1563玩具谜题
08-26
洛谷P1563玩具谜题
预流推进算法的最大流模板
08-13
预流推进算法是不同与增广路思想的另一种求最大流的算法，主要思想是把流从源点经过边，并通过调整顶点的标号来把流压入汇点。 该算法需要维护一个装有活跃顶点的标号，所谓活跃定点则是流出量小于流入量的顶点，则具有超额量，欲使该顶点平衡（流出==流入），则需要把该定点积累的流压出，被压入流的顶点只能是距离标号比该点标号小1且这条边未满流的顶点，若不存在这样的顶点但是改点仍未活跃顶点，则需要把该点的距离标号加1. 这只是利用普通的FIFO队列，为了加速算法，可以维护一个优先队列，以顶点的距离标号为优先级，距离标号大的先出列。
洛谷-4722 【模板】最大流 加强版 / 预流推进
mkopvec的博客
08-12 600
题目描述 给定 n 个点，m 条有向边，给定每条边的容量，求从点 s 到点 t 的最大流。 输入格式 第一行包含四个正整数n、m、s、t，用空格分隔，分别表示点的个数、有向边的个数、源点序号、汇点序号。 接下来m行每行包含三个正整数ui、vi、ciu_i、v_i、c_iui​、vi​、ci​​，用空格分隔，表示第iii条有向边从uiu_iui​​出发，到达viv_ivi​​，容量为cic_ici...
网络流算法//最大流（洛谷3376）//EK算法//dinic算法//最小（大）费用最大流（待续）
我歌月徘徊，我舞影零乱
08-08 438
目录 最大流流算法 （EK算法） 时间复杂度 O(V*（E^2）) Dinic算法 时间复杂度O(（V ^ 2） * E) 在洛谷题解中看到了一句很有启发的话：网络流善于解决各种有要求的匹配 联想到题目是匹配问题，且满足网络流要求的数据范围，可以尝试网络流 V是点的数目，E是边的数目 最大流流算法 （EK算法） 时间复杂度 O(V*（E^2）) 定义： 我们有...
【模板】网络最大流
weixin_34250709的博客
03-04 168
https://www.luogu.org/problem/show?pid=3376 题目描述 如题，给出一个网络图，以及其源点和汇点，求出其网络最大流。 输入输出格式 输入格式： 第一行包含四个正整数N、M、S、T，分别表示点的个数、有向边的个数、源点序号、汇点序号。 接下来M行每行包含三个正整数ui、vi、wi，表示第i条有向边从ui出发，到达vi，边权为wi（即该边最大流...
【模板】最大流详细证明(Dinic+当前弧优化)(附:题- How Many to Be Happy?)
weixin_45819197的博客
04-06 1343
初识网络流-最大流==最小割 什么是网络流 网络流(network-flows)是一种类比水流的解决问题方法，与线性规划密切相关。网络流的理论和应用在不断发展。而我们今天要讲的就是网络流里的一种常见问题——最大流问题。 如下图： 最大流 (Dinic) 进而引出了最大流问题：从源点出发到达汇点，所能流入最大的流量。 增广路 那么我们如何知道S(源点)---->T(汇点)的最大流量呢？ 我们的思路是： 1.先找出一条增广路径，得到这条路径权值的最小值mins，之后该路径边上的流（权值）的值减去min
最大流模板(Edmond-Karp)
he_qingjun的博客
07-28 92
题目：https://www.luogu.com.cn/problem/P3376 学习心得： FF算法，EK算法敲了五遍才真正理解，理解了才发现原来挺简单的。。估计刚开始被网络流这个名头吓住了，，下意识以为很难，，两个星期学完了30多个算法，有一半敲了十遍以上，另一遍估计就敲了五遍左右，，心急以为早点学完打比赛就不会慌了，，然而等到我回头用心看的时候才是真的理解了，，虽然学了好多个算法，但是都是短时间学的，不能灵活的运用，不能根据题目建模，如果用到这些模板算法的时候，建模是最难的一部分，，不会分离出模型，
【洛谷】P3376【模板】网络最大流（Ford-Fulkerson算法）
悲伤无声地涌出来
10-09 587
题目链接 题目描述 如题，给出一个网络图，以及其源点和汇点，求出其网络最大流。 输入输出格式 输入格式： 第一行包含四个正整数N、M、S、T，分别表示点的个数、有向边的个数、源点序号、汇点序号。 接下来M行每行包含三个正整数ui、vi、wi，表示第i条有向边从ui出发，到达vi，边权为wi（即该边最大流量为wi） 输出格式： 一行，包含一个正整数，即为该网络的最大...
最大流模板（dinic）
IDrandom的专栏
08-09 279
//dinic #include&amp;lt;bits/stdc++.h&amp;gt; using namespace std; typedef long long LL; const int INF = 0x3f3f3f3f; const int maxm =120020; const int maxn =1210; struct Edge{ int u,v,w; int nxt; }e[m...
最大网络流的多种解法（洛谷P3376 网络最大流 为例）
NaHCOOO_的博客
03-19 2007
P3376 【模板】网络最大流 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)https://www.luogu.com.cn/problem/P3376 本题是一个最大网络流的模版题，我也是刚刚学完最大网络流，所以就借此机会，练练手。 关于最大网络流的定义这里就不多做介绍（口齿不清）。如果不知道的可以去先看一下大佬的文章 这里给出一个链接(8条消息) 网络流之最大流算法（...
洛谷P1681 最大正方形II怎么写
07-08
洛谷 P1681 最大正方形II 是一个动态规划问题，要求给定一个由 '0' 和 '1' 组成的矩阵，找出其中最大的正方形，并输出其边长。 以下是一个 C++ 编写的解答示例： ```cpp #include #include #include using ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值