HDU1498 枚举+二分图类棋盘问题（最大匹配）

/*HDU1498
题目大意:
给出N*N（100以内）的矩阵，矩阵上填气球的标号（1--50），给定 K
问，哪些编号的气球是不能一次在K次以内拿完的。每次拿取的过程是：选择一行或者一列，再在其中选一个气球，其余的下一次不能再选。

思考：
变种的棋盘问题。
但是要枚举气球编号(分类讨论).
这样，就是找到最大的匹配，即最多能放置互不在同一行或一列的气球。
这个匹配数小于等于K，则气球能全部被拿完。
复杂度为：O(50*200*200)
*/
#include <iostream>
#include <cmath>
#include <algorithm>
#include <string.h>
#include <stdio.h>
#include <set>
#include <stack>
#include <vector>
#define maxn 210
using namespace std;
vector<int> G[maxn];
int match[maxn];
bool visit[maxn];
int ans[55];
int mat[maxn][maxn];
int cnt;
int N,K;
void read(){
    for(int i=1;i<=N;i++){
        for(int j=1;j<=N;j++) scanf("%d",&mat[i][j]);
    }
}
void init(int k){
    for(int i=1;i<=2*N;i++) G[i].clear();
    for(int i=1;i<=N;i++){
        for(int j=1;j<=N;j++){
            if (mat[i][j]==k){
                int a=i,b=j+N;
                G[a].push_back(b);
                G[b].push_back(a);
            }
        }
    }
}
bool dfs(int s)//找到从s点出发的可增广路
{
    for(int i=0;i<G[s].size();i++){
        int v=G[s][i];
        if (!visit[v]){
            visit[v]=true;
            if (match[v]==0 || dfs(match[v])){//说明v对应的项有增广路
                match[v]=s;//修改v的对应项（即互斥点）是s
                return true;
            }
        }
    }
    return false;
}

int hungary(){
    int ans=0;
    memset(match,0,sizeof(match));//清空匹配
    for(int i=1;i<=2*N;i++){
        memset(visit,0,sizeof(visit));//注意清空增广路
        if (dfs(i)) ans++;
    }
    return ans/2;
}
void solve(){
    cnt=0;
    for(int i=1;i<=50;i++){
        init(i);
        if (hungary()>K) ans[cnt++]=i;
    }
}
int main(){
    while(cin>>N>>K){
        if (!N && !K) break;
        read();
        solve();
        if (cnt==0) puts("-1");
        else {
            for(int i=0;i<cnt;i++)
            if (i==cnt-1) printf("%d
",ans[i]);else printf("%d ",ans[i]);
        }
    }
}