匈牙利算法 最大匹配 Matlab实现

二分图匹配

Matlab实现：

% Hungarian_2.m
% 二分图矩阵表示
map=[0 1 1 0 0
     1 1 0 1 1
     0 1 1 0 0
     0 1 1 0 0
     0 0 0 1 1];
m=size(map,1);%行数
n=size(map,2);%列数

s=0;%匹配个数
% cx=zeros(1,m);%集合X
cy=zeros(1,n);%集合Y 记录集合Y中元素所对应的X元素

for i=1:m
        vis = zeros(1,m);
        [flag,cy] = match(i,map,vis,cy);
       	s = s + flag;
end


function [flag,cy] = match(i,map,vis,cy)
n = size(map,2);
for j= 1:n
    if map(i,j)&&~vis(j) % 有边且未访问 此处只考虑i对j
        vis(j) = 1; % 记录状态为访问过
        if cy(j)==0 || match(cy(j),map,vis,cy) % 如果暂无匹配，或者原来匹配的集合X中的元素可以找到新的匹配
            cy(j) = i;
            flag = 1; % 返回匹配成功
            return;
        end
    end
end
flag = 0; % 循环结束，未找到匹配，返回匹配失败
end

参考资料：

1、算法学习笔记(5)：匈牙利算法https://zhuanlan.zhihu.com/p/96229700
2、Matlab实现二分图的最大匹配https://blog.csdn.net/zyxhangiian123456789/article/details/83301541