二分图匹配问题（未完待续）

最新推荐文章于 2024-08-06 20:40:43 发布
最新推荐文章于 2024-08-06 20:40:43 发布
阅读量209 收藏
点赞数
分类专栏： 匹配问题 文章标签： 算法 acm竞赛 图论
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_50966565/article/details/117425689
版权

二分图匹配问题（学习笔记与例题）

文章目录

知识点

一、二分图

1. 定义

二分图又称作二部图，是图论中的一种特殊模型。 设 G =（V,E）是一个无向图，如果顶点V可分割为两个互不相交的子集（A,B），并且图中的每条边（i，j）所关联的两个顶点i和j分别属于这两个不同的顶点集（i in A ，j in B），则称图G为一个二分图。

无向图G为二分图的充分必要条件是：G至少有两个顶点，且其所有回路的长度均为偶数。
在这里插入图片描述

2. 最大匹配

求二分图最大匹配可以用最大流或者匈牙利算法。

给定一个二分图G，在G的一个子图M中，M的边集中的任意两条边都不依附于同一个顶点，则称M是一个匹配。
选择这样的边数最大的子集称为图的最大匹配问题（maximal matching problem)
如果一个匹配中，图中的每个顶点都和图中某条边相关联，则称此匹配为完全匹配，也称作完备匹配。

增广路径

若P是图G中一条连通两个未匹配顶点的路径，并且属M的边和不属M的边（即已匹配和待匹配的边）在P上交替出现，则称P为相对于M的一条增广路径。
（交错路径：给定图G的一个匹配M，如果一条路径的边交替出现在M中和不出现在M中，我们称之为一条M-交错路径。而如果一条M-交错路径，它的两个端点都不与M中的边关联，我们称这条路径叫做M-增广路径。）

由增广路的定义可以推出下述三个结论：
1-P的路径长度必定为奇数，第一条边和最后一条边都不属于M。
2-P经过取反操作可以得到一个更大的匹配M '。
3-M为G的最大匹配当且仅当不存在相对于M的增广路径。
在这里插入图片描述
如图 2，4是属于M的已匹配的边，而1，3，5也是一个匹配且匹配数比原匹配数多1，是一个比原来匹配数更多的一个匹配。

所以，我们寻找最大匹配的任务就相当于我们不断地在已经确定的匹配下，不断找到新的增广路径，因为出现一条增广路径，就意味着目前的匹配中增加一条边。

3. 匈牙利算法

匈牙利算法，就是通过不断寻找增广路径的办法。

理论理解部分可参考这篇文章：匈牙利算法.

一句话概括就是：有机会就上，没机会就创造机会上。（创造增广路径）

模板：

X（n），Y（m）

int g[N][N],used[N],mh[N],n,m;
int Find(int x)
{
    for(int i=1; i<=m; i++)
    {
        if(g[x][i] && !used[i])
        {
            used[i] = 1;
            if(!mh[i] || Find(mh[i]))
            {
                mh[i] = x;
                return 1;
            }
        }
    }
    return 0;
}

int match()
{
    int sum = 0;
    memset(mh,0,sizeof(mh));
    for(int i=1; i<=n; i++)
    {
        memset(used,0,sizeof(used));
        if(Find(i))
            sum++;
    }
    return sum;
}
HK算法（Hopcroft Karp算法）

对匈牙利算法的优化 时间复杂度可达 O（n^0.5*m）。

4. 二分图的判断

常见方法是染色法： 开始对任意一未染色的顶点染色，之后判断其相邻的顶点中，若未染色则将其染上和相邻顶点不同的颜色， 若已经染色且颜色和相邻顶点的颜色相同则说明不是二分图，若颜色不同则继续判断，bfs和dfs可以搞定！

例子可看 B - The Accomodation of Students



实战应用

[kuangbin带你飞]专题十 匹配问题

A - Fire Net

A - Fire Net
额这题我直接用dfs过了，有空再加一下二分图的写法哈哈

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <sstream>
#define N 5000005
#define INF 0x3f3f3f3f
#define mod 1000
using namespace std;
typedef long long ll;
using namespace std;

