poj 2226(最小点覆盖)

构图有点意思。。没想到。。

下面是构图：

把行里面连在一起的坑连起来视为一个点，即一块横木板，编上序号，Sample则转化为：

1 0 2 0
0 3 3 3
4 4 4 0
0 0 5 0

把这些序号加入X集合，再按列做一次则为：

1 0 4 0
0 3 4 5
2 3 4 0
0 0 4 0

 将 每一个 湿地的点 按其分配的行列 号  连边 ，如（i,j）为湿地 ，则用 其分配的  行  ——》 列，这样，我们求的就是 ，最小点覆盖。

// File Name: 2226.cpp
// Author: Missa
// Created Time: 2013/2/10 星期日 23:58:51

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<cstdlib>
#include<map>
using namespace std;

const int maxn = 55;
int r,c;
int mc[maxn][maxn];
int mr[maxn][maxn];
//***************匈牙利***********************
int nx,ny;
bool vis[maxn*maxn];
int ma[maxn*maxn][maxn*maxn];
int link[maxn*maxn];
bool dfs(int x)
{
    for(int y=1;y<=ny;y++)
    {
        if(!vis[y] && ma[x][y])
        {
            vis[y]=1;
            if(link[y]==-1 || dfs(link[y]))
            {
                link[y]=x;
                return true;
            }
        }
    }
    return false;
}
int maxmatch()
{
    int ans=0;
    memset(link,-1,sizeof(link));
    for(int x=1;x<=nx;x++)
    {
        memset(vis,0,sizeof(vis));
        if(dfs(x)) ans++;
    }
    return ans;
}
//********************************************


int main()
{
    while(~scanf("%d%d",&r,&c))
    {
        char s[maxn];
        memset(mr,-1,sizeof(mr));
        memset(mc,-1,sizeof(mc));
        getchar();
        for(int i=1;i<=r;i++)
        {
            gets(s);
            for(int j=1;j<=c;j++)
            {
                if(s[j-1]=='*')
                {
                    mr[i][j]=0;
                    mc[i][j]=0;
                }
            }
        }
        int numr=0,numc=0;
        for(int i=1;i<=r;i++)
        {
            for(int j=1;j<=c;j++)
            {
                if(mr[i][j]==0)
                {
                    numr++;
                    mr[i][j]=numr;
                    for(int k=j;k<=c;k++)
                    {
                        if(mr[i][k]==-1) break;
                        mr[i][k]=numr;
                    }
                }
            }
        }
        for(int j=1;j<=c;j++)
        {
            for(int i=1;i<=r;i++)
            {
                if(mc[i][j]==0)
                {
                    numc++;
                    mc[i][j]=numc;
                    for(int k=i;k<=r;k++)
                    {
                        if(mc[k][j]==-1)break;
                        mc[k][j]=numc;
                    }
                }
            }
        }
        nx=numc,ny=numr;
        memset(ma,0,sizeof(ma));
        for(int i=1;i<=r;i++)
            for(int j=1;j<=c;j++)
                if(mc[i][j] != -1)
                    ma[mc[i][j]][mr[i][j]]=1;
        printf("%d\n",maxmatch());
    }
    return 0;
}