[POJ3020]Antenna Placement（二分图最大匹配，最小边覆盖）

题目链接：http://poj.org/problem?id=3020

题意：要求用最少的1*2的边数覆盖掉图上的所有*

在专题里看这个图很高能所以不敢做，其实是最小边覆盖的裸题，先给*标号，随后按照每个*附近是否有*来构造二分图。

最后结果 最小边覆盖 = 节点数 - 最大匹配 （虽然和最小路径覆盖结果一样，但是是完全不同的两个事情）.

最大匹配是这个图上匹配的1/2，因为建图的时候是双向的。

/*
━━━━━┒ギリギリ♂ eye！
┓┏┓┏┓┃キリキリ♂ mind！
┛┗┛┗┛┃＼○／
┓┏┓┏┓┃ /
┛┗┛┗┛┃ノ)
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┃┃┃┃┃┃
┻┻┻┻┻┻
*/
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
using namespace std;
#define fr first
#define sc second
#define cl clear
#define BUG puts("here!!!")
#define W(a) while(a--)
#define pb(a) push_back(a)
#define Rint(a) scanf("%d", &a)
#define Rs(a) scanf("%s", a)
#define Cin(a) cin >> a
#define FRead() freopen("in", "r", stdin)
#define FWrite() freopen("out", "w", stdout)
#define Rep(i, len) for(int i = 0; i < (len); i++)
#define For(i, a, len) for(int i = (a); i < (len); i++)
#define Cls(a) memset((a), 0, sizeof(a))
#define Clr(a, x) memset((a), (x), sizeof(a))
#define Full(a) memset((a), 0x7f7f7f, sizeof(a))
#define lrt rt << 1
#define rrt rt << 1 | 1
#define pi 3.14159265359
#define RT return
#define lowbit(x) x & (-x)
#define onecnt(x) __builtin_popcount(x)
typedef long long LL;
typedef long double LD;
typedef unsigned long long ULL;
typedef pair<int, int> pii;
typedef pair<string, int> psi;
typedef pair<LL, LL> pll;
typedef map<string, int> msi;
typedef vector<int> vi;
typedef vector<LL> vl;
typedef vector<vl> vvl;
typedef vector<bool> vb;

const int maxn = 1910;
int nu, nv;
int G[maxn][maxn];
int linker[maxn];
bool vis[maxn];

bool dfs(int u) {
    For(v, 1, nv+1) {
        if(G[u][v] && !vis[v]) {
            vis[v] = 1;
            if(linker[v] == -1 || dfs(linker[v])) {
                linker[v] = u;
                return 1;
            }
        }
    }
    return 0;
}

int hungary() {
    int ret = 0;
    Clr(linker, -1);
    For(u, 1, nu+1) {
        Cls(vis);
        if(dfs(u)) ret++;
    }
    return ret;
}

const int dx[5] = {0, 0, 1, -1};
const int dy[5] = {1, -1, 0, 0};
int n, m;
int mp[maxn][maxn];
int ord;
char tmp[maxn];

bool ok(int x, int y) {
    return x >= 0 && y>= 0 && x < n && y < m;
}

int main() {
    // FRead();
    int T;
    Rint(T);
    W(T) {
        ord = 0; Cls(G); Cls(mp);
        Rint(n); Rint(m);
        Rep(i, n) {
            Rs(tmp);
            Rep(j, m) {
                if(tmp[j] == '*') mp[i][j] = ++ord;
            }
        }
        nu = nv = ord;
        Rep(i, n) {
            Rep(j, m) {
                if(mp[i][j]) {
                        Rep(k, 4) {
                        int x = i + dx[k];
                        int y = j + dy[k];
                        if(ok(x, y) && mp[x][y]) {
                            G[mp[i][j]][mp[x][y]] = 1;
                        }
                    }
                }
            }
        }
        printf("%d
", nu - hungary() / 2);
    }
    RT 0;
}