放牧代码和思想
专注自然语言处理、机器学习算法
    This thing called love. Know I would've. Thrown it all away. Wouldn't hesitate.

POJ 2226 Muddy Fields 题解 《挑战程序设计竞赛》

目录

POJ 2226 Muddy Fields

泥地:R*C的地图上*代表泥地,用长度不限宽度为1的木板覆盖泥地,求最小木板数?

3.5借助水流解决问题的网络流 

二分图匹配 

将每块泥地的最左端和最上端连起来构成二分图,其最小顶点覆盖即为所求。因为最小顶点覆盖中的每个顶点构成的边的集合就是全集,又因为每个顶点都代表一条确定的木板,也就是从这些端点引出的木板可以覆盖所有泥地。

#include <iostream>
#include <vector>

using namespace std;

#define MAX_N 50 + 16
////////////////////////最大流开始//////////////////////////////////////
#define MAX_V 2500 * 2
int V;                 // 顶点数
vector<int> G[MAX_V];    // 图的邻接表
int match[MAX_V];      // 所匹配的顶点
bool used[MAX_V];      // DFS中用到的访问标记

// 向图中增加一条连接u和v的边
void add_edge(int u, int v)
{
	G[u].push_back(v);
	G[v].push_back(u);
}

// 通过DFS寻找增广路
bool dfs(int v)
{
	used[v] = true;
	for (vector<int>::iterator it = G[v].begin(); it != G[v].end(); ++it)
	{
		int u = *it, w = match[u];
		if (w < 0 || !used[w] && dfs(w))
		{
			match[v] = u;
			match[u] = v;
			return true;
		}
	}

	return false;
}

// 求解二分图的最大匹配
int bipartite_matching()
{
	int res = 0;
	memset(match, -1, sizeof(match));
	for (int v = 0; v < V; ++v)
	{
		if (match[v] < 0)
		{
			memset(used, 0, sizeof(used));
			if (dfs(v))
			{
				++res;
			}
		}
	}

	return res;
}

void clear()
{
	for (int i = 0; i < V; ++i)
	{
		G[i].clear();
	}
}

///////////////////////////////最大流结束/////////////////////////////////////
char grid[MAX_N][MAX_N];

///////////////////////////SubMain//////////////////////////////////
int main(int argc, char *argv[])
{
#ifndef ONLINE_JUDGE
	freopen("in.txt", "r", stdin);
	freopen("out.txt", "w", stdout);
#endif
	int R, C;
	scanf("%d %d", &R, &C);
	for (int i = 0; i < R; ++i)
	{
		scanf("%s", grid[i]);
	}
	for (int i = 0; i < R; ++i)
	{
		for (int j = 0; j < C; ++j)
		{
			if (grid[i][j] == '*')
			{
				int y = i, x = j;
				while (y > 0 && grid[y - 1][j] == '*')
				{
					--y;
				}
				while (x > 0 && grid[i][x - 1] == '*')
				{
					--x;
				}
				add_edge(y * 50 + j, i * 50 + x + 2500);
			}
		}
	}
	V = MAX_V;
	printf("%d", bipartite_matching());
#ifndef ONLINE_JUDGE
	fclose(stdin);
	fclose(stdout);
	system("out.txt");
#endif
	return 0;
}
///////////////////////////End Sub//////////////////////////////////
13819124 hankcs 2226 Accepted 304K 32MS C++ 1984B 2015-01-26 19:24:09

知识共享许可协议 知识共享署名-非商业性使用-相同方式共享码农场 » POJ 2226 Muddy Fields 题解 《挑战程序设计竞赛》

评论 欢迎留言

  • 昵称 (必填)
  • 邮箱 (必填)
  • 网址

我的作品

HanLP自然语言处理包《自然语言处理入门》