最小点覆盖poj 2226

#include<stdio.h>
#include<iostream>
#include<algorithm>
using namespace std;
const int D=2501;
const int maxn=55;
int x[maxn][maxn],y[maxn][maxn],ma[D*2],vis[D*2],first[D*2];
int X,Y,size,ret,N,M;
char map[55][55];
struct Edge
{
 int v,next;
}edge[D];
void insert(int u,int v)
{
 edge[size].v=v;
 edge[size].next=first[u];
 first[u]=size++;
}
void bulid()
{
 X=1;Y=D;
 memset(x,0,sizeof(x));
 memset(y,0,sizeof(y));
 for(int i=0;i<N;i++)
  for(int j=0;j<M;j++)
   if(map[i][j]=='*')
   {
    if(!x[i][j])
    {
     int k=j;
     while(k<M && map[i][k]=='*') x[i][k++]=X;
     X++;
    }
    if(!y[i][j])
    {
     int k=i;
     while(k<N && map[k][j]=='*') y[k++][j]=Y;
     Y++;
    }
    insert(x[i][j],y[i][j]);
   }
}
int dfs(int u)
{
 for(int e=first[u];e!=-1;e=edge[e].next)
 {
  int v=edge[e].v;
  if( !vis[v])
  {
   vis[v]=1;
   if(ma[v]==-1 || dfs(ma[v]) )
   {
    ma[v]=u;
    return 1;
   }
  }
 }
 return 0;
}
int maxMatch()
{
 int ret=0;
 memset(ma,-1,sizeof(ma));
 for(int u=1;u<X;u++)
 {
  memset(vis,0,sizeof(vis));
  ret+=dfs(u);
 }
 return ret;
}
int main()
{
// freopen("test.txt","r",stdin);
 scanf("%d%d",&N,&M);
 for(int i=0;i<N;i++)
  scanf("%s",map[i]);
 memset(first,-1,sizeof(first));size=0;
 bulid();
 printf("%d\n",maxMatch());
 return 0;
}