POJ 2226 Muddy Fields(匈牙利算法—最小点覆盖)

题目链接：POJ 2226 Muddy Fields

农夫John的养牛场，是一个R 行C 列的矩形，一场大雨后，养牛场低洼的地方都有了积水。John 的牛都很娇贵的，他们吃草的时候，不想把他们的蹄子给弄脏了。为了不让牛儿们把它们的蹄子弄脏，John 决定把有水的地方铺上木板。他的木板是宽度为1，长度没有限制的。
他想用最少数目的木板把所有有水的低洼处给覆盖上，前提是木板不能覆盖草地，但是可以重叠。
Sample:
4 4
*.*.
.***
***.
..*.


把行里面连在一起的坑连起来视为一个点，即一块横木板，编上序号，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

同样加入Y集合，一个坑只能被横着的或者被竖着的木板盖住，将原图的坑的也标上不同的序号，一共九个坑

1 . 2 .
. 3 4 5
67 8 .
. . 9 .

比如7号坑可以被横着的4号木板和竖着的3号木板盖住，把每个点的对应的横木板(4)和竖木板(3)中间连一条边的话,则问题转化为 找尽量少的边把这些点都盖住,根据定理便是求最大匹配数.

以上转自：http://www.cppblog.com/abilitytao/archive/2009/10/21/99124.aspx

#include <iostream>#include <cstdio>#include <cstring>using namespace std;const int MAX_N = 50 + 10;const int MAX_M = MAX_N * MAX_N;bool _map[MAX_M][MAX_M], vis[MAX_M];int link[MAX_M], x_map[MAX_N][MAX_N], y_map[MAX_N][MAX_N], r, c, x_cnt, y_cnt;char G[MAX_N][MAX_N];bool dfs(int u){    for(int i = 1; i <= y_cnt; i++)    {        if(!vis[i] && _map[u][i])        {            vis[i] = 1;            if(link[i] == -1 || dfs(link[i]))            {                link[i] = u;                return true;            }        }    }    return false;}int MaxMatch(){    int num = 0;    memset(link, -1, sizeof(link));    for(int i = 1; i <= x_cnt; i++)    {        memset(vis, 0, sizeof(vis));        if(dfs(i))            num++;    }    return num;}int main(){    while(scanf("%d%d", &r, &c) != EOF)    {        memset(x_map, 0, sizeof(x_map));        memset(y_map, 0, sizeof(y_map));        x_cnt = y_cnt = 1;        for(int i = 0; i < r; i++)            scanf("%s", G[i]);        for(int i = 0; i < r; i++)        {            for(int j = 0; j < c; j++)            {                if(G[i][j] == '*')                {                    if(j != 0 && G[i][j - 1] == '*' && G[i][j] == '*')                        x_map[i][j] = x_cnt;                    else if(j == 0 && G[i][j] == '*')                        x_map[i][j] = ++x_cnt;                    else if(j != 0 && G[i][j - 1] == '.' && G[i][j] == '*')                        x_map[i][j] = ++x_cnt;                    else if(G[i][j] == '.')                        x_map[i][j] = 0;                }            }        }        for(int j = 0; j < c; j++)        {            for(int i = 0; i < r; i++)            {                if(G[i][j] == '*')                {                    if(i != 0 && G[i - 1][j] == '*' && G[i][j] == '*')                        y_map[i][j] = y_cnt;                    else if(i == 0 && G[i][j] == '*')                        y_map[i][j] = ++y_cnt;                    else if(i != 0 && G[i - 1][j] == '.' && G[i][j] == '*')                        y_map[i][j] = ++y_cnt;                    else if(G[i][j] == '.')                        y_map[i][j] = 0;                }            }        }        for(int i = 0; i < r; i++)            for(int j = 0; j < c; j++)                if(G[i][j] == '*')                    if(x_map[i][j] > 0 && y_map[i][j] > 0)                        _map[x_map[i][j]][y_map[i][j]] = 1;        printf("%d\n", MaxMatch());    }    return 0;}