HDU 4185(二分圖匹配)

/**
 第一次做二分圖匹配，剛開始不知道，讀完題便開始貪心，
 後來找了好幾組測試樣例，才知道不能這麼做，看網上別人
 做的，才知道是二分圖匹配，學習了一下。
*/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long LL;
const int maxn = 700;
int kase, n, num, T = 0;
char s[maxn][maxn];
int x[maxn], y[maxn], a[maxn][maxn];
bool vis[maxn], map[maxn][maxn];

bool Find(int i) {
    for(int j = 0; j < num; ++j) {
        if(map[i][j] && !vis[j]) {
            vis[j] = true;
            if(!y[j] || Find(y[j])) {
                y[j] = i;
                x[i] = j;
                return true;
            }
        }
    }
    return false;
}

int solve() {
    memset(x, 0, sizeof(x));
    memset(y, 0, sizeof(y));
    int ans = 0;
    for(int i = 0; i < num; ++i) {
        memset(vis, false, sizeof(vis));
        if(Find(i)) ans++;
    }
    return ans/2;
}

int main() {
    scanf("%d", &kase);
    while(kase--) {
        num = 0;
        memset(map, false, sizeof(map));
        scanf("%d", &n);
        for(int i = 0; i < n; ++i) scanf("%s", s[i]);

        for(int i = 0; i < n; ++i)
            for(int j = 0; j < n; ++j)
                if(s[i][j] == '#') a[i][j] = num++;

        for(int i = 0; i < n; ++i) {
            for(int j = 0; j < n; ++j) {
                if(s[i][j] == '#' && s[i][j+1] == '#' && j+1 < n) {
                    map[a[i][j]][a[i][j+1]] = true;
                    map[a[i][j+1]][a[i][j]] = true;
                }
                if(s[i][j] == '#' && s[i+1][j] == '#' && i+1 < n) {
                    map[a[i][j]][a[i+1][j]] = true;
                    map[a[i+1][j]][a[i][j]] = true;
                }
            }
        }
        int ans = solve();
        printf("Case %d: %d\n", ++T, ans);
    }
}