洛谷P3376 (最大流模板)
1 #include<bits/stdc++.h>
2 #define int long long
3 using namespace std;
4 const int maxn=5005;
5 const int inf=0x3f3f3f3f;
6 int tot,n,m,s,t,x,y,z,d[maxn];
7 int adj[maxn],nex[maxn<<1],to[maxn<<1],w[maxn<<1],flow[maxn<<1];
8 void add(int x,int y,int z){
9 nex[tot]=adj[x];
10 adj[x]=tot;
11 to[tot]=y;
12 w[tot]=z;
13 flow[tot++]=0;
14 }
15 bool bfs(){
16 memset(d,0,sizeof(d));
17 queue<int> q;
18 d[s]=1;
19 q.push(s);
20 while(!q.empty()){
21 int u=q.front();
22 q.pop();
23 for(int i=adj[u];~i;i=nex[i]){
24 int v=to[i];
25 if(!d[v]&&w[i]>flow[i]){
26 d[v]=d[u]+1;
27 q.push(v);
28 if(v==t) return 1;
29 }
30 }
31 }
32 return 0;
33 }
34 int dfs(int u,int fw){
35 if(u==t) return fw;
36 int rest=fw;
37 for(int i=adj[u];~i&&rest;i=nex[i]){
38 int v=to[i];
39 if(d[v]==d[u]+1&&w[i]>flow[i]){
40 int k=dfs(v,min(rest,w[i]-flow[i]));
41 if(!k) d[v]=0;
42 flow[i]+=k;
43 flow[i^1]-=k;
44 rest-=k;
45 }
46 }
47 return fw-rest;
48 }
49 int dinic(){
50 int maxflow=0;
51 while(bfs()){
52 //for(int i=1;i<=n;i++)
53 maxflow+=dfs(s,inf);
54 }
55 return maxflow;
56 }
57 signed main()
58 {
59 cin>>n>>m>>s>>t;
60 memset(adj,-1,sizeof(adj));
61 for(int i=1;i<=m;i++){
62 cin>>x>>y>>z;
63 add(x,y,z);add(y,x,0);
64 }
65 cout<<dinic();
66 return 0;
67 }
注意几个细节:adj初始化为-1,tot下标从0开始;~i
加上当前弧优化
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int maxn=5005;
const int inf=0x3f3f3f3f;
int tot,n,m,s,t,x,y,z,d[maxn],now[maxn];
int adj[maxn],nex[maxn<<1],to[maxn<<1],w[maxn<<1],flow[maxn<<1];
void add(int x,int y,int z){
nex[tot]=adj[x];
adj[x]=tot;
to[tot]=y;
w[tot]=z;
flow[tot++]=0;
}
bool bfs(){
memset(d,0,sizeof(d));
queue<int> q;
d[s]=1;
q.push(s);
for(int i=1;i<=n;i++)
now[i]=adj[i];
while(!q.empty()){
int u=q.front();
q.pop();
for(int i=adj[u];~i;i=nex[i]){
int v=to[i];
if(!d[v]&&w[i]>flow[i]){
d[v]=d[u]+1;
q.push(v);
if(v==t) return 1;
}
}
}
return 0;
}
int dfs(int u,int fw){
if(u==t) return fw;
int rest=fw;
for(int i=now[u];~i&&rest;i=nex[i]){
int v=to[i];
if(d[v]==d[u]+1&&w[i]>flow[i]){
int k=dfs(v,min(rest,w[i]-flow[i]));
if(!k) d[v]=0;
flow[i]+=k;
flow[i^1]-=k;
rest-=k;
}
}
return fw-rest;
}
int dinic(){
int maxflow=0;
while(bfs()){
//for(int i=1;i<=n;i++)
maxflow+=dfs(s,inf);
}
return maxflow;
}
signed main()
{
cin>>n>>m>>s>>t;
memset(adj,-1,sizeof(adj));
for(int i=1;i<=m;i++){
cin>>x>>y>>z;
add(x,y,z);add(y,x,0);
}
cout<<dinic();
return 0;
}
洛谷P3376 (最大流模板)的更多相关文章
- 洛谷 P3376 【【模板】网络最大流】
题目描述 如题,给出一个网络图,以及其源点和汇点,求出其网络最大流. 输入 第一行包含四个正整数N.M.S.T,分别表示点的个数.有向边的个数.源点序号.汇点序号. 接下来M行每行包含三个正整数ui. ...
- 题解 洛谷 P3376 【【模板】网络最大流】
本人很弱,只会Dinic.EK与Ford-Fulkerson...(正在学习ISAP...) 这里讲Dinic... Dinic:与Ford-Fulkerson和的思路相似(话说好像最大流求解都差不多 ...
- [洛谷P3376题解]网络流(最大流)的实现算法讲解与代码
[洛谷P3376题解]网络流(最大流)的实现算法讲解与代码 更坏的阅读体验 定义 对于给定的一个网络,有向图中每个的边权表示可以通过的最大流量.假设出发点S水流无限大,求水流到终点T后的最大流量. 起 ...
- 【洛谷P3369】【模板】普通平衡树题解
[洛谷P3369][模板]普通平衡树题解 题目链接 题意: 您需要写一种数据结构(可参考题目标题),来维护一些数,其中需要提供以下操作:1. 插入x数2. 删除x数(若有多个相同的数,因只删除一个)3 ...
- 洛谷P3376【模板】网络最大流 ISAP
这篇博客写得非常好呀. 传送门 于是我是DCOI这一届第一个网络流写ISAP的人了,之后不用再被YKK她们嘲笑我用Dinic了!就是这样! 感觉ISAP是会比Dinic快,只分一次层,然后不能增广了再 ...
- bzoj 5120 无限之环 & 洛谷 P4003 —— 费用流(多路增广SPFA)
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=5120 https://www.luogu.org/problemnew/show/P4003 ...
- 【洛谷4721】【模板】分治FFT(CDQ分治_NTT)
题目: 洛谷 4721 分析: 我觉得这个 "分治 FFT " 不能算一种特殊的 FFT ,只是 CDQ 分治里套了个用 FFT (或 NTT)计算的过程,二者是并列关系而不是偏正 ...
- 草地排水 洛谷P2740 最大流 入门题目
草地排水 洛谷P2740 最大流入门题目 题意 在农夫约翰的农场上,每逢下雨,贝茜最喜欢的三叶草地就积聚了一潭水.这意味着草地被水淹没了,并且小草要继续生长还要花相当长一段时间.因此,农夫约翰修建了一 ...
- 【洛谷3865】 【模板】ST表(猫树)
传送门 洛谷 Solution 实测跑的比ST表快!!! 这个东西也是\(O(1)\)的,不会可以看我上一篇Blog 代码实现 代码戳这里
- 【最大流ISAP】洛谷P3376模板题
题目描述 如题,给出一个网络图,以及其源点和汇点,求出其网络最大流. 输入输出格式 输入格式: 第一行包含四个正整数N.M.S.T,分别表示点的个数.有向边的个数.源点序号.汇点序号. 接下来M行每行 ...
随机推荐
- 2022-7-9 第五小组 潘堂智 html学习笔记
什么是 HTML? HTML 是用来描述网页的一种语言. HTML 指的是超文本标记语言 (Hyper Text Markup Language) HTML 不是一种编程语言,而是一种标记语言 (ma ...
- docker for windows无法共享硬盘
最近在一次win10安装docker的过程中无法共享D盘.每次一点击共享就出现如下 经过网上搜索后,有反馈是防火墙的问题,结果关掉防火墙还是无法改善,查找日志 因为无法判断是什么server服务,故又 ...
- 海纳百川无所不容,Win10环境下使用Docker容器式部署前后端分离项目Django+Vue.js
原文转载自「刘悦的技术博客」https://v3u.cn/a_id_179 随着现代化产品研发的不断推进,我们会发现,几乎每个产品线都会包含功能各异的服务,而且服务与服务之间存在也会存在着错综复杂的依 ...
- 手搓一个“七夕限定”,用3D Engine 5分钟实现烟花绽放效果
七夕来咯!又到了给重要的人送惊喜的时刻. 今年,除了将心意融入花和礼物,作为程序员,用自己的代码本事手搓一个技术感十足"七夕限定"惊喜,我觉得,这是不亚于车马慢时代手写信的古典主义 ...
- java的stream让我灵光一现
说实话,我是一个到了退役也没有搞明白C++的istream和ostream的. 刚开始的时候我把<iostream>直接拆解成ios和tream 真,果粉暴露 退役之后划水,倒是从java ...
- js--自定义对象
1.直接量对象(JSON) {"name":"zhangsan","age":25} {}代表一个对象,包含多组键值对. 通常key是字符串 ...
- jsp一句话木马总结
一.无回显的命令执行(命令执行后不会在前端页面返回数据) <%Runtime.getRuntime().exec(request.getParameter("i"));%&g ...
- CF662C Binary Table (快速沃尔什变换FWT)
题面 题解 我们会发现,如果单独的一列或一行,它的答案是O1确定的,如果确定了每一行是否变换,那么最后的答案也就简单了许多, 如果确定了行的变换状压下来是x(即x的i位表示第i行是否变换,理解就行), ...
- 动态规划——leetcode5、最长回文子串
1.题目描述: 2.解题方法:动态规划 动态规划解题步骤: 1.确定状态 最后一步:如果s[i,...,j]是回文子串,那么需要满足两个条件 ① s[i] == s[j]: ② s[i+1,...,j ...
- 网卡限速工具之WonderShaper
GreatSQL社区原创内容未经授权不得随意使用,转载请联系小编并注明来源. GreatSQL是MySQL的国产分支版本,使用上与MySQL一致. 什么是WonderShaper 如何安装Wonder ...