"《算法导论》之‘图’":不带权二分图最大匹配(匈牙利算法)

  博文“二分图的最大匹配、完美匹配和匈牙利算法”对二分图相关的几个概念讲的特别形象,特别容易理解。本文介绍部分主要摘自此博文。

  还有其他可参考博文:

  趣写算法系列之--匈牙利算法

  用于二分图匹配的匈牙利算法

  1. 前言

  二分图:简单来说,如果图中点可以被分为两组,并且使得所有边都跨越组的边界,则这就是一个二分图。准确地说:把一个图的顶点划分为两个不相交集 U 和V ,使得每一条边都分别连接U、V中的顶点。如果存在这样的划分,则此图为一个二分图。二分图的一个等价定义是:不含有「含奇数条边的环」的图。图 1 是一个二分图。为了清晰,我们以后都把它画成图 2 的形式。

  匹配:在图论中,一个「匹配」(matching)是一个边的集合,其中任意两条边都没有公共顶点。例如,图 3、图 4 中红色的边就是图 2 的匹配。

      

  我们定义匹配点匹配边未匹配点非匹配边,它们的含义非常显然。例如图 3 中 1、4、5、7 为匹配点,其他顶点为未匹配点;1-5、4-7为匹配边,其他边为非匹配边。

  最大匹配:一个图所有匹配中,所含匹配边数最多的匹配,称为这个图的最大匹配。图 4 是一个最大匹配,它包含 4 条匹配边。

  完美匹配:如果一个图的某个匹配中,所有的顶点都是匹配点,那么它就是一个完美匹配。图 4 是一个完美匹配。显然,完美匹配一定是最大匹配(完美匹配的任何一个点都已经匹配,添加一条新的匹配边一定会与已有的匹配边冲突)。但并非每个图都存在完美匹配。

  举例来说:如下图所示,如果在某一对男孩和女孩之间存在相连的边,就意味着他们彼此喜欢。是否可能让所有男孩和女孩两两配对,使得每对儿都互相喜欢呢?图论中,这就是完美匹配问题。如果换一个说法:最多有多少互相喜欢的男孩/女孩可以配对儿?这就是最大匹配问题。

  

  基本概念讲完了。求解最大匹配问题的一个算法是匈牙利算法,下面讲的概念都为这个算法服务。

  

  交替路:从一个未匹配点出发,依次经过非匹配边、匹配边、非匹配边…形成的路径叫交替路。

  增广路:从一个未匹配点出发,走交替路,如果途径另一个未匹配点(出发的点不算),则这条交替路称为增广路(agumenting path)。例如,图 5 中的一条增广路如图 6 所示(图中的匹配点均用红色标出):

  

  增广路有一个重要特点:非匹配边比匹配边多一条。因此,研究增广路的意义是改进匹配。只要把增广路中的匹配边和非匹配边的身份交换即可。由于中间的匹配节点不存在其他相连的匹配边,所以这样做不会破坏匹配的性质。交换后,图中的匹配边数目比原来多了 1 条。

  我们可以通过不停地找增广路来增加匹配中的匹配边和匹配点。找不到增广路时,达到最大匹配(这是增广路定理)。匈牙利算法正是这么做的。在给出匈牙利算法 DFS 和 BFS 版本的代码之前,先讲一下匈牙利树。

  匈牙利树一般由 BFS 构造(类似于 BFS 树)。从一个未匹配点出发运行 BFS(唯一的限制是,必须走交替路),直到不能再扩展为止。例如,由图 7,可以得到如图 8 的一棵 BFS 树:

         

  这棵树存在一个叶子节点为非匹配点(7 号),但是匈牙利树要求所有叶子节点均为匹配点,因此这不是一棵匈牙利树。如果原图中根本不含 7 号节点,那么从 2 号节点出发就会得到一棵匈牙利树。这种情况如图 9 所示(顺便说一句,图 8 中根节点 2 到非匹配叶子节点 7 显然是一条增广路,沿这条增广路扩充后将得到一个完美匹配)。

  2. 匈牙利算法实现

  匈牙利算法可基于DFS或BFS实现。本文实现的是基于DFS的。

  具体代码头文件定义为:

 1 #ifndef HUNGARIAN_H
 2 #define HUNGARIAN_H
 3
 4 #include "BFS_n_DFS.h" // 详细代码请转到Github.
 5
 6 class Graph_Hun : public Graph_BD
 7 {
 8 public:
 9     Graph_Hun();
10     virtual ~Graph_Hun();
11     virtual bool DFS(int vexID);
12     virtual bool BFS(int vexID);
13     int hungarian();
14
15     void setNumOfLeft(int num);
16     void setNumOfRight(int num);
17
18 private:
19     int matching[MAXNUM];
20     bool checked[MAXNUM];
21     int numOfLeft;
22     int numOfRight;
23 };
24
25 #endif

  实现文件为:

 1 #include "Hungarian.h"
 2 #include <cassert>
 3
 4 Graph_Hun::Graph_Hun()
 5 {
 6     memset(matching, -1, sizeof(matching));
 7     memset(checked, false, sizeof(checked));
 8 }
 9
