题目:

Description

We are supposed to make a budget proposal for this multi-site competition. The budget proposal is a matrix where the rows represent different kinds of expenses and the columns represent different sites. We had a meeting about this, some time ago where we discussed the sums over different kinds of expenses and sums over different sites. There was also some talk about special constraints: someone mentioned that Computer Center would need at least 2000K Rials for food and someone from Sharif Authorities argued they wouldn't use more than 30000K Rials for T-shirts. Anyway, we are sure there was more; we will go and try to find some notes from that meeting.

And, by the way, no one really reads budget proposals anyway, so we'll just have to make sure that it sums up properly and meets all constraints.

Input

The first line of the input contains an integer N, giving the number of test cases. The next line is empty, then, test cases follow: The first line of each test case contains two integers, m and n, giving the number of rows and columns (m <= 200, n <= 20). The second line contains m integers, giving the row sums of the matrix. The third line contains n integers, giving the column sums of the matrix. The fourth line contains an integer c (c < 1000) giving the number of constraints. The next c lines contain the constraints. There is an empty line after each test case.

Each constraint consists of two integers r and q, specifying some entry (or entries) in the matrix (the upper left corner is 1 1 and 0 is interpreted as "ALL", i.e. 4 0 means all entries on the fourth row and 0 0 means the entire matrix), one element from the set {<, =, >} and one integer v, with the obvious interpretation. For instance, the constraint 1 2 > 5 means that the cell in the 1st row and 2nd column must have an entry strictly greater than 5, and the constraint 4 0 = 3 means that all elements in the fourth row should be equal to 3.

Output

For each case output a matrix of non-negative integers meeting the above constraints or the string "IMPOSSIBLE" if no legal solution exists. Put one empty line between matrices.

Sample Input

  1. 2
  2.  
  3. 2 3
  4. 8 10
  5. 5 6 7
  6. 4
  7. 0 2 > 2
  8. 2 1 = 3
  9. 2 3 > 2
  10. 2 3 < 5
  11.  
  12. 2 2
  13. 4 5
  14. 6 7
  15. 1
  16. 1 1 > 10

Sample Output

  1. 2 3 3
  2. 3 3 4
  3.  
  4. IMPOSSIBLE

Source

题解:

算是求有源汇上下界最大流的一道模板题···只不过建图有点麻烦····

先说有源汇上下界最大流最小流的基本求法(引用自农民伯伯-Coding):

有源汇网络的可行流

对于流量有上下界的有源汇网络,原图中存在源点S和汇点T,为了求可行流,先将其转换为无源汇网络。

从T-->S引入一条边,其流量上下界为[0, INF],此时原图变为了无源汇网络,然后按照无源汇网络求解可行流的方法来求解。 

有源汇网络的最大流

要求最大流,先求可行流,通过“有源汇网络的可行流”的求解方法来判断有源汇网络存在可行流。 
若存在可行流,记从S流出的流量sum1,然后将T-->S的边取消,再次从S到T求解网络的最大流,记从S流出的流量sum2. 那么该有源汇网络的最大流为 sum1 + sum2. 
其中,sum1是在网络满足流量下界的条件下,从源点S流出的流量;求出sum1之后,网络中可能还有余量可以继续增广,那么再次求解从S到T的最大流,得到sum2,sum1 + sum2即为最终的最大流。

有源汇网络的最小流

求解有源汇网络最小流分为以下几步: 
(1)对SS到TT求一次最大流,即为f1.(在有源汇的情况下,先把整个网络趋向必须边尽量满足的情况); 
(2)添加一条边T-->S,流量上限为INF,这条边即为P.(构造无源网络) 
(3)对SS到TT再次求最大流,即为f2。(判断此时的网络中存在可行流,同时求出SS-->TT最大流) 
    如果所有必须边都满流,证明存在可行解,原图的最小流为“流经边P的流量”(原图已经构造成无源汇网络,对于S同样满足 入流和==出流和,只有新添加的边流向S,而S的出流就是原图的最小流)。 
注: 最小流求法的正确性不知如何证明,待继续学习。

