1. 程式人生 > >匈牙利演算法——最大匹配

匈牙利演算法——最大匹配

問題簡介:

二分圖是指這樣的圖:圖的頂點分成X和Y兩個集合,同一個集合的頂點不存在邊相連,只有不同集合的頂點才有可能有邊直接相連。二分圖的一個匹配是指求出一個邊的集合,使得集合裡任意兩條邊都沒有公共的頂點,也就是說一個頂點最多隻屬於一條邊。

簡單地看,二分圖的一個匹配給出了一個X集合的某些頂點Y集合的某些頂點的一個一一對映。

二分圖的最大匹配(maximal matching)就是找出邊數最大邊集,也就是求最多有多少個頂點可以被匹配上。

如果一個匹配中,圖中的每個頂點都和圖中某條邊相關聯,則稱此匹配為完全匹配,也稱作完備,完美匹配。

演算法描述:

求最大匹配的一種顯而易見的演算法是:先找出全部匹配,然後保留匹配數最多的。但是這個演算法的時間複雜度為邊數的指數級函式。因此,需要尋求一種更加高效的演算法。下面介紹用增廣路

求最大匹配的方法(稱作匈牙利演算法,匈牙利數學家Edmonds於1965年提出)。

增廣路的定義(也稱增廣軌或交錯軌):

一條增廣路的概念是指找一條已經被匹配上的邊和未被匹配上的邊相互交替的路徑,這條路徑長度為奇數並且未被匹配上的邊比被匹配上的邊1。(因為,我要從左邊匹配到右邊)通過增廣路上未被匹配上的邊改為匹配上的邊,而匹配上的邊全部改為未被匹配上的邊,就可以把匹配數加1

如果這樣的增廣路不能被找到,演算法結束。

這是一種用增廣路求二分圖最大匹配的演算法。它由匈牙利數學家Edmonds於1965年提出,因而得名。 定義 未蓋點:設Vi是圖G的一個頂點,如果Vi 不與任意一條屬於匹配M的邊相關聯,就稱Vi 是一個未蓋點。

交錯路:設P是圖G的一條路,如果P的任意兩條相鄰的邊一定是一條屬於M而另一條不屬於M,就稱P是一條交錯路。

可增廣路:兩個端點都是未蓋點的交錯路叫做可增廣路。

流程圖

l

虛擬碼:

  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. void 匈牙利hungary()  
  18. {  
  19.     for i->1 to n  
  20.     {  
  21.         if (則從i的對應項出有可增廣路)  
  22.             匹配數++;  
  23.     }  
  24.     輸出 匹配數;  
  25. }  

演示:

1

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