图论---匈牙利算法求二分图最大匹配的实现

开始编程前分析设计思路和程序的整体的框架,以及作为数学问题的性质:

程序流程图:

数学原理:

求解二分图最大匹配问题的算法,寻找一个边的子集,使得每个左部点都与右部点相连,并且没有两条边共享相同的端点。使用深度优先搜索来寻找增广路径。当找到一条增广路径时,将路径上的边加入匹配集合中,并更新匹配状态。重复这个过程直到无法找到增广路径为止。

时间复杂度分析:

时间消耗来自于DFS搜索和增广路径的寻找。对于每个左侧未匹配的点u,需要进行一次DFS搜索。在最坏情况下,DFS搜索需要遍历所有的右侧点,因此时间复杂度为O(n)。由于需要进行n次DFS搜索,所以时间复杂度为O(n^2)。

源代码:

cpp 复制代码
#include <iostream>
#include <vector>
#include <cstring>
using namespace std;

const int MAXN = 505; // 最大顶点数
int n, m, k; // n为左部点数,m为右部点数,k为总变数
vector<int> G[MAXN]; // 邻接表表示二分图
int match[MAXN]; // match[i]表示右边第i个点匹配的是左边的哪点,如果没有匹配则为-1
bool vis[MAXN]; // DFS中标记右边各点是否已经被访问过
//进行深度搜索最大匹配边
bool dfs(int u) {
    int i;
    for (i = 0; i < G[u].size(); i++) {
        int v = G[u][i];
        if (!vis[v]) {
            vis[v] = true;
            if (match[v] == -1 || dfs(match[v])) {
                match[v] = u;//对边进行存储
                return true;
            }
        }
    }
    return false;
}
//匈牙利算法
int hungarian() {
    int res = 0;
    memset(match, -1, sizeof(match));
    for (int i = 1; i <= n; i++) {
        memset(vis, false, sizeof(vis));
        if (dfs(i)) res++;
    }
    return res;
}
int main() {
    cin >> n >> m >> k; // 输入左部点数和右部点数
    int u, v;
    int i=0;
    for(;i<k;i++){
        cin >> u >> v ; // 输入边 (u, v),u属于左部,v属于右部
        G[u].push_back(v);
    }
    int max_matches = hungarian(); // 计算最大匹配数
    cout << "最大匹配边数:" << max_matches << endl;
    // 输出最大匹配数
    // 输出最大匹配的所有边
    cout << "最大匹配边为:" << endl;
    for (int v = 1; v <= m; v++) {
        if (match[v] != -1) { // 如果右部点v有匹配
            cout << "Edge: " << match[v] << " - " << v <<endl; // 输出匹配边
        }
    }
    system("pause");
    return 0;
}

测试用例:(两侧顶点数均大于10,且两侧顶点数不相等)

|------------------------------------------------------------------------------------------------------------------------------------------------------|
| |

输出结果:

相关推荐
charlie1145141918 分钟前
嵌入式C++教程——ETL(Embedded Template Library)
开发语言·c++·笔记·学习·嵌入式·etl
陳10308 分钟前
C++:AVL树的模拟实现
开发语言·c++
CSDN_RTKLIB20 分钟前
错进错出得到正确的字节序列
c++
闻缺陷则喜何志丹43 分钟前
【前后缀分解 排序】B4274 [蓝桥杯青少年组省赛 2023] 数字游戏|普及+
c++·蓝桥杯·排序·洛谷·前后缀分解
m0_7369191043 分钟前
C++中的享元模式变体
开发语言·c++·算法
zho_uzhou2 小时前
c++ imgui implot绘图使用示例 visual studio
开发语言·c++·visual studio
dyyx1112 小时前
C++中的过滤器模式
开发语言·c++·算法
机器视觉知识推荐、就业指导2 小时前
用惯了QTimer定时器,如何快速在纯 C++ 项目中替换?
c++
消失的旧时光-19432 小时前
从拷贝到移动:C++ 移动构造与移动赋值是怎么被逼出来的?(附完整示例)
开发语言·c++
2301_817497332 小时前
C++中的装饰器模式高级应用
开发语言·c++·算法