代码:

  1. #include<iostream>
  2. #include<cstdio>
  3. #include<cstdlib>
  4. #include<cmath>
  5. #include<ctime>
  6. #include<cctype>
  7. #include<cstring>
  8. #include<string>
  9. #include<algorithm>
  10. using namespace std;
  11. const int N=;
  12. const int M=;
  13. const int inf=0x3f3f3f3f;
  14. int tot=,first[N],next[M],go[M],rest[M],lev[N],cur[N];
  15. int tr[N],tc[N],n,m,low[N][N],up[N][N],T,src,des,ss,tt,sum1,sum2,k,ans,id[N][N],p;
  16. bool jud[N][N];
  17. inline void comb(int a,int b,int c)
  18. {
  19. next[++tot]=first[a],first[a]=tot,go[tot]=b,rest[tot]=c;
  20. next[++tot]=first[b],first[b]=tot,go[tot]=a,rest[tot]=;
  21. }
  22. inline void comb1(int a,int b,int c)
  23. {
  24. next[tot]=first[a],first[a]=tot,go[tot]=b,rest[tot]=c;
  25. next[++tot]=first[b],first[b]=tot,go[tot]=a,rest[tot]=;
  26. }
  27. inline int R()
  28. {
  29. char c;
  30. int f=;
  31. for(c=getchar();c<''||c>'';c=getchar());
  32. for(;c<=''&&c>='';c=getchar())
  33. f=(f<<)+(f<<)+c-'';
  34. return f;
  35. }
  36. inline void pre()
  37. {
  38. memset(first,,sizeof(first));
  39. memset(low,,sizeof(low));
  40. memset(up,inf,sizeof(up));
  41. memset(tr,,sizeof(tr));
  42. memset(tc,,sizeof(tc));
  43. tot=,src=,des=n+m+,ss=n+m+,tt=n+m+,sum1=,sum2=,ans=;
  44. }
  45. inline bool bfs()
  46. {
  47. for(int i=src;i<=tt;i++) cur[i]=first[i],lev[i]=-;
  48. static int que[N],tail,u,v;
  49. que[tail=]=ss;
  50. lev[ss]=;
  51. for(int head=;head<=tail;head++)
  52. {
  53. u=que[head];
  54. for(int e=first[u];e;e=next[e])
  55. {
  56. if(lev[v=go[e]]==-&&rest[e])
  57. {
  58. lev[v]=lev[u]+;
  59. que[++tail]=v;
  60. if(v==tt) return true;
  61. }
  62. }
  63. }
  64. return false;
  65. }
  66. inline int dinic(int u,int flow)
  67. {
  68. if(u==tt)
  69. return flow;
  70. int res=,delta,v;
  71. for(int &e=cur[u];e;e=next[e])
  72. {
  73. if(lev[v=go[e]]>lev[u]&&rest[e])
  74. {
  75. delta=dinic(v,min(flow-res,rest[e]));
  76. if(delta)
  77. {
  78. rest[e]-=delta;
  79. rest[e^]+=delta;
  80. res+=delta;
  81. if(res==flow) break;
  82. }
  83. }
  84. }
  85. if(flow!=res) lev[u]=-;
  86. return res;
  87. }
  88. inline void maxflow()
  89. {
  90. while(bfs())
  91. ans+=dinic(ss,inf);
  92. }
  93. int main()
  94. {
  95. //freopen("a.in","r",stdin);
  96. T=R();
  97. int a,b,c;
  98. char s[];
  99. while(T--)
  100. {
  101. scanf("\n");
  102. n=R(),m=R();
  103. pre();
  104. for(int i=;i<=n;i++)
  105. {
  106. a=R();
  107. tc[src]+=a,tr[i]+=a;
  108. sum1+=a;
  109. }
  110. for(int i=;i<=m;i++)
  111. {
  112. a=R();
  113. tr[des]+=a,tc[n+i]+=a;
  114. sum2+=a;
  115. }
  116. if(sum1!=sum2)
  117. {
  118. cout<<"IMPOSSIBLE"<<endl;
  119. cout<<endl;
  120. continue;
  121. }
  122. k=R();
  123. while(k--)
  124. {
  125. scanf("%d%d%s%d",&a,&b,s,&c);
  126. if(a==&&b==)
  127. {
  128. for(int i=;i<=n;i++)
  129. for(int j=;j<=m;j++)
  130. {
  131. if(s[]=='=')
  132. {
  133. low[i][j]=max(low[i][j],c);
  134. up[i][j]=min(up[i][j],c);
  135. }
  136. if(s[]=='<')
  137. up[i][j]=min(up[i][j],c-);
  138. if(s[]=='>')
  139. low[i][j]=max(low[i][j],c+);
  140. }
  141. }
  142. else if(a==)
  143. for(int i=;i<=n;i++)
  144. {
  145. if(s[]=='=')
  146. {
  147. low[i][b]=max(low[i][b],c);
  148. up[i][b]=min(up[i][b],c);
  149. }
  150. if(s[]=='<')
  151. up[i][b]=min(up[i][b],c-);
  152. if(s[]=='>')
  153. low[i][b]=max(low[i][b],c+);
  154. }
  155. else if(b==)
  156. for(int i=;i<=m;i++)
  157. {
  158. if(s[]=='=')
  159. {
  160. low[a][i]=max(low[a][i],c);
  161. up[a][i]=min(up[a][i],c);
  162. }
  163. if(s[]=='<')
  164. up[a][i]=min(up[a][i],c-);
  165. if(s[]=='>')
  166. low[a][i]=max(low[a][i],c+);
  167. }
  168. else
  169. {
  170. if(s[]=='=')
  171. {
  172. low[a][b]=max(low[a][b],c);
  173. up[a][b]=min(up[a][b],c);
  174. }
  175. if(s[]=='<')
  176. up[a][b]=min(up[a][b],c-);
  177. if(s[]=='>')
  178. low[a][b]=max(low[a][b],c+);
  179. }
  180. }
  181. bool flag=false;
  182. for(int i=;i<=n;i++)
  183. {
  184. if(flag==true) break;
  185. for(int j=;j<=m;j++)
  186. {
  187. if(low[i][j]>up[i][j])
  188. {
  189. flag=true;
  190. break;
  191. }
  192. else
  193. {
  194. tr[n+j]+=low[i][j];
  195. tc[i]+=low[i][j];
  196. id[i][j]=++tot;
  197. comb1(i,n+j,up[i][j]-low[i][j]);
  198. }
  199. }
  200. }
  201. if(flag==true)
  202. {
  203. cout<<"IMPOSSIBLE"<<endl;
  204. cout<<endl;
  205. continue;
  206. }
  207. int cnt=;
  208. p=++tot;
  209. comb1(des,src,inf);
  210. for(int i=src;i<=des;i++)
  211. {
  212. if(tr[i]>tc[i])
  213. {
  214. comb(ss,i,tr[i]-tc[i]);
  215. cnt+=tr[i]-tc[i];
  216. }
  217. if(tc[i]>tr[i])
  218. comb(i,tt,tc[i]-tr[i]);
  219. }
  220. maxflow();
  221. if(ans!=cnt)
  222. {
  223. cout<<"IMPOSSIBLE"<<endl;
  224. cout<<endl;
  225. continue;
  226. }
  227. rest[p]=rest[p^]=;
  228. maxflow();
  229. for(int i=;i<=n;i++)
  230. {
  231. for(int j=;j<=m;j++)
  232. cout<<rest[id[i][j]^]+low[i][j]<<" ";
  233. cout<<endl;
  234. }
  235. cout<<endl;
  236. }
  237. return ;
  238. }

