匈牙利算法
它是二分图匹配最常见的算法,该算法的核心就是寻找增广路径,它是一种用增广路径求二分图最大匹配的算法。这篇文章讲无权二分图(unweighted
bipartite graph)的最大匹配(maximum matching)和完美匹配( perfect matching),以及用于求解匹配的匈牙利算法(Hungarian Algorithm)
相关知识介绍:
二分图:简单来说,如果图中点可以被分为两组,并且使得所有边都跨越组的边界,即没有边连接组内的点,则这就是一个二分图。
准确地说:把一个图的顶点划分为两个不相交集U和V,使得每一条边都分别连接U、V中的顶点。如果存在这样的划分,则此图为一个二分图。
二分图的一个等价定义是:不含有「含奇数条边的环」的图。图1是一个二
分图。为了清晰,我们以后都把它画成图2的形式。
匹配:
在图论中,一个「 匹配」(matching) 是一-个边的集合,其中任意两条边都没有公共顶点。例如,图3、图4中红色的边就是图2的匹配。
我们定义匹配点、匹配边、未匹配点、非匹配边,它们的含义非常显然。例如图3中1、4、5、7为匹配点,其他顶点为未匹配点; 1-5、 4-7为匹配边,其他边为非匹配边。
最大匹配:
一个图所有匹配中,所含匹配边数最多的匹配,称为这个图的最大匹配。图4是一个最大匹配,它包含4条匹配边。
完美匹配:
如果一个图的某个匹配中,所有的顶点都是匹配点,那么它就是一个完美匹配。图4是一个完美匹配。显然,完美匹配一定是最大匹配(完美匹配的任何一个点都已经匹配,添加一条新的匹配边一定会与已有的匹配边冲突)。但并非每个图都存在完美匹配。
最大匹配数:最大匹配的匹配边的数目;
最小点覆盖数:选取最少的点,使任意一条边至少有一个端点被选择;
最大独立数:选取最多的点,使任意所选两点均不相连;
最小路径覆盖数:对于一个DAG (有向无环图),选取最少条路径,使得
每个顶点属于且仅属于一条路径。路径长可以为0 (即单个点)。
定理1:最大匹配数=最小点覆盖数(这是Konig定理)
定理2:最大匹配数=最大独立数
定理3:最小路径覆盖数=顶点数-最大匹配数
例题:
Problem Description
RPG
girls今天和大家一起去游乐场玩,终于可以坐上梦寐以求的过山车了。可是,过山车的每一排只有两个座位,而且还有条不成文的规矩,就是每个女生必须找个个男生做partner和她同坐。但是,每个女孩都有各自的想法,举个例子把,Rabbit只愿意和XHD或PQK做partner,Grass只愿意和linle或LL做partner,PrincessSnow愿意和水域浪子或伪酷儿做partner。考虑到经费问题,boss刘决定只让找到partner的人去坐过山车,其他的人,嘿嘿,就站在下面看着吧。聪明的Acmer,你可以帮忙算算最多有多少对组合可以坐上过山车吗?
Input
输入数据的第一行是三个整数K , M , N,分别表示可能的组合数目,女生的人数,男生的人数。0<K<=1000 1<=N和M<=500.接下来的K行,每行有两个数,分别表示女生Ai愿意和男生Bj做partner。最后一个0结束输入。
Output
对于每组数据,输出一个整数,表示可以坐上过山车的最多组合数。
Sample Input
6 3 3
1 1
1 2
1 3
2 1
2 3
3 1
0
Sample Output
3
思路:
这就是一个最简单的二分图的最大匹配问题,主要根据此问题分析匈牙利算法。
题解:
#include<bits/stdc++.h>
using namespace std;
#define N 5005
bool used[N];//记录女生是否有男生选中
bool match[N][N];//记录男女之间是否有连线
int girl[N];//记录女生选择那个男生
bool find(int x,int n){
for(int j=1;j<=n;j++){//在n个女生中依次遍历
if(match[x][j] && !used[j]){// x男生与j女生之间有连线,并且j女生未被选择过
used[j] = true;
if(girl[j]==0||find(girl[j],n)){//如果j女生还没有匹配的男生,或者j女生匹配的男生可以在n女生中找到另一个
girl[j] = x;
return true;
}
}
}
return false;
}
int slove(int m,int n){ //匹配过程
int ans = 0;
for(int i=1;i<=m;i++){
memset(used,false,sizeof(used));//每判断一个男生的时候,开始这个男生没有选择的女生,所以该数组要初始化
if(find(i,n)) ans++; //在n个女生中找到与i男生配对的,
}
return ans;
}
int main(){
int k,m,n;
while(cin >> k&&k){ //循环输入
cin >> m >> n;
memset(match,false,sizeof(match));//初始化
memset(girl,0,sizeof(girl));
for(int i=1;i<=k;i++){
int x,y;
cin >> x >> y;
match[x][y] = true;
}
cout << slove(m,n) << endl;
}
return 0;
}