题目链接:
题目分析:
联想到Asteroids这道题,将行和列分别作为二分图的两边的点,把点本身作为二分图的边去处理
那么这个题的限制条件打不穿的墙怎么处理呢,发现这样的话
\[ ...X.. \]
其实左边那段和右边那段在横着考虑的时候是互不影响的,不妨对行和列重新编号
用一下某题解的剪枝讲解的图,样例第一个大概就是这个样子
然后对于每个空地的行编号和列编号连一条边,当空地已经被控制/放满了就是最终答案,跑最大匹配即可
代码:
#include <iostream>
#include <stdio.h>
#include <math.h>
#include <string.h>
#include <time.h>
#include <stdlib.h>
#include <string>
#include <bitset>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <algorithm>
#include <sstream>
#include <stack>
#include <iomanip>
#define N (800 + 10)
using namespace std;
inline int read() {
int cnt = 0, f = 1; char c = getchar();
while (!isdigit(c)) {if (c == '-') f = -f; c = getchar();}
while (isdigit(c)) {cnt = (cnt << 3) + (cnt << 1) + c - '0'; c = getchar();}
return cnt * f;
}
int n, cnt = 0, res = 0, l, r, ans;
char s[N][N];
int mapp[N][N], mapp2[N][N];
bool vis[N];
int nxt[N * 100], first[N * 100], to[N * 100], tot;
int match[N];
void add(int x, int y) {nxt[++tot] = first[x], first[x] = tot, to[tot] = y;}
void build_map() {
for (register int i = 1; i <= n; ++i)
for (register int j = 1; j <= n; ++j) {
if (s[i][j] == 'X') {
mapp[i][j] = 0;
continue;
} else {
if (s[i][j] == '.' && s[i][j - 1] == 'X') ++cnt;
mapp[i][j] = cnt;
}
}
l = cnt, cnt = 0;
for (register int i = 1; i <= n; ++i)
for (register int j = 1; j <= n; ++j) {
if (s[j][i] == 'X') {
mapp2[j][i] = 0;
continue;
} else {
if (s[j][i] == '.' && s[j - 1][i] == 'X') ++cnt;
mapp2[j][i] = cnt;
}
}
r = cnt, cnt = 0;
for (register int i = 1; i <= n; ++i)
for (register int j = 1; j <= n; ++j) {
if (!mapp[i][j]) continue;
add(mapp[i][j], mapp2[i][j]);
}
}
bool find (int u) {
for (register int i = first[u]; i; i = nxt[i]) {
int v = to[i];
if (vis[v]) continue;
else {
vis[v] = true;
if (match[v] == -1 || find(match[v])) {
match[v] = u;
return 1;
}
}
}
return 0;
}
int hungary() {
ans = 0;
for (register int i = 1; i <= l; ++i) {
memset(vis, 0, sizeof(vis));
ans += find(i);
}
return ans;
}
int main() {
// freopen("1.in", "r", stdin);
while (1) {
n = read();
if (!n) break;
memset(first, 0, sizeof (first));
memset(to, 0, sizeof(to));
memset(nxt, 0, sizeof(nxt));
tot = 0, l = 0, r = 0;
for (register int i = 0; i <= n + 1; ++i)
for (register int j = 0; j <= n + 1; ++j) s[i][j] = 'X';
for (register int i = 1; i <= n; ++i) scanf("%s", s[i] + 1);
for (register int i = 0; i <= n + 1; ++i) s[i][n + 1] = 'X';
build_map();
for (register int i = 1; i <= r; ++i) match[i] = -1;
res = hungary();
printf("%d\n", res);
}
return 0;
}