二分图匹配相关算法与题目讲解

来源:互联网 发布:泰国4g网络制式 频段 编辑:程序博客网 时间:2024/05/02 06:13

二分图匹配相关算法讲解


张天翔

blog.csdn.net/hzoi_ztx
ztx97@qq.com

一些基本概念

二分图

二分图(BipartiteGraph)又称作二部图,是图论中的一种特殊模型。

G=(V,E)是一个无向图,如果顶点V可分割为两个互不相交的子集(X,Y),并且图中的每条边(i,j)所关联的两个顶点ij分别属于这两个不同的顶点集(iX,jY),则称图G为一个二分图

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

如果GX的每个顶点都与Y的每个顶点相邻,则称G完全二分图。

注意,二分图是无向图,但是有些问题虽然是有向图的外壳,实际上还是二分图的解决方法,后面会有一些例子。

/////X/////////Y//////////////////////X/X/X/////////

图为一个二分图

匹配

匹配是在图中定义的,二分图中的匹配也是下面这些定义,你可以动手画一画。

G=(V,E)是一个图,ME的一个子集,如果M不含环且任意两边都不相邻,则称MG的一个匹配
G中边数最多的匹配称为G最大匹配

下面设MG的一个匹配。

对于图G=(V,E),在每条边ei上赋一个实数权wi。定义eiMiwi为匹配M的权。G中权最大的匹配称为G最大权匹配。其中,若eiEwi=c,(c为非负常数),则G的最大权匹配就是G的最大匹配。上面的例图显然最大匹配数是4

viVviM中的边相关联,则称viM饱和点,否则称viM非饱和点
如果G中每个顶点都是M饱和点,则称MG完美(备)匹配
显然,XY点集大小需要相等。

举个栗子,我有5只母鸡,4只公鸡,可以看上面那张图。
每只母鸡和公鸡都有自己的好朋友(用边表示),它们可以生下一只小鸡,不是好朋友的则不可以产生后代。我决定让它们实行一夫一妻制。问,如何安排它们才能得到最多的小鸡?拒绝二胎和乱伦 这就是最大匹配的一个问题。。。

PG的一条链。如果P的边交替地一条是M中的边,一条不是M中的边,则称PM交错轨。类似地,我们可以定义G交错圈。易知,G的交错圈一定是偶圈。

一条连接两个不同的M非饱和点的M交错链称为M增广轨

增广轨定理

广

所有匹配算法都是基于增广轨定理,这个定理适用于任意图

覆盖

覆盖同样是在图中定义的,但因为二分图中有其特殊性,直接给出二分图中覆盖的相关概念。
覆盖分为点覆盖边覆盖,分别是用一个点集关联所有边(使所有的边至少有一个端点属于此点集)和用一个边集关联所有点(使所有的点为此边集中至少一条边的端点)。

最小覆盖(最小点覆盖)点最少的点覆盖
König定理:

=

证明:

在最大匹配M中,我们取每一条匹配边上的一个端点(当然不是随便取的,后面会介绍算法求此点集),个数为m.
(1)m个点是足够的。只需要让它们覆盖最大匹配的m条边,则其它边一定被覆盖(如果有边e不被覆盖,把e加入后得到一个更大的匹配)
(2)m个点是必需的。仅考虑形成最大匹配的这m条边,由于它们两两无公共点,因此至少需要m个点才能把它们覆盖。
证毕。

举个栗子,任务安排(黑书P331) 问题。

最小边覆盖边最少的边覆盖

独立集

独立集即一个点集,集合中任两个结点不相邻,则称V独立集。或者说是导出的子图是零图(没有边)的点集。

最大独立集点最多的独立集。
独立数:最大独立集的点。

=|G|

证明:

很显然的把最大匹配两端的点都从顶点集中去掉这个时候剩余的点是独立集,这是|G|2|M|,同时必然可以从每条匹配边的两端取一个点加入独立集并且保持其独立集性质,故为|G||M|

举个栗子,棋盘上的骑士问题(黑书P332)。棋盘上放置最多的骑士(马)不互相攻击。

