HDOJ 3639.Hawk-and-Chicken（Tarjan强连通分量+缩点）

2015-8-10

问题简述：

有一群小孩要投票，有的两人之间存在单向的支持关系，求获得支持最多的人的支持数，并要求把这些人的编号（0~n-1）按升序排列输出来。

原题链接：http://acm.hdu.edu.cn/showproblem.php?pid=3639

解题思路：

我们可以采取反向建图的策略，这样支持数最大的节点只可能在入度为0的点中产生，只需要对所有入度为0的点dfs一遍计算它能到达的所有顶点数即可。但是还有一个问题，当遇到有向环时（也是一个强连通分量），dfs的过程中对于环中的最后一个点会产生错误的判断，显然这个环中的所有点状态完全一样，可以使用Tarjan的强连通分量进行缩点。

源代码：

/*
OJ: HDOJ
ID: forever
TASK: 3639.Hawk-and-Chicken
LANG: C++
NOTE: Tarjan(强连通分量+缩点)
*/
#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <stack>
#define MAX 10010
using namespace std;

struct edge {
    int from,to;
    int next;
}e[30010];
int head[MAX], cnt;

int dfn[MAX], low[MAX], vis[MAX], scc[MAX];
int ans[MAX], id[MAX], in[MAX], ins[MAX];
int time, num, sum;
stack<int> s;

void init() {
    memset( head, -1, sizeof(head) );
    //memset( e, 0, sizeof(e) );
    memset( vis, 0, sizeof(vis) );
    memset( ins, 0, sizeof(ins) );
    cnt = 0;
    time = 0;
}

void addedge(int u, int v) {
    e[cnt].from = u;
    e[cnt].to = v;
    e[cnt].next = head[u];
    head[u] = cnt ++;
}

void tarjan(int u) {
    dfn[u] = low[u] = time ++;
    s.push(u);
    ins[u] = 1;
    vis[u] = 1;
    for( int i = head[u]; i != -1; i = e[i].next ) {
        int v = e[i].to;
        if( !vis[v] ) {
            tarjan(v);
            low[u] = min( low[u], low[v] );
        }
        else if( ins[v] )
            low[u] = min( low[u], dfn[v] );

    }
    if( low[u] == dfn[u] ) {
        int v = -1;
        while( u != v ) {
            v = s.top();
            s.pop();
            id[v] = num;
            scc[num] ++;
            ins[v] = 0;
        }
        num ++;
    }
}

void dfs(int u) {
    vis[u] = 1;
    sum += scc[u];
    for( int i = head[u]; i != -1; i = e[i].next ) {
        int v = e[i].to;
        if( !vis[v] )
            dfs(v);
    }
}

int main()
{
    int t, n, m, x, y, cas = 1;
    scanf( "%d", &t );
    while( t -- ) {
        init();
        scanf( "%d%d", &n, &m );
        for( int i = 0; i < m; i ++ ) {
            scanf( "%d%d", &x, &y );
            addedge( x, y );
        }

        num = 0;
        while(!s.empty()) s.pop();
        memset( scc, 0, sizeof(scc) );
        for( int i = 0; i < n; i ++ )
            if( !vis[i] )
                tarjan(i);

        init();
        memset( in, 0, sizeof(in) );
        for( int i = 0; i < m; i ++ ) {
            int u = id[e[i].from];
            int v = id[e[i].to];
            if( u != v ) in[u] = 1;
            addedge( v, u );
        }

        int maxans = -1;
        memset( ans, -1, sizeof(ans) );
        for( int i = 0; i < num; i ++ ) {
            if( in[i] ) continue;
            sum = 0;
            memset( vis, 0, sizeof(vis) );
            dfs(i);
            ans[i] = sum;
            maxans = max( sum, maxans );
        }

        printf( "Case %d: %d\n", cas ++, maxans - 1 );
        int flag=0;
        for (int i = 0; i < n; i++)
            if ( ans[id[i]] == maxans) {
                if(!flag){
                printf("%d", i);
                flag=1;
                }
                else printf(" %d",i);
            }
        printf("\n");
    }
    return 0;
}