[Luogu P4003] [BZOJ 5120] [THU集训] 无限之环

最新推荐文章于 2019-08-05 20:21:24 发布
最新推荐文章于 2019-08-05 20:21:24 发布
阅读量222 收藏
点赞数
分类专栏： 网络流 费用流 文章标签： 网络流 费用流
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/LPA20020220/article/details/84670965
版权
洛谷传送门
BZOJ传送门

题目描述

曾经有一款流行的游戏，叫做 Infinity Loop，先来简单的介绍一下这个游戏：

游戏在一个 n ∗ m n ∗ m n∗m 的网格状棋盘上进行，其中有些小方格中会有水管，水管可能在

格某些方向的边界的中点有接口，所有水管的粗细都相同，所以如果两个相邻方格的

共边界的中点都有接头，那么可以看作这两个接头互相连接。水管有以下 15 15 15 种形状：

img

游戏开始时，棋盘中水管可能存在漏水的地方。形式化地：如果存在某个接头，没有和其它接头相连接，那么它就是一个漏水的地方。

玩家可以进行一种操作：选定一个含有非直线型水管的方格，将其中的水管绕方格中心顺时针或逆时针旋转 90 90 90 度。

直线型水管是指左图里中间一行的两种水管。现给出一个初始局面，请问最少进行多少次操作可以使棋盘上不存在漏水的地方。

输入输出格式

输入格式：

从文件 infinityloop.in 中读入数据。

第一行两个正整数 n n n， m m m，代表网格的大小。

接下来 n n n 行每行 m m m 个数，每个数是 [ 0 , 15 ] [0,15] [0,15] 中的一个，你可以将其看作一个 4 4 4 位的

二进制数，从低到高每一位分别代表初始局面中这个格子上、右、下、左方向上是否有水管接头。

特别地，如果这个数是 0 0 0 ，则意味着这个位置没有水管。

比如 $3(0011(2)) $代表上和右有接头，也就是一个 L L L 型，而 12 ( 1100 ( 2 ) ) 12(1100(2)) 12(1100(2)) 代表下和左

有接头，也就是将 L L L 型旋转 180 180 180 度。

输出格式：

输出到文件 infinityloop.out 中。

输出共一行，表示最少操作次数。如果无法达成目标，输出 − 1 -1 −1.

输入输出样例

输入样例#1：
2 3
3 14 12
3 11 12
输出样例#1：
2
输入样例#2：
3 2
1 8
5 10
2 4
输出样例#2：
-1
输入样例#3：
3 3
9 11 3
13 15 7
12 14 6
输出样例#3：
16

说明

【样例 1 解释】

样例 1 棋盘如下

旋转方法很显然，先将左上角虚线方格内的水管顺时针转 90 度

img

然后右下角虚线方格内的水管逆时针旋转 90 度，这样就使得水管封闭了

【样例 2 解释】

样例 2 为题目描述中的第一张图片，无法达成目标。

【样例 3 解释】

样例 3 为题目描述中的第二张图片，将除了中心方格以外的每个方格内的水管都转 180 度即可。

img

解题分析

费用流神题啊（听说轮廓线 D P DP DP吊打费用流orz）QAQ…

考虑到可能拼好之后会有很多连通的环，不利于确定源汇， 我们不妨黑白染色， 强制流向相邻的点。 很显然， 我们需要拆成5个点， 分别表示中点， 向上方向， 向下方向， 向左方向， 向右方向。 现在我们会做没有旋转操作的了， 我们可以根据要求从 S S S向中点或从中点向 T T T连流量为水管支线数量的边， 再从中点依次向对应方向或从对应方向向中点连流量为 1 1 1， 费用为 0 0 0的边。 最后我们看是否流量达到所有水管支线数量即可。

考虑如何将旋转的情况考虑进去。先来看正 L L L形的情况， 逆时针旋转90°其实等同于将向右的水管接到向左的水管， 因此我们直接从向右方向的点连向向左方向的点连一条流量为 1 1 1， 花费为 1 1 1的边， 从向上方向向向下方向的点连一条流量为 1 1 1， 花费为 1 1 1的边。 丄形的同理。

然后我们发现， 题目中规定直线型的水管不能旋转的原因就在于其不能通过这样增加花费的方式得到旋转后的结果。

