ZOJ 3229 有上下界最大流
1: /**
2: ZOJ 3229 有上下界的最大流
3: 两次求最大流的过程,非二分
4: 有源汇上下界的最大流问题, 首先连接 sink -> src, [0,INF].
5: 根据net的正负,来建立 Supersrc 与 supersink 之间的边,做一次 maxflow.
6: 若所有的Supersrc 与 Supersink满流,则说明存在可行流.
7: 然后删除 sink -> src之间的边.(cap 置零即可). 从src -> sink 做一次最大流.
8: 两次最大流的和即为整个网络的最大流.
9: */
10:
11: #include<iostream>
12: #include<cmath>
13: #include<memory>
14: #include <string.h>
15: #include <cstdio>
16: #include <vector>
17: using namespace std;
18:
19: #define V 1500 // vertex
20: #define E V *80 // edge
21: #define INF 0x3F3F3F3F // 1061109567
22:
23: int i,j,k;
24: #define REP(i,n) for((i)=0;(i)<(int)(n);(i)++)
25: #define snuke(c,itr) for(__typeof((c).begin()) itr=(c).begin();itr!=(c).end();itr++)
26:
27: struct MaxFlow
28: {
29: struct Edge
30: {
31: int v, w, next; //w for capicity
32: int lb,up;
33: } edge[E];
34:
35: int head[V]; // head[u]表示顶点u第一条邻接边的序号, 若head[u] = -1, u没有邻接边
36: int e; // the index of the edge
37: int src, sink;
38: int net[V]; // 流入此节点的流的下界和 - 流出此节点的流的下界和,对于带上下界的来进行使用
39:
40:
41: void addedge(int u, int v, int w, int lb = 0, int up = INF, int rw = 0)
42: {
43: edge[e].v = v;
44: edge[e].w= w;
45: edge[e].next = head[u];
46: edge[e].lb = lb, edge[e].up = up;
47: head[u] = e++;
48: // reverse edge v -> u
49: edge[e].v = u;
50: edge[e].w = rw;
51: edge[e].lb = lb, edge[e].up = up;
52: edge[e].next = head[v];
53: head[v] = e++;
54: }
55:
56: int ISAP(int VertexNum )
57: {
58: int u, v, max_flow, aug, min_lev;
59: int curedge[V], parent[V], level[V];
60: int count[V], augment[V];
61:
62: memset(level, 0, sizeof(level));
63: memset(count, 0, sizeof(count));
64: REP(i,VertexNum+1) curedge[i] = head[i];
65: max_flow = 0;
66: augment[src] = INF;
67: parent[src] = -1;
68: u = src;
69:
70: while (level[src] < VertexNum)
71: {
72: if (u == sink)
73: {
74: max_flow += augment[sink];
75: aug = augment[sink];
76: for (v = parent[sink]; v != -1; v = parent[v])
77: {
78: i = curedge[v];
79: edge[i].w -= aug;
80: edge[i^1].w += aug;
81: augment[edge[i].v] -= aug;
82: if (edge[i].w == 0) u = v;
83: }
84: }
85: for (i = curedge[u]; i != -1; i = edge[i].next)
86: {
87: v = edge[i].v;
88: if (edge[i].w > 0 && level[u] == (level[v]+1))
89: {
90: augment[v] = min(augment[u], edge[i].w);
91: curedge[u] = i;
92: parent[v] = u;
93: u = v;
94: break;
95: }
96: }
97: if (i == -1)
98: {
99: if (--count[level[u]] == 0) break;
100: curedge[u] = head[u];
101: min_lev = VertexNum;
102: for (i = head[u]; i != -1; i = edge[i].next)
103: if (edge[i].w > 0)
104: min_lev = min(level[edge[i].v], min_lev);
105: level[u] = min_lev + 1;
106: count[level[u]]++;
107: if (u != src ) u = parent[u];
108: }
109: }
110: return max_flow;
111: }
112: // girl 0-m-1, day m,m+n-1, src m+n, sink m+n+1. all m+n+2 point
113: void solve()
114: {
115: int N,M; // n days m girl
116: while(scanf("%d%d", &N,&M) != EOF)
117: {
118: e = 0;
119: memset(head, -1, sizeof(head));
120: memset(net, 0, sizeof(net));
121: int G; src = M+N, sink = M+N+1;
122: for(int i=0; i<M; i++)
123: {
124: scanf("%d", &G);
125: addedge(i,sink,INF-G, G,INF);
126: net[i] -= G;
127: net[sink] += G;
128: }
129: vector<int> CE;
130: for(int i=0; i<N; i++)
131: {
132: int C,D; scanf("%d%d", &C,&D);
133: addedge(src, M+i, D, 0, D);
134: for(int j=0; j<C; j++)
135: {
136: int T,L,R; scanf("%d%d%d", &T,&L,&R);
137: CE.push_back(e);
138: addedge(M+i,T,R-L, L,R);
139: net[M+i] -= L;
140: net[T] += L;
141: }
142: }
143: int spec = e;
144: // 添加从sink -> src 容量为INF的边
145: addedge(sink, src, INF,0,INF);
146: src = M+N+2; sink = M+N+3; // M+N+4 point
147: int rangea = e;
148: for(int i=0; i<M+N+2; i++)
149: {
150: if(net[i] >= 0) addedge(src, i, net[i]);
151: else addedge(i, sink,-net[i]);
152: }
153: double ret = 0;
154: int rangeb = e;
155: // 从super src ->super sink 做一次最大流
156: ret+=ISAP(M+N+4);
157: bool flag = true;
158: // 判断是否满流
159: for(int i= rangea; i<rangeb; i+=2)
160: {
161: if(edge[i].w !=0)
162: {
163: flag = false;
164: break;
165: }
166: }
167: if(flag)
168: {
169: // 修改 src sink,然后把 从sink -> src 的边删除
170: src = M+N; sink = M+N+1;
171: edge[spec].w = 0; edge[spec+1].w = 0;
172: // 从 src->sink 做一次最大流
173: ret += ISAP(M+N+2);
174: int tmp = ret;
175: printf("%d\n", tmp);
176: for(int i=0; i<CE.size(); i++)
177: printf("%d\n", edge[CE[i]+1].lb + edge[CE[i]+1].w);
178: }else printf("-1\n");
179: cout<<endl;
180: }
181: }
182: }sap;
183:
184: int main()
185: {
186: // freopen("1.txt","r",stdin);
187: sap.solve();
188: return 0;
189: }
转载于:https://www.cnblogs.com/sosi/p/3716721.html
ZOJ 3229 有上下界最大流相关推荐
- Shoot the Bullet(ZOJ3229)(有源汇上下界最大流)
描述 ensokyo is a world which exists quietly beside ours, separated by a mystical border. It is a utop ...
- 【有源汇点上下界最小流】[SGU176]Flow construction
题目大意 给出N个点,M条有向边 如果有向边的标号是1的话,就表示该边的上界下界都为容量 如果有向边的标号为0的哈,表示该边的下界为0,上界为容量 现在问,从1到N的最小流是多少,并输出每条边的流量 ...
- 【bzoj2324】[ZJOI2011]营救皮卡丘 最短路-Floyd+有上下界费用流
原文地址:http://www.cnblogs.com/GXZlegend/p/6832504.html 题目描述 皮卡丘被火箭队用邪恶的计谋抢走了!这三个坏家伙还给小智留下了赤果果的挑衅!为了皮卡丘 ...
- ACM-ICPC 2018 沈阳赛区网络预赛 F. Fantastic Graph(有源上下界最大流 模板)
关于有源上下界最大流: https://blog.csdn.net/regina8023/article/details/45815023 #include<cstdio> #includ ...
- LOJ - #116. 有源汇有上下界最大流(有源汇有上下界的最大流)
题目链接:点击查看 题目大意:给出一个 n 个点和 m 条边的有向图,每条边都有一个流量限制 [ lower , upper ],给定源点 s 和汇点 t ,求出源点到汇点的最大流 题目分析:参考博客 ...
- 【BZOJ3876】[Ahoi2014]支线剧情 有上下界费用流
[BZOJ3876][Ahoi2014]支线剧情 Description [故事背景] 宅男JYY非常喜欢玩RPG游戏,比如仙剑,轩辕剑等等.不过JYY喜欢的并不是战斗场景,而是类似电视剧一般的充满恩 ...
- Loj#116-[模板]有源汇有上下界最大流
正题 题目链接:https://loj.ac/p/116 题目大意 nnn个点mmm条边的一张图,每条边有流量上下限制,求源点到汇点的最大流. 解题思路 先别急着求上面那个,考虑一下怎么求无源点汇点的 ...
- LOJ116 有源汇有上下界最大流(上下界网络流)
考虑有源汇上下界可行流:由汇向源连inf边,那么变成无源汇图,按上题做法跑出可行流.此时该inf边的流量即为原图中该可行流的流量.因为可以假装把加上去的那些边的流量放回原图. 此时再从原来的源向原来的 ...
- loj#115. 无源汇有上下界可行流
\(\color{#0066ff}{ 题目描述 }\) 这是一道模板题. \(n\) 个点,\(m\) 条边,每条边 \(e\) 有一个流量下界 \(\text{lower}(e)\) 和流量上界 \ ...
- 【bzoj2406】矩阵 二分+有上下界可行流
题目描述 输入 第一行两个数n.m,表示矩阵的大小. 接下来n行,每行m列,描述矩阵A. 最后一行两个数L,R. 输出 第一行,输出最小的答案: 样例输入 2 2 0 1 2 1 0 1 样例输出 1 ...
最新文章
- slave-pxc后GTID不一致
- Pytorch基础(十)——优化器(SGD,Adagrad,RMSprop,Adam,LBFGS等)
- java中接口私有反方_Java 8:在接口中声明私有和受保护的方法
- 开源的商业意义_为开源项目提供资金具有良好的商业意义
- 第八届“图灵杯”NEUQ-ACM程序设计竞赛个人赛——A题 切蛋糕
- GTK+ 3.5.18 发布,GUI 开发工具包
- python数据驱动读取用例_利用Python如何实现数据驱动的接口自动化测试
- 数据分析进阶-Excel绘制分段折线图
- 2022危险化学品经营单位主要负责人操作证考试题及模拟考试
- 数字转换大写英文字母
- 痛苦的挣扎--msp430g2553我恨你!
- 谁教会老公出轨外面养情人
- 开源电子商城系统:罗列几个电子商城系统,和一个不错的开源电子商城项目:mall,先做技术调研,主要还是学习代码。
- 中国大学 MOOC Android 性能优化:冷启动优化总结
- 中文核心期刊目录(2008年版).
- Go 语言如何自定义 linter(静态检查工具)
- 让OpenAi给我写个JS的set对象的笔记和快速去重方法
- InteractiveDataDisplay.WPF 固定坐标轴
- 详解独角兽应该具备的6个特点,创业者必看
- 杂谈:SDWAN二级运营商的大杀器 1
热门文章
- 想赚钱,你对钱敏感么?
- *第十五周*数据结构实践项目一【验证哈希表及其算法】
- C/C++ Bug记录
- Magento开发文档(二):Magento配置
- .Net中如何操作IIS(原理篇)
- iOS/Android React Native 配置教程
- 基于JAVA语言的selenium测试基础总结
- 544B. Sea and Islands
- 游戏开发中的数学和物理算法(18):缩放
- SharePoint And Ajax Technology(2):Ajax Control Toolkit学习