题目在这里 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;
}