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
",ans);
    }
    return 0;
}

 http://blog.csdn.net/pi9nc/article/details/11848327