POJ 3041 匈牙利算法

大致题意

N * N 的格子上有 K 颗行星，Bessie 有把大威力武器，一次可以消灭某一行或列的行星。求武器的最少使用次数。1 <= N <= 500, 1 <= K <= 10,000

x 或 y 方向的一次攻击看做顶点，坐标为 (x, y) 的点看做 x 指向 y 的一条边，问题就化为 x∈V, y∈U 的二分图匹配问题。求覆盖所有点的最小使用次数，即求覆盖所有边的最小点数。根据 $Konig$ 定理，二分图的最小点集覆盖 = 最大匹配，最后用匈牙利算法求最大匹配即可。

#include <cstdio>
#include <stdlib.h>
#include <algorithm>
#include <vector>
#define min(a,b)    (((a) < (b)) ? (a) : (b))
#define max(a,b)    (((a) > (b)) ? (a) : (b))
#define abs(x)    ((x) < 0 ? -(x) : (x))
#define INF 0x3f3f3f3f
#define delta 0.85
#define eps 1e-3
#define M_PI 3.14159265358979323846
#define MAX_N 505
using namespace std;
int N, K;
vector<int> G[MAX_N * 2];
bool used[MAX_N * 2];
int match[MAX_N * 2];

bool dfs(int v){
	used[v] = 1;
	for(int i = 0; i < G[v].size(); i++){
		int u = G[v][i], w = match[u];
		if(w < 0 || !used[w] && dfs(w)){
			match[v] = u;
			match[u] = v;
			return 1;
		}
	}
	return 0;
}

int main(){
	while(~scanf("%d%d", &N, &K)){
		for(int r, c, i = 0; i < K; i++){
			scanf("%d%d", &r, &c);
			G[r].push_back(c + N);
			G[c + N].push_back(r);
		}
		int res = 0;
		memset(match, -1, sizeof(match));
		for(int v = 1; v <= 2 * N; v++){
			if(match[v] < 0){
				memset(used, 0, sizeof(used));
				if(dfs(v)) ++res;
			}
		}
		printf("%d\n", res);
	}
	return 0;
}