题意:

Blinker最近喜欢上一个奇怪的游戏。
这个游戏在一个 N*M 的棋盘上玩,每个格子有一个数。每次 Blinker 会选择两个相邻
的格子,并使这两个数都加上 1。
现在 Blinker 想知道最少多少次能使棋盘上的数都变成同一个数,如果永远不能变成同
一个数则输出-1。

n,m<=40 a[i,j]<=10^9

思路:假设我们已经询问最终是否都能变成数字D怎么做?

注意到黑白点染色后每次加的点必定是1黑点1白点

设黑格子有k1个,和为s1,白格子有k2个,和为s2

显然k1*d-s1=k2*d-s2

d=(s1-s2)/(k1-k2)

当k1<>k2时D可以直接被计算出,用网络流验证即可,注意当D<max(a[i,j])时显然不可能完成

当k1=k2时注意到每次加必定黑+1,白+1,如果s1<>s2就不可能完成

又因为此时n*m为偶数,任意的D如果符合要求,任意的D'>D必定符合要求

因为总个数为偶数时必定可以设计出方案使所有数都+1

所以可以二分答案,找到最小的能被取到的D,同样用网络流验证

网络流验证:

S——>黑点 流量为D-a[i,j]

白点——>T 流量为D-a[i,j]

黑点——>相邻的白点 流量为oo

看最小割是否等于需要被加上的总和即可,因为每个黑点只向相邻的白点连边所以一个最小割就对应了一个操作的方案

另外二分上界只需要2^50,太大会炸

num数组一定要动态开,直接设完全部会导致标号不连续调了一天全是泪

 const oo=<<;
dx:array[..]of longint=(-,,,);
dy:array[..]of longint=(,,-,);
var head,gap,dis:array[..]of longint;
vet,next,fan:array[..]of longint;
len:array[..]of int64;
num,a:array[..,..]of longint;
n,m,i,j,tot,k1,k2,s,source,src,cas,v,mx:longint;
l,r,mid,last,s1,s2,k:int64; procedure add(a,b:longint;c:int64);
begin
inc(tot);
next[tot]:=head[a];
vet[tot]:=b;
len[tot]:=c;
head[a]:=tot; inc(tot);
next[tot]:=head[b];
vet[tot]:=a;
len[tot]:=;
head[b]:=tot;
end; function min(x,y:int64):int64;
begin
if x<y then exit(x);
exit(y);
end; function dfs(u:longint;aug:int64):int64;
var e,v,val:longint;
flow,t:int64;
begin
if u=src then exit(aug);
e:=head[u]; val:=s-; flow:=;
while e<> do
begin
v:=vet[e];
if len[e]> then
begin
if dis[u]=dis[v]+ then
begin
t:=dfs(v,min(len[e],aug-flow));
len[e]:=len[e]-t;
len[fan[e]]:=len[fan[e]]+t;
flow:=flow+t;
if dis[source]>=s then exit(flow);
if aug=flow then break;
end;
val:=min(val,dis[v]);
end;
e:=next[e];
end;
if flow= then
begin
dec(gap[dis[u]]);
if gap[dis[u]]= then dis[source]:=s;
dis[u]:=val+;
inc(gap[dis[u]]);
end;
exit(flow);
end; function maxflow:int64;
var ans:int64;
begin
fillchar(gap,sizeof(gap),);
fillchar(dis,sizeof(dis),);
gap[]:=s; ans:=;
while dis[source]<s do ans:=ans+dfs(source,oo);
exit(ans);
end; function max(x,y:longint):longint;
begin
if x>y then exit(x);
exit(y);
end; function isok(k:int64):boolean;
begin
if k*k1-s1=maxflow then exit(true);
exit(false);
end; procedure build(k:int64);
var i,j,l,x,y:longint;
begin
fillchar(head,sizeof(head),); tot:=;
for i:= to n do
for j:= to m do
if (i+j) mod = then
begin
for l:= to do
begin
x:=i+dx[l]; y:=j+dy[l];
if (x>)and(x<=n)and(y>)and(y<=m) then add(num[i,j],num[x,y],oo);
end;
end; for i:= to n do
for j:= to m do
if (i+j) mod = then add(source,num[i,j],k-a[i,j])
else add(num[i,j],src,k-a[i,j]);
end; begin
assign(input,'bzoj2756.in'); reset(input);
assign(output,'bzoj2756.out'); rewrite(output);
read(cas);
for i:= to do
if i mod = then fan[i]:=i+
else fan[i]:=i-;
// for i:= to do
// for j:= to do num[i,j]:=(i-)*+j;
for v:= to cas do
begin
read(n,m);
s1:=; s2:=; k1:=; k2:=; mx:=-maxlongint;
for i:= to n do
for j:= to m do
begin
read(a[i,j]);
if (i+j) mod = then begin inc(k1); s1:=s1+a[i,j]; end
else begin inc(k2); s2:=s2+a[i,j]; end;
mx:=max(mx,a[i,j]); num[i,j]:=(i-)*m+j;
end;
s:=n*m+; source:=n*m+; src:=n*m+;
if k1<>k2 then
begin
k:=(s1-s2) div (k1-k2);
if k<mx then writeln(-)
else
begin
build(k);
if isok(k) then writeln(k*k1-s1)
else writeln(-);
end;
end
else
begin
if s1<>s2 then writeln(-)
else
begin
l:=mx; r:=oo; last:=oo;
while l<=r do
begin
mid:=(l+r)>>;
build(mid);
if isok(mid) then begin last:=mid; r:=mid-; end
else l:=mid+;
end;
writeln(last*k1-s1);
end;
end; end; close(input);
close(output);
end.

