ICode9

精准搜索请尝试: 精确搜索
首页 > 编程语言> 文章详细

学习笔记 匈牙利算法

2020-08-12 06:32:17  阅读:326  来源: 互联网

标签:匹配 增广 int 匈牙利 mch 笔记 vis 算法


众所周知,二分图匹配是二分图理论中的基础,而匈牙利算法是一种求解它的基本算法。事实上,匈牙利算法是一个十分简洁的算法,简介到让人感到惊诧。下面就让我们了解一下这个神奇的算法。

先看几个定义:

  • 匹配 匹配是一个的集合,其中任意两条边都没有公共端点。
  • 最大匹配 顾名思义,包含边数最多的匹配
  • 交错路 匹配边和非匹配边依次出现的一条路
  • 增广路 从非匹配边出发,以非匹配边结束的交错路(增广路长度一定是奇数)

你可能已经被这几个定义搞懵了不过没关系,慢慢看下去你就懂了。

匈牙利算法的基本思想,就是不断地扩大匹配的集合。事实上,他就是一个不断尝试让匹配边数加一的过程。初始状态下,所有的边都是非匹配边,而当算法结束时,我们再也不能让匹配边数变多了,最大匹配数也就求出来了,还能附带求出一种方案。

匈牙利算法基于两个性质:

  • 一个匹配\(M\)是最大匹配,当且仅当图中不存在增广路。
  • 从某一点出发,如果某个时刻找不到增广路,那么你将永远找不到从该点出发的最短路。

性质1的必要性是显然的(因为增广路长度是奇数),我们先把这两条性质抛开,先来看看这个算法的过程。

我们首先给二分图定向(名字挺高大上,其实就是划分出左部和右部,规定增广路从左部出发)。枚举每一个左部点(想想这个点会不会已经被匹配),如果找不到增广路就算了,如果找到了增广路,我们就把这条路径上所有边的状态取反,于是我们就得到了一种大小+1的匹配。结合代码理解一下 复杂度\(O(nm)\)

int vis[mxn],mch[mxn];//match
inline bool hungary(int x){
    if(vis[x])return 0;vis[x]=1;
//vis记录左部点右部点皆可
    for(int i=first[x];i;i=nxt[i]) //边是单向边
        if(!mch[to[i]] || hungary(mch[to[i]]))return mch[to[i]]=x,1;
    return 0;
}

//main函数
int ans=0;
for(int i=1;i<=n;++i)memset(vis,0,sizeof(vis)),ans+=hungary(i);

如果你没看懂以上代码,那只能是心理原因 我们回过头看这两个性质,性质1显然是算法的基础,而性质2确保了每个点只搜一次。至此我们已经完全理解了匈牙利算法的内在逻辑。

有没有觉得这两个性质很神奇?如果感兴趣可以看看这篇论文,挺好理解的,性质2尤其好理解,有不懂的名词可以百度一下。再次强调,匹配是边集。

好像还有bfs实现的版本,据说常数更小。

以上。

标签:匹配,增广,int,匈牙利,mch,笔记,vis,算法
来源: https://www.cnblogs.com/happyguy/p/13485224.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有