char m[5][5];
int r[5],c[5],ans;
int n;
void dfs(int id,int sum)
{
    if(id == n*n)
    {
        if(sum > ans)
            ans = sum;
        return;
    }
    int x=id/n, y=id%n;
    if(m[x][y] == 'X')  r[x]=c[y]=0;
    if(m[x][y] == '.' && !r[x] && !c[y])
    {
        m[x][y] = 'Y';
        r[x]=c[y]=1;
        dfs(id+1,sum+1);
        m[x][y] = '.';
        r[x]=c[y]=0;
    }
    dfs(id+1,sum);
}


int main()
{
    while(~scanf("%d",&n) && n)
    {
        ans = 0;
        for(int i=0; i<n; i++)
            scanf("%s",m[i]);
        dfs(0,0);
        printf("%d\n",ans);
    }

    return 0;
}

B - The Accomodation of Students

B - The Accomodation of Students
题意：给了n个点，m个边，问能不能构成二分图，如果可以输出最大匹配数
大致思路：先用染色法判断是不是二分图，然后再用匈牙利算法求最大匹配数。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <vector>
#include <sstream>
#define N 250
#define INF 0x3f3f3f3f
#define mod 1000
using namespace std;
typedef long long ll;
using namespace std;

vector <int> ve[N];
int color[N];
int mat[N],used[N];
int n,m;

bool dfs(int x,int c) // 染色
{
    color[x]=c;
    int s=ve[x].size();
    for(int i=0; i<s; i++)
    {
        if(color[ve[x][i]] == c) // 连接相同颜色则不是二分图
            return 0;
        if(color[ve[x][i]] == 0 && !dfs(ve[x][i],-c))//如果相连点还没染色，搜索相连点
            return 0;
    }
    return 1;
}

bool solve()
{
    memset(color,0,sizeof(color));
    for(int i=1; i<=n; i++)
    {
        if(!color[i])
        {
            if(!dfs(i,1))
                return 0;
        }
    }
    return 1;
}

bool Find(int x)
{
    int s=ve[x].size();
    for(int i=0; i<s; i++)
    {
        int v=ve[x][i];
        if(!used[v])
        {
            used[v]=1;
            if(!mat[v] || Find(mat[v]))
            {
                mat[v] = x;
                return 1;
            }
        }
    }
    return 0;
}

int match()
{
    int ans=0;
    memset(mat,0,sizeof(mat));
    for(int i=1; i<=n; i++)
    {
        memset(used,0,sizeof(used));
        if(Find(i))
            ans++;
    }
    return ans;
}

int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        int u,v;
        for(int i=1; i<=n; i++)
            ve[i].clear();
        for(int i=0; i<m; i++)
        {
            scanf("%d%d",&u,&v);
            ve[u].push_back(v);
            ve[v].push_back(u);
        }
        if(!solve())//看能不能构成二分图
            printf("No\n");
        else
            printf("%d\n",match()/2);
    }

    return 0;
}

C - Courses

C - Courses.
模板题
学生是一个集合，课程是一个集合，用匈牙利算法求出最大匹配，
如果最大匹配等于课程数就YES，否则NO

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <vector>
#include <sstream>
#define N 305
#define INF 0x3f3f3f3f
#define mod 1000
using namespace std;
typedef long long ll;

int g[N][N],used[N],mh[N],n,m;

void Init()
{
    for(int i=1; i<=n ;i++)
        used[i]=mh[i]=0;
    memset(g,0,sizeof(g));
}

int Find(int x)
{
    for(int i=1; i<=n; i++)
    {
        if(g[x][i] && !used[i])
        {
            used[i] = 1;
            if(!mh[i] || Find(mh[i]))
            {
                mh[i] = x;
                return 1;
            }
        }
    }
    return 0;
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&m,&n);
        Init();
        for(int i=1; i<=m; i++)
        {
            int cnt;
            scanf("%d",&cnt);
            while(cnt--)
            {
                int x;
                scanf("%d",&x);
                g[i][x] = 1;
            }
        }
        int sum = 0;
        for(int i=1; i<=m; i++)
        {
            memset(used,0,sizeof(used));
            sum += Find(i);
        }
        if(sum == m)
            printf("YES\n");
        else
            printf("NO\n");
    }

    return 0;
}

