HDU - 4685（完美匹配构图+强连通）

题解:做这题前建议先去做poj1904他的题目已经是完美匹配了直接强连通算，然后这题难在构造完美匹配的图，如何构造呢？

先求原始图的最大匹配cnt数，然后加入m-cnt个王子，这些王子喜欢所有的公主

然后加入n-cnt个公主，这些公主被所有的王子喜欢

然后再求一次完美匹配那么就是按poj1904的方法:

如果王子（u）喜欢这个公主（v）就建一个u->v的有向边，如果公主可以嫁给这个王子就建一个v->u的有向边

图画出来后你会发现公主和王子如果是一个强连通的，那么这个这个王子可以和这个强联通图里面的任意匹配而且这个图中的王子数量和公主数量一样。

求完强联通后如果这个王子喜欢这个公主并且在同一个强连通就可以加进答案。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<stack>
#include<queue>
#include<vector>
using namespace std;
const int inf = 0x3f3f3f3f;
const int mx = 2005;
int dis,dx[mx],dy[mx],vis[mx];
int x[mx],y[mx],head[mx],ans[mx];
int pre[mx],low[mx],sccno[mx];
int cnt,tot,dfn;
stack<int>st;
vector<int>g[mx];
int n,m,nv,nu;
struct node{
    int v,next;
}E[mx*mx];
void init(){
    for(int i = 1; i <= n+m; i++) g[i].clear();
    memset(head,0,sizeof(head));
    memset(x,0,sizeof(x));
    memset(y,0,sizeof(y));
    memset(pre,0,sizeof(pre));
    memset(sccno,0,sizeof(sccno));
    dfn = cnt = tot = 0;
}
void add(int u,int v){
    tot++;
    E[tot].v = v;
    E[tot].next = head[u];
    head[u] = tot;
}
bool search(){
    dis = inf;
    queue<int>q;
    memset(dx,0,sizeof(dx));
    memset(dy,0,sizeof(dy));
    for(int u = 1; u <= nu; u++)
        if(!x[u]){
            dx[u] = 1;
            q.push(u);
        }
    while(!q.empty()){
        int u = q.front();
        q.pop();
        if(dx[u]>dis)  break;
        for(int i = 0; i < g[u].size(); i++){
            int v = g[u][i];
            if(!dy[v]){
                dy[v] = dx[u]+1;
                if(!y[v]) dis = dy[v];
                else{
                    dx[y[v]] = dy[v]+1;
                    q.push(y[v]);
                }
            }
        }
    }
    return dis!=inf;
}
bool find(int u){
    for(int i = 0; i < g[u].size(); i++){
        int v = g[u][i];
        if(!vis[v]&&dy[v] == dx[u]+1){
            vis[v] = 1;
            if(y[v]&&dis == dy[v])  continue;
            if(!y[v]||find(y[v])){
                x[u] = v;
                y[v] = u;
                return true;
            }
        }
    }
    return false;
}
int MaxMatch(){
    int ans = 0;
    while(search()){
        memset(vis,0,sizeof(vis));
        for(int u = 1; u <= nu; u++)
            if(!x[u]&&find(u))  ans++;
    }
    return ans;
}
void tarjan(int u){
    pre[u] = low[u] = ++dfn;
    st.push(u);
    for(int i = head[u]; i; i = E[i].next){
        int v = E[i].v;
        if(!pre[v]){
            tarjan(v);
            low[u] = min(low[v],low[u]);
        }
        else if(!sccno[v])
            low[u] = min(low[v],low[u]);
    }
    if(pre[u] == low[u]){
        cnt++;
        while(1){
            int x = st.top();
            st.pop();
            sccno[x] = cnt;
            if(x == u)  break;
        }
    }
}
void solve(){
    nu = n;
    nv = m;
    int cnt = MaxMatch();
    nu = nv = n+m-cnt;
    for(int u = 1; u <= n; u++)
        for(int v = m+1; v <= nv; v++)
            g[u].push_back(v);
    for(int u = n+1; u <= nu; u++)
        for(int v = 1; v <= m; v++)
            g[u].push_back(v);
    MaxMatch();
    for(int u = 1; u <= nu; u++)
        for(int i = 0; i < g[u].size(); i++){
            int v = g[u][i];
            add(u,v+nu);
        }
    for(int v = 1; v <= nv; v++)
        add(nu+v,y[v]);
    for(int u = 1; u <= 2*nu; u++)
        if(!pre[u])
            tarjan(u);
    for(int u = 1; u <= n; u++){
        int len = 0;
        for(int i = 0; i < g[u].size(); i++){
            int v = g[u][i];
            if(v>m)
                break;
            if(sccno[v+nu] == sccno[u])
                ans[len++] = v;
        }
        sort(ans,ans+len);
        printf("%d",len);
        for(int i = 0; i < len; i++)
            printf(" %d",ans[i]);
        printf("\n");
    }
}
int main(){
    int t;
    scanf("%d",&t);
    for(int casei = 1; casei <= t; casei++){
        scanf("%d%d",&n,&m);
        init();
        for(int u = 1; u <= n; u++){
            int num,v;
            scanf("%d",&num);
            for(int i = 1; i <= num; i++){
                scanf("%d",&v);
                g[u].push_back(v);
            }
        }
        printf("Case #%d:\n",casei);
        solve();
    }
    return 0;
}