hdu-4185.loiol_skimming(简单二分匹配模型)

/*************************************************************************
    > File Name: hdu-4185.oil_skimming.cpp
    > Author: CruelKing
    > Mail: 2016586625@qq.com 
    > Created Time: 2019年09月03日 星期二 09时12分12秒
    本题思路:简单分析过后就可以知道如果一点a被另一个点b匹配,那么和b匹配的点c不可能和点a匹配,因为每个结点都只能匹配四个方向,所以这个匹配的图就可以看做是一个二分图,因此对于每个结点,将能和他匹配的边存入,接着跑一波二分匹配之后我们得到的是无向图情况下的最大匹配,所以除以二就是最后的答案.
 ************************************************************************/

#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;

const int maxn = 600 + 5;
char str[maxn][maxn];

int n;
int linker[maxn * maxn];
bool used[maxn * maxn];
int tot, head[maxn * maxn];
int un[maxn * maxn], cnt;

struct Edge {
    int to, next;
} edge[maxn * maxn * 4 + 5];

void init() {
    memset(head, -1, sizeof head);
    tot = 0;
}

void addedge(int u, int v) {
    edge[tot] = (Edge) {v, head[u]};
    head[u] = tot ++;
}

bool dfs(int u) {
    for(int k = head[u]; ~k; k = edge[k].next) {
        int v = edge[k].to;
        if(!used[v]) {
            used[v] = true;
            if(linker[v] == -1 || dfs(linker[v])) {
                linker[v] = u;
                return true;
            }
        }
    }
    return false;
}

int main() {
    int k, Case = 0;
    scanf("%d", &k);
    while(k --) {
        init();
        cnt = 0;
        scanf("%d", &n);
        for(int i = 0; i < n; i ++) {
            scanf("%s", str[i]);
        }
        for(int i = 0; i < n; i ++) {
            for(int j = 0; j < n; j ++) {
                if(str[i][j] == '#') {
                    un[cnt ++] = (i * n + j + 1);
                    for(int dx = -1; dx <= 1; dx ++) {                    
                        for(int dy = -1; dy <= 1; dy ++) {                        
                            if(abs(dx - dy) == 1) {        
                                if(dx + i >= 0 && dx + i < n && dy + j >= 0 && dy + j < n) {                    
                                    if(str[dx + i][dy + j] == '#') addedge(i * n + j + 1, (i + dx) * n + j + dy + 1);
                                }
                            }
                        }
                    }
                }
            }
        }
        int res = 0;
        memset(linker, -1, sizeof linker);
        for(int i = 0; i < cnt; i ++) {
            memset(used, false, sizeof used);
            if(dfs(un[i])) res ++;
        }
        printf("Case %d: %d
", ++Case, res / 2);
    }
    return 0;
　　}