D - 棋盘游戏

D - 棋盘游戏.
题意：
x,y的匹配就是"车"的摆放位置，而求最多的摆放就是求最大匹配数，求出会改变最大匹配数的点个数（重要点）
思路：
先把把x和y轴分别看作独立的集合用匈牙利算法得出最大匹配数，
然后再逐个去掉已经匹配好的边，
看再次求出的最大匹配数与原来是否一样，一样则不是重要点，不一样则是重要点。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <vector>
#include <sstream>
#define N 105
#define NN 10005
#define INF 0x3f3f3f3f
#define mod 1000
using namespace std;
typedef long long ll;

int g[N][N],used[N],mh[N],n,m;
int X[NN],Y[NN];
int Find(int x)
{
    for(int i=1; i<=m; i++)
    {
        if(g[x][i] && !used[i])
        {
            used[i] = 1;
            if(!mh[i] || Find(mh[i]))
            {
                mh[i] = x;
                return 1;
            }
        }
    }
    return 0;
}

int match()
{
    int sum = 0;
    memset(mh,0,sizeof(mh));
    for(int i=1; i<=n; i++)
    {
        memset(used,0,sizeof(used));
        if(Find(i))
            sum++;
    }
    return sum;
}

int main()
{
    int k,Case=1;
    while(~scanf("%d%d%d",&n,&m,&k))
    {
        memset(g,0,sizeof(g));
        for(int i=1; i<=k; i++)
        {
            scanf("%d%d",&X[i],&Y[i]);
            g[X[i]][Y[i]] = 1;
        }
        int ans = match(), imp=0;
        for(int i=1; i<=k; i++)
        {
            if(g[X[i]][Y[i]])
            {
                g[X[i]][Y[i]] = 0;
                if(ans>match())
                    imp++;
                g[X[i]][Y[i]] = 1;
            }
        }

        printf("Board %d have %d important blanks for %d chessmen.\n",Case++,imp,ans);
    }

    return 0;
}

E - Swap

E - Swap
题意：
题意很简单，让你把一个只含0或1的矩阵进行 行或着列 的交换（如果能换成，那么一行里至少有一列有1，所以交换行列等价，所以简化了题目），使得对角线上全是1。
思路：
在循环中找到需要交换的列，然后再交换匹配，使得其在对角线上。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <vector>
#include <sstream>
#include <utility>
#define N 105
#define NN 10005
#define INF 0x3f3f3f3f
#define mod 1000
using namespace std;
typedef long long ll;

int g[N][N],used[N],mh[N],n;
pair <int,int> a[1005];

int Find(int x)
{
    for(int i=1; i<=n; i++)
    {
        if(g[x][i] && !used[i])
        {
            used[i] = 1;
            if(!mh[i] || Find(mh[i]))
            {
                mh[i] = x;
                return 1;
            }
        }
    }
    return 0;
}

int match()
{
    int sum = 0;
    memset(mh,0,sizeof(mh));
    for(int i=1; i<=n; i++)
    {
        memset(used,0,sizeof(used));
        sum += Find(i);
    }
    return sum;
}

int main()
{
    while(~scanf("%d",&n))
    {
        memset(g,0,sizeof(g));
        for(int i=1; i<=n; i++)
            for(int j=1; j<=n; j++)
            {
                scanf("%d",&g[i][j]);
            }
        if(n > match())
            printf("-1\n");
        else
        {
            int cnt=0;
            for(int i=1; i<=n; i++)
            {
                int p;
                for(int j=1; j<=n; j++)
                {
                    if(mh[j] == i)
                    {
                        p = j;
                        break;
                    }
                }
                if(p != i)
                {
                    a[cnt++] = {i,p};
                    swap(mh[i],mh[p]);
                }
            }
            printf("%d\n",cnt);
            for(int i=0; i<cnt; i++)
                printf("C %d %d\n",a[i].first,a[i].second);
        }

    }

    return 0;
}

