poj 2226 Muddy Fields(二分图最大匹配)

来源:互联网 发布:app store推荐软件 编辑:程序博客网 时间:2024/06/05 07:48

貌似和 hdu 1045 Fire Net 差不多。
按照行列划分为两个点集,对于每一行,连续的 * 为一个木板,看作一个点,给他一个编号。对于每一列,同行进行一样的操作。行的那些编号是一个点集,列的那些编号是一个点集。如果两个 * 的块有交点,则在这两点之间连边,边就相当于木板之间的交点,即泥泞地的交点。这就是个二分图了。
用最少的木板覆盖所有的泥泞地,也就是覆盖上这些交点,即二分图中的边,因为木板长度不限制,所以只要覆盖上所有交点了,泥泞地也就全覆盖了,这就是最小点覆盖了,即二分图的最大匹配。

#include <stdio.h>#include <string.h>char gp[55][55];int number[55][55];const int MAXN = 510;int uN,vN;//u,v的数目,使用前面必须赋值int g[MAXN][MAXN];//邻接矩阵int linker[MAXN];bool used[MAXN];bool dfs(int u){    for(int 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;    memset(linker,-1,sizeof(linker));    for(int u = 0; u < uN; u++)    {        memset(used,false,sizeof(used));        if(dfs(u))res++;    }    return res;}int main(){    int r,c;    scanf("%d %d",&r,&c);    for(int i = 0; i < r; ++i)    {        scanf("%s",gp[i]);        gp[i][c] = '.';    }    ++c;    for(int i = 0; i < c; ++i)        gp[r][i] = '.';    ++r;    int nb = 0;    bool flag = false;    for(int i = 0; i < r; ++i)    {        flag = false;        for(int j = 0; j < c; ++j)        {            if(gp[i][j] == '*')            {                flag = true;                number[i][j] = nb;            }            else if(flag)            {                ++nb;                flag = false;            }        }    }    uN = nb;    nb = 0;    for(int i = 0; i < c; ++i)    {        flag = false;        for(int j = 0; j < r; ++j)        {            if(gp[j][i] == '*')            {                flag = true;                g[number[j][i]][nb] = 1;            }            else if(flag)            {                ++nb;                flag = false;            }        }    }    vN = nb;    int res = hungary();    printf("%d\n",res);    return 0;}
原创粉丝点击