POJ 3041 匈牙利算法 最小覆盖 最大二分匹配

题目在这里 http://poj.org/problem?id=3041

因为最近在准备找工作，所以一直在学习算法，不过我觉得对于一个初学者来说，最难的不是弄懂一个算法，二是如何将题目编程算法的模型，然后编写出程序。个人觉得这个过程应该就是靠多做题，多分析吧，培养自己的灵感来源。

这个题目大概的意思是说， 在一个N*N的矩阵中，有若干个陨石，飞船的大炮可以每次消灭一排或者一列，问如何使得开炮的次数最小，但是却可以消灭掉所有的陨石。

题目分析：

可以将行和列做如下转换：

上面的图中，左图表示行编号，右图表示列编号，如果在(i,j)上有一颗陨石，我们将行i 和 列 j连接。

这样每一条边都表示一个陨石，如果我在第1行开炮，那么和行1相连的所有边就消除了。 所以这样这个问题就转换成了最小覆盖问题。

引入点理论吧，在二分图中：

匹配：

给定一个二分图，在G的一个子图G’中，如果G’的边集中的任意两条边都不依附于同一个顶点，则称G’的边集为G的一个匹配

最大匹配：

在所有的匹配中，边数最多的那个匹配就是二分图的最大匹配了

顶点覆盖：

在顶点集合中，选取一部分顶点，这些顶点能够把所有的边都覆盖了。这些点就是顶点覆盖集

最小顶点覆盖：

在所有的顶点覆盖集中，顶点数最小的那个叫最小顶点集合。

二分图最大匹配的König定理

最大匹配 = 最小顶点覆盖。此处不证明其正确性。

然后我们写出程序如下：

#include <stdio.h>
#include <memory.h>

#define GRID_NUM 505

int Grid[GRID_NUM][GRID_NUM];
int matched[GRID_NUM];
bool visited[GRID_NUM];
int nGrid;

bool dfs(int from)
{
	int i;

	for( i = 0;i < nGrid; i++)
	{
		if(!visited[i] && Grid[from][i])
		{
			visited[i] = true;
			if(matched[i] == -1 || dfs(matched[i]))
			{
				
				matched[i] = from;
				return true;
			}
		}
	}
	return false;
	
}

int main()
{
	int nAstroid, j;
	int f, t;
	while(scanf("%d%d", &nGrid, &nAstroid) != EOF)
	{
		memset(Grid, 0, sizeof(Grid));
		memset(matched, -1, sizeof(matched));
		for(j = 0; j < nAstroid; j++)
		{
			scanf("%d%d", &f, &t);
			Grid[f-1][t-1] = 1;
		}

		int nIncrement = 0;
		for( j = 0; j < nGrid; j++)
		{
			memset(visited, 0, sizeof(visited));
			if(dfs(j))
			{
				nIncrement++;
			}
		}
		printf("%d\n", nIncrement);
	}

	return 0;
}