【ybtoj高效进阶 21254】摆放鞋子(二分图匹配)(网络流)

摆放鞋子

题目链接:ybtoj高效进阶 21254

题目大意

给你一个网格,然后每个点可能是两个类型的其中一种,四个方向的一种,然后你可以无限次操作,每次选两个相邻的点,一个顺时针转 90 度,一个逆时针转 90 度。

然后两个相邻的点能匹配要他们是不同类型的,而且它们的方向要满足四个条件的其中一个。
(四个条件其实就是要按着摆鞋子的方法摆,总之就是左脚右脚要放好)

然后问你能有的最大匹配。

思路

补充一下四个条件:
【ybtoj高效进阶 21254】摆放鞋子(二分图匹配)(网络流)

首先我们要发现它这个转是可以传递下去的,然后如果你把右下左上分别标成 \(0,1,2,3\),你会发现整个图所有点方向的值的和对 \(4\) 的取模结果是不变的。

你考虑从这个地方入手,首先不管这个条件直接用类型匹配,这个二分图匹配。
然后如果匹配没有满,那你剩下的那个随便挑方向,一定会可以使得最后的值是满足的。
但是如果匹配满了,你就要特殊判断一下(二分图的匹配结果拿出来用),如果值是满足的就直接用,否则就要放掉一个匹配来使得条件满足。

然后就好啦。

代码

#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#define INF 0x3f3f3f3f3f3f3f3f

using namespace std;

int n, m, a[105][105], st, ans, S, T;
char c[105][105], w[105][105];
int dx[4] = {0, 1, 0, -1}, dy[4] = {1, 0, -1, 0};

bool ck(int x, int y) {
	if (x < 1 || x > n) return 0;
	if (y < 1 || y > m) return 0;
	return 1;
}

struct node {
	int x, to, nxt, op;
}e[600001];
int le[10005], KK, lee[10005];

void add(int x, int y, int z) {
	e[++KK] = (node){z, y, le[x], KK + 1}; le[x] = KK;
	e[++KK] = (node){0, x, le[y], KK - 1}; le[y] = KK;
}

int get_bh(int x, int y) {
	return (x - 1) * m + y;
} 

int dis[10005];
queue <int> q;

bool bfs() {//网络流
	while (!q.empty()) q.pop();
	memset(dis, 0, sizeof(dis));
	for (int i = 1; i <= T; i++) lee[i] = le[i];
	
	dis[S] = 1;
	q.push(S);
	while (!q.empty()) {
		int now = q.front();
		q.pop();
		
		for (int i = le[now]; i; i = e[i].nxt)
			if (!dis[e[i].to] && e[i].x) {
				dis[e[i].to] = dis[now] + 1;
				q.push(e[i].to);
			}
	}
	
	return dis[T];
}

int dfs(int now, int sum) {
	if (now == T) return sum;
	
	int go = 0;
	for (int &i = lee[now]; i; i = e[i].nxt)
		if (e[i].x && dis[e[i].to] == dis[now] + 1) {
			int this_go = dfs(e[i].to, min(sum - go, e[i].x));
			if (this_go) {
				e[i].x -= this_go;
				e[e[i].op].x += this_go;
				go += this_go;
				if (go == sum) return go;
			}
		}
	
	if (go < sum) dis[now] = -1;
	return go;
}

int dinic() {
	int re = 0;
	while (bfs())
		re += dfs(S, INF);
	return re;
}

int main() {
//	freopen("shoe.in", "r", stdin);
//	freopen("shoe.out", "w", stdout);
	
	scanf("%d %d", &n, &m);
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= m; j++) {
			c[i][j] = getchar();
			while (c[i][j] != 'L' && c[i][j] != 'R') c[i][j] = getchar();
		}
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= m; j++) {
			w[i][j] = getchar();
			while (w[i][j] != 'L' && w[i][j] != 'R' && w[i][j] != 'U' && w[i][j] != 'D') w[i][j] = getchar();
			if (w[i][j] == 'R') a[i][j] = 0;
			if (w[i][j] == 'L') a[i][j] = 2;
			if (w[i][j] == 'D') a[i][j] = 1;
			if (w[i][j] == 'U') a[i][j] = 3;
			st = (st + a[i][j]) % 4;
		}
	
	S = n * m + 1; T = S + 1;
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= m; j++)
			if (c[i][j] == 'L') {
				add(S, get_bh(i, j), 1);
				for (int k = 0; k < 4; k++) {
					int xx = i + dx[k], yy = j + dy[k];
					if (ck(xx, yy) && c[xx][yy] == 'R') {
						add(get_bh(i, j), get_bh(xx, yy), 1);
					}
				}
			}
			else add(get_bh(i, j), T, 1);
	
	ans = dinic();
	
	if (ans * 2 < n * m) {//转移没有满,可以剩下的随便动使得满足条件
		printf("%d", ans);
	}
	else {
		int ck = 0;
		for (int i = 1; i <= n; i++)
			for (int j = 1; j <= m; j++)
				if (c[i][j] == 'L') {
					int x = get_bh(i, j);
					for (int k = le[x]; k; k = e[k].nxt)
						if (!e[k].x && k & 1 && e[k].to != S) {
							if (i == (e[k].to - 1) / m + 1)
								ck = (ck + 2) % 4;
						}
				}
		if (ck == st) printf("%d", ans);//当前已经满足了条件
			else printf("%d", ans - 1);//不能就只能解开一对来满足条件
	}
	
	return 0;
}

上一篇:CSP2021自闭记


下一篇:CES 2022:英飞凌推出全新 AURIX 28nm MCU 和 AIROC BLE SoC