二分图匹配——匈牙利算法

---

二分图

二分图又称作二部图，是图论中的一种特殊模型。 设G=(V,E)是一个无向图，如果顶点V可分割为两个互不相交的子集(A,B)，并且图中的每条边（i，j）所关联的两个顶点i和j分别属于这两个不同的顶点集(i in A,j in B)，则称图G为一个二分图。

简而言之，就是顶点集V可分割为两个互不相交的子集，并且图中每条边依附的两个顶点都分属于这两个互不相交的子集，两个子集内的顶点不相邻。

二分图匹配

给定一个二分图G，在G的一个子图M中，M的边集{E}中的任意两条边都不依附于同一个顶点，则称M是一个匹配。

最大匹配

选择边数最大的子图称为图的最大匹配问题。
如果一个匹配中，图中的每个顶点都和图中某条边相关联，则称此匹配为完全匹配，也称作完备匹配。

增广路径

也称增广轨或交错轨。设M为二分图G已匹配边的集合，若P是图G中一条联通两个未匹配顶点的路径，且属于M的边和不属于M的边在P上交替出现，则称P为相对于M的一条增广路径。P的起点在X部，终点在Y部，反之亦可，路径P的起点终点都是未匹配的点。

增广路径是一条“交错轨”。也就是说, 它的第一条边是目前还没有参与匹配的,第二条边参与了匹配,第三条边没有…最后一条边没有参与匹配,并且起点和终点还没有被选择过，这样交错进行,显然P有奇数条边，因为不属于匹配的边比匹配的边要多一条。

增广路径性质

由增广路的定义可以推出下述三个结论：

1.P的路径长度必定为奇数，第一条边和最后一条边都不属于M，因为两个端点分属两个集合，且未匹配（单独的一条连接两个未匹配点的边显然也增广路径）。
2.P经过取反操作可以得到一个更大的匹配M’。
3.M为G的最大匹配当且仅当不存在相对于M的增广路径。

匈牙利算法

用增广路求最大匹配称作匈牙利算法。
算法轮廓：

1.置M为空。
2.找出一条增广路径P（先找出单个边的增广路径），碰到多个边的增广路径通过取反操作获得更大的匹配M’代替M。
3.重复2操作直到找不出增广路径为止。

找增广路径的算法

我们采用深度优先搜索的办法找一条增广路径：

从X部一个未匹配的顶点Xi开始，找一个未访问的邻接点Yi（Yi一定是Y部顶点）。对于Yi，分两种情况：

如果Yi未匹配，这俩相连，则就算已经找到一条增广路（是单个边的增广路径）
如果Yi已经匹配，则取出yi的匹配顶点Xm(Xm一定是X部顶点)，边(Xm,Yi)目前是匹配的，根据“取反”的想法，要将(Xm,Yi)改为未匹配，(Xi,Yi)设为匹配，能实现这一点的条件是看从Xm为起点能否新找到一条增广路径P’。如果行，则Xi-Yi-P’就是一条以Xi为起点的增广路径。

然后取反，则(Xi,Yi)就匹配上了，总的匹配的边要多一条了。

代码实现

下面的代码求的是最大匹配数量：

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<stdio.h>
 
using namespace std;
 
const int N = 505; 
bool line[N][N];  
int result[N]; 
bool used[N]; 
int k, m, n;
bool found(int x)
{
	for (int i = 1; i <= n; i++)
	{
		if (line[x][i] && !used[i])
		{
			used[i] = true;
			if (result[i] == 0 || found(result[i]))
			{
				result[i] = x;
				return true;
			}
		}
	}
	return false;
}
int main()
{
	int x, y;
	printf("请输入相连边的数量k：\n");
	while (scanf("%d", &k) && k)  
	{
		printf("请输入二分图中x和y中点的数目：\n");
		scanf("%d %d", &m, &n); 
		memset(line, 0, sizeof(line));
		memset(result, 0, sizeof(result));
		for (int i = 0; i < k; i++)
		{
			printf("请输入相连边的两个点：\n");
			scanf("%d %d", &x, &y);
			line[x][y] = 1;
		}
		int sum = 0;
		for (int i = 1; i <= m; i++)
		{
			memset(used, 0, sizeof(used)); 
			if (found(i)) sum++;
		}
		printf("%d\n", sum);
	}
	return 0;
}