学习匈牙利算法总结（求解二分图最大匹配）

学习一下最基础的匈牙利算法。

感谢参考：

http://blog.csdn.net/china8848/article/details/2287769

http://blog.csdn.net/akof1314/article/details/4421262

http://baike.baidu.com/view/501092.htm?fr=aladdin

http://blog.csdn.net/dark_scope/article/details/8880547

    匈牙利算法代码并不复杂，但是其中算法蕴含的思想理解是有难度的，其思想的理解对网络流算法的学习有很大帮助。

    二分图是这样一个图，它的顶点可以分类两个集合X和Y，所有的边关联在两个顶点中，恰好一个属于集合Ｘ，另一个属于集合Ｙ。给定一个二分图G，M为G边集的一个子集，如果M满足当中的任意两条边都不依附于同一个顶点，则称M是一个匹配。图中包含边数最多的匹配称为图的最大匹配。二分图的最大匹配有两种求法，第一种是最大流；第二种就是我现在要讲的匈牙利算法。这个算法说白了就是最大流的算法，但是它跟据二分图匹配这个问题的特点，把最大流算法做了简化，提高了效率。

    概念：

未盖点：设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的对应项出有可增广路)匹配数++;}输出 匹配数;}

    用通俗的话描述匈牙利算法就是：不停的找增广路径,并增加匹配的个数。

    增广路径顾名思义是指一条可以使匹配数变多的路径,在匹配问题中,增广路径的表现形式是一条"交错路径",也就是说这条由图的边组成的路径,它的第一条边是目前还没有参与匹配的,第二条边参与了匹配,第三条边没有..最后一条边没有参与匹配,并且始点和终点还没有被选择过.这样交错进行,显然他有奇数条边.那么对于这样一条路径,我们可以将第一条边改为已匹配,第二条边改为未匹配...以此类推.也就是将所有的边进行"反转",容易发现这样修改以后,匹配仍然是合法的,但是匹配数增加了一对.另外,单独的一条连接两个未匹配点的边显然也是交错轨.可以证明,当不能再找到增广轨时,就得到了一个最大匹配.这也就是匈牙利算法的思路.

#include<stdio.h>#include<string.h>int n1,n2,m,ans;int result[101];//记录V2中的点匹配的点的编号bool state[101];//记录V2中的每个点是否被搜索过bool data[101][101];//邻接矩阵true代表有边相连void init(){int t1,t2;memset(data,0,sizeof(data));memset(result,0,sizeof(result));ans=0;scanf("%d%d%d",&n1,&n2,&m);for(int i=1;i<=m;i++){scanf("%d%d",&t1,&t2);data[t1][t2]=true;}return;}bool find(int a){for(int i=1;i<=n2;i++){if(data[a][i]==1&&!state[i]){//如果节点i与a相邻并且未被查找过state[i]=true;//标记i为已查找过if(result[i]==0//如果i未在前一个匹配M中||find(result[i])){//i在匹配M中，但是从与i相邻的节点出发可以有增广路result[i]=a;//记录查找成功记录return true;//返回查找成功}}}return false;}int main(){init();for(int i=1;i<=n1;i++){memset(state,0,sizeof(state));//清空上次搜索时的标记if(find(i))ans++;//从节点i尝试扩展}printf("%d\n",ans);return 0;}

    程序默认有向图的最大匹配。程序分析：

    result[]表示搜索i的增广路径之后，得到的匹配结果。每次搜索i的增广路径，都要将state[]数组初始化为全0，搜索过后不就用搜索的原理：如果从一个点A出发，没有找到增广路径，那么无论再从别的点出发找到多少增广路径来改变现在的匹配，从A出发都永远找不到增广路径。如果找到增广路径皆大欢喜，返回ture。