边独立集即一个边集,满足边集中的任两边不邻接。
最大边独立集:边最多的边独立集。
边独立数*:最大边独立集的边数。
可以看到:边独立集就是匹配,最大匹配即最大边独立集。一般情况下称匹配。

一个无向图G=(V,E)。取 V 的一个子集U,若对于U中任意两个点 uv,有边(u,v)E,那么称 UG 的一个完全子图。 U 是一个当且仅当 U 不被包含在一个更大的完全子图中。
G最大团指的是定点数最多的一个团。

路径覆盖

最小路径覆盖:是“路径” 覆盖“点”,即用尽量的不相交简单路径覆盖* 有向无环图G *的所有顶点,即每个顶点严格属于一条路径。路径的长度可能为0(单个点)。
可以直观地想象一下,在一个有向无环图中至少放几个人才能走到所有点?最少的人数就是最小路径覆盖数。
容易想到

|G|

应该使得最小路径覆盖中的边数尽量多,同时不让两条边在同一个顶点相交,这样才能得到最小路径覆盖。

有向无环图(DAG)的最小路径覆盖,可以用二分图最大匹配来求解。

=|G|

要将有向无环图G转化为二分图G,具体操作方法为拆点:将每一个顶点i拆成两个顶点xiyi。然后根据原图中边的信息,如果原图存在边ij,那么就在新图中添加边xiyj,这样所有边的方向都是由X部到Y部。于是所转化出的二分图G的最大匹配数是原图G中最小路径覆盖上的边数。即:
=

证明:

(1)如果匹配数为零,那么G中不存在有向边,于是显然有:

=|G|=|G|0=|G|

G的最小路径覆盖数为G
(2)G中暂时不存在匹配边时,路径覆盖数为|G|。此时如果在G中增加一条匹配边xiyj,那么在图P的路径覆盖中就存在边ij,也就是说边ij在某一条路径上,于是路径覆盖数就可以减少一个。如此继续增加匹配边,每增加一条,路径覆盖数就减少一条;直到匹配边不能继续增加时,路径覆盖数也不能再减少了,此时就有了前面的公式。
但是这只是说明了每条匹配边对应于路径覆盖中的一条路径上的某一条有向边,下面来说明一个路径覆盖中的每条有向边对应于一条匹配边。与前面类似,对于路径覆盖中的每条连接两个顶点之间的每条有向边ij,我们可以在匹配图中对应做一条边xiyj, 显然这样做出来图的是一个匹配图(这一点用反证法很容易证明,如果得到的图不是一个匹配图,那么这个图中必定存在这样两条边 xiyjxiyk(jk),那么在路径覆盖图中就存在了两条边ij, ik,那从i出发的路径就不止一条了,这与路径覆盖图是矛盾的;还有另外一种情况就是存在xiyj,xkyj,这种情况也类似可证)
至此,就说明了匹配边与路径覆盖图中连接两顶点之间边的一一对应关系,那么也就说明了前面的公式成立!

铲雪车问题。

二分图的性质

稍稍总结一下二分图的性质

  • 最大匹配数=最小覆盖数
  • DAG最小路径覆盖数=DAG顶点数-二分图最大匹配数
  • 最小边覆盖 = 顶点数 - 最大匹配数 = 最大独立集(最小边覆盖可以看做最小路径覆盖在二分图上的特殊情况)

二分图的判定

要求解二分图问题,首先需要判定其为二分图。
1.

二分图:有两顶点集且图中每条边的的两个顶点分别位于两个顶点集中,每个顶点集中没有边直接相连接

2.

无向图G为二分图的充分必要条件是:G至少有两个顶点,且其所有回路的长度均为偶数。

3.

判断二分图的常见方法是染色法:开始对任意一未染色的顶点染色,之后判断其相邻的顶点中,若未染色则将其染上和相邻顶点不同的颜色,若已经染色且颜色和相邻顶点的颜色相同则说明不是二分图,若颜色不同则继续判断,bfs和dfs可以搞定!(易知:任何无回路的图均是二分图。)

