POJ-2553-The Bottom of a Graph 解题报告

       强连通分量题。题目神描述，要锻炼题意理解能力的话此题为上等。题意：在有向图G中，给这样一个新的定义：对于在G中任何一个点v可达的点w，w都可达v，那么点v是一个sink。G图的bottom子图是由G图中所有的sink点构成，请按照顺序输出G图对应的bottom子图中的所有点编号，如果没有sink点，那么输出一个空行。Orz，题意描述得我好辛苦。

       我的解题思路：首先如果v可达w，那么必须w可达v，v才算是sink点。那么很明显，v可达的所有点都必须可达v，那么v才算是sink点。一个强连通分量中的任意点都是相互可达的，当我们把图缩点成一个DAG之后，只有出度为0的强连通分量里面的点才满足sink点的性质。在一个DAG中至少存在一个出度为0的点，所以题目要求没有sink点时输出空行这个情况是绝对不会发生的。最后，我们按顺序把出度为0的强连通分量里面的点全部按照升序输出就可以了。

       我的解题代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <cmath>
#include <algorithm>
#include <vector>
#include <queue>

using namespace std;

const int N = 5001;

vector<int> e[N];
int dfn[N], low[N];
int stack[N], stop;
bool instack[N];
int belong[N];
int outdeg[N];      //存储强连通分量的出度
int timer, sccn;
int n, m;

void InitRead();

void DataProcess();

void Tarjan(int x);

int main()
{
    while (~scanf("%d", &n))
    {
        if (n == 0) break;
        InitRead();
        DataProcess();
    }
    return 0;
}

void InitRead()
{
    timer = sccn = stop = 0;
    memset(dfn, 0, sizeof(dfn));
    memset(instack, false, sizeof(instack));
    memset(outdeg, 0, sizeof(outdeg));
    scanf("%d", &m);
    for (int i=1; i<=n; ++i)
    {
        e[i].clear();
    }
    int a, b;
    for (int i=0; i<m; ++i)
    {
        scanf("%d %d", &a, &b);
        e[a].push_back(b);
    }
    return;
}

void DataProcess()
{
    for (int i=1; i<=n; ++i)
    {
        if (!dfn[i]) Tarjan(i);
    }
    for (int i=1; i<=n; ++i)    //计算强连通分量出度
    {
        int size = e[i].size();
        for (int j=0; j<size; ++j)
        {
            if (belong[i] != belong[e[i][j]])
            {
                outdeg[belong[i]]++;
            }
        }
    }
    for (int i=1; i<=n; ++i)    //按升序输出出度为0的强连通分量里的点
    {
        if (outdeg[belong[i]] == 0)
        {
            printf("%d", i++);
            for (; i<=n; ++i)
            {
                if (outdeg[belong[i]] == 0)
                {
                    printf(" %d", i);
                }
            }
            break;
        }
    }
    putchar('\n');
    return;
}

void Tarjan(int x)
{
    dfn[x] = low[x] = ++timer;
    stack[stop++] = x;
    instack[x] = true;
    int size = e[x].size();
    int y;
    for (int i=0; i<size; ++i)
    {
        y = e[x][i];
        if (!dfn[y])
        {
            Tarjan(y);
            low[x] = min(low[x], low[y]);
        }
        else if (instack[y])
        {
            low[x] = min(low[x], dfn[y]);
        }
    }
    if (dfn[x] == low[x])
    {
        sccn++;
        do
        {
            y = stack[--stop];
            instack[y] = false;
            belong[y] = sccn;
        } while (x != y);
    }
    return;
}