题目描述

XWW是个影响力很大的人,他有很多的追随者。这些追随者都想要加入XWW教成为XWW的教徒。但是这并不容易,需要通过XWW的考核。
XWW给你出了这么一个难题:XWW给你一个N*N的正实数矩阵A,满足XWW性。
称一个N*N的矩阵满足XWW性当且仅当:(1)A[N][N]=0;(2)矩阵中每行的最后一个元素等于该行前N-1个数的和;(3)矩阵中每列的最后一个元素等于该列前N-1个数的和。
现在你要给A中的数进行取整操作(可以是上取整或者下取整),使得最后的A矩阵仍然满足XWW性。同时XWW还要求A中的元素之和尽量大。

输入

第一行一个整数N,N ≤ 100。
接下来N行每行包含N个绝对值小于等于1000的实数,最多一位小数。

输出

输出一行,即取整后A矩阵的元素之和的最大值。无解输出No。

样例输入

4
3.1 6.8 7.3 17.2
9.6 2.4 0.7 12.7
3.6 1.2 6.5 11.3
16.3 10.4 14.5 0

样例输出

129


题解

有上下界最大流

正常这种题应该是费用流建模,但实际上由于本题的特殊性质,只需要最大流即可解决。

首先一个数最小就是它向下取整,而如果是小数还可以取向上取整。

这里为了方便,设a[i][j]表示某个数的最小取值,b[i][j]表示某个数能否+1,能则为1,不能则为0.

那么容易想到建图:S->第i行最后一列,容量下界为a[i][n],上界为a[i][n]+b[i][n];第i行最后一列->最后一行第j列,容量下界为a[i][j],上界为a[i][j]+b[i][j];最后一行第j列->T,容量下界为a[n][j],上界为a[n][j]+b[n][j]。

然后正常人都会发现这是有上下界最大费用流,而实际上很多题解都是有上下界最大流,为什么?

因为本题的特殊性,一条增广路一定会经过3条边,所以费用为流量*3。因此只需要求出最大流,乘以3即为费用。

至于有上下界最大流的具体求法:对于某条边x->y,容量下界为z,上界为z+w,连x->y,容量为w;SS->y&x->TT,容量为z;加T->S,容量为inf。从SS到TT跑最大流,未满流则无解,满流则记录T->S的容量(即反向边的残量)为ans1;再去掉与SS或TT相连的边,去掉T->S的边,从S到T跑最大流为ans2,ans1+ans2为解。

另外本题规定的a[n][n]=0实际上并无卵用,直接无视就好。

#include <cstdio>
#include <cstring>
#include <queue>
#define N 300
#define M 100000
#define inf 0x3fffffff
using namespace std;
queue<int> q;
int u[N][N] , v[N][N] , head[N] , to[M] , val[M] , next[M] , cnt = 1 , b , e , s , t , dis[N];
void add(int x , int y , int z)
{to[++cnt] = y , val[cnt] = z , next[cnt] = head[x] , head[x] = cnt;to[++cnt] = x , val[cnt] = 0 , next[cnt] = head[y] , head[y] = cnt;
}
bool bfs()
{int x , i;memset(dis , 0 , sizeof(dis));while(!q.empty()) q.pop();dis[s] = 1 , q.push(s);while(!q.empty()){x = q.front() , q.pop();for(i = head[x] ; i ; i = next[i]){if(val[i] && !dis[to[i]]){dis[to[i]] = dis[x] + 1;if(to[i] == t) return 1;q.push(to[i]);}}}return 0;
}
int dinic(int x , int low)
{if(x == t) return low;int temp = low , i , k;for(i = head[x] ; i ; i = next[i]){if(val[i] && dis[to[i]] == dis[x] + 1){k = dinic(to[i] , min(temp , val[i]));if(!k) dis[to[i]] = 0;val[i] -= k , val[i ^ 1] += k;if(!(temp -= k)) break;}}return low - temp;
}
int main()
{int n , i , j , sum = 0 , ans = 0;double tmp;scanf("%d" , &n) , b = 0 , e = 2 * n - 1 , s = 2 * n , t = 2 * n + 1 , add(e , b , inf);for(i = 1 ; i <= n ; i ++ )for(j = 1 ; j <= n ; j ++ )scanf("%lf" , &tmp) , u[i][j] = (int)tmp , v[i][j] = (tmp > u[i][j]) , sum += u[i][j];for(i = 1 ; i < n ; i ++ ) add(b , i , v[i][n]) , add(s , i , u[i][n]) , add(b , t , u[i][n]);for(i = 1 ; i < n ; i ++ ) add(i + n - 1 , e , v[n][i]) , add(s , e , u[n][i]) , add(i + n - 1 , t , u[n][i]);for(i = 1 ; i < n ; i ++ )for(j = 1 ; j < n ; j ++ )add(i , j + n - 1 , v[i][j]) , add(s , j + n - 1 , u[i][j]) , add(i , t , u[i][j]);while(bfs()) sum -= dinic(s , inf);if(sum){printf("No\n");return 0;}ans = val[3];for(i = head[s] ; i ; i = next[i]) val[i] = val[i ^ 1] = 0;for(i = head[t] ; i ; i = next[i]) val[i] = val[i ^ 1] = 0;val[2] = val[3] = 0 , s = b , t = e;while(bfs()) ans += dinic(s , inf);printf("%d\n" , ans * 3);return 0;
}

