匈牙利算法是由匈牙利数学家Edmonds于1965年提出,因而得名。匈牙利算法是基于Hall定理中充分性证明的思想,它是二部图匹配最常见的算法,该算法的核心就是寻找增广路径,它是一种用增广路径求二分图最大匹配的算法。
#include<stdio.h> #include<string.h> #include<stdlib.h> int n1,n2; char map[1005][1005]; //数组开大点 int mapx[1005][1005],mapy[1005][1005]; int ma[1005][1005];//邻接矩阵true代表有边相连 int result[1005],visit[1005]; int x,y; int find(int a){ int i; for(i=1;i<=y;i++){ if(!visit[i]&&ma[a][i]){//如果节点i与a相邻并且未被查找过 visit[i]=1;//标记i为已查找过 if(!result[i]||find(result[i])){//如果i未在前一个匹配M中或者i在匹配M中,但是从与i相邻的节点出发可以有增广路 result[i]=a;//记录查找成功记录 return 1; } } } return 0; } int main(){ int i,j,ans; while(scanf("%d%d",&n1,&n2)!=EOF){ for(i=0;i<n1;i++){ scanf("%s",map[i]); } memset(mapx,0,sizeof(mapx)); memset(mapy,0,sizeof(mapy)); x=0; for(i=0;i<n1;i++){ for(j=0;j<n2;j++){ if(map[i][j]=='*'){ ++x; while(j<n2&&map[i][j]=='*'){ mapx[i][j]=x; j++; } } } } y=0; for(j=0;j<n2;j++){ for(i=0;i<n1;i++){ if(map[i][j]=='*'){ ++y; while(i<n1&&map[i][j]=='*'){ mapy[i][j]=y; i++; } } } } for(i=0;i<n1;i++){ for(j=0;j<n2;j++){ ma[mapx[i][j]][mapy[i][j]]=1; } } ans=0; memset(result,0,sizeof(result)); for(i=1;i<=x;i++){ memset(visit,0,sizeof(visit));//清空上次搜索时的标记 ans+=find(i); //从节点i尝试扩展 } printf("%d\n",ans); } return 0; }
http://blog.csdn.net/pi9nc/article/details/11848327