Kosaraju算法求有向图的极大连通子图。算法首先对逆图做后序遍历，求出后序编码。随后从右到左扫描编码，对于每一个未标记所属分量编号的节点，以它为根做一次DFS，标记沿途所遇到的所有节点为同一个强连通分量。
   算法正确性的证明基于：任意两点s,t在原图中属于强联通分量的话，当且仅当在逆图中两者也属于同一强连通分量

#define  MAXV 100000
  static   int  postR[MAXV],post[MAXV];
  static   int  cnt0,cnt1;

typedef  struct  LinkNod  * Link;
 struct  LinkNod {
     int dest;
     Link next;
} ;

typedef  struct  GraphNod  * Graph;
 struct  GraphNod {
     Link v[MAXV];
     int vnum;
     int id[MAXV];  // 分量标号 
} ;

 int  IsScc(Graph g, int  s, int  t) {
    return g->id[s] == g->id[t];
}

void  DFS(Graph g, int  w) {
     g->id[w] = cnt1;
     for( Link pos = g->v[w]; pos!=NULL; pos=pos->next ) if( g->id[ pos->dest ] == -1 ) DFS(g,pos->dest);
     post[cnt0++] = w; // 后序编码 
}

Graph MakeGraphR(Graph g) {
     Link pos;
     Link tmpNod;
     Graph r = (Graph)malloc(sizeof(struct GraphNod));
     r->vnum = g->vnum;
     for(int i=0;i<r->vnum;i++) r->v[i] = NULL;
     for(int i=0;i<g->vnum;i++){
        pos = g->v[i];
while(pos){  
     tmpNod = (Link)malloc(sizeof(struct LinkNod));
     tmpNod->dest = i; tmpNod->next = r->v[pos->dest]; 
     r->v[ pos->dest ] = tmpNod;
     pos=pos->next; 
  }
     }
     return r;
}

int  FindScc(Graph g) {
   Graph r = MakeGraphR(g);
   memset(r->id,-1,sizeof(r->id));
   cnt0 = cnt1 =0;
   for( int i=0;i<r->vnum;i++ ) if( r->id[i] == -1 ) DFS(r,i);
   for(int i=0;i<g->vnum;i++) postR[i] = post[i];
   memset(g->id,-1,sizeof(g->id)); 
   cnt0 = cnt1 =0;
   for( int i=g->vnum-1;i>-1;i-- ) if( g->id[ postR[i] ] == -1 ){  DFS(g,postR[i]); cnt1++;  }
   printf("分量阵:\n");
   for( int i=0;i<g->vnum;i++ ) printf("%d: %d\n",i,g->id[i]);
   return cnt1;  // 分量数 
}