前面博文对拓扑排序进行了引入
http://blog.csdn.net/derkampf/article/details/57414185,这一篇进行拓扑排序算法的实现。
一.概念
由偏序得到全序的操作称为拓扑排序。在拓扑排序之后,有向图中的每一个顶点都会按照一定的次序进行输出。
二.如何进行拓扑排序
步骤:
1.在有向图中选择一个无前驱顶点(无入度)且输出
2.删除该顶点以及所有以该顶点作为起始顶点的边。
重复1,2操作,最终直至全部顶点全部输出。若最终已经没有无入度的顶点,且有向图中还有顶点,说明AOV-网中有环(这也是拓扑排序的关键)
三.拓扑排序算法的实现
以下图为例,进行实现
1.对之前编写的邻接表代码进行改动,由于之前针对的是无向图,现在是有向图,所以需要更改的有:(1)插入边,(2)删除边,(3)删除顶点
(1)单向插入边
template<class Type>
bool GraphLink<Type>::InsertEdge(const Type &vertex1, const Type &vertex2)
{
int v1 = GetPosVertex(vertex1);
int v2 = GetPosVertex(vertex2);
if(v1 == -1 || v2 == -1)
return false;
//采用单链表的头插方式
//v1-->v2
Edge<Type> *e = new Edge<Type>(v2);
e->link = NodeTable[v1].adj;
NodeTable[v1].adj = e;
numEdges++;
}
(2)单向删除边
template<class Type>
bool GraphLink<Type>::RemoveEdge(const Type &vertex1, const Type &vertex2)
{
int v1 = GetPosVertex(vertex1);
int v2 = GetPosVertex(vertex2);
if(v1 == -1 || v2 == -1)
return false;
//删除v1-->v2
Edge<Type> *p = NodeTable[v1].adj;
Edge<Type> *q = NULL;
while(p != NULL && p->dest != v2)
{
q = p;
p = p->link;
}
if(NULL == p)
return false;
if(q == NULL)//说明删除的是头结点
NodeTable[v1].adj = p->link;
else
q->link = p->link;
delete p;
p == NULL;
numEdges--;
}
(3)删除顶点vertex时,遍历NodeTable,删除所有包含v的边,之后最后一个顶点替换vertex时,将所有边链表上包含的最后一个顶点下标替换成vertex的下标
template<class Type>
bool GraphLink<Type>::RemoveVertex(const Type &vertex)
{
/*思路:
** 1.剔除顶点vertex边链表中所包含顶点中的关于vertex的结点(可以借用RemoveEdage())
** 2.删除所有包含该顶点的边。
** 3.最后一行覆盖所删除行
(1)顶点覆盖顶点
1)p保存最后一行的单链表,删除行的顶点指向最后一行的边链表,tmp保存numVertices-1
2)覆盖顶点
3)与最后一行所关联的结点,将他们所包含的最后一个结点的下标更改成所删除的下标
(2)边链表覆盖边链表(指向tmp)
** 4.减少顶点数
*/
int v = GetPosVertex(vertex);
if(v == -1)
return false;
//1.
Edge<Type> *p = NodeTable[v].adj;
while(p != NULL)
{
RemoveEdge(vertex, NodeTable[p->dest].data);
p = NodeTable[v].adj;
}
//2.
int num = numVertices-1;
while(num >= 0)
{
p = NodeTable[num].adj;
if(NULL != p && p->dest != v)
p = p->link;
if(NULL != p && p->dest == v)
RemoveEdge(NodeTable[num].data, NodeTable[v].data);
num--;
}
//3.
//1)
p = NodeTable[numVertices-1].adj;
NodeTable[v].adj = p;
int tmp = numVertices-1;
//2)
NodeTable[v].data = NodeTable[numVertices-1].data;
//3)
num = numVertices - 2;
while(num >= 0)
{
p = NodeTable[num].adj;
if(NULL != p && p->dest != tmp)
p = p->link;
if(NULL != p && p->dest == tmp)
p->dest = v;
num--;
}
//4.
numVertices--;
return true;
}
2.进行拓扑排序
为了找到入度为0的顶点,我们需要借助一个数组模拟栈结构,步骤如下:
(1)建立一个count数组用来保存每个顶点的入度 ,top==-1,先全部初始化为0,然后根据邻接表将入度填充在数组中;
(2)当遇到入度为0的顶点(下标i)时,count[i] = top;top = i;
用来模拟入栈;
(3)将所有入度为0的点入栈完毕后,在这里如果一开始有向图就是一个闭环,则top依然为-1,这也是输出“有向图有环”的依据。
(4)进行打印”入度为0的顶点”的操作。即模拟出栈的操作v=top;top=count[top]
,然后打印下标为v的顶点,此时top回退到上一个入度为0的顶点(即模拟的栈顶)。
(5)进行模拟出栈操作后,将以出栈元素v为起始点的顶点元素(即第一个邻接顶点)的入度进行-1操作,如果-1之后入度为0,则该元素为新的栈顶元素,寻找下一个邻接顶点重复上述动作直至到达边链表的尾部
(6)如果有向图中有环,那么在进行(4)(5)时,top的值从未改变,那么会一直回退到top==-1的状态,即在循环时,到达(3)的操作,会输出”有环”,退出程序。
下面是具体实现
bool GraphLink<Type>::TopologicalSort()
{
//(1)
int top = -1;
Edge<Type> *p;
for(int i=0; i<numVertices; ++i)
{
p = NodeTable[i].adj;
while(NULL != p)
{
count[p->dest]++;
p = p->link;
}
}
//(2)
for(int i=0; i<numVertices; ++i)
{
if(count[i] == 0)
{
count[i] = top;
top = i;
}
}
int v,w;
for(int i=0; i<numVertices; ++i)
{
//(3)
if(top == -1)
{
cout<<"无向图有环"<<endl;
return false;
}
//(4)
v = top;
top = count[top];
cout<<NodeTable[v].data<<"-->";
w = GetFirstNeighbor(NodeTable[v].data);
while(w != -1)
{
//(5)
if(--count[w] == 0)
{
count[w] = top;
top = w;
}
w = GetNextNeighbor(NodeTable[v].data, NodeTable[w].data);
}
}
cout<<"over."<<endl;
return true;
}
源码展示:
/*图之拓扑排序算法实现
**GraphLink.h
**2016.2.20
*/
#pragma once
#include<iostream>
using namespace std;
#define DEFAULT_VERTEX_SIZE 10
template<class Type> class GraphLink;
template<class Type>
class Edge
{
friend class GraphLink<Type>;
public:
Edge(int num = -1) : dest(num),link(NULL)
{}
~Edge()
{}
private:
int dest;
Edge *link;
};
template<class Type>
class Vertex
{
friend class GraphLink<Type>;
public:
Vertex():data(),adj(NULL)
{}
~Vertex()
{}
private:
Type data;
Edge<Type> *adj;
};
template<class Type>
class GraphLink
{
public:
GraphLink(int sz = DEFAULT_VERTEX_SIZE)
{
maxVertices = sz > DEFAULT_VERTEX_SIZE ? sz : DEFAULT_VERTEX_SIZE;
numEdges = numVertices = 0;
NodeTable = new Vertex<Type>[maxVertices];
count = new int[maxVertices];
for(int i=0; i<maxVertices; ++i)
count[i] = 0;
}
~GraphLink()
{
int n = numVertices;
for(int i = 0; i < n; ++i)
RemoveVertex(NodeTable[i].data);
delete []NodeTable;
delete []count;
}
public:
bool InsertVertex(const Type &v); //插入顶点v
bool InsertEdge(const Type &vertex1, const Type &vertex2);//插入vertex1-->vertex2边
int NumberOfVertice()const; //获取顶点总数
int NumberOfEdge()const; //获取边总数
int GetFirstNeighbor(const Type &vertex)const; //获取vertex的第一个邻接顶点
int GetNextNeighbor(const Type &vertex1, const Type &vertex2)const;//获取vertex1的邻接顶点vertex2的下一个邻接顶点
bool RemoveVertex(const Type &vertex); //删除顶点vertex
bool RemoveEdge(const Type &vertex1, const Type &vertex2);//删除vertex1和vertex2构成的边
bool TopologicalSort();
public:
void ShowGraph()const
{
for (int i = 0; i < numVertices; ++i)
{
cout<<NodeTable[i].data<<"-->";
Edge<Type> *e = NodeTable[i].adj;
while(e != NULL)
{
cout<<e->dest<<"-->";
e = e->link;
}
cout<<"Nul"<<endl;
}
}
private:
int GetPosVertex(const Type v)const
{
for (int i = 0; i < numVertices; ++i)
{
if(NodeTable[i].data == v)
return i;
}
return -1;
}
Vertex<Type> *NodeTable;//顶点表
int maxVertices;
int numVertices;
int numEdges;
int *count;
};
template<class Type>
bool GraphLink<Type>::InsertVertex(const Type &v)
{
if(numVertices > maxVertices)
return false;
NodeTable[numVertices++].data = v;
return true;
}
template<class Type>
bool GraphLink<Type>::InsertEdge(const Type &vertex1, const Type &vertex2)
{
int v1 = GetPosVertex(vertex1);
int v2 = GetPosVertex(vertex2);
if(v1 == -1 || v2 == -1)
return false;
//采用单链表的头插方式
//v1-->v2
Edge<Type> *e = new Edge<Type>(v2);
e->link = NodeTable[v1].adj;
NodeTable[v1].adj = e;
numEdges++;
}
template<class Type>
int GraphLink<Type>::NumberOfVertice()const
{return numVertices;}
template<class Type>
int GraphLink<Type>::NumberOfEdge()const
{return numEdges;}
template<class Type>
int GraphLink<Type>::GetFirstNeighbor(const Type &vertex)const
{
int v = GetPosVertex(vertex);
if(v == -1)
return -1;
if(NodeTable[v].adj != NULL)
return NodeTable[v].adj->dest;
return -1;
}
template<class Type>
int GraphLink<Type>::GetNextNeighbor(const Type &vertex1, const Type &vertex2)const
{
int v1 = GetPosVertex(vertex1);
int v2 = GetPosVertex(vertex2);
if(v1 == -1 || v2 == -1)
return -1;
Edge<Type> *p = NodeTable[v1].adj;
while(p != NULL && p->dest != v2)
p = p->link;
if(NULL == p)
return -1;
if(p->link != NULL)
return p->link->dest;
return -1;
}
template<class Type>
bool GraphLink<Type>::RemoveEdge(const Type &vertex1, const Type &vertex2)
{
int v1 = GetPosVertex(vertex1);
int v2 = GetPosVertex(vertex2);
if(v1 == -1 || v2 == -1)
return false;
//删除v1-->v2
Edge<Type> *p = NodeTable[v1].adj;
Edge<Type> *q = NULL;
while(p != NULL && p->dest != v2)
{
q = p;
p = p->link;
}
if(NULL == p)
return false;
if(q == NULL)//说明删除的是头结点
NodeTable[v1].adj = p->link;
else
q->link = p->link;
delete p;
p == NULL;
numEdges--;
}
template<class Type>
bool GraphLink<Type>::RemoveVertex(const Type &vertex)
{
/*思路:
** 1.剔除顶点vertex边链表中所包含顶点中的关于vertex的结点(可以借用RemoveEdage())
** 2.删除所有包含该顶点的边。
** 3.最后一行覆盖所删除行
(1)顶点覆盖顶点
1)p保存最后一行的单链表,删除行的顶点指向最后一行的边链表,tmp保存numVertices-1
2)覆盖顶点
3)与最后一行所关联的结点,将他们所包含的最后一个结点的下标更改成所删除的下标
(2)边链表覆盖边链表(指向tmp)
** 4.减少顶点数
*/
int v = GetPosVertex(vertex);
if(v == -1)
return false;
//1.
Edge<Type> *p = NodeTable[v].adj;
while(p != NULL)
{
RemoveEdge(vertex, NodeTable[p->dest].data);
p = NodeTable[v].adj;
}
//2.
int num = numVertices-1;
while(num >= 0)
{
p = NodeTable[num].adj;
if(NULL != p && p->dest != v)
p = p->link;
if(NULL != p && p->dest == v)
RemoveEdge(NodeTable[num].data, NodeTable[v].data);
num--;
}
//3.
//1)
p = NodeTable[numVertices-1].adj;
NodeTable[v].adj = p;
int tmp = numVertices-1;
//2)
NodeTable[v].data = NodeTable[numVertices-1].data;
//3)
num = numVertices - 2;
while(num >= 0)
{
p = NodeTable[num].adj;
if(NULL != p && p->dest != tmp)
p = p->link;
if(NULL != p && p->dest == tmp)
p->dest = v;
num--;
}
//4.
numVertices--;
return true;
}
template<class Type>
bool GraphLink<Type>::TopologicalSort()
{
//(1)建立一个count数组用来保存每个顶点的入度 ,top==-1,先全部初始化为0,然后根据邻接表将入度填充在数组中;
//(2)当遇到入度为0的顶点(下标i)时,`count[i] = top;top = i;` 用来模拟入栈;
//(3)将所有入度为0的点入栈完毕后,在这里如果一开始有向图就是一个闭环,则top依然为-1,这也是输出“有向图有环”的依据。
//(4)进行打印"入度为0的顶点"的操作。即模拟出栈的操作`v=top;top=count[top]`,然后打印下标为v的顶点,此时top回退到上一个入度为0的顶点(即模拟的栈顶)。
//(5)进行模拟出栈操作后,将以出栈元素v为起始点的顶点元素(即第一个邻接顶点)的入度进行-1操作,如果-1之后入度为0,则该元素为新的栈顶元素,寻找下一个邻接顶点重复上述动作直至到达边链表的尾部
//(6)如果有向图中有环,那么在进行(4)(5)时,top的值从未改变,那么会一直回退到top==-1的状态,即在循环时,到达(3)的操作,会输出"有环",退出程序。
//(1)
int top = -1;
Edge<Type> *p;
for(int i=0; i<numVertices; ++i)
{
p = NodeTable[i].adj;
while(NULL != p)
{
count[p->dest]++;
p = p->link;
}
}
//(2)
for(int i=0; i<numVertices; ++i)
{
if(count[i] == 0)
{
count[i] = top;
top = i;
}
}
int v,w;
for(int i=0; i<numVertices; ++i)
{
//(3)
if(top == -1)
{
cout<<"无向图有环"<<endl;
return false;
}
//(4)
v = top;
top = count[top];
cout<<NodeTable[v].data<<"-->";
w = GetFirstNeighbor(NodeTable[v].data);
while(w != -1)
{
//(5)
if(--count[w] == 0)
{
count[w] = top;
top = w;
}
w = GetNextNeighbor(NodeTable[v].data, NodeTable[w].data);
}
}
cout<<"over."<<endl;
return true;
}
测试
//test.cpp
/*Test.cpp
**2016.2.16
*/
#include"GraphLink.h"
#include<vld.h>
int main()
{
GraphLink<char> gl;
gl.ShowGraph();
gl.InsertVertex('A');
gl.InsertVertex('B');
gl.InsertVertex('C');
gl.InsertVertex('D');
gl.InsertVertex('E');
/*gl.InsertEdge('A', 'B');
gl.InsertEdge('A', 'C');
gl.InsertEdge('B', 'E');
gl.InsertEdge('C', 'D');
gl.InsertEdge('D', 'B');
gl.InsertEdge('E', 'D');*/
gl.InsertVertex('F');
gl.InsertEdge('A', 'B');
gl.InsertEdge('A', 'C');
gl.InsertEdge('A', 'D');
gl.InsertEdge('C', 'E');
gl.InsertEdge('C', 'B');
gl.InsertEdge('D', 'E');
gl.InsertEdge('E', 'G');
gl.InsertEdge('F', 'D');
gl.InsertEdge('F', 'E');
gl.ShowGraph();
gl.TopologicalSort();
}
测试结果
可以看到最后拓扑排序的结果是
F–>A–>C–>B–>D–>E–>over
在回忆一下原来的有向无环图: