POJ 2226二分图最大匹配

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

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
int n1,n2;
char map[][];  //数组开大点
int mapx[][],mapy[][];
int ma[][];//邻接矩阵true代表有边相连
int result[],visit[];
int x,y;

int find(int a){
    int i;
    for(i=;i<=y;i++){
        if(!visit[i]&&ma[a][i]){//如果节点i与a相邻并且未被查找过
            visit[i]=;//标记i为已查找过
            if(!result[i]||find(result[i])){//如果i未在前一个匹配M中或者i在匹配M中，但是从与i相邻的节点出发可以有增广路
                result[i]=a;//记录查找成功记录
                return ;
            }
        }
    }
    return ;
}

int main(){
    int i,j,ans;
    while(scanf("%d%d",&n1,&n2)!=EOF){
        for(i=;i<n1;i++){
            scanf("%s",map[i]);
        }
        memset(mapx,,sizeof(mapx));
        memset(mapy,,sizeof(mapy));
        x=;
        for(i=;i<n1;i++){
            for(j=;j<n2;j++){
                if(map[i][j]=='*'){
                    ++x;
                    while(j<n2&&map[i][j]=='*'){
                        mapx[i][j]=x;
                        j++;
                    }
                }
            }
        }
        y=;
        for(j=;j<n2;j++){
            for(i=;i<n1;i++){
                if(map[i][j]=='*'){
                    ++y;
                    while(i<n1&&map[i][j]=='*'){
                        mapy[i][j]=y;
                        i++;
                    }
                }
            }
        }

        for(i=;i<n1;i++){
            for(j=;j<n2;j++){
                ma[mapx[i][j]][mapy[i][j]]=;
            }
        }

        ans=;

        memset(result,,sizeof(result));
        for(i=;i<=x;i++){
            memset(visit,,sizeof(visit));//清空上次搜索时的标记
            ans+=find(i); //从节点i尝试扩展
        }
        printf("%d\n",ans);
    }
    return ;
}

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