UVA 10888 - Warehouse(二分图完美匹配)

UVA 10888 - Warehouse

option=com_onlinejudge&Itemid=8&page=show_problem&category=562&problem=1829&mosmsg=Submission+received+with+ID+14222079" target="_blank" style="">题目链接

题意:就是推箱子游戏,问最少要几步

思路:每一个箱子和目标位置建边。权值为负权值,然后进行二分图完美匹配就可以,注意不能到达的位置权值应该置为最小

代码:

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
using namespace std;

const int MAXNODE = 45;

typedef int Type;
const Type INF = 0x3f3f3f3f;

struct KM {
	int n;
	Type g[MAXNODE][MAXNODE];
	Type Lx[MAXNODE], Ly[MAXNODE], slack[MAXNODE];
	int left[MAXNODE];
	bool S[MAXNODE], T[MAXNODE];

	void init(int n) {
		this->n = n;
	}

	void add_Edge(int u, int v, Type val) {
		g[u][v] = val;
	}

	bool dfs(int i) {
		S[i] = true;
		for (int j = 0; j < n; j++) {
			if (T[j]) continue;
			Type tmp = Lx[i] + Ly[j] - g[i][j];
			if (!tmp) {
				T[j] = true;
				if (left[j] == -1 || dfs(left[j])) {
					left[j] = i;
					return true;
				}
			} else slack[j] = min(slack[j], tmp);
		}
		return false;
	}

	void update() {
		Type a = INF;
		for (int i = 0; i < n; i++)
			if (!T[i]) a = min(a, slack[i]);
		for (int i = 0; i < n; i++) {
			if (S[i]) Lx[i] -= a;
			if (T[i]) Ly[i] += a;
		}
	}

	int km() {
		for (int i = 0; i < n; i++) {
			left[i] = -1;
			Lx[i] = -INF; Ly[i] = 0;
			for (int j = 0; j < n; j++)
				Lx[i] = max(Lx[i], g[i][j]);
		}
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < n; j++) slack[j] = INF;
			while (1) {
				for (int j = 0; j < n; j++) S[j] = T[j] = false;
				if (dfs(i)) break;
				else update();
			}
		}
		int ans = 0;
		for (int i = 0; i < n; i++)
			ans -= g[left[i]][i];
		return ans;
	}
} gao;

#define MP(a,b) make_pair(a,b)
const int d[4][2] = {{0, 1}, {0, -1}, {1, 0}, {-1, 0}};

typedef pair<int, int> pii;

const int N = 45;

int t, n, m;
char g[N][N];
int idb[N][N], idx[N][N], bn, xn, vis[N][N];
pii tob[N], tox[N];

int dis(pii a, pii b) {
	return abs(a.first - b.first) + abs(a.second - b.second);
}

void bfs(int xs, int ys) {
	queue<pii> Q;
	Q.push(MP(xs, ys));
	memset(vis, INF, sizeof(vis));
	vis[xs][ys] = 0;
	while (!Q.empty()) {
		pii u = Q.front();
		if (g[u.first][u.second] == 'X')
			gao.add_Edge(idb[xs][ys], idx[u.first][u.second], -vis[u.first][u.second]);
		Q.pop();
		for (int i = 0; i < 4; i++) {
			int x = u.first + d[i][0];
			int y = u.second + d[i][1];
			if (x < 0 || x >= n || y < 0 || y >= m || g[x][y] == '#') continue;
			if (vis[x][y] <= vis[u.first][u.second] + 1) continue;
			vis[x][y] = vis[u.first][u.second] + 1;
			Q.push(MP(x, y));
		}
	}
}

int main() {
	scanf("%d", &t);
	while (t--) {
		scanf("%d%d", &n, &m);
		bn = xn = 0;
		memset(gao.g, -INF, sizeof(gao.g));
		for (int i = 0; i < n; i++) {
			scanf("%s", g[i]);
			for (int j = 0; j < m; j++) {
				if (g[i][j] == 'B') {
					tob[bn] = MP(i, j);
					idb[i][j] = bn++;
				}
				if (g[i][j] == 'X') {
					tox[xn] = MP(i, j);
					idx[i][j] = xn++;
				}
			}
		}
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < m; j++) {
				if (g[i][j] == 'B') {
					bfs(i, j);
				}
			}
		}
		gao.init(bn);
		printf("%d\n", gao.km());
	}
	return 0;
}
时间: 2024-10-10 06:14:54

UVA 10888 - Warehouse(二分图完美匹配)的相关文章

UVA 1411 - Ants(二分图完美匹配)

