poj3041 Asteroids 最小点覆盖 二分图 匈牙利算法BFS

                                 POJ 3041

大致题意：

题意比较好理解，在一个nxn的矩阵中，有几个点，你可以任意选择一个点，然后将它所在行后者所在列的所有点都清空。问你最少选几个点，才能将所有点清空掉？

大致思路：

将x的值和y的值进行划分，再根据点的坐标进行连接，形成一个二分图。那么这道题的本质就是，寻求最小的一个点集合，使得所有的边都与这个集合的点有关系。

最小点覆盖问题~~~根据konig定理：最小点覆盖数=最大匹配数。

代码：

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<vector>
#include<queue>
using namespace std;
int vis[1505],pre[1505];
int ml[1505],mr[1505];
vector<int>V[1505];
int n,k;
int MinCov()
{
    int ans=0;
    memset(vis,0,sizeof(vis));
    memset(ml,-1,sizeof(ml));
    memset(mr,-1,sizeof(mr));
    for(int i=1; i<=n; i++)
    {
        if(ml[i]==-1)
        {
           queue<int>Q;
           Q.push(i);
           pre[i]=-1;
           int flag=0;
           while(!Q.empty()&&!flag)
           {
               int st=Q.front();
               Q.pop();
               for(int j=0;j<V[st].size()&&!flag;j++)
               {
                   int to=V[st][j];
                   if(vis[to]!=i)
                   {
                       vis[to]=i;
                       Q.push(mr[to]);
                       if(mr[to]>=0)
                       {
                           pre[mr[to]]=st;
                       }
                       else
                       {
                           flag=1;
                           int d=st,e=to;
                           while(d!=-1)
                           {
                               int temp=ml[d];
                               ml[d]=e;
                               mr[e]=d;
                               d=pre[d];
                               e=temp;
                           }
                       }
                   }
               }
           }
          if(flag)ans++;
        }
    }
    return ans;
}
int main()
{

    while(~scanf("%d%d",&n,&k))
    {
        for(int i=1; i<=n; i++)V[i].clear();
        int x,y;
        for(int i=0; i<k; i++)
        {
            scanf("%d%d",&x,&y);
            V[x].push_back(y);
        }
        printf("%d\n",MinCov());
    }
}