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>
#include<string>
#define Min(a,b) a<b?a:b
#define Max(a,b) a>b?a:b
#define CL(a,num) memset(a,num,sizeof(a));
#define eps 1e-12
#define inf 0x7fffffff

//freopen("data.txt","r",stdin);
const double pi = acos(-1.0);
typedef __int64 ll;
const int maxn = 510 ;
using namespace std;
int n , m;
int mat[maxn][maxn],result[maxn],vis[maxn] ;
void init()
{
memset(mat,0,sizeof(mat));
memset(result,0,sizeof(result));
}
int find(int a)
{
int i;
for(i=1;i<=n;i++)
{
if(mat[a][i]&&!vis[i])
{
vis[i]=1;
if(result[i]==0||find(result[i]))
{
result[i]=a;
return 1;
}
}
}
return 0;
}
int main()
{
int i ,x,y ;
while(scanf("%d%d",&n,&m)!=EOF)
{
init() ;
for(i = 0 ; i < m;i++)
{
scanf("%d%d",&x,&y);
mat[x][y] = 1;

}
int ans = 0 ;
for(i = 1 ; i<=n;i++ )
{
CL(vis,0);
if(find(i))ans ++ ;
}
printf("%d\n",ans) ;

}
}