Asteroids【最小点覆盖问题】

$$Asteroids$$

题目链接：ssl 1341 / poj 3041

题目

$Bessie$想以$N$ $\times$ $N$网格$（1<=N<=500）$的形状在危险的小行星场中导航她的飞船。网格包含$K$个小行星$（1<=K<=10,000）$，它们方便地位于网格的晶格点处。

幸运的是，贝茜拥有强大的武器，可以一击就蒸发掉网格中任何给定行或列中的所有小行星，这种武器非常昂贵，因此她希望谨慎使用。字段中，找到贝西需要射击以消除所有小行星的最小射击次数。

输入

第$1$行：两个整数$N$和$K$，以单个空格分隔。
第$2～K+1$行：每行包含两个空格分隔的整数$R$和$C（1<=R，C<=N）$，分别表示小行星的行和列坐标。

输出

第$1$行：整数，表示$Bessie$必须拍摄的最小次数。

样例输入

3 4
1 1
1 3
2 2
3 2

样例输出

2

样例解释

下图表示数据，其中$“X”$是小行星，$“.”$是空白空间：
$X.X$
$.X.$
$.X.$

$Bessie$可以在第$1$行射击以摧毁$（1,1）$和$（1,3）$处的小行星，然后她可以击落第$2$列以摧毁$（2,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;
}