匈牙利算法——最小点覆盖及行列建图法例题

P1524柯南开锁

时间限制 : 10000 MS   空间限制 : 65536 KB

问题描述

柯南决定深入OIBH组织内部, 一探虚实.他经过深思熟虑, 决定从OIBH组织大门进入...........
OIBH组织的大门有一个很神奇的锁.锁是由M*N个格子组成, 其中某些格子凸起(灰色的格子). 每一次操作可以把某一行或某一列的格子给按下去.

如果柯南能在组织限定的次数内将所有格子都按下去, 那么他就能够进入总部. 但是OIBH组织不是吃素的, 他们的限定次数恰是最少次数
请您帮助柯南计算出开给定的锁所需的最少次数.

输入格式

第一行 两个不超过100的正整数N, M表示矩阵的长和宽
以下N行 每行M个数 非0即1 1为凸起方格

输出格式

一个整数 所需最少次数

样例输入

4 4
0000
0101
0000
0100

样例输出

2

本题选自NKOJ1524。

图论题首先考虑如何建图，对于像本题这样，在矩阵之中且某一操作影响一整行或一整列的题目，应多考虑使用行列建图法。把行看成一种类型的点，把列看成另一种类型的点，以此构造二分图。行与列有交点（如本题的突出格子），便把对应的“行点”与“列点”连边，于是本题要解决的问题就变为了求把所有的边都覆盖完的“行点”或“列点”的最少点数，即最小点覆盖问题。

有结论：最小点覆盖=最大匹配数。

于是只需把建好的二分图的最小匹配数求出即可。

以下为代码：

#include<cstdio>
using namespace std;
char map[105][105];
int n,m,ans,road[105],Link[105];
bool Find(int x,int t)
{
	for(int i=1;i<=m;i++)
		if(map[x][i]=='1'&&road[i]!=t)
		{
			road[i]=t;
			if(!Link[i]||Find(Link[i],t))
			{
				Link[i]=x;
				return true;
			}
		}
	return false;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%s",&map[i][1]);
	for(int i=1;i<=n;i++) ans+=Find(i,i);
	printf("%d\n",ans);
	return 0;
}