『题解』洛谷P3376 【模板】网络最大流
Problem Portal
Portal1:Luogu
Description
如题,给出一个网络图,以及其源点和汇点,求出其网络最大流。
Input
第一行包含四个正整数\(N,M,S,T\),分别表示点的个数,有向边的个数,源点序号,汇点序号。
接下来\(M\)行每行包含三个正整数\(u_i,v_i,w_i\),表示第\(i\)条有向边从\(w_i\)出发,到达\(v_i\),边权为\(w_i\)(即该边最大流量为\(w_i\))。
Output
一行,包含一个正整数,即为该网络的最大流。
Sample Input
4 5 4 3
4 2 30
4 3 20
2 3 20
2 1 30
1 3 40
Sample Output
50
Hint
数据规模:
对于\(30\%\)的数据:\(N \leq 10,M \leq 25\);
对于\(70\%\)的数据:\(N \leq 200,M \leq 1000\);
对于\(100\%\)的数据:\(N \leq 10000,M \leq 100000\)。
样例说明:
题目中存在\(3\)条路径:
\(4 \to 2 \to 3\),该路线可通过\(20\)的流量
\(4 \to 3\),可通过\(20\)的流量
\(4 \to 2 \to 1 \to 3\),可通过\(10\)的流量(边\(4 \to 2\)之前已经耗费了\(20\)的流量)
故流量总计\(20+20+10=50\),输出\(50\)。
Solution
模板题,求最大流。
Code
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
const int INF=0x3f3f3f3f, MAXN=10005, MAXM=200005;
int n, m, s, t, u, v, val, cnt, head[MAXN], dis[MAXN], pre[MAXN];
queue<int> Q;
struct node {
int u, v, val, flow, nxt;
} edge[MAXM];
inline void addedge(int u, int v, int w) {//前向星存图
edge[++cnt].u=u; edge[cnt].v=v; edge[cnt].val=w; edge[cnt].nxt=head[u]; head[u]=cnt;
}
inline bool Emonds_Karp() {
memset(pre, 0, sizeof(pre));
memset(dis, 0, sizeof(dis));//初始化
dis[s]=INF;
pre[s]=0;
Q.push(s);
while(!Q.empty()) {//类似于SPFA
int x=Q.front();
Q.pop();
for (int i=head[x]; i; i=edge[i].nxt) {
int v=edge[i].v;
if (!dis[v] && edge[i].val>edge[i].flow) {
Q.push(v);
dis[v]=min(dis[x], edge[i].val-edge[i].flow);
pre[v]=i;
}
if (dis[t]) break;
}
}
return dis[t];
}
int main() {
scanf("%d%d%d%d",&n, &m, &s, &t);
cnt=1;
for(int i=1; i<=m; i++) {
scanf("%d%d%d",&u, &v, &val);
addedge(u, v, val);
addedge(v, u, 0);//前向星存图
}
int ans=0;
while (Emonds_Karp()) {//Emonds Karp算法
for (int i=t; i!=s; i=edge[pre[i]].u) {
edge[pre[i]].flow+=dis[t];
edge[pre[i] xor 1].val-=dis[t];
}
ans+=dis[t];
}
printf("%d\n",ans);
return 0;
}
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
const int INF=0x3f3f3f3f, MAXN=10005, MAXM=200005;
struct node {
int nxt, to, val;
} edge[MAXM];
queue<int> Q;
int n, m, s, t, u, v, cnt, val, dis[MAXN], head[MAXN];
inline void addedge(int u, int v, int w) {//前向星存图
edge[++cnt].val=w; edge[cnt].to=v; edge[cnt].nxt=head[u]; head[u]=cnt;
}
inline bool bfs() {
memset(dis, -1, sizeof(dis));
dis[s]=0;
Q.push(s);
while (!Q.empty()) {
int x=Q.front();
Q.pop();
for (int i=head[x]; ~i; i=edge[i].nxt) {
int v=edge[i].to;
if (edge[i].val && dis[v]==-1) {
dis[v]=dis[x]+1;
Q.push(v);
}
}
}
if (~dis[t]) return 1;
return 0;
}
inline int dfs(int u, int flow) {
if (u==t) return flow;
int ret=flow;
for (int i=head[u]; ~i; i=edge[i].nxt) {
if (ret<=0) break;
int v=edge[i].to;
if (edge[i].val && dis[v]==dis[u]+1) {
int x=dfs(v, min(edge[i].val, ret));
ret-=x;
edge[i].val-=x;
edge[i xor 1].val+=x;
}
}
return flow-ret;
}
inline int Dinic() {//Dinic算法
int ret=0;
while (bfs()) ret+=dfs(s, INF);
return ret;
}
int main() {
scanf("%d%d%d%d",&n, &m, &s, &t);
memset(head, -1, sizeof(head));
cnt=1;
for (int i=1; i<=m; i++) {
scanf("%d%d%d",&u, &v, &val);
addedge(u, v, val);
addedge(v, u, 0);//双边
}
printf("%d\n",Dinic());
return 0;
}
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
const int INF=0x3f3f3f3f, MAXN=10005, MAXM=200005;
queue<int> Q;
int n, m, s, t, u, v, val, cnt, ans, head[MAXN], head1[MAXN], deep[MAXN], pre[MAXN], a[MAXN];
struct node {
int v, w, nxt;
} edge[MAXM];
inline void addedge(int u, int v, int w) {//前向星存图
edge[cnt].v=v; edge[cnt].w=w; edge[cnt].nxt=head[u]; head[u]=cnt++;
}
inline void bfs(int t) {
for (int i=1; i<=n; i++)
head1[i]=head[i];
for (int i=1; i<=n; i++)
deep[i]=n;
deep[t]=0;
Q.push(t);
while (!Q.empty()) {
int u=Q.front();
Q.pop();
for (int i=head[u]; ~i; i=edge[i].nxt)
if (deep[edge[i].v]==n && edge[i^1].w) {
deep[edge[i].v]=deep[u]+1;
Q.push(edge[i].v);
}
}
}
int calc(int s,int t) {//计算
int ans=INF, u=t;
while (u!=s) {
ans=min (ans,edge[pre[u]].w);
u=edge[pre[u] xor 1].v;
}
u=t;
while (u!=s) {
edge[pre[u]].w-=ans;
edge[pre[u] xor 1].w+=ans;
u=edge[pre[u] xor 1].v;
}
return ans;
}
inline void ISAP(int s, int t) {//ISAP算法
int u=s;
bfs(t);
for (int i=1; i<=n; i++)
a[deep[i]]++;
while (deep[s]<n) {
if (u==t) {
ans+=calc(s, t);
u=s;
}
bool flag=0;
for (int &i=head1[u]; ~i; i=edge[i].nxt)
if (deep[u]==deep[edge[i].v]+1 && edge[i].w) {
flag=1;
u=edge[i].v;
pre[edge[i].v]=i;
break;
}
if (!flag) {
int Min=n-1;
for (int i=head[u]; ~i; i=edge[i].nxt)
if (edge[i].w) Min=min(Min, deep[edge[i].v]);
if ((--a[deep[u]])==0) break;
a[deep[u]=Min+1]++;
head1[u]=head[u];
if (u!=s) u=edge[pre[u] xor 1].v;
}
}
}
int main() {
memset(head, -1, sizeof (head));
scanf("%d%d%d%d",&n, &m, &s, &t);
for (int i=1; i<=m; i++) {
scanf("%d%d%d",&u, &v, &val);
addedge(u, v, val);
addedge(v, u, 0);//双边
}
ISAP(s, t);
printf ("%d\n",ans);
return 0;
}
『题解』洛谷P3376 【模板】网络最大流的更多相关文章
- 『题解』洛谷P1063 能量项链
原文地址 Problem Portal Portal1:Luogu Portal2:LibreOJ Portal3:Vijos Description 在\(Mars\)星球上,每个\(Mars\)人 ...
- 『题解』洛谷P3384 【模板】树链剖分
Problem Portal Portal1: Luogu Description 如题,已知一棵包含\(N\)个结点的树(连通且无环),每个节点上包含一个数值,需要支持以下操作: 操作\(1\): ...
- 『题解』洛谷P1993 小K的农场
更好的阅读体验 Portal Portal1: Luogu Description 小\(K\)在\(\mathrm MC\)里面建立很多很多的农场,总共\(n\)个,以至于他自己都忘记了每个农场中种 ...
- 『题解』洛谷P2296 寻找道路
更好的阅读体验 Portal Portal1: Luogu Portal2: LibreOJ Description 在有向图\(\mathrm G\)中,每条边的长度均为\(1\),现给定起点和终点 ...
- 『题解』洛谷P1351 联合权值
更好的阅读体验 Portal Portal1: Luogu Portal2: LibreOJ Description 无向连通图\(\mathrm G\)有\(n\)个点,\(n - 1\)条边.点从 ...
- 『题解』洛谷P2170 选学霸
更好的阅读体验 Portal Portal1: Luogu Description 老师想从\(N\)名学生中选\(M\)人当学霸,但有\(K\)对人实力相当,如果实力相当的人中,一部分被选上,另一部 ...
- 『题解』洛谷P1083 借教室
更好的阅读体验 Portal Portal1: Luogu Portal2: LibreOJ Portal3: Vijos Description 在大学期间,经常需要租借教室.大到院系举办活动,小到 ...
- 『题解』洛谷P1314 聪明的质监员
更好的阅读体验 Portal Portal1: Luogu Portal2: LibreOJ Portal3: Vijos Description 小T是一名质量监督员,最近负责检验一批矿产的质量.这 ...
- 『题解』洛谷P2357 守墓人
Portal Portal1: Luogu Description 在一个荒凉的墓地上有一个令人尊敬的守墓人,他看守的墓地从来没有被盗过, 所以人们很放心的把自己的先人的墓安顿在他那守墓人能看好这片墓 ...
随机推荐
- 《Windows内核分析》专题-索引目录
该篇博客整理了<Windows内核分析>专题的各篇博文,方便查找. 一.保护模式 二.进程与线程 [Windows内核分析]KPCR结构体介绍 (CPU控制区 Processor Cont ...
- 关于人工智能和python
人工智能的话题在近几年可谓是相当火热,前几天看快本时其中有一个环节就是关于人工智能的,智能家电.智能机器人.智能工具等等,在我的印象里,提到人工智能就会出现 Python,然后我便在网上查找了相关信息 ...
- 本人亲测-SSM环境搭建(使用eclipse作为示例,过程挺全的,可作为参考)
本人亲测-SSM环境搭建(使用eclipse作为示例,过程挺全的,可作为参考) 本人亲测-SSM环境搭建(使用eclipse作为示例,过程挺全的,可作为参考) 本人亲测-SSM环境搭建(使用eclip ...
- C#刷遍Leetcode面试题系列连载(4) No.633 - 平方数之和
上篇文章中一道数学问题 - 自除数,今天我们接着分析 LeetCode 中的另一道数学题吧~ 今天要给大家分析的面试题是 LeetCode 上第 633 号问题, Leetcode 633 - 平方数 ...
- 想转行做程序员,目前想学WEB前端,想问该自学还是报培训班
首先我们还是先看一下WEB前端目前的工资情况吧,我在IT招聘网站拉勾网来进行搜索1-3年WEB前端工作经验大专学历的条件来看. 深圳: 可以看出目前深圳的平均的工资都在10K以上,因为大城市给的机会多 ...
- shark恒破解笔记3-EAX决定胜负
PEID查壳 od载入 输入假的注册码 查找出错字符串 往上查找是否有关键跳转和关键call 可以看到此处有个je跳转 实现了跳转,并且跳过了我们注册成功的地址 网上查找这个跳转的关键call,这个c ...
- [Luogu3070][USACO13JAN]岛游记Island Travels
题目描述 Farmer John has taken the cows to a vacation out on the ocean! The cows are living on N (1 < ...
- 实用---生命游戏 Java
本程序由四个类组成: 其中Init_data,用于初始化各个活细胞的状态judge_state,用于判断下一代的细胞状态,并进行更新.set_color,用于给GUI界面中各个细胞涂色set_fram ...
- 【Python秘籍】ASCII码与字符的转换
如何在python中显示ASCII码呢?其实你只需要记住两个函数即可:ord()和 chr(),这两个函数都是python内置的函数,不需要引入任何的包,直接就可以使用. 一.显示ASCII码 显示A ...
- Python编程系列---Python中装饰器的几种形式及万能装饰器
根据函数是否传参 是否有返回值 ,可以分析出装饰器的四种形式: 形式一:无参无返回值 def outer(func): def wrapper(): print("装饰器功能1" ...