首页 > 代码库 > POJ 2226二分图最大匹配

POJ 2226二分图最大匹配

    匈牙利算法是由匈牙利数学家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

POJ 2226二分图最大匹配