To be continued …

总结

Royen_
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
概率论经典问题之匹配问题
pdcxs007的专栏
10-23 3万+
一个屋子里面有N个人，每个人有一顶帽子。假如所有人把帽子扔到屋子中央，然后每个人都随机选一顶帽子。 a) 没有人捡到自己帽子的概率； b) 有k(k⩽N)k(k \leqslant N)个人捡到自己帽子的概率。
二分图与匹配问题
weixin_42374938的博客
09-01 1217
一、二分图的定义 给定图G=(V,E)。如果顶点集V可分为两个互不相交的非空子集X和Y，并且图中的每条边（i，j）的两个端点i和j分别属于这两个不同的顶点子集，那么就称图G为一个二分图，记为G=(X,E,Y)。如果在子集X中的任一顶点i和子集Y中的任一顶点j之间都存在一条边,那么就称图G为一个完全二分图。如下图所示： 在网络科学研究中，二分图也称为二分网络、从属网络和二模网络。 二、二分图的判断方法 在网络中任意选取一点，把该节点以及与该节点距离为偶数的所有其他节点的集合记为X，把与该节点的距离为奇数的节
二分图匹配问题
__surtr的博客
09-17 359
二分图匹配问题 定义 ：无向图中n个节点分成A，B两个非空集合， A∩B = ∅，并且同一集合内点互不相连， 称此无向图为二分图，A中的点称左部点， B中的点称右部点。 二分图的判定：当且仅当图中不存在长度为奇数的环。 染色法： void dfs(int x,int color){ v[x] = color; for(int i = head[x]; i; i = Next[i]){ in...
BZOJ1305 [CQOI2009]dance跳舞【二分图匹配 / 网络流】
qq_43710979的博客
10-11 129
[CQOI2009]dance跳舞[CQOI2009]dance跳舞[CQOI2009]dance跳舞 Description: 一次舞会有n个男孩和n个女孩。每首曲子开始时，所有男孩和女孩恰好配成n对跳交谊舞。每个男孩都不会和同一个女孩跳两首（或更多）舞曲。有一些男孩女孩相互喜欢，而其他相互不喜欢（不会“单向喜欢”）。每个男孩最多只愿意和k个不喜欢的女孩跳舞，而每个女孩也最多只愿意和k个不喜欢...
linux服务器性能调优文档（未完待续）
07-23
教程名称：linux服务器性能调优文档(未完待续)课程目录：linux服务器性能调优01linux服务器性能调优02linux服务器性能调优03linux服务器性能调优04linux服务器性能调优05linux服务器性能调优06linux服务器性能调优07...
2015年高中语文优秀作文青春未完待续…素材
08-06
10. 青春的延续：尽管青春已逝，但它的精神和影响却可以"未完待续"，持续影响着我们的生活和决策。 总结，本文通过对青春的描绘，探讨了时间、成长、分别、记忆、诺言等主题，展现了青春的美好、脆弱和短暂，同时也...
Visual Studio问题集锦by低调年华,未完待续.docx
11-02
【Visual Studio问题集锦】 在使用Visual Studio进行软件开发时，可能会遇到各种各样的问题，以下是一些常见的问题及其解决方案： 1. **安装问题**：VS2013的安装需要网络接入，否则可能导致部分功能无法正常使用...
Visual Studio问题集锦by低调年华,未完待续.pdf
10-27
**Visual Studio问题集锦** 在使用Visual Studio进行开发时，我们可能会遇到各种各样的问题，以下是一些常见的问题及其解决方案： 1. **错误MSB 3073：已退出，代码1** 这个错误通常出现在项目的生成事件中。解决...
quasar_clone:未完待续
03-17
Matthias Mann的续集库，为 。 原始文档如下。马蒂亚斯·曼（Matthias Mann）的续集图书馆这个小的库允许您以一种简单的顺序方式来编写游戏AI或动画代码，从而阻塞诸如walkTo()类的命令，而所有这些都不会阻塞您的AI...
代码随想录算法训练营第七天（一）| 454.四数相加II 383. 赎金信
xiaowutongxue666的博客
08-06 511
给你四个整数数组nums1nums2nums3和nums4，数组长度都是n，请你计算有多少个元组2两个元组如下：1。
P1145 约瑟夫
zty20120913的博客
08-06 303
n个人站成一圈，从某个人开始数数，每次数到m的人就被杀掉，然后下一个人重新开始数，直到最后只剩一个人。现在有一圈人，k个好人站在一起，k个坏人站在一起。从第一个好人开始数数。你要确定一个最小的m，使得在第一个好人被杀死前，k个坏人先被杀死。
5力扣：最长回文子串
m0_75015083的博客
08-02 323
从后到前遍历是为了让每一次遍历时所需用到的子串提前被计算过。
VRF 算法流程及代码实现
smilejiasmile的博客
08-03 162
Verifiable Random Functions (VRFs) 是一种密码学原语，能够为给定输入生成一个随机输出，并提供一个可以验证的证明，证明输出是由特定输入生成的。VRFs 具有许多应用，包括区块链、分布式系统和安全协议。使用secp256k1曲线可以产生确定性签名，适用于 VRF 实现。通过使用 HMAC-SHA256 生成确定性随机数k，可以确保签名过程对相同的输入始终产生相同的结果，从而满足 VRF 的要求。
day021二叉树最后章节
最新发布
m0_62956971的博客
08-06 444
思路: 二叉搜索树的中序遍历就是有序升序数组，所以想要将它在转换成二叉树，我想到的方法就是通过递归来做，通过中间节点找到二叉搜索树的根节点，然后根据left和right双指针逐层递归构造左右子树。思路:最开始没有判断当当前节点小于最小值时，当前节点的右侧节点的情况，和当前节点大于最大值时，当前节点的左节点的情况，所以出错了，因此要在代码中判断一下这两种情况。思路: 累加树 由题意可知，在二叉搜索树中的递归顺序是右中左的顺序，创建一个新的函数，进行右中左的遍历，在遍历中要将节点的值记录下来，用来做累加。
算法第十七天：leetcode242.有效的字母异位词
m0_75068978的博客
08-02 547
如下表所示，您可以直接复制下面网址进入力扣学习，在观看下面的内容之前您一定要先做一遍哦，以便让你印象更加深刻！
CSP - J day7
2301_81608637的博客
08-06 144
答案具有单调性：对于正确答案 ans ，小于 ans 的都一定符合条件，大于 ans 的一定不符合条件，因此可以二分答案。先dfs到叶节点，再逐步递归出来，过程中节点数增加，一旦达到 mid ，ans 加一，同时向上贡献节点数变为零。
埃拉托斯特尼筛法，厄拉多塞筛算法，C#实现
GoodCooking的博客
08-06 251
埃拉托斯特尼筛法，厄拉多塞筛算法，C#实现
代码随想录算法训练营第三十三天 | 动态规划理论基础、509. 斐波那契数、70. 爬楼梯、746. 使用最小花费爬楼梯
weixin_43581681的博客
08-01 323
1. 动态规划基础；2. 背包问题；3. 打家劫舍；4. 股票问题；5. 子序列问题。1. 确定 dp[ i ] 数组及下标含义；2. 确定递推公式；3.确定dp数组如何初始化；4. 确定遍历顺序（背包问题中很重要）；5. 举例推导dp数组。
"PMPBOK第六版自整理（未完待续）：项目与运营引论
关于PMP（项目管理专业人士）和PMBOK（项目管理知识体系）的相关介绍和概述，以下是针对这些内容进行整理和总结的一些要点。 PMP（项目管理专业人士）是由世界上最知名的国际管理学会之一——美国项目管理学会...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值