二分图最大匹配及常用建图方法.docVIP

算法———艺术 二分图匹配剖析 很多人说，算法是一种艺术。但是对于初学者的我，对算法认识不是很深刻，但偶尔也能感受到他强大的魅力与活力。 这让我追求算法的脚步不能停止。下面我通过分析匈牙利算法以及常用建图方式，与大家一起欣赏算法的美。 匈牙利算法 匈牙利算法是用来解决最大二分图匹配问题的，所谓二分图即 “一组点集可以分为两部分，且每部分内各点互不相连，两部分的点之间可以有边”。所谓最大二分图匹配即 ”对于二分图的所有边，寻找一个子集，这个子集满足两个条件， 1：任意两条边都不依赖于同一个点。 2：让这个子集里的边在满足条件一的情况下尽量多。 首先可以想到的是，我们可以通过搜索，找出所有的这样的满足上面条件的边集，然后从所有的边集中选出边数最多的那个集合，但是我们可以感觉到这个算法的时间复杂度是边数的指数级函数，因此我们有必要寻找更加高效的方法。目前比较有效的方法有匈牙利算法和通过添加汇点和源点的网络流算法，对于点的个数都在200 到300 之间的数据，我们是采取匈牙利算法的，因为匈牙利算法实现起来要比网络流简单些。下面具体说说匈牙利算法： 介绍匈牙利之前，先说说“增广轨”。 定义： 若P是图G中一条连通两个未匹配顶点的路径，并且属最大匹配边集M的边和不属M的边(即已匹配和待匹配的边)在P上交替出现，则称P为相对于M的一条增广轨 定义总是抽象的下面通过图来理解它。 图中的线段（2-3, 3-1, 1-4）便是上面所说的p路径， 我们假定边（1,3）是以匹配的边，（2,3）（1,4）是未匹配 的边，则边（4,1）边（1,3）和 边（3,2） 在路径p上交 替的出现啦，那么p就是相对于M的一条增广轨，这样我 们就可以用边1,4 和 边2,3来替换边1,3 那么以匹配的边 集数量就可以加1,。 匈牙利算法就是同过不断的寻找增广轨实现的。很明显如果二分图的两部分点分别为n和m，那么最大匹配的数目应该小于等于MIN（n,m）; 因此我们可以枚举任第一部分（的二部分也可以）里的每一个点，我们从每个点出发寻找增广轨，最后吧第一部分的点找完以后，就找到了最大匹配的数目，当然我们也可以通过记录找出这些边。 下面给出关于二分图最大匹配的两个定理 1:最大匹配数 + 最大独立集 = n + m 2:二分图的最小覆盖数 = 最大匹配数 3:最小路径覆盖 = 最大独立集 最大独立集是指求一个二分图中最大的一个点集，该点集内的点互不相连。 最小顶点覆盖是指 在二分图中，用最少的点，让所有的边至少和一个点有关联。 最小路径覆盖是指一个不含圈的有向图G中，G的一个路径覆盖是一个其结点不相交的 路径集合P，图中的每一个结点仅包含于P中的某一条路径。路径可以从任意结点 开始和结束，且长度也为任意值，包括0 下面给出二分图最大匹配的伪代码： /* use[i] 数组时标记第二部分点中的第j个点是否使用过。 match[i] 与第二部分中的点i 匹配的点是 match[i] */ Int GetAugmentPath ( number ) //通过number这个点寻找增广轨，找到返回1找不到返回0； P ? related[number] . next; // 与 number 相连的点； While（p != NULL）{ If not used[p] then // p 点没有被使用过 If match[p] == 0 or GetAugmentPath(match[p]) // p点没有和另一部分的点配对 // 或 和p配对的那个点可以找到其它点和他配对（找到增广轨 ） Return 1; p ? p . next // 与 number 相连的下几个点 } Return 0; } end GetAugmentPath; 我们只需要在主程序中对某一部分点集的每个点进行增广轨的寻找； Int main{ ans ? 0 For I = 1 ? n // 枚举第一部分的n个点； For j = 1 ? m use[j] ? false // 把第二部分的m个点表示为没有使用过。 if GetAugmentPath (i) then ans ? ans + 1 // 找到增广轨就给边数加一； } 程序非常好写，也非常好懂。 1