双DFS求强连通分量

用双DFS实现强连通分量是3个方法中比较好理解的。注意事项：1.读入文件时，我用的是fstream，流的方式读入！2.如果是数据集中没有别的多余符号，可以直接读入。若是由逗号等符号分隔开，我用的是CString中的readstring来按行读入，需引入头文件.在按行截取，用到find()函数，TrimLeft()和TrimRight()函数！3.按行读入数据后建立正反两个邻接表！按头插法建立效率更高些！4.建立一个Graph类和Readfile类，其中Graph类中包含图的各种信息。Readfile类中有读文件的方法，使得调用这个方法后会产生一个Graph类对象，包含数据集的所有信息，即调用这个函数后我们把数据集转换成了一个Graph对象！5.其中一个较难点是输出结点的出栈顺序，我用了一个time变量的引用来记录每个结点的出栈时间！具体代码如下：这是读文件的代码，结果是产生一个图对象！void Readdata::readfile(Graph &G)//产生了正反两个邻接表和标志访问数组{ifstream infile("F:\\Email-EuAll.txt");//ifstream infile("F:\\ceshi.txt");infile >> G.nodenum >> G.edgenum;int i,j,k;G.visited = new bool[G.nodenum];//bool *visited=new bool[G.nodenum];G.dfn = new int[G.nodenum];for(i = 0;i < G.nodenum;i ++)//构造访问标志数组{G.visited[i] = false;}G.Adj = new node[G.nodenum];//建立正反两个邻接表的顶点向量组G.VerAdj = new node[G.nodenum];for(j = 0;j < G.nodenum;j ++){G.Adj[j].data = j;G.Adj[j].firstarc = NULL;G.VerAdj[j].data = j;G.VerAdj[j].firstarc = NULL;}for(k= 0;k < G.edgenum;k ++){int node1,node2;infile >> node1 >>node2;arcnode* adjvex = new arcnode;//构造正邻接表的边结点adjvex->data = node2;adjvex->nextarc = NULL; if(G.Adj[node1].firstarc == NULL){G.Adj[node1].firstarc = adjvex;}else{adjvex->nextarc = G.Adj[node1].firstarc;G.Adj[node1].firstarc = adjvex;}arcnode* veradjvex = new arcnode;//构造反邻接表的边结点veradjvex->data = node1;veradjvex->nextarc = NULL;if(G.VerAdj[node2].firstarc == NULL){G.VerAdj[node2].firstarc = veradjvex;}else{veradjvex->nextarc = G.VerAdj[node2].firstarc;G.VerAdj[node2].firstarc = veradjvex;}}}这是第一次深度优先遍历，用&time来记录每个顶点的出栈时间：void dfsfir(Graph &G,int i,int &time){G.visited[i] = true;arcnode* w; for(w = G.Adj[i].firstarc;w != NULL;w = w->nextarc){if(!G.visited[w->data]) dfsfir(G,w->data,time);}G.dfn[time++] = i;}这是第二次深度优先遍历，输出各个强连通分量按行隔开！void dfssec(Graph &G,int i){G.visited[i] = true;//cout<<i<<" ";arcnode* w;for(w = G.VerAdj[i].firstarc; w != NULL; w = w->nextarc){if(! G.visited[w->data]) dfssec(G,w->data);}}这是主程序，通过调用读文件函数生成一个图对象，在这个对象上分别进行两次DFS即可!int main(){double start = GetTickCount();int time = 0;int count = 0; //记录强连通分量的个数Readdata R;Graph G;R.readfile(G);//readfile(G);int i,j,k;for(i = 0;i < G.nodenum;i ++){if(!G.visited[i]){dfsfir(G,i,time);}}for(j = 0;j < G.nodenum;j ++){G.visited[j] = false;}for(k = G.nodenum-1;k >= 0;k --){if(!G.visited[G.dfn[k]]){dfssec(G,G.dfn[k]);//cout<<endl;count ++;}}double end = GetTickCount();cout<<"强连通分量的个数为："<<count<<endl;cout<<"运行时间为："<<end - start<<"毫秒！"<<endl;system("pause");return 0;}