[POJ3041] Asteroids（最小点覆盖-匈牙利算法）

传送门

题意：

给一个N*N的矩阵，有些格子有障碍，要求我们消除这些障碍，问每次消除一行或一列的障碍，最少要几次。

 

解析：

把每一行与每一列当做二分图两边的点。

某格子有障碍，则对应行与列连边。

选出最少的点，使得所有边被覆盖。

最小点覆盖。

——代码

 #include <cstdio>
 #include <cstring>
 #define M(x, a) memset(a, x, sizeof(a))

 using namespace std;

 const int MAXN = ;
 int n, k, cnt, ans;
 int head[MAXN], next[MAXN * MAXN], to[MAXN * MAXN], belong[MAXN];
 bool vis[MAXN];

 inline void add(int x, int y)
 {
     to[cnt] = y;
     next[cnt] = head[x];
     head[x] = cnt++;
 }

 inline bool find(int u)
 {
     int i, v;
     for(i = head[u]; i != -; i = next[i])
     {
         v = to[i];
         if(!vis[v])
         {
             vis[v] = ;
             if(!belong[v] || find(belong[v]))
             {
                 belong[v] = u;
                 return ;
             }
         }
     }
     return ;
 }

 int main()
 {
     int i, x, y;
     scanf("%d %d", &n, &k);
     M(-, head);
     for(i = ; i <= k; i++)
     {
         scanf("%d %d", &x, &y);
         add(x, y);
     }
     for(i = ; i <= n; i++)
     {
         M(, vis);
         if(find(i)) ans++;
     }
     printf("%d", ans);
     return ;
 }