图论二分图匈牙利算法求最大匹配

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

 

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

 

交错路：设P是图G的一条路，如果P的任意两条相邻的边一定是一条属于M而另一条不属于M，就称P是一条交错路。

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

 

 

流程图

 

伪代码：

[cpp]  view plain copy

1. bool 寻找从k出发的对应项出的可增广路  
2. {  
3.     while (从邻接表中列举k能关联到顶点j)  
4.     {  
5.         if (j不在增广路上)  
6.         {  
7.             把j加入增广路;  
8.             if (j是未盖点 或者 从j的对应项出发有可增广路)  
9.             {  
10.                 修改j的对应项为k;  
11.                 则从k的对应项出有可增广路,返回true;  
12.             }  
13.         }  
14.     }  
15.     则从k的对应项出没有可增广路,返回false;  
16. }  
17.    
18. void 匈牙利hungary()  
19. {  
20.     for i->1 to n  
21.     {  
22.         if (则从i的对应项出有可增广路)  
23.             匹配数++;  
24.     }  
25.     输出 匹配数;  
26. }  

 

演示：

 

C实现（作者BYVoid）

[cpp]  view plain copy

1. #include <stdio.h>  
2. #include <string.h>  
3. #define MAX 102  
4.    
5. long n,n1,match;  
6. long adjl[MAX][MAX];  
7. long mat[MAX];  
8. bool used[MAX];  
9.    
10. FILE *fi,*fo;  
11.    
12. void readfile()  
13. {  
14.     fi=fopen("flyer.in","r");  
15.     fo=fopen("flyer.out","w");  
16.     fscanf(fi,"%ld%ld",&n,&n1);  
17.     long a,b;  
18.     while (fscanf(fi,"%ld%ld",&a,&b)!=EOF)  
19.         adjl[a][ ++adjl[a][0] ]=b;  
20.     match=0;  
21. }  
22.    
23. bool crosspath(long k)  
24. {  
25.     for (long i=1;i<=adjl[k][0];i++)  
26.     {  
27.         long j=adjl[k][i];  
28.         if (!used[j])  
29.         {  
30.             used[j]=true;  
31.             if (mat[j]==0 || crosspath(mat[j]))  
32.             {  
33.                 mat[j]=k;  
34.                 return true;  
35.             }  
36.         }  
37.     }  
38.     return false;  
39. }  
40.    
41. void hungary()  
42. {  
43.     for (long i=1;i<=n1;i++)  
44.     {  
45.         if (crosspath(i))  
46.             match++;  
47.         memset(used,0,sizeof(used));  
48.     }  
49. }  
50.    
51. void print()  
52. {  
53.     fprintf(fo,"%ld",match);  
54.     fclose(fi);  
55.     fclose(fo);  
56. }  
57.    
58. int main()  
59. {  
60.     readfile();  
61.     hungary();  
62.     print();  
63.     return 0;  
64. }  

Pascal实现（作者魂牛）

[delphi]  view plain copy

1. var  
2.   a:array[1..1000,1..1000] of boolean;  
3.   b:array[1..1000] of longint;  
4.   c:array[1..1000] of boolean;  
5.   n,k,i,x,y,ans,m:longint;  
6.    
7. function path(x:longint):boolean;  
8. var  
9.   i:longint;  
10. begin  
11.   for i:=1 to n do  
12.   if a[x,i] and not c[i] then  
13.   begin  
14.     c[i]:=true;  
15.     if (b[i]=0) or path(b[i]) then  
16.     begin  
17.       b[i]:=x;  
18.       exit(true);  
19.     end;  
20.   end;  
21.   exit(false);  
22. end;  
23.    
24. procedure hungary;  
25. var  
26.   i:longint;  
27. begin  
28.   fillchar(b,sizeof(b),0);  
29.   for i:=1 to m do  
30.   begin  
31.     fillchar(c,sizeof(c),0);  
32.     if path(i) then inc(ans);  
33.   end;  
34. end;  
35.    
36. begin  
37.   fillchar(a,sizeof(a),0);  
38.   readln(m,n,k);  
39.   for i:=1 to k do  
40.   begin  
41.     readln(x,y);  
42.     a[x,y]:=true;  
43.   end;  
44.   ans:=0;  
45.   hungary;  
46.   writeln(ans);  
47. end.