POJ 3020 -Antenna Placement-二分图匹配

题意：一个N*M的矩阵里有K个观测点，你必须放置天线覆盖所有观测点。每个雷达只能天线两个观测点，这两点必须相邻。计算最少天线数。

做法：将所有相邻的观测点连起来，建图。跑一遍匈牙利算法就计算出了最大的覆盖数，除以二就是天线数。还要加上落单的观测点，每个都需要一个天线。

/*--------------------------------------------------------------------------------------*/
//        Helica's header 
//        Second Edition
//        2015.11.7
//
#include <algorithm>
#include <iostream>
#include <cstring>
#include <ctype.h>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <string>
#include <queue>
#include <stack>
#include <cmath>
#include <set>
#include <map>

//debug function for a N*M array 
#define debug_map(N,M,G) printf("
");for(int i=0;i<(N);i++)
{for(int j=0;j<(M);j++){
printf("%d",G[i][j]);}printf("
");}                
//debug function for int,float,double,etc.
#define debug_var(X) cout<<#X"="<<X<<endl;
/*--------------------------------------------------------------------------------------*/
using namespace std;

const int maxn = 100;
int N,M,T;
int G[10*maxn][10*maxn],Map[maxn][maxn];
int dx[4]={1,-1,0,0},dy[4]={0,0,1,-1};
int linker[10*maxn];
bool used[10*maxn];
int uN,vN;

bool dfs(int u)
{
    for(int v=1;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()
{
    memset(linker,-1,sizeof linker);
    int res = 0;
    for(int u=1;u<=uN;u++)
    {
        memset(used,false,sizeof used);
        if(dfs(u)) res++;    
    }    
    return res;
}

int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&N,&M);
        char c;
        uN = 0;
        for(int i=0;i<N;i++)
        {
            getchar();
            for(int j=0;j<M;j++)
            {
                scanf("%c",&c);
                if(c == '*') Map[i][j] = ++uN;
                else Map[i][j] = 0;
            }
        }
        memset(G,0,sizeof G);
        vN = uN;
        for(int i=0;i<N;i++)
        {
            for(int j=0;j<M;j++) if(Map[i][j])
            {
                int u = Map[i][j];
                for(int p=0;p<4;p++)
                {
                    int nx = i+dx[p],ny = j+dy[p];
                    if(nx >=0 && nx < N && ny >=0 && ny < M)
                    {
                        if(int v = Map[nx][ny])
                           {
                            //printf("u:%d v:%d
",u,v);
                            G[u][v] = G[v][u] = 1;
                        }    
                    }
                }
            }
        }
        
        int ans = hungary();
        //printf("ans:%d uN:%d
",ans,uN);
        printf("%d
",(uN-ans)+ans/2);
    }    
}