A s t e r o i d s Asteroids Asteroids
题目链接:ssl 1341 / poj 3041
题目
B e s s i e Bessie Bessie想以 N N N × × × N N N网格 ( 1 < = N < = 500 ) (1 <= N <= 500) (1<=N<=500)的形状在危险的小行星场中导航她的飞船。网格包含 K K K个小行星 ( 1 < = K < = 10 , 000 ) (1 <= K <= 10,000) (1<=K<=10,000),它们方便地位于网格的晶格点处。
幸运的是,贝茜拥有强大的武器,可以一击就蒸发掉网格中任何给定行或列中的所有小行星,这种武器非常昂贵,因此她希望谨慎使用。字段中,找到贝西需要射击以消除所有小行星的最小射击次数。
输入
第
1
1
1行:两个整数
N
N
N和
K
K
K,以单个空格分隔。
第
2
~
K
+
1
2~K + 1
2~K+1行:每行包含两个空格分隔的整数
R
R
R和
C
(
1
<
=
R
,
C
<
=
N
)
C(1 <= R,C <= N)
C(1<=R,C<=N),分别表示小行星的行和列坐标。
输出
第 1 1 1行:整数,表示 B e s s i e Bessie Bessie必须拍摄的最小次数。
样例输入
3 4
1 1
1 3
2 2
3 2
样例输出
2
样例解释
下图表示数据,其中
“
X
”
“X”
“X”是小行星,
“
.
”
“.”
“.”是空白空间:
X
.
X
X.\ X
X. X
.
X
.
.\ X.
. X.
.
X
.
.\ X.
. X.
B
e
s
s
i
e
Bessie
Bessie可以在第
1
1
1行射击以摧毁
(
1
,
1
)
(1,1)
(1,1)和
(
1
,
3
)
(1,3)
(1,3)处的小行星,然后她可以击落第
2
2
2列以摧毁
(
2
,
2
)
(2,2)
(2,2)和
(
3
,
2
)
(3,2)
(3,2)处的小行星。
因此,用两个子弹就可以解决所有敌人。
思路
这道题其实是最小点覆盖问题。
然后,通过某种神奇的证明,我们可以得出最小点覆盖问题的答案就是最大匹配。
我也不会证,我是这样理解的:
在选了那几个点之后,所有的边都必定有一端是那几个点的其中一个。
也就是说,无论你选那一条边,都会碰到那几个点的其中一个。
我们要让碰到的点不同,线尽可能的多(最大匹配的要求),就只能每的每一条边碰到的那几个点都不相同。
那么也就是说,最小点覆盖就是最大匹配了。
代码
#include<cstring>
#include<cstdio>
#include<queue>
using namespace std;
struct note {
int to, next, op, now;
}e[3000001];
int n, m, o, x, y, le[10001], k, s, t, ans, dis[10001];
queue<int>q;
bool bfs() {
while (!q.empty()) q.pop();
memset(dis, 0x7f, sizeof(dis));
dis[s] = 0;
q.push(s);
while (!q.empty()) {
int now = q.front();
q.pop();
for (int i = le[now]; i; i = e[i].next)
if (dis[e[i].to] > dis[now] + 1 && e[i].now) {
dis[e[i].to] = dis[now] + 1;
if (e[i].to == t) return 1;
q.push(e[i].to);
}
}
return 0;
}
int dfs(int now, int an) {
if (now == t) return an;
int go = 0;
for (int i = le[now]; i; i = e[i].next)
if (dis[e[i].to] == dis[now] + 1 && e[i].now) {
int line_go = dfs(e[i].to, min(e[i].now, an - go));
if (!line_go) dis[e[i].to] = -1;
e[i].now -= line_go;
e[e[i].op].now += line_go;
go += line_go;
if (go == an) break;
}
return go;
}
int main() {
scanf("%d %d", &n, &o);//读入
m = n;
s = n + m + 1;//设置源点汇点
t = n + m + 2;
for (int i = 1; i <= o; i++) {
scanf("%d %d", &x, &y);//读入
y += n;
e[++k] = (note){y, le[x], k + 1, 1}; le[x] = k;//连边
e[++k] = (note){x, le[y], k - 1, 0}; le[y] = k;
}
for (int i = 1; i <= n; i++) {//连源点
e[++k] = (note){i, le[s], k + 1, 1}; le[s] = k;
e[++k] = (note){s, le[i], k - 1, 0}; le[i] = k;
}
for (int i = n + 1; i <= n + m; i++) {//连汇点
e[++k] = (note){t, le[i], k + 1, 1}; le[i] = k;
e[++k] = (note){i, le[t], k - 1, 0}; le[t] = k;
}
while (bfs())//dinic算法
ans += dfs(s, 2147483647);
printf("%d", ans);//输出
return 0;
}