二分图常见模型

最小点覆盖

最小的点集使得其相连的边能覆盖所有边。

König定理：最小点覆盖=最大匹配

证明很简单，在得到最大匹配之后，一条边必然只有两种情况，两点都是匹配点，只有一个是匹配点。又因为对于一条匹配边来说，其两点上不可能同时连有非匹配点（即第二种情况的边），那么只需要把匹配边上两点的其中一点（含有第二种情况的点）包含到点集里，其另一点就可以通过选择的这个点来覆盖了（因为匹配边之间都是独立的）。

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

矩阵上每一个点只需要选中一次，那么对应从行部连到列部的一条边，把所有矩阵上的点消除就需要用最少的点来覆盖这些边，所以匈牙利跑一边即可。

类似的一道题：
UVA 11419 SAM I AM



只不过要输出方案，其他方法不好处理，那么可以用S T 数组记录访问到的点，那么求出最大匹配之后，再从X部的非匹配点出发进行匹配，此时不影响答案，但是过程中访问到的点是必须要选的，不然这个非匹配点无法覆盖，另外不用考虑访问过程中访问到Y部的非匹配点，因为这种情况根本不存在。然后对于Y部非匹配点同样不会从X部匹配点出发，所以这样找出X中没有访问过，Y中访问到的店即使答案。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<iomanip>
#include<ctime>
#include<climits>
#include<cctype>
#include<algorithm>
#ifdef WIN32
#define AUTO "%I64d"
#else
#define AUTO "%lld"
#endif
using namespace std;
#define smax(x,tmp) x=max((x),(tmp))
#define smin(x,tmp) x=min((x),(tmp))
#define maxx(x1,x2,x3) max(max(x1,x2),x3)
#define minn(x1,x2,x3) min(min(x1,x2),x3)
const int INF=0x3f3f3f3f;
const int maxn = 1005;
struct Edge
{
    int to,next;
}edge[maxn*maxn];
int head[maxn];
int maxedge;
inline void addedge(int u,int v)
{
    edge[++maxedge] = (Edge) { v,head[u] };
    head[u] = maxedge;
}
int r,c,n;
bool init()
{
    if(!~scanf("%d%d%d",&r,&c,&n) || !(r|c|n)) return false;
    memset(head,-1,sizeof(head));
    maxedge = -1;
    for(int i=1;i<=n;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        addedge(x,y);
    }
    return true;
}
bool S[maxn],T[maxn];
int matchx[maxn],matchy[maxn];
bool dfs(int u)
{
    S[u]=true;
    for(int i=head[u];~i;i=edge[i].next)
    {
        int v = edge[i].to;
        if(T[v]) continue;
        T[v]=true;
        if(!matchy[v] || dfs(matchy[v]))
        {
            matchx[u]=v;
            matchy[v]=u;
            return true;
        }
    }
    return false;
}
int hungary()
{
    int ret=0;
    memset(matchx,0,sizeof(matchx));
    memset(matchy,0,sizeof(matchy));
    for(int i=1;i<=r;i++)
    {
        memset(S,0,sizeof(S));
        memset(T,0,sizeof(T));
        if(dfs(i)) ret++;
    }
    return ret;
}
void work()
{
    int ans = hungary();
    printf("%d",ans);
    memset(S,0,sizeof(S));
    memset(T,0,sizeof(T));
    for(int i=1;i<=r;i++) <