一些二分图求解问题

二分图最大匹配求解

问题

给定一个二分图G,求G的最大匹配数

匈牙利算法

算法思想:

根据一个匹配是最大匹配当且仅当没有增广路,求最大匹配就是找增广轨,直到找不到增广轨,就找到了最大匹配。遍历每个点,查找增广路,若找到增广路,则修改匹配集和匹配数,否则,终止算法,返回最大匹配数。

增广路径必须满足的性质

  1. 有奇数条边。
  2. 起点在二分图的左半边,终点在右半边。
  3. 路径上的点一定是一个在左半边,一个在右半边,交替出现。(其实二分图的性质就决定了这一点,因为二分图同一边的点之间没有边相连,不要忘记哦。)
  4. 整条路径上没有重复的点。
  5. 起点和终点都是目前还没有配对的点,而其它所有点都是已经配好对的。
  6. 路径上的所有第奇数条边都不在原匹配中,所有第偶数条边都出现在原匹配中。
  7. 最后,也是最重要的一条,把增广路径上的所有第奇数条边加入到原匹配中去,并把增广路径中的所有第偶数条边从原匹配中删除(这个操作称为增广路径的取反),则新的匹配数就比原匹配数增加了1个(奇数=偶数+1)。而每次查找得到的增广路径的长度都是在上一次查找到的增广路径的基础上延伸的,这样每次更新匹配数都是增加1。

匈牙利算法步骤

G=(X,,Y)是一个二分图,其中,X={x1,x2,...xm},Y={y1,y2,...yn}。令MG中的任一个匹配
1. 置M为空
2. 从G中找出一个未匹配点v(增广路性质5要求的),如果没有则算法结束,否则,以v为起点,查找增广路(邻接点是为未匹配点,则返回寻找完成,若v的邻接点u是匹配点,则从u开始查找,直至查找到有未匹配点终止)即满足增广路的性质,如果没有找到增广路,则算法终止
3. 找出一条增广路径P,通过异或(取反)操作获得更大的匹配M代替M(方便要输出增广矩阵以及进一步查找),匹配数加1(性质7得到)
4. 重复(2)(3)操作直到找不出增广路径为止

彻底理解增广路

  1. 总是从X集的未匹配点出发,寻找匹配点或者未匹配点,如查找到未匹配点则该增广路终止,否则以该点的增广路不存在。
  2. 每次查找增广路都是在之前形成的匹配(上面步骤3中异或后的匹配)的基础上进行延伸的,也就是查找匹配点总是在匹配M中,其实就是用起点和终点两个未匹配点将得到匹配的边尽可能的连接起来的增广路,这样增广路长度就延长了,当然也可以是直接就是以两个个未匹配点的边(就直接添加进匹配中)。总而言之,每次扩充匹配不是通过延伸增广路径就是新增增广路径(当然长度为1)。

时间空间复杂度

时间复杂度
邻接矩阵:最坏为O(n3) 邻接表:O(mn)
空间复杂度
邻接矩阵:O(n2)邻接表:O(m+n)

匈牙利算法实现

匈牙利算法只需要以每个节点为起点找一次增广路即可求得最大匹配,寻找增广路的复杂度为O(E),总的复杂度为O(VE)
下面的实现是查找从X到Y的匹配,X中每一个点最多只能被遍历一次用于查找增广路(当已经是匹配点是就不遍历了)

