hdu 3639 Hawk-and-Chicken(求被指（直接+间接）最多的点有哪些and值)

题目链接

题意：
n个人搞投票，m个关系a->b，表示a投票给b。
a投票给b，b投票给c，表明a也投票给c。
获得票数最多的那些人当老大，现在问你票数最多是多少，以及有哪些人。
也就是求图中点被指（直接+间接）的最大值+集合。

思路：
强连通缩点，得到dag，然后反向建边，现在只要走入度为0的点就行了，因为入度为0肯定是最优的，如果你还有入度，那入度连接的那个点肯定走的至少比你多一个（本身），现在问题来了，我一开始用拓扑排序求点的个数，类似与求dag最长路那种思路，设个dp数组，然后遍历图的过程中累加，这种思路是错的，为什么？

比如实说上面那张图，拓扑下来4加上的值是（1+2）和（1+3），1被算了两次，重复了，所以不能拓扑，看了别人的题解是对每个入度为0的点dfs，什么，那不应该是n*n的复杂度吗（但跑不满吧），看了一下n，5000，t，50，好吧，应该没问题，每次dfs到一个团，加上的值是他团里面的元素，对每个入度为0的点dfs，维护个最大值，也就是答案要求的获得最多的票数，然后遍历1~n（输入点时，我都+1），如果他的团的值是最大值就算入集合（一个强连通分量中每个人获得的票数是相同的），然后就没了。

#include <bits/stdc++.h>

#define eb emplace_back
#define mp make_pair
#define mt make_tuple
#define fi first
#define se second
#define pb push_back
#define all(x) (x).begin(), (x).end()
#define rall(x) (x).rbegin(), (x).rend()
#define forn(i, n) for (int i = 0; i < (int)(n); ++i)
#define for1(i, n) for (int i = 1; i <= (int)(n); ++i)
#define ford(i, a, b) for (int i = (int)(a); i >= (int)b; --i)
#define fore(i, a, b) for (int i = (int)(a); i <= (int)(b); ++i)
#define rep(i, l, r) for (int i = (l); i <= (r); i++)
#define per(i, r, l) for (int i = (r); i >= (l); i--)
#define ms(x, y) memset(x, y, sizeof(x))
#define SZ(x) ((int)(x).size())

using namespace std;

typedef pair<int, int> pii;
typedef vector<int> vi;
typedef vector<pii> vpi;
typedef vector<vi> vvi;
typedef long long i64;
typedef vector<i64> vi64;
typedef vector<vi64> vvi64;
typedef pair<i64, i64> pi64;
typedef double ld;

template<class T> bool uin(T &a, T b) { return a > b ? (a = b, true) : false; }
template<class T> bool uax(T &a, T b) { return a < b ? (a = b, true) : false; }

const int maxn = 5000+100;   //点数
int scc, top, tot;
vector<int> G[maxn];
vector<int> dag[maxn];
int low[maxn], dfn[maxn], belong[maxn];
int stk[maxn], vis[maxn];
int in[maxn], maxx, dis[maxn], sum, flag[maxn], sz[maxn];
void init(int n) {
    for (int i = 0; i <= n; ++i) {
        G[i].clear();
        dag[i].clear();
        low[i] = 0;
        dfn[i] = 0;
        stk[i] = 0;
        vis[i] = 0;
        in[i] = 0;
        dis[i] = 0;
        sz[i] = 0;
    }
    maxx = 0;
    scc = top = tot = 0;
}
void tarjan(int x) {
    stk[top++] = x;
    low[x] = dfn[x] = ++tot;
    vis[x] = 1;
    for (int to : G[x]) {
        if (!dfn[to]) {
            tarjan(to);
            low[x] = min(low[x], low[to]);
        } else if (vis[to]) {
            low[x] = min(low[x], dfn[to]);
        }
    }
    if (low[x] == dfn[x]) {
        ++scc;
        int temp;
        do {
            temp = stk[--top];
            belong[temp] = scc;
            vis[temp] = 0;
        } while (temp != x);
    }
}
int tc, n, m;
void dfs(int x) {
    sum += sz[x];
    flag[x] = 1;
    for (int to : dag[x]) {
        if (!flag[to]) {
            dfs(to);
        }
    }
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.precision(10);
    cout << fixed;
#ifdef LOCAL_DEFINE
    freopen("input.txt", "r", stdin);
#endif

    cin >> tc;
    int kase = 1;
    while (tc--) {
        cin >> n >> m;
        init(n);
        forn(i, m) {
            int u, v;
            cin >> u >> v;
            ++u; ++v;
            G[u].eb(v);
        }
        for1(i, n) if (!dfn[i]) tarjan(i);
        int ans = -1;
        for (int i = 1; i <= n; ++i) {
            ++sz[belong[i]];
            for (int j : G[i]) {
                if (belong[i] == belong[j]) continue;
                dag[belong[j]].eb(belong[i]);
                ++in[belong[i]];
            }
        }
        for1(i, scc) {
            dis[i] = sz[i]-1;
            if (!in[i]) {
                ms(flag, 0);
                sum = 0;
                dfs(i);
                dis[i] += sum-sz[i];
            }
            maxx = max(maxx, dis[i]);
        }
        vi res;
        for1(i, n) if (maxx == dis[belong[i]]) res.eb(i-1);
        cout << "Case " << kase++ << ": " << maxx << '\n';
        for (int i = 0; i < SZ(res); ++i) {
            if (i == SZ(res)-1) cout << res[i] << '\n';
            else cout << res[i] << ' ';
        }
    }

#ifdef LOCAL_DEFINE
    cerr << "Time elapsed: " << 1.0 * clock() / CLOCKS_PER_SEC << " s.\n";
#endif
    return 0;
}