【BZOJ2756】奇怪的游戏(二分,最小割)的更多相关文章

  1. BZOJ-2756 奇怪的游戏 黑白染色+最大流+当前弧优化+二分判断+分类讨论

    这个题的数据,太卡了,TLE了两晚上,各种调试优化,各种蛋疼. 2756: [SCOI2012]奇怪的游戏 Time Limit: 40 Sec Memory Limit: 128 MB Submit ...

  2. 最大密集子图(01分数规划+二分+最小割)POJ3155

    题意:给出一副连通图,求出一个子图令g=sigma(E)/sigma(V); h[g]=sigma(E)-g*sigma(V):设G是最优值 则当h[g]>0:g<G h[g]<0, ...

  3. P5038 [SCOI2012]奇怪的游戏 二分+网络流

    $ \color{#0066ff}{ 题目描述 }$ Blinker最近喜欢上一个奇怪的游戏. 这个游戏在一个 \(N \times M\) 的棋盘上玩,每个格子有一个数.每次\(Blinker\)会 ...

  4. 【bzoj2756 奇怪的游戏】

    Time Limit: 40 Sec  Memory Limit: 128 MBSubmit: 4403  Solved: 1226[Submit][Status][Discuss] Descript ...

  5. bzoj 2756 [SCOI2012]奇怪的游戏 二分+网络流

    2756:[SCOI2012]奇怪的游戏 Time Limit: 40 Sec  Memory Limit: 128 MBSubmit: 4926  Solved: 1362[Submit][Stat ...

  6. 【二分 最小割】cf808F. Card Game

    Digital collectible card games have become very popular recently. So Vova decided to try one of thes ...

  7. poj 3155 二分+最小割求实型最小割(最大密集子图)

    /* 最大密集子图子图裸题 解法:设源点s和汇点t 根据胡波涛的<最小割模型在信息学中的应用> s-每个点,权值为原边权和m, 每个点-t,权值为m+2*g-degree[i], 原来的边 ...

  8. BZOJ.2756.[SCOI2012]奇怪的游戏(二分 黑白染色 最大流ISAP)

    题目链接 \(Description\) \(Solution\) 这种题当然要黑白染色.. 两种颜色的格子数可能相同,也可能差1.记\(n1/n2\)为黑/白格子数,\(s1/s2\)为黑/白格子权 ...

  9. 洛谷$P5038\ [SCOI2012]$奇怪的游戏 二分+网络流

    正解:二分+网络流 解题报告: 传送门$QwQ$ 这种什么,"同时增加",长得就挺网络流的$QwQ$?然后看到问至少加多少次,于是考虑加个二分呗?于是就大体确定了做题方向,用的网络 ...

  10. [bzoj2756]奇怪的游戏

    对棋盘黑白染色后,若n和m都是奇数(即白色和黑色点数不同),可以直接算得答案(根据白-黑不变):若n和m不都是奇数,二分答案(二分的上限要大一点,开$2^50$),最后都要用用网络流来判定.考虑判定, ...

随机推荐

  1. hdu2874 LCA

    题意:现在有 n 个点与 m 条边的无向无环图,但是图不一定完全连通,边有各自的边权,给出多组询问,查询两点之间的路径权值和,或者输出两点不连通. 一开始有最短路的想法,但是由于询问有 1e6 组,做 ...

  2. C#多线程网摘 1

    1.每个窗体都有自己的都在不同的线程上运行,如果需要在窗体之间交互,就需要在线程之间交互. 2.当线程Sleep时,系统就退出执行队列一段时间,当睡眠结束时,系统会产生一个时钟中断,从而使线程回到执行 ...

  3. (转) 一张图解AlphaGo原理及弱点

    一张图解AlphaGo原理及弱点 2016-03-23 郑宇,张钧波 CKDD 作者简介: 郑宇,博士, Editor-in-Chief of ACM Transactions on Intellig ...

  4. Winform中调用js函数

    var wb = new WebBrowser(); wb.AllowNavigation = true; wb.ScriptErrorsSuppressed = false; wb.Navigate ...

  5. maven的SNAPSHOT版本和正式版本不同

    转载文章: http://www.huangbowen.net/blog/2016/01/29/understand-official-version-and-snapshot-version-in- ...

  6. NET-SNMP配置

    配置/etc/snmp/snmpd.conf such as below : ============================================== com2sec notCon ...

  7. MVC缓存OutputCacheAttribute 类提高网站效率(转)

    原文转自:http://www.cnblogs.com/iamlilinfeng/p/4419362.html 命名空间:  System.Web.Mvc 程序集:  System.Web.Mvc(在 ...

  8. ROC曲线

    1.混淆矩阵(confusion matrix)     针对预测值和真实值之间的关系,我们可以将样本分为四个部分,分别是:     真正例(True Positive,TP):预测值和真实值都为1 ...

  9. [git]git开发流程

    git开发正确做法: 本地要有一个分支A和远端的分支保持对应 然后本地新开分支B开发,提交记录 如果需要将代码推送到远端的话,就切换回A,首先在A分支上pull同步远端的代码(pull还是fetch+ ...

  10. Random:产生随机数的类

    /* * Random:产生随机数的类 * * 构造方法: * public Random():没有给种子,用的是默认种子,是当前时间的毫秒值 * public Random(long seed):给 ...