UVA 1411 - Ants 题目链接 题意:给定一些黑点白点,要求一个黑点连接一个白点,而且全部线段都不相交 思路:二分图完美匹配,权值存负的欧几里得距离,这种话,相交肯定比不相交权值小,所以做一次完美匹配就能够了 代码: #include <cstdio> #include <cstring> #include <cmath> #include <algorithm> using namespace std; const int MAXNODE = 1

UVA 1045 - The Great Wall Game(二分图完美匹配)

UVA 1045 - The Great Wall Game 题目链接 题意:给定一个n*n的棋盘,有n个棋子在上面,现在要移动棋子,每一步代价是1,现在要把棋子移动到一行,一列,或者在主副对角线上,问最小代价 思路:二分图完美匹配,枚举每种情况,建边,边权为曼哈顿距离,然后km算法做完美匹配算出值即可,由于要求最小值所以边权传负数,这样做出来的值的负就是答案 代码: #include <cstdio> #include <cstring> #include <cmath&g

UVA 11383 - Golden Tiger Claw(二分图完美匹配扩展)

UVA 11383 - Golden Tiger Claw 题目链接 题意:给定每列和每行的和,给定一个矩阵,要求每个格子(x, y)的值小于row(i) + col(j),求一种方案,并且所有行列之和的和最小 思路:A二分图完美匹配的扩展,行列建二分图,权值为矩阵相应位置的值,做一次KM算法后,所有顶标之和就是最小的 代码: #include <cstdio> #include <cstring> #include <cmath> #include <algor

POJ 1904 King&#39;s Quest强连通分量+二分图完美匹配

题目描述: Description Once upon a time there lived a king and he had N sons. And there were N beautiful girls in the kingdom and the king knew about each of his sons which of those girls he did like. The sons of the king were young and light-headed, so i

POJ - 1486 Sorting Slides 二分图 完美匹配

题目大意:有n张透明的矩形纸张散乱在桌面上,每张纸张上面都有一个数字. 现在给出每个矩形的左下角和右上角坐标和每个数字所在的位置,问能否找出每个矩形唯一的对应数字 解题思路:分析该题可得到,二分图匹配的结果肯定是完美匹配,匹配的结果肯定为n. 接着就要判断每个点是否唯一匹配了 判断能否唯一匹配,就要在完美匹配的情况下,删除该点的那条匹配,如果还能再找出一个完美匹配,那么就表示该点表示不唯一 #include<cstdio> #include<vector> #include<

KM算法(二分图完美匹配)

前言 其实这是一个很玄学的算法,我也不是很懂.但是慢慢地写,就慢慢地懂了. 例子 我们先看到这样一个例子:现有n个男生和n个女生,每一个女生对于一个男生都有一个期望值.现在想知道期望值最大是多少. 算法 那么显然,我们需要一种算法来完成这个任务,KM腾空出世! 在此之前,我们需要了解一些东西: 前置知识 期望值:女生连的边中边权最大的值,称之为女生的期望值 l(x)+l(y)>=w(i,j) 当l(x)+l(y)=w(i,j)时,是M图的完美匹配! 过程 原则 能换就换(思想很像匈牙利求二分图最

Hall定理 二分图完美匹配

充分性证明就先咕了,因为楼主太弱了,有一部分没看懂 霍尔定理内容 二分图G中的两部分顶点组成的集合分别为X, Y(假设有\(\lvert X \rvert \leq \lvert Y \rvert\)).G中有一组无公共点的边,一端恰好为组成X的点(也就是存在完美匹配)的充分必要条件是:X中的任意k个点至少与Y中的k个点相邻,即对于X中的一个点集W ,令N(W)为W的所有邻居, 霍尔定理即对于任意W,\(\lvert W\rvert \leq \lvert N(W)\rvert\) 证明 1.必

KM算法【带权二分图完美匹配】

先orz litble--KM算法 为什么要用KM算法 因为有的题丧心病狂卡费用流 KM算法相比于费用流来说,具有更高的效率. 算法流程 我们给每一个点设一个期望值[可行顶标] 对于左边的点来说,就是期望能匹配到多大权值的右边的点 对于右边的点来说,就是期望能在左边的点的期望之上还能产生多少贡献 两个点能匹配,当且仅当它们的期望值之和为这条边的权值 一开始初始化所有左点的期望是其出边的最大值,因为最理想情况下当然是每个点都匹配自己能匹配最大的那个 右点期望为0 然后我们逐个匹配,当一个点匹配失败

UVA 1349 - Optimal Bus Route Design(KM完美匹配)

UVA 1349 - Optimal Bus Route Design 题目链接 题意:给定一些有向带权边,求出把这些边构造成一个个环,总权值最小 思路:由于环入度出度为1,所以可以把每个点拆成入度点和出度点,然后建图做一次二分图完美匹配即可,注意这题坑点,有重边 代码: #include <cstdio> #include <cstring> #include <cmath> #include <algorithm> using namespace std