匈牙利算法

匈牙利算法

全栈程序员-站长 • 2021年12月15日上午11:00 • 未分类 • 阅读 52

大家好，又见面了，我是全栈君，今天给大家准备了Idea注册码。

此处内容已经被作者隐藏，请输入验证码查看内容

请关注本站微信公众号，回复“验证码”，获取验证码。在微信里搜索“全栈程序员社区”或者“www_javaforall_cn”或者微信扫描右侧二维码都可以关注本站微信公众号。

本文转自大牛博客：http://www.byvoid.com/blog/hungary/

这是一种用增广路求二分图最大匹配的算法。它由匈牙利数学家Edmonds于1965年提出，因而得名。定义未盖点：设Vi是图G的一个顶点，假设Vi 不与随意一条属于匹配M的边相关联，就称Vi 是一个未盖点。

匈牙利算法

交错路：设P是图G的一条路，假设P的随意两条相邻的边一定是一条属于M而还有一条不属于M，就称P是一条交错路。

可增广路：两个端点都是未盖点的交错路叫做可增广路。

匈牙利算法

流程图

伪代码：

bool 寻找从k出发的相应项出的可增广路 { while (从邻接表中列举k能关联到顶点j) { if (j不在增广路上) { 把j增加增广路; if (j是未盖点或者从j的相应项出发有可增广路) { 改动j的相应项为k; 则从k的相应项出有可增广路,返回true; } } } 则从k的相应项出没有可增广路,返回false; } void 匈牙利hungary() { for i->1 to n { if (则从i的相应项出有可增广路) 匹配数++; } 输出匹配数; }

演示：

匈牙利算法

C实现（作者BYVoid）

#include <stdio.h> #include <string.h> #define MAX 102 long n,n1,match; long adjl[MAX][MAX]; long mat[MAX]; bool used[MAX]; FILE *fi,*fo; void readfile() { fi=fopen(“flyer.in”,”r”); fo=fopen(“flyer.out”,”w”); fscanf(fi,”%ld%ld”,&n,&n1); long a,b; while (fscanf(fi,”%ld%ld”,&a,&b)!=EOF) adjl[a][ ++adjl[a][0] ]=b; match=0; } bool crosspath(long k) { for (long i=1;i<=adjl[k][0];i++) { long j=adjl[k][i]; if (!used[j]) { used[j]=true; if (mat[j]==0 || crosspath(mat[j])) { mat[j]=k; return true; } } } return false; } void hungary() { for (long i=1;i<=n1;i++) { if (crosspath(i)) match++; memset(used,0,sizeof(used)); } } void print() { fprintf(fo,”%ld”,match); fclose(fi); fclose(fo); } int main() { readfile(); hungary(); print(); return 0; }

Pascal实现（作者魂牛）

var a:array[1..1000,1..1000] of boolean; b:array[1..1000] of longint; c:array[1..1000] of boolean; n,k,i,x,y,ans,m:longint; function path(x:longint):boolean; var i:longint; begin for i:=1 to n do if a[x,i] and not c[i] then begin c[i]:=true; if (b[i]=0) or path(b[i]) then begin b[i]:=x; exit(true); end; end; exit(false); end; procedure hungary; var i:longint; begin fillchar(b,sizeof(b),0); for i:=1 to m do begin fillchar(c,sizeof(c),0); if path(i) then inc(ans); end; end; begin fillchar(a,sizeof(a),0); readln(m,n,k); for i:=1 to k do begin readln(x,y); a[x,y]:=true; end; ans:=0; hungary; writeln(ans); end.

ps:近期由于比赛，须要用到二分图，图论这东西好多要好好学的。

发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/117991.html原文链接：https://javaforall.net

相关推荐

icem网格划分如何给内部面网格,ICEM CFD处理混合网格划分中低质量的问题

fbx文件导入3dmax_3d中z轴的值没办法输入

机器学习模型中的损失函数loss function

host地址_hostname在哪个目录

常用数据库排名及分类介绍[通俗易懂]

外企入职第一封英文邮件_投外企要英文简历吗

发表回复