题目地址:

https://www.acwing.com/problem/content/2190/

给定一个包含nnn个点mmm条边的有向图,每条边都有一个流量下界和流量上界。求一种可行方案使得在所有点满足流量平衡条件的前提下,所有边满足流量限制。

输入格式:
第一行包含两个整数nnn和mmm。接下来mmm行,每行包含四个整数a,b,c,da,b,c,da,b,c,d表示点aaa和bbb之间存在一条有向边,该边的流量下界为ccc,流量上界为ddd。点编号从111到nnn。

输出格式:
如果存在可行方案,则第一行输出YES,接下来mmm行,每行输出一个整数,其中第iii行的整数表示输入的第iii条边的流量。如果不存在可行方案,直接输出一行NO。如果可行方案不唯一,则输出任意一种方案即可。

数据范围:
1≤n≤2001≤n≤2001≤n≤200
1≤m≤102001≤m≤102001≤m≤10200
1≤a,b≤n1≤a,b≤n1≤a,b≤n
0≤c≤d≤100000≤c≤d≤100000≤c≤d≤10000

无源汇上下界可行流又称为循环流,像是一定量的水在水管里循环流动。由于每条边有容量下限,所以容易想到构建新的网络,将每条边的容量设为是原图的容量上下限之差,这样在新网络里任一一个可行流必然满足容量限制,但是流量守恒就不一定满足了。我们可以采用这样的方法,如果入边的总容量下限ccc大于出边的总容量下限c′c'c′,那么可以从源点SSS补充差值给这个点(即从源点到这个点开一条容量是c−c′c-c'c−c′的边);反之,则从该点把多出来的差补给汇点TTT(即从该点到汇点开一条容量是c′−cc'-cc′−c的边)。题目问是否存在满足容量限制的可行流,相当于问新网络里的最大流是否可以达到源点的出边容量之和(这里有个细节,为什么我们只需要看新网络里最大流是源点出边的满流就行了呢,不需要看其是否是汇点入边的满流吗?答案是,确实不用看。原因在于,设每个点uuu的入边容量下界之和减去出边容量下界之和等于A[u]A[u]A[u],考虑∑u≠S,TA[u]\sum_{u\ne S,T} A[u]∑u​=S,T​A[u],原图每条边u→vu\to vu→v的容量下界clc_lcl​,都会使得A[u]A[u]A[u]减少clc_lcl​,也会使得A[v]A[v]A[v]增加clc_lcl​,所以∑u≠S,TA[u]=0\sum_{u\ne S,T} A[u]=0∑u​=S,T​A[u]=0,而SSS出边满流就等于∑u≠S,T∧A[u]>0A[u]\sum_{u\ne S,T\land A[u]>0}A[u]∑u​=S,T∧A[u]>0​A[u],其等于TTT入边满流−∑u≠S,T∧A[u]<0A[u]-\sum_{u\ne S,T\land A[u]<0}A[u]−∑u​=S,T∧A[u]<0​A[u])。可以粗略的看一下对应关系,对于新网络里的任意一个满流,去掉源点和汇点(和它们连接的边)之后,将每个边里的流加一下其在原图里的下限,就得到的一个原图的可行流;对于原图的一个可行流,将每个边里的流减去其在原图里的下限,再加上源点和汇点按上述所说方式构图,就能得到新网络里的一个满流。所以问题就转化为在新网络里求最大流,可以用Dinic算法。代码如下:

#include <iostream>
#include <cstring>
using namespace std;const int N = 210, M = (10200 + N) * 2, INF = 1e8;
int n, m, S, T;
int h[N], e[M], ne[M], f[M], l[M], idx;
int q[N], d[N], cur[N], A[N];// 构造新网络的残留网络。每条边的容量是上限 - 下限
void add(int a, int b, int c, int d) {e[idx] = b, ne[idx] = h[a], f[idx] = d - c, l[idx] = c, h[a] = idx++;e[idx] = a, ne[idx] = h[b], f[idx] = 0, h[b] = idx++;
}bool bfs() {memset(d, -1, sizeof d);int hh = 0, tt = 0;q[tt++] = 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 v = e[i];if (d[v] == -1 && f[i]) {d[v] = d[t] + 1;if (v == T) return true;cur[v] = h[v];q[tt++] = v;}}}return false;
}int dfs(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 v = e[i];if (d[v] == d[u] + 1 && f[i]) {int t = dfs(v, min(limit - flow, f[i]));if (!t) d[v] = -1;f[i] -= t, f[i ^ 1] += t, flow += t;}}return flow;
}int dinic() {int r = 0, flow;while (bfs()) while(flow = dfs(S, INF)) r += flow;return r;
}int main() {cin >> n >> m;S = 0, T = n + 1;memset(h, -1, sizeof h);for (int i = 0; i < m; i++) {int a, b, c, d;cin >> a >> b >> c >> d;add(a, b, c, d);A[a] -= c, A[b] += c;}int tot = 0;// 求一下满流的流量for (int i = 1; i <= n; i++)if (A[i] > 0) add(S, i, 0, A[i]), tot += A[i];else if (A[i] < 0) add(i, T, 0, -A[i]);// 如果最大流不是源点出边全满的流,则说明无解,否则输出流量if (dinic() < tot) cout << "NO" << endl;else {cout << "YES" << endl;for (int i = 0; i < m * 2; i += 2) // 流量是残留网络的反向边的容量,再加下限cout << f[i ^ 1] + l[i] << endl;}return 0;
}

