POJ 3020 Antenna Placement(二分图匹配）

题目链接

无向二分图的最小路径覆盖 = 顶点数 – 最大二分匹配数/2

看题解过的。看了网上的题解，说是拆点，然后把所有的城市分为两个集合，然后用匈牙利算法计算出最大二分图匹配。

#include <cstdio>
#include <cstring>
using namespace std;
char str[51][51];
int p[501][501],used[501],linker[501],o[51][51];
int x[4] = {0,0,1,-1};
int y[4] = {1,-1,0,0};
int num;
int dfs(int x)
{
    int i;
    for(i = 1;i <= num;i ++)
    {
        if(p[x][i]&&!used[i])
        {
            used[i] = 1;
            if(linker[i] == -1||dfs(linker[i]))
            {
                linker[i] = x;
                return 1;
            }
        }
    }
    return 0;
}
int main()
{
    int t,i,j,k,ans,n,m,r,c;
    scanf("%d",&t);
    while(t--)
    {
        memset(p,0,sizeof(p));
        memset(o,0,sizeof(o));
        memset(linker,-1,sizeof(linker));
        scanf("%d%d",&n,&m);
        num = 1;
        for(i = 0;i < n;i ++)
        scanf("%s",str[i]);
        for(i = 0;i < n;i ++)
        {
            for(j = 0;j < m;j ++)
            {
                if(str[i][j] == '*')
                o[i][j] = num ++;
            }
        }
        num --;
        for(i = 0;i < n;i ++)
        {
            for(j = 0;j < m;j ++)
            {
                if(str[i][j] == '*')
                {
                    for(k = 0;k <= 3;k ++)
                    {
                        r = i + x[k];
                        c = j + y[k];
                        if(r >= 0&&r < n&&c >= 0&&c < m&&o[r][c])
                        {
                            p[o[r][c]][o[i][j]] = 1;
                        }
                    }
                }
            }
        }
        ans = 0;
        for(i = 1;i <= num;i ++)
        {
            memset(used,0,sizeof(used));
            if(dfs(i))
            ans ++;
        }
        printf("%d\n",num-ans/2);
    }
    return 0;
}