POJ 2226 Muddy Fields(二分匹配-hungary)

Description
一个n*m的草地,草地上有一些泥坑,现在可以用宽度为1,长度任意的长木板去覆盖泥坑,但是木板不能覆盖草地,问最少需要多少长木板
Input
第一行为两个整数n和m表示草地尺寸,之后为一n*m矩阵表示这块草地(’.’表示草地,’*’表示泥坑)
Output
问最少需要多少块木板才能覆盖所有泥坑
Sample Input
这里写图片描述
Sample Output
4
Solution
二分匹配,两排点,第一排点每个点代表一段横向连续的泥坑,第二排点每个点代表一段纵向连续的泥坑,对于每个泥坑,在其在第一排点中的编号i与其在第二排点中的编号j之间建边,对于新建的图问题转化为用最少的点覆盖所有的边,即求图的最小点覆盖,即也为最大匹配数
Code

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define maxn 2555
int uN,vN;//u,v数目
int g[maxn][maxn];//编号是0~n-1的 
int linker[maxn];
bool used[maxn];
bool dfs(int u)
{
    int v;
    for(v=0;v<vN;v++)
        if(g[u][v]&&!used[v])
        {
            used[v]=true;
            if(linker[v]==-1||dfs(linker[v]))
            {
                linker[v]=u;
                return true;
            }    
        }  
    return false;  
}    
int hungary()
{
    int res=0;
    int u;
    memset(linker,-1,sizeof(linker));
    for(u=0;u<uN;u++)
    {
        memset(used,0,sizeof(used));
        if(dfs(u))  res++;
    } 
    return res;   
}  
int main()
{
    int n,m,M[55][55][2];
    char Map[55][55];
    while(~scanf("%d%d",&n,&m)) 
    {
        for(int i=0;i<n;i++)
            scanf("%s",Map[i]);
        memset(g,0,sizeof(g));
        uN=vN=0;
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
                if(Map[i][j]=='*')
                {
                    while(Map[i][j]=='*')
                        M[i][j++][0]=uN;
                    uN++;
                }
        for(int j=0;j<m;j++)
            for(int i=0;i<n;i++)
                if(Map[i][j]=='*')
                {
                    while(Map[i][j]=='*')
                        M[i++][j][1]=vN;
                    vN++;
                }
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
                if(Map[i][j]=='*')
                    g[M[i][j][0]][M[i][j][1]]=1;
        int ans=hungary(); 
        printf("%d\n",ans);
    }
    return 0;
} 

你可能感兴趣的:(POJ 2226 Muddy Fields(二分匹配-hungary))