时间复杂度O(n2m)O(n^2m)O(n2m),空间O(n)O(n)O(n)。

【ACWing】2188. 无源汇上下界可行流相关推荐

  1. bzoj 2406 矩阵——有源汇上下界可行流

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2406 二分答案.把 b 的 n 个行作为一排, m 个列作为一排,每行和每列之间连上下界为 ...

  2. 一篇网络流 基本模型超全总结(最大流 费用流 多源汇最大流 上下界可行流) 思路+代码模板

    文章目录 一.网络流与最大流 二.网络流三个基本性质 三.重要定义定理 四.最大流算法 <Ⅰ> Edmonds-Karp算法(EK算法) 1.EK算法 2.算法思想: 3.代码模板 4.模 ...

  3. BZOJ 2406 LuoguP4194 矩阵 有上下界可行流

    分析: 这道题乍一看--卧槽这都什么玩意-- 然后发现给了个A矩阵,要求一个可行的B矩阵,使得矩阵C=A-B的每一行的和的绝对值和每一列的和的绝对值的最大值最小-- 好拗口啊-- 什么最大值最小之类的 ...

  4. 网络流__4 上下界可行流

    网络流 __4 上下界可行流 对于容量有上下界规范的网络流问题 无源汇上下界可行流 n个点m条边的有向图,每条边有一个流量下界和流量上界规范,求是否存在一个可行流 设原网络为(G,F),变换后网络为( ...

  5. Shoot the Bullet(ZOJ3229)(有源汇上下界最大流)

    描述 ensokyo is a world which exists quietly beside ours, separated by a mystical border. It is a utop ...

  6. loj#115. 无源汇有上下界可行流

    \(\color{#0066ff}{ 题目描述 }\) 这是一道模板题. \(n\) 个点,\(m\) 条边,每条边 \(e\) 有一个流量下界 \(\text{lower}(e)\) 和流量上界 \ ...

  7. 无源汇有上下界可行流(网络流进阶)

    无源汇有上下界可行流(也就是循环流) 模型:一个网络,求出一个流,使得每条边的流量必须>=Li且<=Hi, 每个点必须满足总流入量=总流出量(流量守恒)(这个流的特点是循环往复,无始无终) ...

  8. 有上下界网络流 ---- Zoj3229 Shoot the Bullet|东方文花帖|【模板】有源汇上下界最大流

    学习资料+题目链接 题目大意: 模板讲解 #include <bits/stdc++.h> #define mid ((l + r) >> 1) #define Lson rt ...

  9. HDU 4940 Destroy Transportation system(无源汇上下界网络流)

    Problem Description Tom is a commander, his task is destroying his enemy's transportation system. Le ...

最新文章

  1. laravel mysql sum查询并排行_必看!PHP常见面试题——MySQL篇(二)
  2. 美国康普SYSTIMAX iPatch智能配线系统介绍
  3. vs release 调试 路径设置
  4. QQ 5.0侧滑HorizontalScrollView以及自定义ViewGroup
  5. syslog数据接收并处理
  6. 使用SampleGrabber过滤器扑获图像
  7. 设计模式 (一) 初次体验
  8. HDU-2084(简单DP)
  9. vSAN ReadyNode™中可以(也不能)更改的内容
  10. android markdown软件,安卓上的Markdown写作软件——Markdrop
  11. 面试者为何从来得不到反馈?
  12. python安装教程-Python安装包+安装教程
  13. 俄亥俄州立大学计算机科学转学成功,俄亥俄州立大学转学要求
  14. 用 VC++ 调用最小的浏览器内核创建一个简易的浏览器
  15. FW:nbsp;司马相如琴挑文君真相:劫色劫…
  16. win10进入bios步骤
  17. matlab1到100求和for_一个简单的MATLAB程序(1到100求和)
  18. note-Linux大棚命令百篇1 Shell基础及工具篇
  19. 【建站】手把手教你域名注册(腾讯云)
  20. Android解析服务器Json数据实例

热门文章

  1. 打造一支有服务能力的销售队伍
  2. emqttd 2.2安装和测试使用
  3. linux sftp拷贝文件夹,Linux下sftp命令传输文件的例子
  4. 会议OA之我的会议(查询)
  5. java时间段分割_任意一个起止时间段(如:20160101-20161009),用java将这个时间段拆分成一个个按自然周组成的时间段...
  6. html表格图片垂直居中 css,利用Display: table;实现img图片垂直居中
  7. 基于灰色模型的端口短信预测和垃圾短信治理研究
  8. 用Bibtex导出GB/T 7714等格式引用的方法
  9. ARM GICv3中断控制器
  10. type、dtype和astype