kosaraju+缩点

1、对图G进行深度优先搜索，按递归调用完成的先后顺序对各顶点编号。
2、将原图的每条边反向——直接转置邻接矩阵即可。得到新的一个图G‘。
3、从编号最大的顶点开始对G’进行深度优先搜索。如果搜索的过程中没有访问遍G‘的所有顶点，则从未被访问过的顶点中再选取编号最大的顶点，并从此顶点开始继续做深度优先搜索。（即这里的深度优先搜索是按顶点编号降序进行的）
4、在最后得到的G’的深度优先森林中，每棵树上的顶点组成G的一个强连通分支。

#include <iostream>
#include <fstream>
#include <memory.h>
using namespace std;

#define init(l) memset(l,0,sizeof(l))
#define N 101
struct edge
{
 int dest;
 edge *next;
};
edge *G1[N],*G2[N],*G[N];
int used[N],path[N],part[N],mark[N][N];
int m,n;
void readdata();
void dfs1(int s);
void dfs2(int s);
void kosaraju();
void addedge(edge *G[],int i, int j);
void shrink();

int main()
{
 readdata();
 kosaraju();
 shrink();
 work();
 return 0;
}

void addedge(edge *T[],int i, int j)
{
 edge *L;
 L=new edge;
 L->dest=j;
 L->next=T[i];
 T[i]=L;
}

void shrink()
{
 int i,j,k;
 edge *L;
 init(mark);
 for (k=1;k<=n;k++)
  for (L=G1[k],i=part[k];L!=NULL;L=L->next)
  {
  j=part[L->dest];
  if (i!=j&&!mark[i][j])
  {
  mark[i][j]=1;
  addedge(G,i,j);
  }
  }
}

void dfs1(int s)
{
 if(!used[s])
 {
  used[s]=1;
  for (edge *l=G1[s];l!=NULL;l=l->next)
  dfs1(l->dest);
  path[0]++;
  path[path[0]]=s;
 }
}

void dfs2(int s)
{
 if (!used[s])
 {
  used[s]=1;
  for (edge *l=G2[s];l!=NULL;l=l->next)
  dfs2(l->dest);
  part[s]=part[0];
 }
}

void kosaraju()
{
 int i;
 init(used);
 path[0]=part[0]=0;
 for (i=1;i<=n;i++)
  dfs1(i);
 init(used);
 for (i=n;i>=1;i--)
  if (!used[path[i]])
  {
  part[0]++;
  dfs2(path[i]);
  }
}

void readdata()
{
// ifstream cin("data.txt");
 int i,j;
 cin>>n;
 for (i=1;i<=n;i++)
  G1[i]=G2[i]=G[i]=NULL;
 for (i=1;i<=n;i++)
  while (cin>>j&&j)
  {
  addedge(G1,i,j);
  addedge(G2,j,i);
  }
}