10 Graph_Hun::~Graph_Hun()
11 {
12
13 }
14
15 bool Graph_Hun::DFS(int vexID)
16 {
17     assert(vexID >= 0 && vexID < MAXNUM);
18
19     vector<int> adjVexes = adj(vexID);
20     int sz = adjVexes.size();
21     for (int i = 0; i < sz; i++)
22     {
23         int adjVex = adjVexes[i];
24         if (!checked[adjVex])
25         {
26             checked[adjVex] = true;
27             if (matching[adjVex] == -1 || DFS(matching[adjVex]))
28             {
29                 matching[adjVex] = vexID;
30                 matching[vexID] = adjVex;
31                 return true;
32             }
33         }
34     }
35     return false;
36 }
37
38 bool Graph_Hun::BFS(int vexID)
39 {
40     return true;
41 }
42
43 int Graph_Hun::hungarian()
44 {
45     int maxMatching = 0;
46     for (int i = 0; i < numOfLeft; i++)
47     {
48         if (matching[i] == -1)
49         {
50             memset(checked, false, sizeof(checked));
51             if (DFS(i))
52                 maxMatching++;
53         }
54     }
55     return maxMatching;
56 }
57
58 void Graph_Hun::setNumOfLeft(int num)
59 {
60     numOfLeft = num;
61 }
62
63 void Graph_Hun::setNumOfRight(int num)
64 {
65     numOfRight = num;
66 }

Hungarian.cpp

  下边将用一个例子来说明利用匈牙利算法求二分图最大匹配的过程。作为例子的二分图如下:

  

  图2-1 一个二分图(其中0~1表示二分图左边,6~12表示二分图右边)

  详细的二分图最大匹配过程如下图所示:

  

  图2-2 详细的二分图最大匹配过程

  在子图1中,0直接连接到6;在子图2中,1直接连接到7;在子图3中,2试图连接到6(因为是DFS搜索),但6已经连接到0,这时候就要找增广路径了,找到的路径如下图:

  

  图2-3 增广路径“2->6->0->7->1->10”

这时将已匹配和未匹配的路径反转就能够得到图2-2中的子图3了。在图2-2的子图4中,3、4分别直接连接到8、9,没有寻找增广路径的过程。而对5来说,并未能找到增广路径,因此5无法匹配。最终的匹配结果如下图:

  

  图2-4 二分图最终匹配结果

  

  详细代码请参考自Github.

  

时间: 2024-12-10 15:20:42

"《算法导论》之‘图’":不带权二分图最大匹配(匈牙利算法)的相关文章

POJ1274:The Perfect Stall(二分图最大匹配 匈牙利算法)

The Perfect Stall Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 17895   Accepted: 8143 Description Farmer John completed his new barn just last week, complete with all the latest milking technology. Unfortunately, due to engineering pr

Ural1109_Conference(二分图最大匹配/匈牙利算法/网络最大流)

解题报告 二分图第一题. 题目描述: 为了参加即将召开的会议,A国派出M位代表,B国派出N位代表,(N,M<=1000) 会议召开前,选出K队代表,每对代表必须一个是A国的,一个是B国的; 要求每一个代表要与另一方的一个代表联系,除了可以直接联系,也可以电话联系,求电话联系最少 思路: 电话联系最少就要使直接联系最大,又是一一匹配关系,就是二分图的最大匹配. 下面是匈牙利算法. #include <cstdio> #include <cstring> #include <

UESTC 919 SOUND OF DESTINY --二分图最大匹配+匈牙利算法

