题解 P5038 [SCOI2012]奇怪的游戏
题解
做这题之前,做了一道叫星际战争的题,很容易想到二分 \(+\) 网络流,那么二分啥呢?
我们先推一下式子,因为是对相邻格子加数,那么可以联想到黑白染色类问题。
设有黑色格子 \(B\) 个,其格子中初始数的和为 \(b\),白色格子同理,个数为 \(W\) 个,初始权值和为 \(w\) 个,最后变成的同一个数为 \(num\)。
可以得出 \(B×num-b=W×num-w\) 化简得 \(num=\frac{b-w}{B-W}\)。
首先对于化简式,其必要条件是 \(B\neq W\),那么我们分两种情况讨论。
对于 \(B\neq W\) 的情况,那么我们可以直接算出 \(num\) ,但是我们要检验其是否合法。
对于 \(B=W\) 的情况,首先我们要保证 \(b=w\) 否则直接无解。然后,对于一个合法的 \(num\),\(num+1\) 也是合法的,因为 \(B=W\) 所以在 \(n,m\) 中一定有一个是偶数。于是我们就可以二分求解,找到临界的 \(num\) 就是答案。
这个 \(check()\) 我们可以用最大流求解,对每个黑格子,我们分别向源点 \(s\),权值为 \(num-w_i\),它能到的白格子,权值为 \(INF\) 连边,对每个白格子,我们向汇点 \(t\) 连边,权值为 \(num-w_i\)。
对这张图跑最大流,如果最大流等于 \(B×num-b\) 那么说明这张图跑满了,说明答案正确。
Code:
\(AC\kern 0.4em CODE:\)
#include<bits/stdc++.h>
#define ri register signed
#define p(i) ++i
using namespace std;
namespace IO{
char buf[1<<21],*p1=buf,*p2=buf;
#define gc() p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++
inline int read() {
ri x=0,f=1;char ch=gc();
while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=gc();}
while(ch>='0'&&ch<='9') {x=(x<<1)+(x<<3)+(ch^48);ch=gc();}
return x*f;
}
}
using IO::read;
namespace nanfeng{
#define int long long//懒人必备
#define cmax(x,y) ((x)>(y)?(x):(y))
#define cmin(x,y) ((x)>(y)?(y):(x))
#define FI FILE *IN
#define FO FILE *OUT
#undef bool
static const int N=45;
int ch[N][N],id[N][N],n,m,T,tot=1,B,b,W,w,mx;
namespace NetworkFlows{
#define jud(i,j) ((i)&&(j)&&(i)<=n&&(j)<=m)
static const int INF=1e18+7;
int first[N*N],dep[N*N],cur[N*N],que[N*N],t=2,sw,s=1,et;
signed dx[5]={0,1,0,-1,0},dy[5]={0,0,1,0,-1};
struct edge{int v,nxt,w;}e[N*N*5];
inline void add(int u,int v,int w) {
e[t].v=v,e[t].w=w;
e[t].nxt=first[u];
first[u]=t++;
}
inline void init() {et=n*m+2,t=2,sw=0;memset(first,0,sizeof(first));}
inline void build(int w) {
for (ri i(1);i<=n;p(i)) {
for (ri j(1);j<=m;p(j)) {
if (!((i+j)%2)) {
sw+=w-ch[i][j];
add(s,id[i][j],w-ch[i][j]);
add(id[i][j],s,0);
for (ri d(1);d<=4;p(d)) {
if (jud(i+dx[d],j+dy[d])) add(id[i][j],id[i+dx[d]][j+dy[d]],INF),add(id[i+dx[d]][j+dy[d]],id[i][j],0);
}
} else add(id[i][j],et,w-ch[i][j]),add(et,id[i][j],0);
}
}
}
inline bool bfs(int s,int t) {
memset(dep,0,sizeof(dep));
ri hd=1,tl=0;
dep[que[p(tl)]=s]=1;
cur[s]=first[s];
while(hd<=tl) {
s=que[hd++];
for (ri i(first[s]),v;i;i=e[i].nxt) {
if (e[i].w&&!dep[v=e[i].v]) {
dep[v]=dep[s]+1;
cur[que[p(tl)]=v]=first[v];
if (v==t) return 1;
}
}
}
return 0;
}
int dfs(int x,int flow) {
if (x==et||!flow) return flow;
int rst=flow;
for (ri i(cur[x]),v;i&&rst;i=e[i].nxt) {
if (e[i].w&&dep[v=e[i].v]==dep[x]+1) {
register int k;
if (!(k=dfs(v,cmin(e[i].w,rst)))) dep[v]=0;
e[i].w-=k,e[i^1].w+=k,rst-=k;
}
cur[x]=i;
}
return flow-rst;
}
inline int dinic() {
int res=0;
while(bfs(s,et)) res+=dfs(s,INF);
return res;
}
}
inline bool check(int w) {
NetworkFlows::init();
NetworkFlows::build(w);
return NetworkFlows::dinic()==NetworkFlows::sw;
}
inline void init() {mx=B=b=W=w=0,tot=1;}
inline int main() {
// FI=freopen("nanfeng.in","r",stdin);
// FO=freopen("nanfeng.out","w",stdout);
T=read();
for (ri z(1);z<=T;p(z)) {
init();
n=read(),m=read();
for (ri i(1);i<=n;p(i)) {
for (ri j(1);j<=m;p(j)) {
ch[i][j]=read(),id[i][j]=p(tot);
mx=cmax(mx,ch[i][j]);
if (!((i+j)%2)) p(B),b+=ch[i][j];
else p(W),w+=ch[i][j];
}
}
if (B!=W) {
int num=(b-w)/(B-W);
if (num>=mx&&check(num)) printf("%lld\n",B*num-b);
else puts("-1");
} else {
if (b!=w) {puts("-1");continue;}
int l=mx,r=1e11,res=-1;
while(l<=r) {
int mid((l+r)>>1);
if (check(mid)) r=mid-1,res=mid;
else l=mid+1;
}
printf("%lld\n",res==-1?-1:B*res-b);
}
}
return 0;
}
#undef int
}
int main() {return nanfeng::main();}
题解 P5038 [SCOI2012]奇怪的游戏的更多相关文章
- [题目] Luogu P5038 [SCOI2012]奇怪的游戏
学习资料 -----1----- -----2----- P5038 [SCOI2012]奇怪的游戏 一道甚神但没用到高深模型的题 思路 没思路,看题解吧 代码 #include <iostre ...
- P5038 [SCOI2012]奇怪的游戏 二分+网络流
$ \color{#0066ff}{ 题目描述 }$ Blinker最近喜欢上一个奇怪的游戏. 这个游戏在一个 \(N \times M\) 的棋盘上玩,每个格子有一个数.每次\(Blinker\)会 ...
- P5038 [SCOI2012]奇怪的游戏
题目链接 题意分析 首先我们需要求的是统一以后的值\(x\) 并且一般的棋盘操作我们都需要黑白染色 那么对于棋盘格子是偶数的情况的话 答案是存在单调性的 因为如果统一之后 两两搭配还是可以再加一个的 ...
- 洛谷$P5038\ [SCOI2012]$奇怪的游戏 二分+网络流
正解:二分+网络流 解题报告: 传送门$QwQ$ 这种什么,"同时增加",长得就挺网络流的$QwQ$?然后看到问至少加多少次,于是考虑加个二分呗?于是就大体确定了做题方向,用的网络 ...
- BZOJ 2756: [SCOI2012]奇怪的游戏 网络流/二分
2756: [SCOI2012]奇怪的游戏 Time Limit: 40 Sec Memory Limit: 128 MBSubmit: 1594 Solved: 396[Submit][Stat ...
- bzoj2756: [SCOI2012]奇怪的游戏(网络流+分情况)
2756: [SCOI2012]奇怪的游戏 题目:传送门 题解: 发现做不出来的大难题一点一个网络流 %大佬 首先黑白染色(原来是套路...)染色之后就可以保证每次操作都一定会使黑白各一个各自的值加1 ...
- BZOJ 2756: [SCOI2012]奇怪的游戏 [最大流 二分]
2756: [SCOI2012]奇怪的游戏 Time Limit: 40 Sec Memory Limit: 128 MBSubmit: 3352 Solved: 919[Submit][Stat ...
- Bzoj2756 [SCOI2012]奇怪的游戏
2756: [SCOI2012]奇怪的游戏 Time Limit: 40 Sec Memory Limit: 128 MBSubmit: 3220 Solved: 886 Description ...
- bzoj 2756 [SCOI2012]奇怪的游戏 二分+网络流
2756:[SCOI2012]奇怪的游戏 Time Limit: 40 Sec Memory Limit: 128 MBSubmit: 4926 Solved: 1362[Submit][Stat ...
随机推荐
- 从源码角度谈谈MySQL "Too many open files"错误的根本原因
"Too many open files"是一个比较常见的错误,不仅仅是在 MySQL 中.只要是在 Linux 中启动的进程,都有可能遇到这个错误. 究其原因,是进程打开的文件描 ...
- C语言:获取汉字的编码
#include <stdio.h> #include <locale.h> #include <wchar.h> int main() { setlocale(L ...
- xmind8 Mac序列号
1.首先去官网下载xmind8的安装包:XMind for Mac 也可以去我的百度网盘下载: 链接:https://pan.baidu.com/s/1eY52YsSaPmr-YFhB62Cli ...
- 【剑指offer】50.数组中重复出现的数字
50.数组中重复出现的数字 知识点:数组:Set的不可重复性 题目描述 在一个长度为n的数组里的所有数字都在0到n-1的范围内. 数组中某些数字是重复的,但不知道有几个数字是重复的.也不知道每个数字重 ...
- Pytest单元测试框架之简单操作示例
前言: Pytest是第三方单元格测试框架,更加简单,灵活,而且提供了更多丰富的扩展: Pytest与UnitTest框架的区别 UnitTest测试用例执行顺序是依照ascii码执行,而Pytest ...
- JavaScript实现拖放效果
JavaScript实现拖放效果 笔者实现该效果也是套用别人的轮子的.传送门 然后厚颜无耻的贴别人的readme~,笔者为了方便查阅就直接贴了,有不想移步的可以看这篇.不过还是最好请到原作者的GitH ...
- 分享一个自己画div的技巧
分享一个自己画div的技巧 笔者是小白,前端不是很懂.现在想总结下自己画div布局的小技巧和思路. 先对着设计图把div给好好框选出来 我个人觉得这一步是很重要的,要先分析大局,再细节处理.一定要先决 ...
- Redis双写一致性与缓存更新策略
一.双写一致性 双写一致性,也就是说 Redis 和 mysql 数据同步 双写一致性数据同步的方案有: 1.先更新数据库,再更新缓存 这个方案一般不用: 因为当有两个请求AB先后更新数据库后,A应该 ...
- Lniux上的tomcat部署web项目路径问题
一.项目路径及部署.启动位置 1.在tomcat下部署应用,部署路径:一般直接部署在tomcat/webapps/ROOT下即可.2.默认访问路径:tomcat的默认访问路径为webapps/ROOT ...
- TypeError: attrib() got an unexpected keyword argument 'convert'
使用pyinstaller -F aaa.py时,报错 TypeError: attrib() got an unexpected keyword argument 'convert' 没有exe生成 ...