poj3041 Asteroids 二分匹配之最小顶点覆盖

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

在N*N的平面上有K颗小行星，现在你要摧毁他们，你的每一发子弹可以摧毁同一行，或者是同一列上的小行星，现在问你最少要多少子弹才能摧毁所有的小行星？
构图：如果在i行j列上有一颗小行星 则graph[i][j]=1,再求最大匹配即可。
这一题用到的结论是 ：最小顶点覆盖 = 最大匹配(最小覆盖要求用最少的点（Ｘ集合或Ｙ集合的都行）让每条边都至少和其中一个点关联)
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
int nx, ny; // X的點數目、Y的點數目
int mx[505], my[505]; // X各點的配對對象、Y各點的配對對象
bool vy[505]; // 紀錄Graph Traversal拜訪過的點
bool adj[505][505]; // 精簡過的adjacency matrix

// 以DFS建立一棵交錯樹
bool DFS(int x)
{
for (int y=0; y<ny; ++y)
if (adj[x][y] && !vy[y])
{
vy[y] = true;

// 找到擴充路徑
if (my[y] == -1 || DFS(my[y]))
{
mx[x] = y; my[y] = x;
return true;
}
}
return false;
}

int bipartite_matching()
{
// 全部的點初始化為未匹配點。
memset(mx, -1, sizeof(mx));
memset(my, -1, sizeof(my));

// 依序把X中的每一個點作為擴充路徑的端點，
// 並嘗試尋找擴充路徑。
int c = 0;
for (int x=0; x<nx; ++x)
// if (mx[x] == -1) // x為未匹配點，這行可精簡。
{
// 開始Graph Traversal
memset(vy, false, sizeof(vy));
if (DFS(x)) c++;
}
return c;
}
main()
{
int a,b,i,n,k;
while(scanf("%d%d",&n,&k)!=EOF)
{
nx=ny=n;
memset(adj, 0, sizeof(adj));
for(i=1;i<=k;i++)
{
scanf("%d%d",&a,&b);
adj[a-1][b-1]=1;
}
printf("%d\n",bipartite_matching());
}
system("pause");
}