转载于:https://www.cnblogs.com/GXZlegend/p/7072443.html

【bzoj3698】XWW的难题 有上下界最大流相关推荐

  1. bzoj3698: XWW的难题 有上下界的网络流

    bzoj3698: XWW的难题 Description XWW是个影响力很大的人,他有很多的追随者.这些追随者都想要加入XWW教成为XWW的教徒.但是这并不容易,需要通过XWW的考核. XWW给你出 ...

  2. [bzoj3698]XWW的难题——有上下界的最大流

    题目大意: XWW是个影响力很大的人,他有很多的追随者.这些追随者都想要加入XWW教成为XWW的教徒.但是这并不容易,需要通过XWW的考核. XWW给你出了这么一个难题:XWW给你一个N*N的正实数矩 ...

  3. BZOJ3698 XWW的难题(上下界网络流+gap优化的细节处理)

    XWW给你一个N*N的正实数矩阵A,满足XWW性. 称一个N*N的矩阵满足XWW性当且仅当:(1)A[N][N]=0:(2)矩阵中每行的最后一个元素等于该行前N-1个数的和:(3)矩阵中每列的最后一个 ...

  4. 【BZOJ3698】XWW的难题 有上下界的最大流

    [BZOJ3698]XWW的难题 Description XWW是个影响力很大的人,他有很多的追随者.这些追随者都想要加入XWW教成为XWW的教徒.但是这并不容易,需要通过XWW的考核. XWW给你出 ...

  5. BZOJ 3698: XWW的难题|有上下界的网络流之最大流

    先构建可行流的模型 s−>s->每一行 流量为(a[i][n],a[i][n]+1)(a[i][n],a[i][n]+1) 每一列−>t->t 流量为(a[n][i],a[n] ...

  6. BZOJ 3698: XWW的难题(有源汇上下界最大流)

    题面 XWW是个影响力很大的人,他有很多的追随者.这些追随者都想要加入XWW教成为XWW的教徒.但是这并不容易,需要通过XWW的考核. XWW给你出了这么一个难题:XWW给你一个N*N的正实数矩阵A, ...

  7. 3698: XWW的难题 有源汇上下界最大流

    有源汇上下界最大流,行列建图,下界为ai,ja_{i,j},上界为ai,j+1a_{i,j}+1,跑的飞起. 不要忘记判断上取整下取整相同的情况. #include<iostream> # ...

  8. bzoj3698 XWW的难题

    bzoj3698 Description XWW是个影响力很大的人,他有很多的追随者.这些追随者都想要加入XWW教成为XWW的教徒.但是这并不容易,需要通过XWW的考核. XWW给你出了这么一个难题: ...

  9. [bzoj3698]XWW的难题

    传送门 有源汇上下界最大流,建模很经典,将每行和每列看做一个点,然后去按照有源汇上下界建图就行了,因为这个题,我又发现了自己网络流板子里的一个错误写法,bfs的时候得增广完再退出,改了我半个小时... ...

最新文章

  1. BZOJ-2049 [SDOI2008]洞穴勘测
  2. python中xlsxwriter_python xlsxwriter使用方法汇总
  3. [zz]很详细,涵盖了多数场景!推荐 - python 的日志logging模块学习
  4. paste linux 相同字符,Linux命令之字符串处理命令paste命令使用实例
  5. 宝塔linux修改默认编码,宝塔linux面板防护CC设置(示例代码)
  6. Android中通过ActionBar为标题栏添加搜索以及分享视窗
  7. 做游戏,学编程(C语言) 10 僵尸危机
  8. 美国空手道俱乐部的 Python 库原来长这样!
  9. springboot 前缀_SpringBoot2.x 给Controller的RequestMapping添加统一前缀
  10. C++中for循环的5种语法
  11. 计算机网络基础ios指令,IOS快捷指令制作真正的贴吧每天全自动签到
  12. 这45个场景,正在被区块链抽筋扒皮…
  13. 洛阳地铁一号线无人驾驶_无人驾驶超级拖拉机现身洛阳
  14. css 剪辑图片_css剪裁GIF背景图片动画特效
  15. LaTeX 编辑器对比
  16. 相遇,是天意; 相守,是人意
  17. 回文树(模板)+ 例题
  18. Libuv源码分析 —— 9. DNS
  19. mysql优化 wordpress_WordPress优化Mysql InnoDB转MYISAM
  20. 广州市海珠区投资促进局领导一行莅临袋鼠云参观交流

热门文章

  1. 西游记中孙悟空的等级地位
  2. linux运行python脚本_将Python脚本文件包装成可执行文件
  3. 腾讯,百度,高德地图兴趣点(POI)的获取以及查询,逆解析解析
  4. 怎样用 16.7 个小时做 40 小时的工作
  5. 【优化覆盖】基于matlab粒子群算法优化无人机编队布局求解车载网络通信覆盖优化问题【含Matlab源码 2021期】
  6. 本地同城地方门户网这样运营轻轻松松拿到5W广告费
  7. Verilgo实现的FPGA奇偶校验
  8. 播放数万 | 一文学会Cytoscape网络图绘制 - 最新教程
  9. 文件上传绕过的一些姿势
  10. 微众银行七年,从拓荒者到深耕者