poj 2553 强连通

题意：给出一个有向图，定义：若节点v所有能到达的点{wi}，都能反过来到达v，那么称节点v是sink。题目要求所有的sink点。

思路：强连通缩点找出出度为零的点，输出即可。

这题主要问题是读题，了解题意之后就好做了，然后在数组开小了导致WA？挺莫名其妙的。。

代码：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
using namespace std;

#define MAXN 10500
#define MAXM 20000000

struct Edge
{
      int v, next;
}edge[MAXM];    //边结点数组

int first[MAXN], stack[MAXN], DFN[MAXN], Low[MAXN], Belong[MAXN];
int indegree[MAXN],outdegree[MAXN];
// first[]头结点数组,stack[]为栈,DFN[]为深搜次序数组，Belong[]为每个结点所对应的强连通分量标号数组
// Low[u]为u结点或者u的子树结点所能追溯到的最早栈中结点的次序号
int instack[MAXN];  // instack[]为是否在栈中的标记数组
int n, m, cnt, scnt, top, tot;

void init()
{
    cnt = 0;
    scnt = top = tot = 0;    //初始化连通分量标号，次序计数器，栈顶指针为0
    for(int i=0;i<=n+100;i++)
    {
        first[i]=-1;
        outdegree[i]=0;
        DFN[i]=0;
    }
}

void read_graph(int u, int v) //构建邻接表
{
     edge[tot].v = v;
     edge[tot].next = first[u];
     first[u] = tot++;
}
void Tarjan(int v)       //Tarjan算法求有向图的强连通分量
{
     int min, t;
     DFN[v] = Low[v] = ++tot;    //cnt为时间戳
     instack[v] = 1;    //标记在栈中
     stack[top++] = v;      //入栈
     for(int e = first[v]; e != -1; e = edge[e].next)
     {   //枚举v的每一条边
           int j = edge[e].v;   //v所邻接的边
           if(!DFN[j])
           {   //未被访问
               Tarjan(j);    //继续向下找
               if(Low[v] > Low[j]) Low[v] = Low[j];  // 更新结点v所能到达的最小次数层
           }
           else if(instack[j] && DFN[j] < Low[v])
           {   //如果j结点在栈内,
               Low[v] = DFN[j];
           }
     }
     if(DFN[v] == Low[v])
     {     //如果节点v是强连通分量的根
           scnt++;   //连通分量标号加1
           do
           {
               t = stack[--top];   //退栈
               instack[t] = 0;   //标记不在栈中
               Belong[t] = scnt;   //出栈结点t属于cnt标号的强连通分量
           }while(t != v);  //直到将v从栈中退出
     }
}

void solve()
{
     for(int i = 1; i <= n; i++)   //枚举每个结点，搜索连通分量
        if(!DFN[i])  //未被访问
           Tarjan(i);  //则找i结点的连通分量
}
int e1[MAXN];int e2[MAXN];
int main()
{
    while(scanf("%d",&n),n)
    {
        init();
        scanf("%d",&m);
        for(int i=1;i<=m;i++)
        {
            int v,w;
            scanf("%d%d",&v,&w);
            e1[tot]=v;
            e2[tot]=w;
            read_graph(v, w);
        }
        int num=tot;
        solve();     //求强连通分量
        for(int i=0;i<num;i++)
        {
            if(Belong[e1[i]]!=Belong[e2[i]])
                outdegree[Belong[e1[i]]]++;
        }
        bool ff=false;
        for(int i=1;i<=n;i++)
        {
            if(!outdegree[Belong[i]])
            {
                if(ff==false)
                {
                    printf("%d",i);
                    ff=true;
                }
                else
                {
                    printf(" %d",i);
                }
            }
        }
        printf("\n");
    }
    return 0;
}