HDU-1498 50 years, 50 colors（最小点覆盖）

题意

一个 $n \times n$ 的矩阵，有很多不同颜色的气球，每次操作可以选择某一种颜色，再选择某一行（列），使得这一行（列）的这种颜色的气球撞裂。求有多少种颜色的气球无法在 $k$ 次之内全部撞裂。
$1 \leq n \leq 100$
$1 \leq k \leq n$

思路

先引出一个“最小点覆盖”的概念。选取一个点集 $S$ 使得二分图边集 $E$ 满足 $\forall (u,v)\in E$ 有 $u\in S$ 或 $v\in S$ ，则称 $|S|$ 为最小点覆盖。
有一个基本结论：最小点覆盖 $=$ 最大匹配
证明可以采用归纳法。当对任意一张二分图（包括空二分图）加入一个节点时，两者值均不变。而加入一条边 $(u,v)$ 时分两种情况：
$1.$ 最大匹配加一，说明存在一种最大匹配，满足 $u,v$ 均不为匹配点。那就需要选择 $u,v$ 中的其中一点覆盖这条边，最小点覆盖也加一。
$2.$ 最大匹配不变，说明原来的最大匹配， $u$ 或 $v$ 是匹配点，也就是说，新加的这条边将被覆盖，最小点覆盖不变。
得证。
有了上述定理，这道题的建图就比较明显了。转化所求为最少撞多少次可以撞灭所有气球。先枚举颜色，对于一个位于 $(x,y)$ 的气球，可以在 $x$ 行上撞，也可以在 $y$ 列上撞。所以对于每一行建立一个集合，每一列建立另一个集合，连边的条件是交点上有一个这种颜色的气球，然后只用覆盖所有的边就行了，求出最大匹配，就是最小点覆盖。

代码

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define FOR(i,x,y) for(int i=(x);i<=(y);i++)
#define DOR(i,x,y) for(int i=(x);i>=(y);i--)
#define N 103
typedef long long LL;
using namespace std;
template<const int maxn,const int maxm>struct Linked_list
{
    int head[maxn],to[maxm],nxt[maxm],tot;
    void clear(){memset(head,-1,sizeof(head));tot=0;}
    void add(int u,int v){to[++tot]=v,nxt[tot]=head[u];head[u]=tot;}
    #define EOR(i,G,u) for(int i=G.head[u];~i;i=G.nxt[i])
};
Linked_list<N,N*N>G;
int mark[N],mc[N],mp[N][N];
int n,k;

bool match(int u,int stmp)
{
    EOR(i,G,u)
    {
        int v=G.to[i];
        if(mark[v]==stmp)continue;
        mark[v]=stmp;
        if(!mc[v]||match(mc[v],stmp))
        {
            mc[v]=u;
            return true;
        }
    }
    return false;
}
bool solve(int col)
{
    int ans=0;bool flag=0;
    G.clear();
    FOR(i,1,n)FOR(j,1,n)if(mp[i][j]==col)
    {
        G.add(i,j);
        flag=1;
    }
    if(!flag)return false;
    memset(mark,0,sizeof(mark));
    memset(mc,0,sizeof(mc));
    FOR(i,1,n)if(match(i,i))ans++;
    return ans>k;
}

int main()
{
    while(scanf("%d%d",&n,&k),n||k)
    {
        bool flag=0;
        FOR(i,1,n)FOR(j,1,n)scanf("%d",&mp[i][j]);
        FOR(i,1,50)if(solve(i))
        {
            if(flag)printf(" ");
            flag=1;
            printf("%d",i);
        }
        if(!flag)printf("-1");
        printf("\n");
    }
    return 0;
}