图之有向无环图的应用--拓扑排序算法的实现

前面博文对拓扑排序进行了引入
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
在回忆一下原来的有向无环图:
这里写图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值