FZOJ1638求强连通分量

FZOJ1638求强连通分量
描述
输入一个图，输出该图中的最大强连通分量。
输入
第一行：n和m（n<=10000,m<=100000,n为节点个数，m为边的条数）
接下来m行，每行两个数：a，b，表示a指向b的边（a,b为非负整数）；
输出
输出最大强连通分量的节点，按照节点编号从小到大输出，如果有多个强连通分量节点数相同，则输出节点编号字典序较小的。
样例输入
6 8
1 3
3 5
5 6
1 2
4 1
2 4
4 6
3 4
样例输出
1 2 3 4

强连通分量是由若干个环组成的。所以，当有环形成时（也就是搜索的下一个点已在栈中），我们将这一条路径的low值统一，即这条路径上的点属于同一个强连通分量。
如果遍历完整个搜索树后某个点的dfn值等于low值，则它是该搜索子树的根。这时，它以上（包括它自己）一直到栈顶的所有元素组成一个强连通分量。

#include <cstdio>
#include<algorithm>
using namespace std;
int n,m,a,b;
struct node{int v,first;}no[10010];
struct edge{int v,next;}e[100010];
int cnt=0;//边的计数器 
void addedge(int a,int b){
    cnt++;
    e[cnt].v=b;
    e[cnt].next=no[a].first;
    no[a].first=cnt;
}

int low[10010],dfn[10010],scc[10010],st[10010],ins[10010],cscc=1,index=1,top=1;
//low数组是该节点所在强连通子图上根节点的dfn值
//dfn数组是发现该节点的时间，用递加的index表示
//scc用来记录每个节点所在强连通子图的标号，及cscc
//ins记录是否在栈中 st即为栈 top栈顶 
void tarjan(int rt){
    low[rt]=dfn[rt]=index++;//刚发现该节点时low与dfn相等 
    st[top++]=rt;//入栈 
    ins[rt]=1;//记录他在栈中 
    int v;
    for(int i=no[rt].first;i!=0;i=e[i].next){//找到从他出去的点 
        v=e[i].v;
        if(!dfn[v]){//如果还没有被发现 
            tarjan(v);//带入递归里 
            low[rt]= low[rt]<low[v]?low[rt]:low[v];//如果有环，这颗树的最后会发现一个low比low【rt】小的根节点
                                            //以小的那个为low 

        }else if(ins[v]==1){//如果v已经被发现过了，那我们又发现了一个环 
            low[rt]=low[rt]<dfn[v]?low[rt]:dfn[v];
        }
    }
    if(low[rt]==dfn[rt]){//如果low等于dfn，要么是一个点的强连通子图，要么是一个环的根节点 
        do{
            v=st[--top];//把新发现的强连通子图出栈 
            ins[v]=0;//标记这些节点不在栈中 
            scc[v]=cscc;//记录这个强连通子图  v是一个节点 cscc是他们的标记 
        }while(v!=rt);
        cscc++;
    }
}
int cc[10010],j,i,k=0,max_o=0,max_x,ans[10010];
void out(int n)
{
    for(i=1;i<=cscc-1;i++)//循环cscc-1个强连通子图 
    {
       for(j=1;j<=n;j++)//j写成了i 心好累 ，循环每一个节点 
       {if(scc[j]==i)cc[i]++;}}// 记录每一个 强连通子图包含节点的大小 
    for(i=1;i<cscc;i++)
    {if(cc[i]>max_o){max_o=cc[i];max_x=i;}}//把最大的那个挑出来 max_x是强连通子图的标号 
    for(i=1;i<=n;i++){if(scc[i]==max_x)ans[++k]=i;}//==用了= 又很心酸     i是节点 如果i在最大的强连通子图中开个数组把他们记录下来 
    sort(ans+1,ans+k+1);//整理一下 顺序输出 
    for(i=1;i<=k;i++)printf("%d ",ans[i]-1);//因为输入时为了方便全部加一，现在就要减一了 
}
int main()
{
    freopen("tarjan.in","r",stdin);
    freopen("tarjan.out","w",stdout);
    scanf("%d %d",&n,&m);
    for(int i=1;i<=m;i++){
        scanf("%d %d",&a,&b);//读入 
        addedge(a+1,b+1);//从题目上看非负数 可以有0，为了方便全部加一，输出时再减一 
    }
    for(int i=1;i<=n;i++)//因为不知道输入哪个节点 所以全部试一试 
        if(!dfn[i])//如果已经发现过了就不必在入栈了 
            tarjan(i);
    out(n);
    return 0;
}