分别使用Kahn和DFS实现拓扑排序

1，先了解什么是偏序？

偏序就是图中存在无先后顺序的顶点对。全序即找不到这样的顶点对。

2，什么是拓扑排序？

把图由偏序变成全序的操作即是拓扑排序。前提是有向无环图，拓扑才能成功。

3，Kahn算法思想：先把入度为0的点入栈，每次从中出栈一个顶点，删除从它出发的边，相应顶点的入度减少1，若有入度为0的点，则入栈，重复至栈为空。

4，DFS思想：前提是有向无环图，因为最先完成DFS的顶点恰好拓扑排序中最后一个完成的，最后完成DFS正好是拓扑排序中入度为0的顶点。

下边是代码：

#include "stdafx.h" #include <fstream> #include <iostream> #include <stack> using namespace std; stack <int> s; typedef int DataType; #define Max_Vertex_Num 20 #define path "F:\\test2.txt" typedef struct ADG{ int vexnum,arcnum; int *indegree; int *visited; struct vnode *Adjlist;//邻接表 int *dfs;//记录dfs完成的顺序 }Graph;//有向无环图的结构体定义 struct vnode{ DataType data; struct arcnode* firstarc; vnode(DataType x) { data = x; firstarc = NULL; } vnode() { firstarc = NULL; } }; struct arcnode{ DataType data; struct arcnode* nextarc; arcnode(DataType x) { data = x; nextarc = NULL; } arcnode() { nextarc = NULL; } }; void Readfile(Graph &g)//读文件，建立邻接表，访问数组，入度数组 { ifstream infile(path); cout<<"读入第一行数字，即顶点树和边数！"<<endl infile="">> g.vexnum >> g.arcnum; cout<<g.vexnum<<' '<<g.arcnum<<endl; g.indegree = new int[g.vexnum]; g.visited = new int[g.vexnum]; g.Adjlist = new vnode[g.vexnum]; g.dfs = new int[g.vexnum]; memset(g.dfs,0,sizeof(int)*(g.vexnum)); memset(g.indegree,0,sizeof(int)*(g.vexnum)); memset(g.visited,0,sizeof(int)*(g.vexnum)); cout<<"建立邻接表！"<<endl; int i = 0,j = 0; for(i;i < g.vexnum;i ++) { g.Adjlist[i].data = i; } int head,tail; for(j;j < g.arcnum;j ++) { infile >> head >> tail;//读入每一条边，弧头为head，弧尾为tail arcnode *newarc = new arcnode(tail);//采用头插法，避免扫描了 newarc->nextarc = g.Adjlist[head].firstarc; g.Adjlist[head].firstarc = newarc; g.indegree[tail]++; } } bool KahnTopusort(Graph g)//用栈实现 { int i =0; int count = 0;//记录出栈的次数 int key; arcnode *p; DataType key1; for(i;i < g.vexnum;i ++) { if (!g.indegree[i]) { s.push(i);//入度为0的顶点入栈 } } while (!s.empty()) { key = s.top();//取栈顶元素 s.pop();//出栈 cout<<key<<"-->"; count++; for (p = g.Adjlist[key].firstarc;p;p = p->nextarc) { key1 = p->data; g.indegree[key1]--; if (!g.indegree[key1]) { s.push(key1); } } } if(count < g.vexnum)return false; else return true; } void DFS(Graph &g,int v,int &time)//用&time来记录出栈顺序 { g.visited[v] = 1; arcnode *p; for(p = g.Adjlist[v].firstarc;p;p = p->nextarc) { if (!g.visited[p->data]) { DFS(g,p->data,time); } } g.dfs[time++] = v; } void DFSTopuSort(Graph &g,int &time) { int j = 0; for(j;j < g.vexnum;j ++) { if(!g.visited[j]) { DFS(g,j,time); } } int i = g.vexnum - 1; for(i;i >= 0;i --) { cout<<g.dfs[i]<<"-->"; } } int _tmain(int argc, _TCHAR* argv[]) { Graph A; Readfile(A); //if (KahnTopusort(A)) //{ // cout<<"完成拓扑排序！"<<endl; //} //else // cout<<"图有环！"<<endl; int time = 0; DFSTopuSort(A,time); system("pause"); return 0; } </endl></int></stack></iostream></fstream>

test2.txt中内容为：

12 16
0 1
0 2
1 2
0 3
2 4
3 4
10 5
4 6
2 6
2 7
5 7
8 9
8 10
8 11
9 11
0 11

用Kahn算法求出来的结果为：

用DFS算法求出来的结果为：

两个结果是一样的。

Kahn的时间复杂度为O(n+e);

DFS的时间复杂度为O(n+e),空间复杂度都为O(n);

应用背景：在有关工程进度/次序的规划问题中，拓扑排序有着大量应用。

少了个释放空间的操作：

void deletegraph(Graph &g)//怎么delete邻接表 { delete []g.indegree; g.indegree = NULL; delete []g.visited; g.visited = NULL; delete []g.dfs; g.dfs = NULL; int i = 0; arcnode *s,*p; for (i;i < g.vexnum;i++) { p = g.Adjlist[i].firstarc; while(p) { s = p->nextarc; delete p; p = s; } } delete []g.Adjlist; g.Adjlist = NULL; }

若数据集为下边这种格式：第一行为顶点数，边数；其他行的意思分别为：弧尾，出度，弧头结点列

10 12
0 2 1 2
1 2 3 4
2 2 3 5
3 1 7
4 1 8
5 1 6
6 1 8
7 2 8 9
8 0
9 0

则建立邻接表的代码为：void Readfile(Graph &g)//读文件，建立邻接表，访问数组，入度数组 { ifstream infile(path); cout<<"读入第一行数字，即顶点树和边数！"<<endl infile="">> g.vexnum >> g.arcnum; cout<<g.vexnum<<' '<<g.arcnum<<endl; g.indegree = new int[g.vexnum]; g.visited = new int[g.vexnum]; g.Adjlist = new vnode[g.vexnum]; g.dfs = new int[g.vexnum]; memset(g.dfs,0,sizeof(int)*(g.vexnum)); memset(g.indegree,0,sizeof(int)*(g.vexnum)); memset(g.visited,0,sizeof(int)*(g.vexnum)); cout<<"建立邻接表！"<<endl; int i = 0,j = 0; for(i;i < g.vexnum;i ++) { g.Adjlist[i].data = i; } int head,count; for (i = 0; i < g.vexnum; i++) { infile >> head >> count; //cout << head << ' '<< count << endl; int *tail = new int[count]; for (int t = 0; t < count; t++) { infile >> tail[t]; //cout << tail[t] << endl; arcnode *newarc = new arcnode(tail[t]); newarc->nextarc = g.Adjlist[head].firstarc; g.Adjlist[head].firstarc = newarc; g.indegree[tail[t]]++; //cout << g.indegree[tail[t]] <<"---------------"<<endl; } delete []tail; } } </endl>