首页 > 其他分享 >二分图与匹配 I :二分图的最大匹配

二分图与匹配 I :二分图的最大匹配

时间:2023-08-08 15:44:13浏览次数:26  
标签:二分 匹配 最大 增广 int ans matching

引入:什么是二分图,什么是匹配

口头语言描述:一个图,你把他的点集划为两个集合,让每个集合之间的点没有连边,就是一个二分图。

二分图的一个等价定义是:不含有奇数条边的环的图就是一个二分图。

证明:显然,观察每一条路径,都是从一个点集走向另外一个点集,则一个环必定得走偶数次。

匹配,则是一个边集。在二分图中选择一些边玩连连看游戏,如果选择的边两两没有公共点,这个边集就是原图的一个匹配。

为了方便后面的描述,我们再引入一些概念:

  • 匹配点:被匹配边连接的点;
  • 匹配边:匹配里的边;
  • 非匹配点:没被匹配边连接的点;
  • 非匹配边:不在匹配里的边;
  • 最大匹配:一个图所有匹配中,所含匹配边数最多的匹配,称为这个图的最大匹配。
  • 完美匹配:如果一个图的某个匹配中,所有的顶点都是匹配点,那么它就是一个完美匹配。

显然,完美匹配一定是最大匹配,但并非每个图都存在完美匹配。

我们有时候希望匹配最大,那该怎么办呢?

解决:匈牙利算法

为了理解这个算法,我们给出几个概念:

  • 交替路:从一个未匹配点出发,依次经过非匹配边,匹配边,非匹配边……的路径;
  • 增广路:从一个未匹配点出发,走交替路,如果途径另一个未匹配点(出发的点不算),则这条交替路称为增广路。

增广路有一个重要特点:非匹配边比匹配边多一条,我们可以通过不停地找增广路来增加匹配中的匹配边和匹配点。找不到增广路时,达到最大匹配(这是增广路定理)。这就是匈牙利算法。

  1. 从左边第 1 个顶点开始,挑选未匹配点进行搜索,寻找增广路;
    1. 如果经过一个未匹配点,说明寻找成功。更新路径信息,匹配边数 +1,停止搜索;
    2. 如果一直没有找到增广路,则不再从这个点开始搜索。事实上,我们可以永久性地把它从图中删去,而不影响结果。
  2. 由于找到增广路之后需要沿着路径更新匹配,所以我们需要一个结构来记录路径上的点。

代码使用 BFS 探路,复杂度为 \(O(|V||E|)\)。

int Hungarian() {
    int ans = 0;
    memset(matching, -1, sizeof matching), memset(check, -1, sizeof check);
    for (int i = 1; i <= leftnode; ++ i) {
        if (matching[i] == -1) { // 点 i 还未被匹配过 
            // 清空队列,开始对 i 进行探路
            while (!q.empty()) q.pop();
            q.push(i);
            prev[i] = -1;
            bool flag = false; 
            
            // 开始探路 
            while(!q.empty() && !flag) {
                auto u = Q.front();
                for(int i = head[u] ; i; i = G[u].nxt) {
                    int v = G[i].to;
                    if(check[v] != i) {
                        check[v] = i;
                        Q.push(matching[v]);
                        if(matching[v] >= 0) prev[matching[v]] = u; // 找到匹配点
                        else { // 找到未匹配点,交替路变为增广路
                            flag = true;
                            int d = u, e = v;
                            while(d != -1) {
                            	// 一些存边操作 
                                int t = matching[d];
                                matching[d] = e, matching[e] = d;
                                d = prev[d], e = t;
                            }
                        }
                    }
                }
                q.pop();
            }
            if (matching[i] != -1) ++ ans; // 如果到最后 匹配成功,那么ans ++  
        }
    }
    return ans;
}

标签:二分,匹配,最大,增广,int,ans,matching
From: https://www.cnblogs.com/sdltf/p/17614280.html