// 由于时间紧凑,代码摘自http://dsqiu.iteye.com/blog/1689505#define maxn 10//表示x集合和y集合中顶点的最大个数!  int nx,ny;//x集合和y集合中顶点的个数  int edge[maxn][maxn];//edge[i][j]为1表示ij可以匹配  int cx[maxn],cy[maxn];//用来记录x集合中匹配的y元素是哪个!  int visited[maxn];//用来记录该顶点是否被访问过!  int path(int u){      int v;      for(v=0;v<ny;v++)      {          if(edge[u][v]&&!visited[v])          {              visited[v]=1;             if(cy[v]==-1||path(cy[v]))//如果y集合中的v元素没有匹配或者是v已经匹配,但是从cy[v]中能够找到一条增广路              {                  cx[u]=v; //找到增广路,修改匹配M                  cy[v]=u;                  return 1;              }          }      }      return 0;  }  int maxmatch()  {      int res=0;      memset(cx,0xff,sizeof(cx));//初始值为-1表示两个集合中都没有匹配的元素!      memset(cy,0xff,sizeof(cy));      for(int i=0;i<=nx;i++)      {          if(cx[i]==-1)   //还没被匹配,执行内部代码          {              memset(visited,0,sizeof(visitited));  //重置标记为为访问              res+=path(i);   //以 i 为起点开始查找增广路,返回true ,匹配数+1          }      }      return res;  }

Hopcroft-Karp算法

太长,不会,日后请自行了解

二分图最小覆盖集和最大独立集的构造

二分图有下面两个性质:

定理1:最小覆盖数 = 最大匹配数定理2:最大独立集S 与 最小覆盖集T 互补

算法实现步骤:

  1. 做最大匹配,没有匹配的空闲点∈S
  2. 如果u∈S那么u的邻点必然属于T
  3. 如果一对匹配的点中有一个属于T那么另外一个属于S
  4. 还不能确定的,把左子图的放入S,右子图放入T算法结束

二分图最小路径覆盖求解

有向无环图最小不相交路径覆盖

定义:用最少的不相交路径覆盖所有顶点。算法:把原图中的每个点i拆成Xi和Yi,如果有一条有向边i->j,那么就加边Xi-Yj。这样就得到了一个二分图,下面的任务就是求最大匹配,然后最小路径覆盖=原图的节点数-新图最大匹配。证明见前文。

有向无环图最小可相交路径覆盖

定义:用最小的可相交路径覆盖所有顶点。算法:先用floyd求出原图的传递闭包(即连通性),即如果i到j有路,那么就加边i->j。然后就转化成了最小不相交路径覆盖问题。

带权二分图最佳匹配求解

问题

对于二分图的每条边都有一个权(非负),要求一种完备匹配方案,使得所有匹配边的权和最大,记做最佳(优)完美(备)匹配。(特殊的,当所有边的权为1时,就是最大完美(备)匹配问题)

二分图最大匹配和二分图完美匹配的区别

二分最大匹配:使得匹配中的边数量不小于任何其他的匹配。
二分完备匹配:使得图中所有点出现在该匹配中。

二分图的带权匹配与二分图的最佳匹配的区别

二分图的带权匹配:求出一个匹配集合,使得集合中边的权值之和最大或最小。
二分图的最佳匹配:匹配必须为完备匹配,在此基础上,才要求匹配的边权值之和最大或最小。
二分图的带权匹配与最佳匹配不等价,也不互相包含。

Kuhn-Munkers算法求二分图的最佳匹配
Kuhn-Munkers算法可以实现为O(N^3)的时间复杂度。
这算法太神///会,你们暂时用不到/////////过,可以自行了解。
推荐网址
http://philoscience.iteye.com/blog/1754498
http://www.linuxidc.com/Linux/2013-02/79458p2.htm

二分图最大团求解

例题:已知班级有g个女孩和b个男孩,所有女生之间都相互认识,所有男生之间也相互认识,给出m对关系表示哪个女孩与哪个男孩认识。现在要选择一些学生来组成一个团,使得里面所有人都认识,求此团最大人数。

定理:

原图的最大团=补图的最大独立集

原图的最大独立集=补图的最大团。

由于这个题的补图显然是一个二分图,而二分图的补图的最大独立集可以由匈牙利算法求的,所以该题的最大团问题可以转化成补图的最大独立集来做。

推荐网址
http://www.cnblogs.com/pushing-my-way/archive/2012/08/09/2629700.html
http://www.cnblogs.com/zhj5chengfeng/archive/2013/07/29/3224092.html

1 0
原创粉丝点击