算法复习——有源汇上下界可行流(bzoj2396)的更多相关文章

  1. POJ2396 Budget [有源汇上下界可行流]

    POJ2396 Budget 题意:n*m的非负整数矩阵,给出每行每列的和,以及一些约束关系x,y,>=<,val,表示格子(x,y)的值与val的关系,0代表整行/列都有这个关系,求判断 ...

  2. 有源汇上下界可行流(POJ2396)

    题意:给出一个n*m的矩阵的每行和及每列和,还有一些格子的限制,求一组合法方案. 源点向行,汇点向列,连一条上下界均为和的边. 对于某格的限制,从它所在行向所在列连其上下界的边. 求有源汇上下界可行流 ...

  3. 计蒜客 31447 - Fantastic Graph - [有源汇上下界可行流][2018ICPC沈阳网络预赛F题]

    题目链接:https://nanti.jisuanke.com/t/31447 "Oh, There is a bipartite graph.""Make it Fan ...

  4. poj2396 Budget(有源汇上下界可行流)

    [题目链接] http://poj.org/problem?id=2396 [题意] 知道一个矩阵的行列和,且知道一些格子的限制条件,问一个可行的方案. [思路] 设行为X点,列为Y点,构图:连边(s ...

  5. poj2396有源汇上下界可行流

    题意:给一些约束条件,要求算能否有可行流,ps:刚开始输入的是每一列和,那么就建一条上下界相同的边,这样满流的时候就一定能保证流量相同了,还有0是该列(行)对另一行每个点都要满足约束条件 解法:先按无 ...

  6. ZOJ1994有源汇上下界可行流

    http://fastvj.rainng.com/contest/236779#problem/G Description: n 行 m 列 给你行和 与 列和 然后有Q个限制,表示特定单元格元素大小 ...

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

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2406 这题,首先把题目那个式子的绝对值拆成两个限制,就成了网络流的上下界: 有上下界可行流原 ...

  8. bzoj千题计划158:bzoj2406: 矩阵(有源汇上下界可行流)

    http://www.lydsy.com/JudgeOnline/problem.php?id=2406 设矩阵C=A-B 最小化 C 一行或一列和的最大值 整体考虑一行或者一列的和 二分最大值 这样 ...

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

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

