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

8 篇文章 0 订阅
1 篇文章 0 订阅

      今天介绍 匈牙利算法 : 匈牙利算法,是基于Hall定理中充分性证明的思想,它是部图匹配最常见的算法,该算法的核心就是寻找增广路径,由匈牙利数学家Edmonds于1965年提出,因而得名。

      先介绍一下增广路径:若P是图G中一条连通两个未匹配顶点的路径,并且属于M的边和不属于M的边(即已匹配和待匹配的边)在P上交替出现,则称P为相对于M的一条增广路径。文字难以理解,看图:

      首先 假设 图G中已经两两匹配了4个点 即 M 中 又 (1,A)(4,C)两条边;

      然后路径p是一条连接未匹配点2和D的一条路径 其中属于M的边和不属于M的边交错:(2-A-1-C-4-D),这样就称p为M的一条增广路径;

       增广路径的性质:

               1-P的路径长度必定为奇数,第一条边和最后一条边都不属于M。

               2-不断寻找增广路可以得到一个更大的匹配M’,直到找不到更多的增广路。

               3-M为G的最大匹配当且仅当不存在M的增广路径。

               4-最大匹配数M+最大独立数N=总的结点数。

      我们将P中不属于M的路径加入M,将原来属于M的路径去除。M的匹配数增加到了 3,如果我们一只找M的增广路径然后取反,直到M没有增广路,这时我们就找到了M的最大匹配数,这也时 匈牙利算法的核心所在;

      匈牙利算法 :首先向M中增加一条路径,然后一只寻找M的增广路径直到M中不存在增广路径为止,这时我们就完成了最大匹配的寻找。

      用一个问题来展现匈牙利算法:

      假设有一群人生了不同的病(一个人只能生一种病),这里有多种可以治病的药丸,并且一种药丸只能治一种病且仅有一个,问最多能救治多少人;

      输入 :在第一行中给出 n ,m, k 其中 n是人数,m 是 药的种类数,k是 人和药共有多少种匹配,然后 跟着 k 行,每行 一个 整数:人的标号,一个字符 :药的种类。

      输出:在一行中输出最多可以救治多少人。

      代码:

//二分图最大匹配-匈牙利算法
#include<stdio.h>
#include<string.h>
#define MAXN 100       //定义最大人数
#define TYPE 6        //定义药的种类数 
int map[MAXN][MAXN]; //用来储存人和药的关系  
int vis[MAXN];      //用来标记药的被使用人
int used[MAXN];	   //用在递归的时候标记修改的药的种类 
char array[6] = {'a','b','c','d','e','f'};
int GetPos(char x){ //获取药的索引 
	int i, position;
	for (i=0;i<TYPE;i++)
	{
		if (array[i] == x)
		return i; // 找到的话就返回索引 
	}
	return -1; // 找不到就返回-1 
}
void init(){ // 初始化函数 
	memset(map,-1,sizeof(map));
	memset(vis,-1,sizeof(vis));
	memset(used,-1,sizeof(used));
}
int MaxMatch(int n){ // n是人的代号 
	int i;
	for (i=0;i<TYPE;i++)
	{
		if ( used[i]==-1 && map[n][i]==1) //寻找增广路 
		{
			used[i] = 1;   
			if (vis[i] == -1||MaxMatch(vis[i])==1)
			{
				vis[i] = n;
				return 1;
			}
		}
	}
	return -1;
}
int main (){
	int n, m, i, k, cnt = 0, b, index;//n是患病的人数,m是药的种类数,i是循环变量,k是人和药对应关系的数量,cnt记录最大匹配数,b是某人的代号,index记录药的位置 
	char a; //a是药的品种 
	scanf("%d%d%d",&n,&m,&k);
	init();
	for (i=0;i<k;i++)
	{
		scanf("%d %c",&b,&a);
		index = GetPos(a);
		map[b][index] = 1;
	}
	for (i=1;i<=n;i++)
	{
		memset(used,-1,sizeof(used));//每次都要初始化一遍 每次寻找增广路都与之前无关 
		if (MaxMatch(i)==1)
		cnt++;
	}
	printf("%d\n",cnt);
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值