poj 3041 Asteroids （二分图 匈牙利算法）

http://poj.org/problem?id=3041

题意：

给你N*N的矩阵，里面有的方格里有小行星，你需要用激光射掉它。。。激光可以射掉一行 或者一列的小行星，问最小需要发射多少次

 

这道题 很久 以前就做过了  ，现在有做了 一下 ，对 匈牙利 有 个 更好的了解。。。

 

转自 别处 ：

匈牙利算法是寻找最大匹配的优秀算法，那么与这个看上去一点也不像二分图的题来说有什么用处呢？让我们来做一个尝试：把样例数据里面的横坐标作为二分图的一部，纵坐标作为二分图的另一部，坐标为(x, y)的小行星表示为从横坐标x到纵坐标y的一段弧，就有了下图：

 

 

可以看出，原问题变成了下面这个问题：给定一个二分图G = (V, E)，定义一个点如果被覆盖，那么称所有与这个点相邻的弧被覆盖。求出最少需要覆盖多少个点才能覆盖所有的边。如上两个图，原问题中在x = 1和y = 2两处使用超级武器，等价于在右图中覆盖左边的点1和右边的点2。我们称这个问题的答案为最小覆盖数P，显然在样例中P = 2。

现在我们给出结论，若将二分图最大匹配的边集设为M，则P = | M |。也就是最大匹配数。证明如下：

首先证明P ≥ | M | 。对于M中的每一条边，它们都不相邻。那么至少需要覆盖 | M | 个点才能覆盖M中的所有边。故有P ≥ | M | 。

现在给出一种令P = | M |的覆盖方法。如果当前图中对于所有的v ∈ V，都有v ∈ M，显然有P = | M |。否则，任取一条边v0，使v0 ∈ V且v0 ∉ M。那么显然v0的两个端点至少有一个被匹配。

如 果有一端被匹配，设这个端点为a0。首先覆盖a0，然后观察与a0匹配的点a1的状况。如果存在与a1相邻的边v1，满足v1 ∈ V且v1 ∉ M，则v1一定不与没有匹配过的点相邻，否则，v0，v(a0, a1)及v1就形成了增广路。现在继续覆盖与v1相邻的另外一个点a2，然后依次操作，直到新连接到的匹配过的点不与V \ M中的边相邻为止。

如果两端都被匹配，就向两个方向进行同样的操作。从v0的一系列操作结束之后，考察还没有被覆盖的边，继续进行操作。这样最后所有的边一定会被全部覆盖。

如 果说有没有被覆盖到的边，这条边不可能同时与两个匹配过的点相邻。如果这样，在刚才的过程中一定会检索到这条边。所以说它必然是与一个匹配过的点ax相 邻，与一个没有匹配过的点ar相邻，并且ax的匹配点ay一定被覆盖过。那么，从ay向上找寻刚才的过程遍历过的边，与最后的v(ax, ar)一定会形成一条增广路，矛盾。

观察到每条M中的边都有且仅有一个端点被覆盖，所以说P = | M |。

由于P ≥ | M | 且存在P = | M |的状况，故P = | M | ，证明完毕。也就是说，这个题只需要用匈牙利算法求出最大匹配数就可以了。

 

 

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cmath>
 4 #include<iostream>
 5 #include<algorithm>
 6 #include< set>
 7 #include<map>
 8 #include<queue>
 9 #include<vector>
10 #include< string>
11  #define Min(a,b) a<b?a:b
12  #define Max(a,b) a>b?a:b
13  #define CL(a,num) memset(a,num,sizeof(a));
14  #define eps  1e-12
15  #define inf   0x7fffffff
16 
17  // freopen("data.txt","r",stdin);
18  const  double pi  = acos(- 1.0);
19 typedef   __int64  ll;
20  const  int maxn =  510 ;
21  using  namespace std;
22  int n , m;
23  int mat[maxn][maxn],result[maxn],vis[maxn] ;
24  void init()
25 {
26     memset(mat, 0, sizeof(mat));
27     memset(result, 0, sizeof(result));
28 }
29  int find( int a)
30 {
31      int i;
32      for(i= 1;i<=n;i++)
33     {
34          if(mat[a][i]&&!vis[i])
35         {
36             vis[i]= 1;
37              if(result[i]== 0||find(result[i]))
38             {
39                 result[i]=a;
40                  return  1;
41             }
42         }
43     }
44      return  0;
45 }
46  int main()
47 {
48      int  i ,x,y ;
49      while(scanf( " %d%d ",&n,&m)!=EOF)
50     {
51         init() ;
52          for(i =  0 ; i < m;i++)
53         {
54             scanf( " %d%d ",&x,&y);
55             mat[x][y] =  1;
56 
57         }
58          int ans =  0 ;
59          for(i =  1 ; i<=n;i++ )
60         {
61             CL(vis, 0);
62              if(find(i))ans ++ ;
63         }
64         printf( " %d\n ",ans) ;
65 
66     }
67 }