UVA 11419 SAM I AM（二分图最小点覆盖）

题目：http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2414
题目大意：给你一个 r*c 的方格图，有n个小方格上有小怪，你可以从边界上发射子弹，子弹会把你发射过去的那一行或者列的全部小怪都清掉。问你最少需要多少发子弹，可以把怪全清掉，然后输出具体方案。

解题思路：二分图的最小点覆盖。做二分图专题，在加上前面几题的操练，应该回去想把 行 当做左边，列 当做右边来建图，然后如果某行某列有小怪就相应的连边。那么问题就是你要选最少的点，使这些点能把所有的边都至少有一个点在你选的集合里。

        首先，何为最小点覆盖，其实就是：二分图，让你选出最少的点，使得所有的边都至少有一个端点在你选的那个集合里。

        先做二分图最大匹配，这样就不存在增广路了。左边每个未匹配点出发到右边肯定是一个匹配点，或者 它本身就是一个匹配点。
先来考虑前面一种，即未匹配点，从这个点出发是一棵匈牙利树，要使这上面的每条边都至少有一个端点被选中，我们选那些呢？画棵树看看（树一定为奇数个层），很明显，我们选那些在偶数层的点，把那些点全选了，那么树上的边肯定满足我们的条件，而且选的点还最少。这就相当于，我们走交替路，左边出发，我们选的路径上右边的点能把边都碰到，而且数目比其他至少-1。
        然后第二种情况，即匹配点。如果这个点在之前的所有匈牙利树里面都没有出现过，那么就选其中一个点，都可以。
那么这中选法，选出来的点是几个？它等于最大匹配数。这就是König定理。为什么？很简单，看看我们刚才的选法，就是每条匹配边选一个点嘛！然后其实如果对上面为什么那么选是最少的点数还有疑问，那么这里其实再一次反证了，假设有 n 条匹配边，我们最少也需要 n 个点吧，而它现在就是 n 个点。
        上面只是鄙人的一些短见，如有错误，还望大家指正。
        所以，综上：
        König定理：最小点覆盖数 = 最大匹配数。

        解方案为：从左侧所有未盖点出发拓展匈牙利树，把树上所有的点都标记，所求答案即为，左侧未标记的点和右侧标记的点。

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int MAXN = 1111;

bool w[MAXN][MAXN],vis_x[MAXN],vis_y[MAXN];
int s[MAXN],t[MAXN],left[MAXN];

int match(int i,int c)
{
    s[i] = 1;
    for(int j = 1;j <= c;j++)
        if(w[i][j] && !t[j])
        {
            t[j] = 1;
            if(!left[j] || match(left[j],c))
            {
                left[j] = i;
                return 1;
            }
        }
    return 0;
}

void solve(int r,int c)
{
    memset(left,0,sizeof(left));
    memset(vis_x,0,sizeof(vis_x));
    memset(vis_y,0,sizeof(vis_y));
    int ans = 0;
    for(int i = 1;i <= r;i++)
        for(int j = 1;j <= c;j++)
            if(w[i][j] && !vis_x[i] && !vis_y[j])
            {
                vis_x[i] = vis_y[j] = 1;
                left[j] = i;
                ans++;
            }
    for(int i = 1;i <= r;i++)
    {
        memset(s,0,sizeof(s));
        memset(t,0,sizeof(t));
        if(!vis_x[i] && match(i,c))
            ans++;
    }
    printf("%d",ans);
}

int right[MAXN];

void print(int r,int c)
{
    memset(s,0,sizeof(s));
    memset(t,0,sizeof(t));
    memset(right,0,sizeof(right));
    for(int i = 1;i <= c;i++)
        if(left[i])
        {
            right[left[i]] = 1;
        }
    for(int i = 1;i <= r;i++)
        if(!right[i])
        {
            match(i,c);
        }
    for(int i = 1;i <= r;i++)
        if(!s[i]) printf(" r%d",i);
    for(int i = 1;i <= c;i++)
        if(t[i]) printf(" c%d",i);
    puts("");
}

int main()
{
    int r,c,n;
    while(~scanf("%d%d%d",&r,&c,&n) && r)
    {
        memset(w,0,sizeof(w));
        while(n--)
        {
            int a,b;
            scanf("%d%d",&a,&b);
            w[a][b] = 1;
        }
        solve(r,c);
        print(r,c);
    }
    return 0;
}

/*
4 4 3
1 1
1 4
3 2

4 4 2
1 1
2 2

0 0 0

*/