BZOJ 2668: [cqoi2012]交换棋子
2017-01-10 15:01
218 查看
2668: [cqoi2012]交换棋子
Time Limit: 3 Sec Memory Limit: 128 MBSubmit: 1112 Solved: 409
[Submit][Status][Discuss]
Description
有一个n行m列的黑白棋盘,你每次可以交换两个相邻格子(相邻是指有公共边或公共顶点)中的棋子,最终达到目标状态。要求第i行第j列的格子只能参与mi,j次交换。Input
第一行包含两个整数n,m(1<=n, m<=20)。以下n行为初始状态,每行为一个包含m个字符的01串,其中0表示黑色棋子,1表示白色棋子。以下n行为目标状态,格式同初始状态。以下n行每行为一个包含m个0~9数字的字符串,表示每个格子参与交换的次数上限。Output
输出仅一行,为最小交换总次数。如果无解,输出-1。Sample Input
3 3110
000
001
000
110
100
222
222
222
Sample Output
4HINT
Source
[Submit][Status][Discuss]对于每一个格子,拆成三个点,分别为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 = 25, siz = 500005, inf = 1e9, mv[4][2] = {{0,1},{1,0},{1,1},{1,-1}}; char a[mxn][mxn], b[mxn][mxn], c[mxn][mxn]; int n, m, d[mxn][mxn], id[mxn][mxn][3], 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 >= 1 && x <= n && y >= 1 && 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] = 0; vl[tot] = -w; hd[v] = tot++; } inline bool spfa(void) { register int head = 0, tail = 0; rep(i, s, t)dis[i] = inf, inq[i] = 0; dis[que[tail++] = s] = 0, inq[s] = 1, pre[s] = -1; while (head != tail) { int u = que[head++], v, i; inq[u] = 0; 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^1; if (!inq[v])inq[que[tail++] = v] = 1; } } return dis[t] < inf; } inline int maxFlow(void) { int ret = 0; while (spfa()) { int flow = inf, i; for (i = pre[t]; ~i; i = pre[to[i]])if (flow > fl[i^1])flow = fl[i^1]; for (i = pre[t]; ~i; i = pre[to[i]])fl[i] += flow, fl[i^1] -= flow; ret += flow, cost += flow * dis[t]; } return ret; } signed main(void) { scanf("%d%d", &n, &m); s = 0, t = n * m * 3 + 1; memset(hd, -1, sizeof(hd)); rep(i, 1, n)scanf("%s", a[i] + 1); rep(i, 1, n)scanf("%s", b[i] + 1); rep(i, 1, n)scanf("%s", c[i] + 1); rep(i, 1, n)rep(j, 1, m)d[i][j] = c[i][j] - '0'; rep(i, 1, n)rep(j, 1, m)rep(k, 0, 2)id[i][j][k] = ++ID; rep(i, 1, n)rep(j, 1, m) { if (a[i][j] == '0' && b[i][j] == '1') add(s, id[i][j][1], 1, 0), ++cnt0, add(id[i][j][0], id[i][j][1], d[i][j] >> 1, 0), add(id[i][j][1], id[i][j][2], (d[i][j] + 1) >> 1, 0); if (a[i][j] == '1' && b[i][j] == '0') add(id[i][j][1], t, 1, 0), ++cnt1, add(id[i][j][0], id[i][j][1], (d[i][j] + 1) >> 1, 0), add(id[i][j][1], id[i][j][2], d[i][j] >> 1, 0); if (a[i][j] == b[i][j]) add(id[i][j][0], id[i][j][1], d[i][j] >> 1, 0), add(id[i][j][1], id[i][j][2], d[i][j] >> 1, 0); } rep(i, 1, n)rep(j, 1, m)rep(k, 0, 3)if (legal(i + mv[k][0], j + mv[k][1])) add(id[i][j][2], id[i + mv[k][0]][j + mv[k][1]][0], inf, 1), add(id[i + mv[k][0]][j + mv[k][1]][2], id[i][j][0], inf, 1); printf("%d\n", cnt0 == cnt1 && cnt0 == maxFlow() ? cost : -1); }
@Author: YouSiki
相关文章推荐
- [CQOI2012][bzoj2668] 交换棋子 [费用流]
- 【BZOJ2668】【cqoi2012】交换棋子 费用流
- BZOJ2668 [cqoi2012]交换棋子 【费用流】
- 【bzoj2668】【cqoi2012】【交换棋子】【费用流】
- BZOJ 2668: [cqoi2012]交换棋子(费用流)
- bzoj2668 [cqoi2012]交换棋子
- bzoj 2668: [cqoi2012]交换棋子
- bzoj2668 [cqoi2012]交换棋子
- BZOJ 2668 [cqoi2012]交换棋子 | 最小费用最大流
- bzoj2668 [cqoi2012]交换棋子(费用流)
- [bzoj2668] [cqoi2012]交换棋子
- 【BZOJ】【2668】【CQOI2012】交换棋子
- [BZOJ2668][CQOI2012]交换棋子
- BZOJ 2668 CQOI 2012 交换棋子 费用流
- BZOJ2668: [cqoi2012]交换棋子
- [BZOJ2668][CQOI2012]交换棋子(费用流)
- 【bzoj2668】[cqoi2012]交换棋子 费用流
- 【BZOJ2668】[cqoi2012]交换棋子 费用流
- 【bzoj 2668】: [cqoi2012]交换棋子
- 2668: [cqoi2012]交换棋子