poj 3041 Asteroids

来源：互联网发布：年度总结数据分析编辑：程序博客网时间：2024/05/16 17:10

匈牙利算法求最大匹配。

把x和y看作两个集合，如果 (x, y)有一个asteroid，就在x和y之间连一条边，这样就有了一个二分图。

这里因为不需要求出具体的匹配，所以只需不断寻找交错链（也就是增广路、交错轨），每找到一条，结果加1（因为找到一条交错链意味着匹配数+1）

代码：

/*Poj: 3041 Asteroids*/#include <iostream>#include <cstdio>#include <cstring>#define Max 505using namespace std;int graph[Max][Max];int N, K;bool checked[Max];int xM[Max], yM[Max];int MaxMatch(void);bool SearchPath(int u);int main(){    //freopen("data.in", "rb", stdin);while(scanf("%d%d", &N, &K) != EOF) {    for(int i = 0; i < K; i++) {            int r, c;    scanf("%d%d", &r, &c);    graph[r-1][c-1] = 1;    }    printf("%d\n", MaxMatch());    }return 0;}int MaxMatch(void){int matches = 0;for(int i = 0; i < Max; i++)xM[i] = yM[i] = -1;for(int u = 0; u < N; u++) {if(xM[u] == -1) {memset(checked, false, sizeof(checked));if(SearchPath(u))matches++;}}return matches;}bool SearchPath(int u){for(int v = 0; v < N; v++) {    //寻找增广路if(graph[u][v] == 1 && !checked[v]) {checked[v] = true;if(yM[v] == -1 || SearchPath(yM[v])) {yM[v] = u;xM[u] = v;return true;}}}return false;}