NOCOW-匈牙利算法.pdfVIP

匈牙利算法 求最大匹配的一种显而易见的算法是：先找出全部匹配，然后保留匹配数最 多的。但是这个算法的复杂度为边数的指数级函数。因此，需要寻求一种更 加高效的算法。 增广路的定义(也称增广轨或交错轨)：若P 是图G 中一条连通两个未匹配顶 点的路径，并且属M 的边和不属M 的边(即已匹配和待匹配的边)在P 上交替 出现，则称P 为相对于M 的一条增广路径。（M 为一个匹配） 由增广路的定义可以推出下述三个结论： 1. P 的路径长度必定为奇数，第一条边和最后一条边都不属于M。 2. P 经过取反操作可以得到一个更大的匹配M’。 3. M 为G 的最大匹配当且仅当不存在相对于M 的增广路径。 用增广路求最大匹配(称作匈牙利算法，匈牙利数学家Edmonds 于1965 年 提出) 算法轮廓： 1. 置M 为空 2. 找出一条增广路径P，通过取反操作获得更大的匹配M’代替M 3. 重复(2)操作直到找不出增广路径为止 程序清单： parsical const maxm=200; maxn=200; var i,j,k,m,n:longint; g:array[1..maxm,1..maxn]of boolean ; y:array[1..maxn]of boolean ; lk:array[1..maxn]of longint; function find(x:longint):boolean ; var i:longint; begin for i:=1 to n do if g[x,i] and (not y[i]) then begin y[i]:=true; if (lk[i]=0)or find (lk[i]) then begin lk[i]:=x; lk[x]:=i; find:=true; exit; end; end; find:=false ; end; C++ #includestdio.h #includestring.h bool g [201][201]; int n,m,ans; bool b [201]; int link[201]; bool init () { int _x,_y; memset (g,0,sizeof(g)); memset (link,0,sizeof(link)); ans=0; if(scanf(%d%d,n,m)==EOF)return false ; for(int i=1;i=n;i++) { scanf(%d,_x); for(int j=0;j_x;j++) { scanf(%d,_y); g[ i ][_y]=true; } } return true; } bool find (int a) { for(int i=1;i=m;i++) { if(g[a][ i ]==1!b [ i ]) { b [ i ]=true; if(link[ i ]==0 ||find(link[ i ])) {