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);
}
}