例题5.27 我是SAM UVa11419

1.题目描述：点击打开链接

2.解题思路：本题就是经典的二分图最小覆盖问题。即选择尽量少的点，使得每条边至少有一个端点被选中。可以证明，最小覆盖数=最大匹配数。建模的方法为：将每一行看做一个X结点，每一列看做一个Y结点，每个目标对应了一条边，那么如果可以成功求出最小覆盖，也就是所有的边都被选中了，即所有的目标都被覆盖了。

那么现在的问题是：如何求一个最小覆盖呢？可以从匈牙利树出发。从X结点中的所有未盖点出发，扩展匈牙利树。标记树中的所有点。则X中未标记的点和Y中已经标记的点组成的就是最小覆盖。

3.代码：

#include<iostream>
#include<algorithm>
#include<cassert>
#include<string>
#include<sstream>
#include<set>
#include<bitset>
#include<vector>
#include<stack>
#include<map>
#include<queue>
#include<deque>
#include<cstdlib>
#include<cstdio>
#include<cstring>
//#include<cmath>
#include<ctime>
#include<cctype>
#include<list>
#include<complex>
#include<functional>
using namespace std;

#define me(s) memset(s,0,sizeof(s))
#define rep(i,n) for(int i=0;i<(n);i++)
#define pb push_back
typedef long long ll;
typedef pair <int,int> P;


const int N=1000+5;

struct BPM //求二分图的最小覆盖
{
    int n,m;
    vector<int>g[N];
    int Left[N];
    bool T[N];
    int Right[N];
    bool S[N];

    void init(int n,int m)
    {
        this->n=n;
        this->m=m;
        for(int i=0;i<n;i++)g[i].clear();
    }
    void addedge(int u,int v)
    {
        g[u].push_back(v);
    }
    bool match(int u)
    {
        S[u]=true;
        for(int i=0;i<g[u].size();i++)
        {
            int v=g[u][i];
            if(!T[v])
            {
                T[v]=true;
                if(Left[v]==-1||match(Left[v]))
                {
                    Left[v]=u;
                    Right[u]=v;
                    return true;
                }
            }
        }
        return false;
    }

    int solve()
    {
        memset(Left,-1,sizeof(Left));
        memset(Right,-1,sizeof(Right));
        int ans=0;
        for(int u=0;u<n;u++)
        {
            memset(S,0,sizeof(S));
            memset(T,0,sizeof(T));
            if(match(u))ans++;
        }
        return ans;
    }

    int mincover(vector<int>&X,vector<int>&Y)
    {
        int ans=solve();
        memset(S,0,sizeof(S));
        memset(T,0,sizeof(T));
        for(int u=0;u<n;u++)
            if(Right[u]==-1)match(u); //从X结点中的未盖点出发继续扩展
        for(int u=0;u<n;u++)
            if(!S[u])X.push_back(u);//X中的未标记的点放入X数组
        for(int v=0;v<m;v++)
            if(T[v])Y.push_back(v);//Y中已标记的点放入Y数组
        return ans;
    }
};

BPM solver;

int R,C,n;

int main()
{
    int kase=0;
    while(~scanf("%d%d%d",&R,&C,&n)&&(R||C||n))
    {
        solver.init(R,C);
        for(int i=0;i<n;i++)
        {
            int r,c;
            scanf("%d%d",&r,&c);r--,c--;
            solver.addedge(r,c);
        }
        vector<int>X,Y;
        int ans=solver.mincover(X,Y);
        printf("%d",ans);
        for(int i=0;i<X.size();i++)
            printf(" r%d",X[i]+1);
        for(int i=0;i<Y.size();i++)
            printf(" c%d",Y[i]+1);
        puts("");
    }
}