随机推荐

  1. LINUX 安装JDK (rpm格式和tar.gz格式)

    谷歌博客地址:http://tsaiquinn.blogspot.com/2014/10/linux-jdk-rpmtargz.html JDK rpm方式: 我使用的是SecureCRT,先下载了然 ...

  2. Python 元组、字典、集合操作总结

    元组 a=('a',) a=('a','b') 特点 有序 不可变,不可以修改元组的值,无法为元组增加或者删除元素 元组的创建 a=('a',) a=('a','b') tuple('abcd') 转 ...

  3. UVA821 PageHopping (Floyd)

    求所有点直接的平均最短距离,保存一下出现过的点,题目保证是所有点连通,Floyd求出最短路以后两个for统计一下. #include<bits/stdc++.h> using namesp ...

  4. 利用python实现整数转换为任意进制字符串

    假设你想将一个整数转换为一个二进制和十六进制字符串.例如,将整数 10 转换为十进制字符串表示为 10 ,或将其字符串表示为二进制 1010 . 实现 以 2 到 16 之间的任何基数为参数: def ...

  5. C#中Lock关键字的使用

    C# 中的 Lock 语句通过隐式使用 Monitor 来提供同步功能.lock 关键字在块的开始处调用 Enter,而在块的结尾处调用 Exit. 通常,应避免锁定 public 类型,否则实例将超 ...

  6. robotframework接口测试实例

    *** Settings *** Library Collections Library RequestsLibrary *** Test Cases *** test Create Session ...

  7. shell脚本调试打印日志问题

    shell脚本调试打印日志问题 1. 需求 我们在编写脚本的时候,有时候需要做调试,便于我们定位问题,有时候等脚本上线之后,我们需要保留脚本执行过程中的记录.便于我们在出问题的时候,定位问题. 2. ...

  8. baidumap demo(二)

    接口说明 百度地图API提供的搜索服务包括:POI检索,多关键字检索,公交方案检索,驾车路线检索,步行路线检索,地理编码,反地理编码,公交详情检索,在线建议查询,短串分享. 所有检索请求接口均为异步接 ...

  9. vue建项目并使用

    今天来回顾下vue项目的建立和使用,好久不用感觉不会用了. 下面两个都要全局安装 首先安装git,地址  https://gitforwindows.org/ 安装node, 地址 https://n ...

  10. 有趣的this以及apply,call,bind方法

    看this指向谁,要看执行时而非定义时(箭头函数除外).函数没有绑定在对象上调用,非'strict'模式下,this指向window,否则为undefined 改变this指向的方法 1. apply ...