UVA11419 我是SAM (最小点覆盖）

题意：给一个r*c的矩阵，某些格子中可能有一些怪物，可以在一行或一列防止一枚大炮，大炮会扫光整行/列的怪，问最少需要多少炮？输出炮的位置

思路：首先把每行看成x集合里的r个点，每列看成y集合里的c个点，然后就会发现对于每个目标的x，y值其实就可以抽象成了这个x点和这个y点通过目标连成了一条边，而同一行的x点是通用的，同一列的y点是通用的，这样不就变成了用最少的点使得所有的边至少有一端被覆盖了么， 最小点覆盖，好了，二分图；重点是打印解， 这个比较麻烦， 首先我们先从x集合里找到没有被覆盖的点， 然后从这些点开始标记并再次进行一遍匈牙利找到一条已经匹配过的增光路径， 然后把沿路的x集合里的点和y集合里的点都标记上， 之后一直重复把x集合里没有被覆盖的点都跑一边， 这样最后x集合里没有被标记的点和y集合里被标记的点就是我们要的答案了， 原因在于：已经求出了最小点覆盖了， 按理来说我们不管别的， 直接把x集合拿来是不是就能把目标全部干掉， 但是x集合点的个数是大于等于最小点覆盖个数的， 所以必然有重复， 重复怎么找呢， 因为没有匹配过的在y集合里一定会有y与之匹配，且这个y也一定已经匹配了x集合里的另外一个， 我们就利用这个性质找到这些y， 这样就可以去掉重复了， 至于x集合中没有标记的和刚才的道理一样， 也是要重复掉y里面没有匹配和匹配的。

感想：被这个题QiangJIan了一个晚上。。。痛苦。。。。就是想不明白打印解那个地方。。。最后还是看了大神代码才会。。。。明明自己已经刷了一套匹配专题了。。。结果一个假期过去。。。全忘了。。。心累。。。。

代码：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
const int maxed = 1000 + 10;
struct E
{
    int v, bef;
}e[maxed * maxed];
int r, c, n, ans, head[maxed], d[2][maxed];
bool vis_r[maxed], vis_c[maxed], vis[maxed];
int main()
{
    void add_(int x, int y);
    bool slove(int u);
    while (scanf("%d%d%d", &r, &c, &n) != EOF) {
        if (!r && !c && !n)
            break;
        ans = 1;
        memset(head, -1, sizeof(head));
        int a, b;
        for (int i = 1; i <= n; ++i) {
            scanf("%d%d", &a, &b);
            add_(a, b);
        }
        memset(d, -1, sizeof(d));
        int cnt = 0;
        for (int i = 1; i <= r; ++i)
            if (d[1][i] == -1) {
                memset(vis, false, sizeof(vis));
                if (slove(i))
                    cnt += 1;
            }
        printf("%d", cnt);
        memset(vis_r, false, sizeof(vis_r));
        memset(vis_c, false, sizeof(vis_c));
        for (int i = 1; i <= r; ++i)
            if (d[1][i] == -1) {
                memset(vis, false, sizeof(vis));
                slove(i);
            }
        for (int i = 1; i <= r; ++i)
            if (d[1][i] != -1 && !vis_r[i])
                printf(" r%d", i);
        for (int i = 1; i <= c; ++i)
            if (d[0][i] != -1 && vis_c[i])
                printf(" c%d", i);
        printf("\n");
    }
    return 0;
}
void add_(int x, int y)
{
    e[ans].v = y;
    e[ans].bef = head[x];
    head[x] = ans++;
}
bool slove(int u)
{
    for (int i = head[u]; i != -1; i = e[i].bef) {
        int v = e[i].v;
        if (d[0][v] != -1)
            vis_c[v] = true;
        if (d[1][u] != -1)
            vis_r[u] = true;
        if (!vis[v]) {
            vis[v] = true;
            if (d[0][v] == -1 || slove(d[0][v])) {
                d[0][v] = u;
                d[1][u] = v;
                return true;
            }
        }
    }
    return false;
}