相关文章

  • 代码随想录算法训练营第十三天| 239. 滑动窗口最大值 347.前 K 个高频元素 总结
    239.滑动窗口最大值 (一刷至少需要理解思路)    卡哥建议:之前讲的都是栈的应用,这次该是队列的应用了。本题算比较有难度的,需要自己去构造单调队列,建议先看视频来理解。    题目链接/文章讲解/视频讲解:https://programmercarl.com/0239.%E6%BB%91%E5%8A%A8%E7%AA%......
  • DataFrame 对某列求和、平均值、计数、最大值、最小值
    importpandasaspd#创建示例DataFramedata={'A':[1,2,3,4,5],'B':[10,20,30,40,50]}df=pd.DataFrame(data)#对列'B'求和column_sum=df['B'].sum()print("SumofcolumnB:",column_sum......
  • Codeforces 1857D:Strong Vertices 与图论无关的出度最大统计
    1857D.StrongVerticesDescription:给定两个长度均为\(n\)的数组\(a\)和\(b\)(编号\(1\)~\(n\)),如果\(a_u-a_v\geqb_u-b_v\)\((u\neqv)\),那么从\(u\)到\(v\)建立一条有向边。"Strong"定义为:一个点\(V\)可以经过有向图中合法的通路到达其他所有的点。请求解出"......
  • js replace方法 (字典表匹配替换字符)
    常规replace使用consttext1='abcdefg'consttext2=text1.replace('bc','00')//text2='a00defg'字典表匹配替换constreplacements={'<':'<','>':'>',......
  • 网工这行最大的内卷,就是养了一批“HCIE”
    下午好,我是老杨。这段时间,我个人工作上的事儿挺多,一直在处理。手底下的一个重要项目组,想招个新人,面了几批都不满意,难搞。前两天刚定下人选,现下和你聊上几句。上回说过,我的部门招人,HR过了一遍的简历,100%都是有IE证书的。但我面的这么多人里,虽然都有IE,但人才质量可谓是天差地别。十多......
  • ubuntu系统升级软件sudo apt upgrade后GPU崩溃报错,显示驱动版本不匹配——ubuntu系统
     ubuntu系统升级软件(sudoaptupgrade)后,GPU崩溃报错,查看系统日志: Aug206:25:02lcwtrsyslogd:[originsoftware="rsyslogd"swVersion="8.32.0"x-pid="2059"x-info="http://www.rsyslog.com"]rsyslogdwasHUPedAug207:17:01lcwtCRON......
  • 二分图相关定理
    最长反链:一张有向无环图的最长反链为一个集合\(S\subseteqV\),满足对于\(S\)中的任意两个不同的点\(u,v\inS(u\nev)\),\(u\)不能到达\(v\),\(v\)也不能到达\(u\),且\(S\)的大小尽量大最小不可重链覆盖:在DAG中选出若干条链,经过每个点一次,且链数尽量少最小点覆盖:......
  • 7-18 二分法求多项式单根 (20分)
    7-18 二分法求多项式单根 (20分)二分法求函数根的原理为:如果连续函数f(x)在区间[a,b]的两个端点取值异号,即f(a)f(b)<0,则它在这个区间内至少存在1个根r,即f(r)=0。二分法的步骤为:检查区间长度,如果小于给定阈值,则停止,输出区间中点(a+b)/2;否则如果f(a)f(b)<0,则计算中点的值f((a+b)/2)......
  • Pythonre.compile:用于优化正则表达式匹配的工具
    https://blog.csdn.net/www_xuhss_com/article/details/130858409?spm=1001.2101.3001.6650.2&utm_medium=distribute.pc_relevant.none-task-blog-2%7Edefault%7EYuanLiJiHua%7EPosition-2-130858409-blog-86527810.235%5Ev38%5Epc_relevant_sort_base3&depth_1-utm_......
  • 使用SQL 统计表最大值,最小值,平均值,中位数值
    问:有一个postgres数据表,表中有update_time,create_time关于时间戳的字段,选择某个时间段,计算出update_time减去create_time的值做为耗时时间的字段duration_time,统计出在这段时间内耗时时间的平均值,中位数值,最大值,和最小值,请问如何通过sql语句实现此统计?答:您好!您可以使用以下......