代码如下：

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cctype>
#include <algorithm>
#include <cmath>
#include <queue>
#define R register
#define IN inline
#define W while
#define gc getchar()
#define INF 1e8
#define get(i, j) (((i - 1) * m + j) * 5)
#define DEBUG(i, j, k, l) printf("%d %d %d %d\n", i, j, k, l)
#define MX 10050
#define m(i, j) (get(i, j) + 1)
#define u(i, j) (get(i, j) + 2)
#define d(i, j) (get(i, j) + 3)
#define l(i, j) (get(i, j) + 4)
#define r(i, j) (get(i, j) + 5)
template <class C> IN C max(C a, C b) {return a > b ? a : b;}
template <class C> IN C min(C a, C b) {return a < b ? a : b;}
template <class C>
IN void in(C &x)
{
	x = 0; R char c = gc;
	for (; !isdigit(c); c = gc);
	for (;  isdigit(c); c = gc)
	x = (x << 1) + (x << 3) + c - 48;
}
int n, m, sum, S, T, cnt, flow, ans;
int head[MX], pre[MX], dis[MX], del[MX];
int mp[2050][2050];
bool inq[MX];
struct Edge {int to, fl, len, nex;} edge[200500];
IN void add(R int from, R int to, R int fl, R int cst, R bool inv)
{//如果中点连向T， 就把方向反过来连
	if (inv) std::swap(from, to);
	edge[++cnt] = {to, fl, cst, head[from]}, head[from] = cnt;
	edge[++cnt] = {from, 0, -cst, head[to]}, head[to] = cnt;
}
namespace MCMF
{
	std::queue <int> q;
	IN bool SPFA()
	{
		std::memset(dis, 63, sizeof(dis));
		dis[S] = 0; del[S] = INF; q.push(S); R int now;
		W (!q.empty())
		{
			now = q.front(); q.pop();
			for (R int i = head[now]; ~i; i = edge[i].nex)
			{
				if (dis[edge[i].to] > dis[now] + edge[i].len && edge[i].fl)
				{
					dis[edge[i].to] = dis[now] + edge[i].len;
					del[edge[i].to] = min(del[now], edge[i].fl);
					pre[edge[i].to] = i;
					if (!inq[edge[i].to]) inq[edge[i].to] = true, q.push(edge[i].to);
				}
			}
			inq[now] = false;
		}
		return dis[T] ^ dis[0];
	}
	IN void update()
	{
		R int now = T, pr;
		flow += del[T], ans += del[T] * dis[T];
		W (now ^ S)
		{
			pr = pre[now];
			edge[pr].fl -= del[T], edge[pr ^ 1].fl += del[T];
			now = edge[pr ^ 1].to;
		}
	}
	IN void solve()
	{
		W (SPFA()) update();
		printf("%d", flow * 2 == sum ? ans : -1);
	}
}
/*
0001 up
0010 rig
0100 down
1000 lef
1 mid
2 up
3 down
4 lef
5 rig
*/
int main(void)
{
	int typ, ct;
	in(n), in(m); S = n * m * 5 + 10, T = n * m * 5 + 11;
	std::memset(head, cnt = -1, sizeof(head));
	for (R int i = 1; i <= n; ++i) for (R int j = 1; j <= m; ++j) in(mp[i][j]);
	for (R int i = 1; i <= n; ++i) for (R int j = 1; j <= m; ++j)
	{
		typ = (i + j) & 1; ct = __builtin_popcount(mp[i][j]);
		if (!mp[i][j]) continue;
		if (typ) add(m(i, j), T, ct, 0, 0), sum += ct;
		else
		{
			add(S, m(i, j), ct, 0, 0), sum += ct;
			if (mp[i - 1][j]) add(u(i, j), d(i - 1, j), 1, 0, 0);
			if (mp[i + 1][j]) add(d(i, j), u(i + 1, j), 1, 0, 0);
			if (mp[i][j - 1]) add(l(i, j), r(i, j - 1), 1, 0, 0);
			if (mp[i][j + 1]) add(r(i, j), l(i, j + 1), 1, 0, 0);
		}
		switch(mp[i][j])
		{
			case 1://u
			{
				add(m(i, j), u(i, j), 1, 0, typ);
				add(u(i, j), l(i, j), 1, 1, typ);
				add(u(i, j), r(i, j), 1, 1, typ);
				add(u(i, j), d(i, j), 1, 2, typ);
				break;
			}
			case 2://r
			{
				add(m(i, j), r(i, j), 1, 0, typ);
				add(r(i, j), u(i, j), 1, 1, typ);
				add(r(i, j), d(i, j), 1, 1, typ);
				add(r(i, j), l(i, j), 1, 2, typ);
				break;
			}
			case 3://ur
			{
				add(m(i, j), r(i, j), 1, 0, typ);
				add(m(i, j), u(i, j), 1, 0, typ);
				add(u(i, j), d(i, j), 1, 1, typ);
				add(r(i, j), l(i, j), 1, 1, typ);
				break;
			}
			case 4://d
			{
				add(m(i, j), d(i, j), 1, 0, typ);
				add(d(i, j), l(i, j), 1, 1, typ);
				add(d(i, j), r(i, j), 1, 1, typ);
				add(d(i, j), u(i, j), 1, 2, typ);
				break;
			}
			case 5://ud, cannot rotate
			{
				add(m(i, j), u(i, j), 1, 0, typ);
				add(m(i, j), d(i, j), 1, 0, typ);
				break;
			}
			case 6://dr
			{
				add(m(i, j), r(i, j), 1, 0, typ);
				add(m(i, j), d(i, j), 1, 0, typ);
				add(r(i, j), l(i, j), 1, 1, typ);
				add(d(i, j), u(i, j), 1, 1, typ);
				break;
			}
			case 7://urd
			{
				add(m(i, j), u(i, j), 1, 0, typ);
				add(m(i, j), r(i, j), 1, 0, typ);
				add(m(i, j), d(i, j), 1, 0, typ);
				add(u(i, j), l(i, j), 1, 1, typ);
				add(d(i, j), l(i, j), 1, 1, typ);
				add(r(i, j), l(i, j), 1, 2, typ);
				break;
			}
			case 8://l
			{
				add(m(i, j), l(i, j), 1, 0, typ);
				add(l(i, j), u(i, j), 1, 1, typ);
				add(l(i, j), d(i, j), 1, 1, typ);
				add(l(i, j), r(i, j), 1, 2, typ);
				break;
			}
			case 9://ul
			{
				add(m(i, j), u(i, j), 1, 0, typ);
				add(m(i, j), l(i, j), 1, 0, typ);
				add(u(i, j), d(i, j), 1, 1, typ);
				add(l(i, j), r(i, j), 1, 1, typ);
				break;
			}
			case 10://lr, cannot rotate
			{
				add(m(i, j), l(i, j), 1, 0, typ);
				add(m(i, j), r(i, j), 1, 0, typ);
				break;
			}
			case 11://lru
			{
				add(m(i, j), l(i, j), 1, 0, typ);
				add(m(i, j), u(i, j), 1, 0, typ);
				add(m(i, j), r(i, j), 1, 0, typ);
				add(l(i, j), d(i, j), 1, 1, typ);
				add(r(i, j), d(i, j), 1, 1, typ);
				add(u(i, j), d(i, j), 1, 2, typ);
				break;
			}
			case 12://ld
			{
				add(m(i, j), l(i, j), 1, 0, typ);
				add(m(i, j), d(i, j), 1, 0, typ);
				add(l(i, j), r(i, j), 1, 1, typ);
				add(d(i, j), u(i, j), 1, 1, typ);
				break;
			}
			case 13://udl
			{
				add(m(i, j), l(i, j), 1, 0, typ);
				add(m(i, j), u(i, j), 1, 0, typ);
				add(m(i, j), d(i, j), 1, 0, typ);
				add(u(i, j), r(i, j), 1, 1, typ);
				add(d(i, j), r(i, j), 1, 1, typ);
				add(l(i, j), r(i, j), 1, 2, typ);
				break;
			}
			case 14://lrd
			{
				add(m(i, j), l(i, j), 1, 0, typ);
				add(m(i, j), r(i, j), 1, 0, typ);
				add(m(i, j), d(i, j), 1, 0, typ);
				add(l(i, j), u(i, j), 1, 1, typ);
				add(r(i, j), u(i, j), 1, 1, typ);
				add(d(i, j), u(i, j), 1, 2, typ);
				break;
			}
			case 15:
			{
				add(m(i, j), l(i, j), 1, 0, typ);
				add(m(i, j), r(i, j), 1, 0, typ);
				add(m(i, j), u(i, j), 1, 0, typ);
				add(m(i, j), d(i, j), 1, 0, typ);
				break;
			}
		}
	}
	MCMF::solve();
}
LPA20020220
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
【清华集训】2017 无限之环
晚来天欲雪
08-11 370
无限之环 无限之环 题目 链接 描述 输入格式 输出格式 子任务 解题法 评价：神奇的网络流神题 Round 1: 黑白染色 题目 链接 洛谷 描述 曾经有一款流行的游戏，叫做 Infinity Loop，先来简单的介绍一下这个游戏： 游戏在一个 n ∗ m 的网格状棋盘上进行，其中有些小方格中会有水管，水管可能在 格某些方向的边界的中点有接口，所有...
BZOJ5120无限之环-费用流
Fizzmy的博客
01-05 355
传送门 题意：看原题吧，想不出比原题更好的描述了- - Solution： 一开始思考过用网络流，但是想不出如何建图，最后还是去看了题解QwQ，建图思路很妙啊，我们先把每个点拆成四个小点，分别对应上，下，左，右，然后对应每种水管在点内分别建图（细节大家可以结合代码思考一下），由于这是一个二分图（拆点之前），所以说我们对每个点进行黑白染色，这样便可以确定每个点是连向T还是从S连过来，点
[LOJ2321][清华集训2017]无限之环-最小费用最大流
「嫌われ者のフィロソフィ」
01-06 1353
无限之环题目描述曾经有一款流行的游戏，叫做 Infinity Loop，先来简单的介绍一下这个游戏:游戏在一个 n×mn \times m 的网格状棋盘上进行，其中有些小方格中会有水管，水管可能在方格某些方向的边界的中点有接口，所有水管的粗细都相同，所以如果两个相邻方格的公共边界的中点都有接头，那么可以看作这两个接头互相连接。水管有以下 1515 种形状： 游戏开始时，棋盘中水管可能存在漏水的地方
LOJ2321「清华集训 2017」无限之环
ShadyPi的IT笔记
06-30 495
原题链接：https://www.lydsy.com/JudgeOnline/problem.php?id=5120 无限之环 Description 曾经有一款流行的游戏，叫做InfinityLoop，先来简单的介绍一下这个游戏: 游戏在一个n×m的网格状棋盘上进行，其中有些小方格中会有水管，水管可能在方格某些方向的边界的中点有接口，所有水管的粗细都相同，所以如果两个相邻方格的公...
【清华集训2017】【LOJ2331】【洛谷P4003】无限之环（费用流）
zxyoi_dreamer的博客（不定期诈尸）
04-04 293
LOJ传送门 洛谷传送门 解析： 考虑将所有点交替黑白染色，源点向黑点连边，白点向汇点连边。每个点在四周建立四个点，表示这个位置是否被连接，然后所有挨着的点连边，比如iii上面的点的下标点和iii的上标点连边，iii左边点的右标点和iii的左标点连边。 现在我们需要考虑怎么连接iii和iii四周的点。 显然初始状态接触的点直接连接费用为000的边。 现在对于可以旋转的点分情况考虑： 1.单插头类...
Luogu P2278 [HNOI2003]操作系统
01-20
题面 原来是道大水题，但是它的题面有点意思，于是我就手残把它加进了解题计划中。 题面描述 对于操作系统，我们只拥有一个CPU，只能同时处理一个任务。现在有很多任务需要操作系统解决，它们将按照产生时间输入。...
luogu P4568代码
01-27
洛谷P4568题解 其实就是个简单的NOIP..
vscode-luogu:解决VSCode中落谷问题
05-30
vscode-luogu轻松的让你在VSCode 上提交洛谷代码使用说明查看题目按下F1或Ctrl+Shift+P ，查找Luogu: 查看题目命令，点击后输入题号完成搜索点击界面左侧洛谷图片，然后点击搜索图标，输入题号完成搜索完成功能查看...
luogu-cli：罗古在线裁判的cli工具
03-04
如果luogu-cli在运行时产生了类似下面的错误： /usr/bin/env: “node \r ”: 没有那个文件或目录 请尝试下面的操作： 罗古里 洛谷命令行工具 支持以下功能： 译文查看 译文提交 译文搜索 奔奔发送 奔奔查看 云剪切...
Luogu-Painter:洛谷「冬日绘版」自动 Paint 姬
05-29
Luogu-Painter Usage 在洛谷「冬日绘版」中自动 Painting！ Configuration Clone 项目到本地。 npm install。 使用 scripts/loadPic.py 生成图片的 json 格式，将其放到项目 pictures 目录下。 需要安装 PIL 库：pip...
洛谷SP2885 WORDRING - Word Rings（0/1分数规划）（SPFA）
A_Bright_CH的博客
10-26 433
题目 洛谷SP2885 WORDRING - Word Rings 题解 0/1分数规划+spfa判正环 先套0/1分数规划的公式： 再一变 接下来就是一个判正环的问题了。如果有正环则上式&gt;=0，即L有更大取值；否则L无法达到这么优。 提一下建边，真的要做一个什么KMP，AC/后缀自动机，还是马拉车？都是乱来！这样想，它说的是后两个字符，所以不妨直接用两个字符表示点的编号...
洛谷P4003 无限之环（infinityloop）（网络流，费用流）
a83229442的博客
01-12 435
洛谷题目传送门 题目 题目描述 曾经有一款流行的游戏，叫做 Infinity Loop，先来简单的介绍一下这个游戏： 游戏在一个 n ∗ m 的网格状棋盘上进行，其中有些小方格中会有水管，水管可能在格子某些方向的边界的中点有接口，所有水管的粗细都相同，所以如果两个相邻方格的共边界的中点都有接头，那么可以看作这两个接头互相连接。水管有以下 15 种形状： 游戏开始时，棋盘中水管可能存在漏...
【BZOJ】5120: [2017国家集训队测试]无限之环-最小费用流最大流
远行客
12-23 203
巧妙的费用流
BZOJ5120 [2017国家集训队测试]无限之环
yzyyylx的博客
08-30 277
题面 题意 给出一个网格图，每个格子内可能有一个水管，水管有15种形状，也就是说它向每个方向都有可能有分支，请问至少旋转（顺时针或逆时针90度）几个非直线型的水管，使没有水管漏水，即每个水管的每个分支都与其他的某个分支相连。 做法 算法是费用流，难点在于建图： 因为要求每个分支都与其他某个分支相连，因此就要使它满流，不难想到根据每个点横纵坐标和的奇偶性将它与源点或汇点相连，然后将...
【洛谷P4003】无限之环【费用流】
stoorz
08-05 225
题目大意： 题目链接：https://www.luogu.org/problem/P4003 曾经有一款流行的游戏，叫做 Infinity Loop，先来简单的介绍一下这个游戏： 游戏在一个n×mn\times mn×m的网格状棋盘上进行，其中有些小方格中会有水管，水管可能在 格某些方向的边界的中点有接口，所有水管的粗细都相同，所以如果两个相邻方格的 共边界的中点都有接头，那么可以看作这两个接头互...
HDU-IMNU集训之不及AC送我情~：Climbing Worm-新学到一个函数
That's all
01-27 384
题目大致意思为： 一只小虫要从井底爬出 井深为n，上升一秒走u，休息一秒下降d. 可以这样考虑： 每两秒走（u-d)； 无论走了多少个（u-d),而最后一秒只走了u. 可以说最后一秒可以一步到达。 #include #include #include #include #include//ceil()向“大”去整 #include//strstr() using namespa
[BZOJ1596][Usaco2008 Jan]电话网络（贪心||树形dp）
zyf2000
04-26 1260
我们目送着那消失的飞空机云 因为太过於炫丽而想逃开　不知何时变的如此的软弱
[Luogu P4043] [BZOJ 3876] [AHOI2014 JSOI2014]支线剧情
LPA的博客
12-01 472
洛谷传送门 BZOJ传送门 题目背景 宅男JYY非常喜欢玩RPG游戏，比如仙剑，轩辕剑等等。不过JYY喜欢的并不是战斗场景，而是类似电视剧一般的充满恩怨情仇的剧情。这些游戏往往都有很多的支线剧情，现在JYY想花费最少的时间看完所有的支线剧情。 题目描述 JYY现在所玩的RPG游戏中，一共有NNN个剧情点，由111到NNN编号，第iii个剧情点可以根据JYY的不同的选择，而经过不同的支线剧情，前往K...
Luogu P1088 火星人
最新发布
06-07
解题思路：首先将题目转化为背包问题，将火星人分为两组，使得两组的能力值之差最小，可以看作是在一组火星人中选择若干个人，使得他们的能力值之和尽可能接近所有火星人的能力值之和的一半。 然后，使用动态规划...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值