BZOJ 2668: [cqoi2012]交换棋子
2668: [cqoi2012]交换棋子
Time Limit: 3 Sec Memory Limit: 128 MB
Submit: 1112 Solved: 409
[Submit][Status][Discuss]
Description
Input
Output
Sample Input
110
000
001
000
110
100
222
222
222
Sample Output
HINT
Source
对于每一个格子,拆成三个点,分别为x,y,z。
其中,x为格子间转移边的入点,z为出点,y为与S和T的连接点。
S向每个原图中黑色,新图中白色的y点连一条容量为1,费用为0的边,流过这条边表示把该棋子换走替代其他棋子。
T向每个原图中白色,新图中黑色的y点连一条容量为1,费用为0的边,流经这条边表示用其他棋子把该棋子换走。
所有原图中黑色,新图中白色的点,x向y连容量为$\frac{limit_{i,j}}{2}$,费用为0的边,表示至多接受这么多替换;y向z连容量为$\frac{limit_{i,j}+1}{2}$,费用为0的边,表示之多提供这么多替换。所有原图中黑色,新图中白色的点,和这个反着连。新图原图中一样的,容量则都是$\frac{limit_{i,j}}{2}$。
按照八联通,所有相邻点对,连接其对应为x,z即可,容量无限,费用为1,表示交换一次的代价为1。
#include <bits/stdc++.h>
#define rep(a,b,c) for(register int a=b;a<=c;++a)
const int mxn = , siz = , inf = 1e9, mv[][] = {{,},{,},{,},{,-}};
char a[mxn][mxn], b[mxn][mxn], c[mxn][mxn];
int n, m, d[mxn][mxn], id[mxn][mxn][], ID, cnt0, cnt1;
int s, t, tot, hd[siz], to[siz], nt[siz], fl[siz], vl[siz], dis[siz], pre[siz], que[siz], inq[siz], cost;
inline bool legal(int x, int y) { return x >= && x <= n && y >= && y <= m; }
inline void add(int u, int v, int f, int w) {
nt[tot] = hd[u]; to[tot] = v; fl[tot] = f; vl[tot] = +w; hd[u] = tot++;
nt[tot] = hd[v]; to[tot] = u; fl[tot] = ; vl[tot] = -w; hd[v] = tot++;
}
inline bool spfa(void) {
register int head = , tail = ;
rep(i, s, t)dis[i] = inf, inq[i] = ;
dis[que[tail++] = s] = , inq[s] = , pre[s] = -;
while (head != tail) {
int u = que[head++], v, i; inq[u] = ;
for (i = hd[u]; ~i; i = nt[i])if (fl[i] && dis[v = to[i]] > dis[u] + vl[i]) {
dis[v] = dis[u] + vl[i], pre[v] = i^; if (!inq[v])inq[que[tail++] = v] = ;
}
}
return dis[t] < inf;
}
inline int maxFlow(void) {
int ret = ;
while (spfa()) {
int flow = inf, i;
for (i = pre[t]; ~i; i = pre[to[i]])if (flow > fl[i^])flow = fl[i^];
for (i = pre[t]; ~i; i = pre[to[i]])fl[i] += flow, fl[i^] -= flow;
ret += flow, cost += flow * dis[t];
}
return ret;
}
signed main(void) {
scanf("%d%d", &n, &m);
s = , t = n * m * + ;
memset(hd, -, sizeof(hd));
rep(i, , n)scanf("%s", a[i] + );
rep(i, , n)scanf("%s", b[i] + );
rep(i, , n)scanf("%s", c[i] + );
rep(i, , n)rep(j, , m)d[i][j] = c[i][j] - '';
rep(i, , n)rep(j, , m)rep(k, , )id[i][j][k] = ++ID;
rep(i, , n)rep(j, , m) {
if (a[i][j] == '' && b[i][j] == '')
add(s, id[i][j][], , ), ++cnt0,
add(id[i][j][], id[i][j][], d[i][j] >> , ),
add(id[i][j][], id[i][j][], (d[i][j] + ) >> , );
if (a[i][j] == '' && b[i][j] == '')
add(id[i][j][], t, , ), ++cnt1,
add(id[i][j][], id[i][j][], (d[i][j] + ) >> , ),
add(id[i][j][], id[i][j][], d[i][j] >> , );
if (a[i][j] == b[i][j])
add(id[i][j][], id[i][j][], d[i][j] >> , ),
add(id[i][j][], id[i][j][], d[i][j] >> , );
}
rep(i, , n)rep(j, , m)rep(k, , )if (legal(i + mv[k][], j + mv[k][]))
add(id[i][j][], id[i + mv[k][]][j + mv[k][]][], inf, ),
add(id[i + mv[k][]][j + mv[k][]][], id[i][j][], inf, );
printf("%d\n", cnt0 == cnt1 && cnt0 == maxFlow() ? cost : -);
}
@Author: YouSiki
BZOJ 2668: [cqoi2012]交换棋子的更多相关文章
- BZOJ 2668 [cqoi2012]交换棋子 | 最小费用最大流
传送门 BZOJ 2668 题解 同时分别限制流入和流出次数,所以把一个点拆成三个:入点in(x).中间点mi(x).出点ou(x). 如果一个格子x在初始状态是黑点,则连(S, mi(x), 1, ...
- BZOJ.2668.[CQOI2012]交换棋子(费用流zkw)
题目链接 首先黑白棋子的交换等价于黑棋子在白格子图上移动,都到达指定位置. 在这假设我们知道这题用网络流做. 那么黑棋到指定位置就是一条路径,考虑怎么用流模拟出这条路径. 我们发现除了路径的起点和终点 ...
- 2668: [cqoi2012]交换棋子
Description 有一个n行m列的黑白棋盘,你每次可以交换两个相邻格子(相邻是指有公共边或公共顶点)中的棋子,最终达到目标状态.要求第i行第j列的格子只能参与mi,j次交换. Input 第一行 ...
- BZOJ2668: [cqoi2012]交换棋子
题解: 可以戳这里:http://www.cnblogs.com/zig-zag/archive/2013/04/21/3033485.html 其实自己yy一下就知道这样建图的正确性了. 感觉太神奇 ...
- [cqoi2012]交换棋子
2668: [cqoi2012]交换棋子 Time Limit: 3 Sec Memory Limit: 128 MBSubmit: 1334 Solved: 518[Submit][Stat ...
- 【BZOJ2668】[cqoi2012]交换棋子 费用流
[BZOJ2668][cqoi2012]交换棋子 Description 有一个n行m列的黑白棋盘,你每次可以交换两个相邻格子(相邻是指有公共边或公共顶点)中的棋子,最终达到目标状态.要求第i行第j列 ...
- 洛谷 P3159(BZOJ 2668)[CQOI2012]交换棋子
有一个\(n\)行\(m\)列的黑白棋盘,你每次可以交换两个相邻格子(相邻是指有公共边或公共顶点)中的棋子,最终达到目标状态.要求第\(i\)行第\(j\)列的格子只能参与\(m[i][j]\)次交换 ...
- BZOJ2668:[CQOI2012]交换棋子——题解
http://www.lydsy.com/JudgeOnline/problem.php?id=2668 https://www.luogu.org/problemnew/show/P3159#sub ...
- P3159 [CQOI2012]交换棋子
思路 相当神奇的费用流拆点模型 最开始我想到把交换黑色棋子看成一个流流动的过程,流从一个节点流向另一个节点就是交换两个节点,然后把一个位置拆成两个点限制流量,然后就有了这样的建图方法 S向所有初始是黑 ...
随机推荐
- windows 平台安装 ffmpeg
一.从https://ffmpeg.zeranoe.com/builds/中下载ffmpeg的static版本: 二.将下载下来的“ffmpeg-4.0.2-win64-static.zip”解压到任 ...
- nodejs 搭建自己的简易缓存cache管理模块
http://www.infoq.com/cn/articles/built-cache-management-module-in-nodejs/ 为什么要搭建自己的缓存管理模块? 这个问题其实也是在 ...
- 关于SQL while 循环嵌套 外部循环数据无法进入内部循环
下面一般是,作为SQL新手第一次写循环嵌套的办法,但是大家会发现一个问题,那就是变量@i总是不能进入第二个循环. declare @i int ,@j int, @k int set @j = 1 - ...
- Java调用XML的方法:DocumentBuilderFactory
(1)首先得到:得到 DOM 解析器的工厂实例 DocumentBuilderFactory domfac=DocumentBuilderFactory.newInstance(); (2)然后从 D ...
- 在 Linux 下执行安装 bin 文件
原文: http://www.linuxidc.com/Linux/2014-03/98541.htm
- Daily Scrum 1 --团队项目所需时间估计以及任务分配
考虑到所有的任务不可能逐一细化分配给成员,我们将需要完成的任务进行了大致的分配.任务所需要的具体实现可以参看<学霸网站NABC> 所需要的总时间一共为44h. 我们会在以后的每日任务中进行 ...
- 什么是REST编程
参考:什么是REST编程:http://www.ruanyifeng.com/blog/2011/09/restful.html 一.REST是Representational State Trans ...
- 07慕课网《进击Node.js基础(一)》HTTP小爬虫
获取HTML页面 var http = require('http') var url='http://www.imooc.com/learn/348' http.get(url,function(r ...
- P4论文粗读笔记(一)
一 文章名称:SNAP: Stateful Network-Wide Abstractions for Packet Processing 数据包处理的带状态网络概念 发表时间:2016 期刊来源:S ...
- python学习笔记05:贪吃蛇游戏代码
贪吃蛇游戏截图: 首先安装pygame,可以使用pip安装pygame: pip install pygame 运行以下代码即可: #!/usr/bin/env python import pygam ...