BZOJ 3698 XWW的难题:有上下界的最大流
题意
给你一个 $ n*n $ 的正实数矩阵 $ A $ ,满足XWW性。
称一个 $ n*n $ 的矩阵满足XWW性当且仅当:
- $ A[n][n] = 0 $
- 矩阵中每行的最后一个元素等于该行前 $ n-1 $ 个数的和(除最后一行)
- 矩阵中每列的最后一个元素等于该列前 $ n-1 $ 个数的和(除最后一列)
现在你要给 $ A $ 中的数进行取整操作(可以是上取整或者下取整),使得最后的 $ A $ 矩阵仍然满足XWW性。
问你 $ A $ 中元素之和最大为多少。如果无解,输出"No"。
题解
考虑将每一行和每一列看做一个点。
首先从源点向每一行 $ R(i) $ 连一条上下界分别为 $ (\lfloor A[i][n] \rfloor, \lceil A[i][n] \rceil) $ 的边,从每一列 $ C(i) $ 向汇点连一条上下界分别为 $ (\lfloor A[n][i] \rfloor, \lceil A[n][i] \rceil) $ 的边。
然后对于每一个 $ A[i][j] $ 来说,连一条从 $ R(i) $ 到 $ C(i) $ 的上下界为 $ (\lfloor A[i][j] \rfloor, \lceil A[i][j] \rceil) $ 边。
这样就保证了最大流一定满足了后两个条件。
然后跑有上下界的有源汇最大流就好。
AC Code
#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <vector>
#include <queue>
#define MAX_N 205
#define INF 1000000000
#define r(x) (x)
#define c(x) (n-1+(x))
using namespace std;
struct Edge
{
int dst,cap,rev;
Edge(int _dst,int _cap,int _rev) { dst=_dst,cap=_cap,rev=_rev; }
Edge(){}
};
int n,s,t,S,T,tot,dif=0;
int a[MAX_N];
int it[MAX_N];
int lv[MAX_N];
double w[MAX_N][MAX_N];
vector<Edge> edge[MAX_N];
queue<int> q;
void read()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
scanf("%lf",&w[i][j]);
}
}
}
inline void add(int s,int t,int c)
{
edge[s].push_back(Edge(t,c,edge[t].size()));
edge[t].push_back(Edge(s,0,edge[s].size()-1));
}
void build()
{
s=(n<<1)-1,t=s+1,S=t+1,T=S+1,tot=T;
for(int i=1;i<n;i++)
{
add(s,r(i),ceil(w[i][n])-floor(w[i][n]));
add(c(i),t,ceil(w[n][i])-floor(w[n][i]));
a[s]-=floor(w[i][n]),a[r(i)]+=floor(w[i][n]);
a[c(i)]-=floor(w[n][i]),a[t]+=floor(w[n][i]);
}
for(int i=1;i<n;i++)
{
for(int j=1;j<n;j++)
{
add(r(i),c(j),ceil(w[i][j])-floor(w[i][j]));
a[r(i)]-=floor(w[i][j]),a[c(j)]+=floor(w[i][j]);
}
}
add(t,s,INF);
for(int i=1;i<=(n<<1);i++)
{
if(a[i]>0) dif+=a[i],add(S,i,a[i]);
else if(a[i]<0) add(i,T,-a[i]);
}
}
void bfs(int s)
{
memset(lv+1,0,sizeof(int)*tot);
q.push(s),lv[s]=1;
while(!q.empty())
{
int x=q.front(); q.pop();
for(int i=0;i<edge[x].size();i++)
{
Edge temp=edge[x][i];
if(temp.cap>0 && !lv[temp.dst])
{
lv[temp.dst]=lv[x]+1;
q.push(temp.dst);
}
}
}
}
int dfs(int x,int t,int f)
{
if(x==t) return f;
for(int &i=it[x];i<edge[x].size();i++)
{
Edge &temp=edge[x][i];
if(temp.cap>0 && lv[x]<lv[temp.dst])
{
int d=dfs(temp.dst,t,min(f,temp.cap));
if(d>0)
{
temp.cap-=d;
edge[temp.dst][temp.rev].cap+=d;
return d;
}
}
}
return 0;
}
int max_flow(int s,int t)
{
int ans=0,f;
while(true)
{
bfs(s);
if(!lv[t]) return ans;
memset(it+1,0,sizeof(int)*tot);
while((f=dfs(s,t,INF))>0) ans+=f;
}
}
void work()
{
build();
int now=max_flow(S,T);
if(now!=dif)
{
printf("No\n");
return;
}
printf("%d\n",max_flow(s,t)*3);
}
int main()
{
read();
work();
}
BZOJ 3698 XWW的难题:有上下界的最大流的更多相关文章
- 【BZOJ3698】XWW的难题 有上下界的最大流
[BZOJ3698]XWW的难题 Description XWW是个影响力很大的人,他有很多的追随者.这些追随者都想要加入XWW教成为XWW的教徒.但是这并不容易,需要通过XWW的考核.XWW给你出了 ...
- BZOJ 3698: XWW的难题 [有源汇上下界最大流]
3698: XWW的难题 题意:(1)A[N][N]=0:(2)矩阵中每行的最后一个元素等于该行前N-1个数的和:(3)矩阵中每列的最后一个元素等于该列前N-1个数的和.给A中的数进行取整操作(可以是 ...
- 【bzoj3698】XWW的难题 有上下界最大流
题目描述 XWW是个影响力很大的人,他有很多的追随者.这些追随者都想要加入XWW教成为XWW的教徒.但是这并不容易,需要通过XWW的考核.XWW给你出了这么一个难题:XWW给你一个N*N的正实数矩阵A ...
- BZOJ 3698: XWW的难题(有源汇上下界最大流)
题面 XWW是个影响力很大的人,他有很多的追随者.这些追随者都想要加入XWW教成为XWW的教徒.但是这并不容易,需要通过XWW的考核. XWW给你出了这么一个难题:XWW给你一个N*N的正实数矩阵A, ...
- bzoj 3698 XWW的难题(有源汇的上下界最大流)
[题意] 对每个格子确定上下取整,使得满足1.A[n][n]=0 2.每行列前n-1个之和为第n个 3.格子之和尽量大. [思路] 设格子(i,j)上下取整分别为up(i,j)down(i,j),构图 ...
- BZOJ.3698.XWW的难题(有源汇上下界最大流ISAP)
题目链接 按套路行列作为两部分,连边 \(S->row->column->T\). S向代表行的元素连边cap(A[i][n])(容量上下界为上下取整),代表列的元素向T连边cap( ...
- BZOJ 3698: XWW的难题
Description XWW是个影响力很大的人,他有很多的追随者.这些追随者都想要加入XWW教成为XWW的教徒.但是这并不容易,需要通过XWW的考核. XWW给你出了这么一个难题:XWW给你一个N ...
- bzoj 2406 二分+有源有汇上下界网络流可行流判定
弱爆了,典型的行列建模方式,居然想不到,题做少了,总结少了...... 二分答案mid s----------------------->i行-----------------------> ...
- 3698: XWW的难题[有源汇上下界最大流]
3698: XWW的难题 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 354 Solved: 178[Submit][Status][Discus ...
随机推荐
- 【Python之路】第二十篇--MySQL(二)
视图 视图是一个虚拟表(非真实存在),其本质是[根据SQL语句获取动态的数据集,并为其命名], 用户使用时只需使用[名称]即可获取结果集,并可以将其当作表来使用. 1.创建视图 --格式:CREATE ...
- pro-select-like
DELIMITER | drop procedure if exists pro_query; CREATE PROCEDURE pro_query ( cname VARCHAR() ) BEGIN ...
- MongoDB资料汇总(转)
原文:MongoDB资料汇总 上一篇Redis资料汇总专题很受大家欢迎,这里将MongoDB的系列资料也进行了简单整理.希望能对大家有用. 最后更新时间:2013-04-22 1.MongoDB是什么 ...
- npm更新和nodejs更新
npm更新和nodejs更新 更新你已经安装的NPM库,这个很简单,只需要运行. npm update -g 更新Nodejs自身.一直依赖我都是下载最新版的源码,然后make install,及其繁 ...
- java 多线程 day05 线程范围内的数据共享
import java.util.HashMap;import java.util.Map;import java.util.Random;/** * Created by chengtao on 1 ...
- wtforms Form实例化流程(源码解析)
class LoginForm(Form): #首先执行后得到的结果是UnboundField()对象 name=simple.StringField( label='用户名', validato ...
- django-ORM复习补充
建表 class Author(models.Model): name = models.CharField(max_length=32) age = models.IntegerField() # ...
- day6 字典的增减查删
字典的使用 一.字典的特性 字典是无序的,数据关联性强,键值对,唯一一个映射的数据类型 字典的键必须是可哈希的(不可变的数据类型:字符串,数字,布尔值,元祖),并且是唯一的 不可哈希的(可变的数据类型 ...
- 第七课 GDB调试 (下)
1序言: 通过前面一节第六课 GDB调试 (下)文章,可以掌握理解了gdb调试:怎么启动.运行,打断点.查看变量.甚至改变变量等的知识,今天来大概讲解下调试bug的类型. 2知识点: 2.1 就像之前 ...
- JavaScript:判断当前浏览器是否为微信浏览器
<script type="text/javascript"> function CheckWeChatBrowser() { var ua = navigator.u ...