二分图最大匹配的匈牙利算法模板题. 由题目易知,需求二分图的最大匹配数,采取匈牙利算法,并采用邻接表来存储边,用邻接矩阵会超时,因为邻接表复杂度O(nm),而邻接矩阵最坏情况下复杂度可达O(n^3). 代码: #include <iostream> #include <cstdio> #include <cstring> #include <cmath> #include <algorithm> #include <vector> u

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

这篇文章讲无权二分图(unweighted bipartite graph)的最大匹配(maximum matching)和完美匹配(perfect matching),以及用于求解匹配的匈牙利算法(Hungarian Algorithm):不讲带权二分图的最佳匹配. 二分图:简单来说,如果图中点可以被分为两组,并且使得所有边都跨越组的边界,则这就是一个二分图.准确地说:把一个图的顶点划分为两个不相交集 U  和 V ,使得每一条边都分别连接U . V  中的顶点.如果存在这样的划分,则此图为一

HDU - 1045 Fire Net (二分图最大匹配-匈牙利算法)

(点击此处查看原题) 匈牙利算法简介 个人认为这个算法是一种贪心+暴力的算法,对于二分图的两部X和Y,记x为X部一点,y为Y部一点,我们枚举X的每个点x,如果Y部存在匹配的点y并且y没有被其他的x匹配,那就直接匹配:如果Y中已经没有可以和x匹配的点(包括可以匹配的点已经被其他的x匹配),那就让已经匹配的y的原配x'寻找其他可以匹配的y’,并将y和x匹配,最后,统计出匹配的对数 (详细了解的话,可以看看这位的博客:https://blog.csdn.net/sunny_hun/article/de

poj - 3041 Asteroids (二分图最大匹配+匈牙利算法)

http://poj.org/problem?id=3041 在n*n的网格中有K颗小行星,小行星i的位置是(Ri,Ci),现在有一个强有力的武器能够用一发光速将一整行或一整列的小行星轰为灰烬,想要利用这个武器摧毁所有的小行星最少需要几发光束. 主要是构图,将每一行当成一个点,构成集合1,每一列也当成一个点,构成集合2,每一个障碍物的位置坐标将集合1和集合2的点连接起来,也就是将每一个障碍物作为连接节点的边,这样可以得出本题是一个最小点覆盖的问题==二分图的最大匹配. 就可以通过匈牙利算法求解.

二分图最大匹配---匈牙利算法BFS 实现

二分图指的是这样一种图,其所有顶点可以分成两个集合X和Y,其中X或Y中任意两个在同一集合中的点都不相连,所有的边关联在两个顶点中,恰好一个属于集合X,另一个属于集合Y.给定一个二分图G,M为G边集的一个子集,如果M满足当中的任意两条边都不依附于同一个顶点,则称M是一个匹配.图中包含边数最多的匹配称为图的最大匹配. 二分图的最大匹配有两种求法,第一种是最大流:第二种就是我现在要讲的匈牙利算法.这个算法说白了就是最大流的算法,但是它跟据二分图匹配这个问题的特点,把最大流算法做了简化,提高了效率. 增

51Nod 2006 飞行员配对(二分图最大匹配)-匈牙利算法

2006 飞行员配对(二分图最大匹配) 题目来源: 网络流24题 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题  收藏  关注 第二次世界大战时期,英国皇家空军从沦陷国征募了大量外籍飞行员.由皇家空军派出的每一架飞机都需要配备在航行技能和语言上能互相配合的2名飞行员,其中1名是英国飞行员,另1名是外籍飞行员.在众多的飞行员中,每一名外籍飞行员都可以与其他若干名英国飞行员很好地配合.如何选择配对飞行的飞行员才能使一次派出最多的飞机.对于给定的外籍飞行员与英国飞行员的

51N偶读2006 飞行员配对(二分图最大匹配)-匈牙利算法

2006 飞行员配对(二分图最大匹配) 题目来源: 网络流24题 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题  收藏  关注 第二次世界大战时期,英国皇家空军从沦陷国征募了大量外籍飞行员.由皇家空军派出的每一架飞机都需要配备在航行技能和语言上能互相配合的2名飞行员,其中1名是英国飞行员,另1名是外籍飞行员.在众多的飞行员中,每一名外籍飞行员都可以与其他若干名英国飞行员很好地配合.如何选择配对飞行的飞行员才能使一次派出最多的飞机.对于给定的外籍飞行员与英国飞行员的