POJ ~ 1904 ~ King's Quest （强连通缩点，思维）

题意

先输入N表示有n个王子和n个公主，接下来N行，每行先输入K表示，第i个王子有K个喜欢的公主，然后输入这K个公主的编号。每个王子只能和喜欢的妹子结婚，大臣给出一个匹配表，每个王子都和一个妹子结婚，但是国王不满意，他要求大臣给他另一个表，每个王子可以和几个公主结婚，按序号升序输出妹子的编号。这个表应满足所有的王子最终都有妹子和他结婚，即必须保证王子与这些对象中的任意一个结婚，都不会影响让剩余的王子中突然有一个人没婚可结了。

思路

王子编号为[0，n），公主编号为[n，n+n)。
王子u喜欢公主v，建一条u->(v+n)的边。公主u可以和王子v结婚建一条(u+n)->v的边。
求强连通并缩点，在同一个强连通分量中的王子和公主均可以结婚。
因为这样建图保证了，在某个强连通分量中王子和公主个数相等。因为是强连通分量，所以王子任意选择公主都没问题。

//#include <bits/stdc++.h>
#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <stack>
#include <algorithm>
using namespace std;
const int MAXN = 1e5+5;
struct Edge
{
    int from, to;
    Edge (int from, int to): from(from), to(to) {}
};
struct SCC
{
    int n, m;
    int DFN[MAXN], LOW[MAXN], sccno[MAXN], dfs_clock, scc_cnt;
    vector<Edge> edges;
    vector<int> G[MAXN];
    stack<int> S;
    void init(int n)
    {
        this->n = n, m = 0;
        edges.clear();
        for (int i = 0; i <= n; i++) G[i].clear();
    }
    void AddEdge (int from, int to)
    { 
        edges.push_back(Edge(from, to));
        m = edges.size();
        G[from].push_back(m-1);
    }
    void dfs(int u)
    {
        DFN[u] = LOW[u] = ++dfs_clock;
        S.push(u);
        for (int i = 0; i < G[u].size(); i++)
        {
            Edge e = edges[G[u][i]];
            int v = e.to;
            if (!DFN[v])
            {
                dfs(v);
                LOW[u] = min(LOW[u], LOW[v]);
            }
            else if (!sccno[v])
                LOW[u] = min(LOW[u], DFN[v]);
        }
        if (LOW[u] == DFN[u])
        {
            scc_cnt++;
            while (1)
            {
                int x = S.top(); S.pop();
                sccno[x] = scc_cnt;
                if (x == u) break;
            }
        }
    }
    void find_scc()
    {
        dfs_clock = scc_cnt  = 0;
        memset(DFN, 0, sizeof(DFN)), memset(sccno, 0, sizeof(sccno));
        for (int i = 0; i < n; i++)
            if (!DFN[i]) dfs(i);
    }

}gao;

int n;
int main()
{
    scanf("%d", &n);
    gao.init(n*2);
    for (int u = 0; u < n; u++)
    {
        int k; scanf("%d", &k);
        while (k--)
        {
            int v; scanf("%d", &v);
            v--;
            gao.AddEdge(u, v+n);
        }
    }
    for (int v = 0; v < n; v++)
    {
        int u; scanf("%d", &u);
        u--;
        gao.AddEdge(u+n, v); 
    }
    gao.find_scc();
    for (int u = 0; u < n; u++)
    {
        vector<int> ans;
        for (int i = 0; i < gao.G[u].size(); i++)
        {
            Edge e = gao.edges[gao.G[u][i]];
            int v = e.to;
            if (gao.sccno[u] == gao.sccno[v])
                ans.push_back(v-n+1);
        }
        sort(ans.begin(), ans.end());
        printf("%d ", ans.size());
        for (int i = 0; i < ans.size(); i++)
            printf("%d%c", ans[i], i==ans.size()-1?'\n':' ');
    }
    return 0;
}
/*
4
2 1 2
2 1 2
2 2 3
2 3 4
1 2 3 4
*/