POJ 2226二分图最大匹配

#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;
}