[POJ 3041][USACO 2005] Asteroids Hungary算法求最小点覆盖

题目传送门：【POJ 3041】

题目大意：贝茜想驾驶她的飞船穿过危险的小行星群。小行星群是一个 N * N 的网格 (1 ≤ N ≤ 500)，在网格内有 K 个小行星 ( 1 ≤ K ≤ 10000)。幸运地是贝茜有一个很强大的武器，一次可以消除所有在一行或一列中的小行星。这种武器很贵，所以她希望尽量地少用。给出所有的小行星的位置，算出贝茜最少需要多少次射击就能消除所有的小行星。

输入第一行为 N 和 K，之后的 K 行，每行两个整数 R，C，代表一个小行星的位置。输出为一个整数，代表使用武器的最少次数。

题目分析：
这道题是我拿来当模板练手的。
由题，由于每次我们可以消除整行或整列的小行星，根据题目性质，所以我们可以将其转化为二分图，求出其最小点覆盖即可。又因为最小点覆盖等于最大匹配，所以直接跑一次关于二分图最大匹配的算法即可。

下面附上代码：

[cpp] view plain copy

print ?

1. #include<cstdio>  
2. #include<cstring>  
3. #include<iostream>  
4. #include<algorithm>  
5. using namespace std;  
6. const int MXN = 505,MXK = 10005;  
7.   
8. struct Edge{  
9.     int to,next;  
10. };  
11. Edge edge[MXK * 2];  
12. int maps[MXN][MXN];  
13. int head[MXK],now = 0,match[MXN * 2],n,k,ans = 0;  
14. bool vis[MXN * 2];  
15.   
16. void adde(int u,int v){  
17.     edge[++now].to = v;  
18.     edge[now].next = head[u];  
19.     head[u] = now;  
20. }  
21.   
22. bool hungary(int u){                    //最小点覆盖 = 最大匹配   
23.     for (int i = head[u];i;i = edge[i].next){  
24.         int v = edge[i].to;  
25.         if (!vis[v]){  
26.             vis[v] = true;  
27.             if (match[v] == 0 || hungary(match[v])){  
28.                 match[v] = u;  
29.                 return true;  
30.             }  
31.         }  
32.     }  
33.     return false;  
34. }  
35.   
36. int main(){  
37.     int T,a,b;  
38.     cin>>n>>k;  
39.     for (int i = 1;i <= k;i++){  
40.         scanf(”%d%d”,&a,&b);  
41.         adde(a,b + 500);  
42.         adde(b + 500,a);  
43.     }  
44.     for (int i = 1;i <= n;i++){  
45.         memset(vis,0,sizeof(vis));  
46.         if (hungary(i)) ++ans;  
47.     }  
48.     cout<<ans<<endl;  
49.     return 0;  
50. }