【最小边覆盖 && 无向图】POJ - 3020 Antenna Placement

Problem Description

感兴趣的点用*表示
一个n*m的矩阵， 一个雷达可覆盖相邻的两个位置，一个感兴趣的点可被多个雷达覆盖, 至少需要多少个雷达，才能把所有的感兴趣的点覆盖掉。

思路

对于不存在孤立点的图，最大匹配数 + 最小边覆盖数 = 顶点数
紧接着就是建图：
相邻的感兴趣的点可以到达，我们就给它们之间连接一条边。所以我们对于每个感兴趣的点都给它们一个编号。

#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
int id[50][50], n, m, match[2500], vis[2500];
int x[5] = {1, -1, 0, 0};
int y[5] = {0, 0, 1, -1};
char s[50][50];
vector<int> Map[2500];
bool dfs(int u)//匈牙利算法
{
    int mm = Map[u].size();
    for(int i = 0; i < mm; i++)
    {
        int to = Map[u][i];
        if(!vis[to])
        {
            vis[to] = 1;
            if(!match[to] || dfs(match[to]))
            {
                match[to] = u;//只适用于双向图
                match[u] = to;
                return 1;
            }
        }
    }
    return 0;
}
bool judge(int cx, int cy)
{
    if(cx >= 0 && cx < n && cy >= 0 && cy < m && s[cx][cy] == '*') return 1;
    else return 0;
}
void add(int u, int v)
{
    Map[u].push_back(v);
}
int main()
{
    int T, i, j, k;
    scanf("%d", &T);
    while(T--)
    {
        scanf("%d %d", &n, &m);
        int len = 1;
        for(i = 0; i < n; i++)
        {
            scanf("%s", s[i]);
            for(j = 0; j < m; j++)
            {
                if(s[i][j] == '*')//给每个感兴趣的点标号
                    id[i][j] = len++;
            }
        }
        for(i = 0; i <= len; i++)//初始化
            Map[i].clear();
        for(i = 0; i < n; i++)//建图
        {
            for(j = 0; j < m; j++)
            {
                if(s[i][j] == '*')
                {
                    for(k = 0; k < 4; k++)
                    {
                        int cx = x[k] + i, cy = y[k] + j;
                        if(judge(cx, cy))
                        {
                            add(id[i][j], id[cx][cy]);
                        }
                    }
                }
            }
        }
        memset(match, 0, sizeof(match));
        int ans = 0;
        for(int i = 1; i <= len; i++)//求最大匹配数
        {
            if(!match[i])
            {
                memset(vis, 0, sizeof(vis));
                if(dfs(i)) ans++;
            }
        }
        printf("%d\n",len - ans - 1);
    }
    return 0;
}