第8章_图论

图是又一种非线性数据结构。在图结构中，数据元素之间的关系是多对多的，即如果任选一个结点作为初始结点，则图中任意一个结点有多个前驱结点和多个后继结点。

图

图是由结点集合及结点间的关系集合组成的一种数据结构。

图的实现

领接矩阵存储结构下图操作的实现

例：带权图为

SeqList.h

typedef struct {    DataType list[MaxSize];    int size;}SeqList;//顺序表操作实现//1.初始化void ListInitiate(SeqList*L){    L->size = 0;//定义初始数据元素个数}//2.求当前数据元素个数int ListLength(SeqList L){    return L.size;//返回顺序表L的当前数据元素个数}//3.插入数据元素int ListInsert(SeqList*L,int i,DataType x){    /*    在顺序表L的第i(0<=i<=size)个位置前插入数据元素值x    插入成功返回1，插入失败返回0    */    int j;    if(L->size >= MaxSize)    {        printf("顺序表已满无法插入！\n");        return 0;    }    else if(i < 0 || i > L->size)    {        printf("参数i不合法！\n");        return 0;    }    else    {        /*为插入做准备*/        for(j=L->size;j > i; j--)        {            L->list[j] = L->list[j-1];        }        L->list[i] = x; /*插入元素x*/        L->size++;      /*元素个数加1*/        return 1;    }}//4.删除数据元素int ListDelete(SeqList*L,int i,DataType*x){    /*    删除顺序表L中位置为i(0<=i<=size-1)的数据元素并存放到x中    删除成功返回1，删除失败返回0  */    int j;    if(L->size <= 0)    {        printf("顺序表已空无数据元素可删！\n");        return 0;    }    else if(i < 0 || i > L->size - 1)    {        printf("参数i不合法\n");        return 0;    }    else    {        *x = L->list[i];    /* 保存删除的元素到x中*/        for(j = i;j < L->size-1; j++)        {            L->list[j] = L->list[j+1];        }        L->size--;          /*数据元素个数减1*/        return 1;    }}//5.取数据元素int ListGet(SeqList L,int i,DataType* x){    /* 取顺序表L中第i个数据元素存于x中，成功返回1，失败返回0*/    if(i < 0 || i > L.size - 1)    {        printf("参数i不合法！\n");        return 0;    }    else    {        *x=L.list[i];        return 1;    }}

AdjMGraph.h

#include "SeqList.h"        //包含顺序表头文件typedef struct {    SeqList Vertices;                       //存放结点的顺序表    int edge[MaxVertices][MaxVertices];     //存放边的邻接矩阵    int numOfEdges;                         //边的条数} AdjMGraph;                                //图的结构体定义void Initiate(AdjMGraph* G,int n)   //初始化{    int i,j;    for(i = 0;i < n; i++)    {        for(j = 0;j < n; j++)        {            if(i == j) G->edge[i][j] = 0;            else                G->edge[i][j] = MaxWeight;        }    }    G->numOfEdges = 0;              //连的条数置为0    ListInitiate(&G->Vertices);     //顺序表初始化}//在图G中插入结vertexvoid InsertVertex(AdjMGraph*G,DataType vertex){    ListInsert(&G->Vertices,G->Vertices.size,vertex);       //在顺序表尾插入}//在图G中插入边<v1,v2>,边<v1,v2>的权为weightvoid InsertEdge(AdjMGraph*G,int v1,int v2,int weight){    if(v1 < 0 || v1 > G->Vertices.size || v2 < 0 || v2 > G->Vertices.size)    {        printf("参数v1或v2越界出错！\n");        exit(1);    }    G->edge[v1][v2] = weight;    G->numOfEdges ++;}//在图G中删除边<v1,v2>void DeleteEdge(AdjMGraph*G,int v1,int v2){    if(v1 < 0 || v1 > G->Vertices.size || v2 < 0 || v2 > G->Vertices.size || v1 == v2)    {        printf("参数v1或v2越界出错！\n");        exit(1);    }    if(G->edge[v1][v2] == MaxWeight || v1 == v2)    {        printf("该边不存在！\n");        exit(0);    }    G->edge[v1][v2] = MaxWeight;    G->numOfEdges--;}//删除结点Vvoid DeleteVerten(AdjMGraph*G,int V){    int n = ListLength(G->Vertices),i,j;    DataType x;    for(i = 0;i < n; i++)    {        for(j = 0;j < n; j++)        {            if((i == V || j == V) && G->edge[i][j] > 0 && G->edge[i][j] < MaxWeight)                G->numOfEdges--;            //计算被删除边        }    }    for(i = V;i < n; i++)        for(j = 0;j < n; j++)            G->edge[i][j] = G->edge[i + 1][j];    for(i = 0;i < n; i++)        for(j = V; j < n; j++)                  //删除结第V列            G->edge[i][j] = G->edge[i][j + 1];    ListDelete(&G->Vertices,V,&x);              //删除结V}//在图G中寻找序号为v的结点的第一个邻接结点//如果这样的邻 接结存在，返回该邻接结的序号；否则，返回-1int GetFirstVex(AdjMGraph G,int v){    int col;    if(v < 0 || v > G.Vertices.size)    {        printf("参数v1越界出错！\n");        exit(1);    }    for(col = 0;col < G.Vertices.size; col++)        if(G.edge[v][col] > 0 && G.edge[v][col] < MaxWeight) return col;    return -1;}//在图G中寻找v1结点的邻接结点v2的下一个邻接结点//如果这样的邻接结点存在，返回该邻接结点的序号；否则，返回-1//v1和v2都是相应结点的序号int GetNextVex(AdjMGraph G,int v1,int v2){    int col;    if(v1 < 0 || v1 > G.Vertices.size || v2 < 0 || v2 > G.Vertices.size)    {        printf("参数v1或v2越界出错！\n");        exit(1);    }    for(col = v2 + 1;col < G.Vertices.size; col++)        if(G.edge[v1][col] > 0 && G.edge[v1][col] < MaxWeight) return col;    return -1;}

AdjMGraphCreate.h

typedef struct {    int row;        //行下标    int col;        //列下标    int weight;     //权值} RowColWeight;     //边信息结构体定义//在图G中插入n个结点信息V和e条边信息Evoid CreatGraph(AdjMGraph*G,DataType V[],int n,RowColWeight E[],int e){    int i,k;    Initiate(G,n);    for(i = 0;i < n; i++)        InsertVertex(G,V[i]);                           //结点插入    for(k = 0;k < e; k++)        InsertEdge(G,E[k].row,E[k].col,E[k].weight);    //边插入}

main.c

#include <stdio.h>#include <stdlib.h>typedef char DataType;#define MaxSize 100#define MaxVertices 10#define MaxEdges 100#define MaxWeight 10000#include "AdjMGraph.h"#include "AdjMGraphCreate.h"void main(void){    AdjMGraph g1;    DataType a[] = {'A','B','C','D','E'};    RowColWeight rcw[] = {        {0,1,10},        {0,4,20},        {0,3,30},        {2,1,40},        {3,2,50},    };    int n = 5,e = 5;    int i,j;    CreatGraph(&g1,a,n,rcw,e);    DeleteEdge(&g1,0,4);        //删除<0,4>    DeleteVerten(&g1,3);        //删除结点3    printf("结点集合为：");    for(i = 0;i < g1.Vertices.size; i++)        printf("%c ",g1.Vertices.list[i]);    printf("\n");    printf("权值集合为：\n");    for(i = 0;i < g1.Vertices.size; i++)    {        for(j = 0;j < g1.Vertices.size; j++)            printf("%5d ",g1.edge[i][j]);        printf("\n");    }}

邻接表的存储结构

typedef struct Node{    int dest;               //邻接边的弧头结点序号    struct Node * next;     } Edge;                     //邻接边单链表的结点结构体typedef struct {    DataType data;          //结点数据元素    int score;              //邻接边的弧尾结点序号    Edge * adj;             //邻接边的头指针} AdjLHeight;               //数组的数据元素类型结构体typedef struct {    AdjLHeight a[MaxVertices];  //邻接表数组    int numOfVerts;             //结点个数    int numOfEdges;             //边个数} AdjLGraph;                    //邻接表结构体//初始化操作函数void AdjInitiate(AdjLGraph* G){    int i;    G->numOfVerts = 0;    G->numOfEdges = 0;    for(i = 0;i < MaxVertices; i++)    {        G->a[i].score = i;          //置邻接边的弧头结点        G->a[i].adj = NULL;         //置邻接边单链表头指针初值    }}//撤消操作函数void AdjDestroy(AdjLGraph*G){    int i;    Edge*p,*q;    for(i = 0;i < G->numOfVerts; i++)    {        p = G->a[i].adj;        while(p != NULL)        {            q = p->next;            free(p);            p = q;        }    }}//插入结点操作函数//在图G中的第i(0<=i<MaxVertices)个位置插入结点数据元素vertexvoid InsertVertex(AdjLGraph*G,int i,DataType vertex){    if(i >= 0 && i < MaxVertices)    {        G->a[i].data = vertex;          //存储结点数据元素vertex        G->numOfVerts++;                //个数加1    }    else        printf("结点越界");}//插入边操作函数//在图G中加入边<v1,v2>的信息void InsertEdge(AdjLGraph*G,int v1,int v2){    Edge*p;    if(v1 < 0 || v1 >= G->numOfVerts || v2 < 0 || v2 >= G->numOfVerts)    {        printf("参数v1或v2越界出错！");        exit(0);    }    p = (Edge*)malloc(sizeof(Edge));        //申请邻接边单链表结点空间    p->dest = v2;                           //置邻接边弧头序号    p->next = G->a[v1].adj;                 //新结点插入单链表的表头    G->numOfEdges++;                        //边个数加1}//删除边操作函数//删除图G中的边<v1,v2>信息void DeleteEdge(AdjLGraph*G,int v1,int v2){    Edge*curr,*pre;    if(v1 < 0 || v1 >= G->numOfVerts || v2 < 0 || v2 >= G->numOfVerts)    {        printf("参数v1或v2越界出错！");        exit(0);    }    pre = NULL;    curr = G->a[v1].adj;    while(curr != NULL && curr->dest != v2)    {        //在v1结点的邻接边单链表中查找v2结点        pre = curr;        curr = curr->next;    }    //删除表示邻接边<v1,v2>的结点    if(curr != NULL && curr->dest == v2 && pre == NULL)    {        //当邻 接边<v1,v2>的结点是单链表的第一个结点时        G->a[v1].adj = curr->next;        free(curr);        G->numOfEdges--;    }    else if(curr != NULL && curr->dest == v2 && pre != NULL)    {        //当邻 接边<v1,v2>的结点不是单链表的第一个结点时        pre->next = curr->next;        free(curr);        G->numOfEdges--;    }    else    {        //当邻接边<v1,v2>的结点不存在时        printf("边<v1,v2>不存在！");    }}//取第1个邻接结点int GetFirstVex(AdjLGraph G,int v){    //取图G中结点v的第一个邻接结点    //取到时返回该邻接结点的对应序号，否则返回-1    Edge*p;    if(v < 0 || v > G.numOfVerts)    {        printf("参数v1或v2越界出错！");        exit(0);    }    p = G.a[v].adj;    if(p != NULL) return p->dest;   //返回该邻接结点的对应序号    else        return -1;                  //返回-1}//取下一个邻接结点//取图G中结点v1的邻 接结点v2的下一个邻接结点//取到时返回该邻 接结瓣对应序号，否则返回-1int GetNextVex(AdjLGraph G,int v1,const int v2){    Edge * p;    if(v1 < 0 || v1 > G.numOfVerts || v2 < 0 || v2 > G.numOfVerts)    {        printf("参数v1或v2越界出错！");        exit(0);    }    p = G.a[v1].adj;    while(p != NULL)    {        if(p->dest != v2)        {            p = p->next;            continue;        }        else break;    }    p = p->next;                        //p指向邻接结点v2的下一个邻 接结点    if(p != NULL) return p->dest;       //返回该邻接结点的对应序号    else return -1;                     //返回-1}typedef struct {    int row;        //行下标    int col;        //列下标    int weight;     //权值} RowCol;       //边信息结构体定义//创建图//创建有n个结点e条边的图G//结点信息存放在数组v中，边信息存放在数组d中void CreateGraph(AdjLGraph*G,DataType v[],int n,RowCol d[],int e){    int i,k;    AdjInitiate(G);    for(i = 0;i < n; i++) InsertVertex(G,i,v[i]);           //插入结点    for(k = 0;k < e; k++) InsertEdge(G,d[k].row,d[k].col);  //插入边}

图的遍历

图的深度和广度优先遍历算法

与树的遍历操作相似，图遍历操作的定义是访问图中的每一个结点且每个结点只被访问一次。图的遍历方法主要有两种：一种是深度优先搜索遍历，另一种是广度优先搜索遍历。图的深度优先搜索遍历类同于树的先根遍历，图的广度优先搜索遍历类同于树的层序遍历。

深度优先遍历
AdjMGraphTraverse.h

#include "SeqCQueue.h"//连通图G以v为初始结点的访问操作为Visit()的深度优先遍历//数组visited标记了相应结点是否已访问过，0表示未访问，1表示已访问void DepthFSearch(AdjMGraph G,int v,int visited[],void Visit(DataType item)){    int w;    Visit(G.Vertices.list[v]);          //访问结点    visited[v] = 1;                     //置已访问标记    w = GetFirstVex(G,v);               //取第一个邻接结点    while(w != -1)    {        if(!visited[w])            DepthFSearch(G,w,visited,Visit);    //递归        w = GetNextVex(G,v,w);                  //取下一个邻接结点    }}//非连通图G的访问操作为Visit()的深度优先遍历void DepthFirstSearch(AdjMGraph G,void Visit(DataType item)){    int i;    int * visited = (int*)malloc(sizeof(int)*G.Vertices.size);    for(i = 0;i < G.Vertices.size; i++)        visited[i] = 0;    for(i = 0;i < G.Vertices.size; i++)        if(!visited[i])            DepthFSearch(G,i,visited,Visit);    free(visited);}//广度优先遍历//连通图G以v为初始结点访问操作为Visit()的广度优先遍历//数组visited标记了相应结点是否已访问过，0表示未访问，1表示已访问void BroadFSearch(AdjMGraph G,int v,int visited[],void Visit(DataType item)){    DataType u,w;    SeqCQueue queue;    Visit(G.Vertices.list[v]);      //访问结点    visited[v] = 1;    QueueInitiate(&queue);          //队列初始化    QueueAppend(&queue,v);          //初始结点v入入    while(QueueNotEmpty(queue))     //队列非空时    {        QueueDelete(&queue,&u);     //出队列        w = GetFirstVex(G,u);       //取结点u的第一个邻接结点        while(w != -1)        {            if(!visited[w])            {                Visit(G.Vertices.list[w]);      //访问结点w                visited[w] = 1;                 //置已访问标记                QueueAppend(&queue,w);          //结点w入队列            }            w = GetNextVex(G,u,w);              //取下一个邻接结点        }    }}//非连通图G访问操作为Visit()的广度优先遍历void BroadFirstSearch(AdjMGraph G,void Visit(DataType item)){    int i;    int * visited = (int*)malloc(sizeof(int)*G.Vertices.size);    for(i = 0;i < G.Vertices.size; i++)        visited[i] = 0;    for(i = 0;i < G.Vertices.size; i++)        if(!visited[i]) BroadFSearch(G,i,visited,Visit);    free(visited);}

main.c

#include <stdio.h>#include <stdlib.h>#include <malloc.h>typedef char DataType;#define MaxSize 100#define MaxVertices 10#define MaxEdges 100#define MaxWeight 10000#define MaxQueueSize 100#include "AdjMGraph.h"#include "AdjMGraphCreate.h"#include "AdjMGraphTraverse.h"void Visit(DataType item) //定义访问操作的函数{    printf("%c ",item);}void main(void){    AdjMGraph g1;    DataType a[] = {'A','B','C','D','E'};    RowColWeight rcw[] = {        {0,1,10},        {0,4,20},        {0,3,30},        {2,1,40},        {3,2,50},    };    int n = 5,e = 5;    CreatGraph(&g1,a,n,rcw,e);    printf("深度优先搜索序列为：");    DepthFirstSearch(g1,Visit);    printf("\n广度优先搜索序列：");    BroadFirstSearch(g1,Visit);    printf("\n");}

最小生成树

如果无向连通图是一个带权图，那么它的所有生成树中必有一棵边的权值总和为最小的生成树，称这棵生成树为最小代价生成树，简称最小生成树。

普里姆算法

Prim.h

typedef struct {    VerT vertex;    int weight;} MinSpanTree;//用普里姆算法建立带权图G的最小生成树closeVertexvoid Prim(AdjMGraph G,MinSpanTree closeVertex[]){    VerT x;    int n = G.Vertices.size,minCost;    int * lowCost = (int*)malloc(sizeof(int)*n);    int i,j,k;    for(i = 1;i < n; i++)               //初始化        lowCost[i] = G.edge[0][i];    //从结点0出发构造最小生成树    ListGet(G.Vertices,0,&x);           //取结点0    closeVertex[0].vertex = x;          //保存结点0    lowCost[0] = -1;                    //标记结点0    for(i = 1;i < n; i++)    {        //寻找当前最小权值的边所对应的弧头结点k        minCost = MaxWeight;        //MaxWeight为定义的最大权值        for(j = 1;j < n; j++)        {            if(lowCost[j] < minCost && lowCost[j] > 0)            {                minCost = lowCost[j];                k = j;            }        }        ListGet(G.Vertices,k,&x);       //取弧头结点k        closeVertex[i].vertex = x;      //保存弧头结点k        closeVertex[i].weight = minCost;    //保存相应的权值        lowCost[k] = -1;                //标记结点k        //根据加入集合U的结点k修改lowCost中的数值        for(j = 1;j < n; j++)        {            if(G.edge[k][j] < lowCost[j])                lowCost[j] = G.edge[k][j];        }    }}

main.c

#include <stdio.h>#include <stdlib.h>#include <malloc.h>typedef char DataType;typedef char VerT;#define MaxSize 100#define MaxVertices 10#define MaxWeight 10000#define N 7#include "AdjMGraph.h"#include "AdjMGraphCreate.h"#include "Prim.h"void main(void){    AdjMGraph g;    DataType a[] = {'A','B','C','D','E','F','G'};    RowColWeight rcw[] = {        {0,1,50},        {1,0,50},        {0,2,66},        {2,0,60},        {1,3,65},        {3,1,65},        {1,4,40},        {4,1,40},        {2,3,52},        {3,2,52},        {2,6,45},        {6,2,45},        {3,4,50},        {4,3,50},        {3,5,30},        {5,3,30},        {3,6,42},        {6,3,42},        {4,5,70},        {5,4,70},    };    int n = 7,e = 20,i;    MinSpanTree closeVertex[7];    CreatGraph(&g,a,n,rcw,e);    Prim(g,closeVertex);    //输出prim函数得到的最小生成树的结点序列和权值序列    printf("初始结点 = %c\n",closeVertex[0].vertex);    for(i = 1;i < n; i++)        printf("结点 = %c 边的权值 = %d\n",closeVertex[i].vertex,closeVertex[i].weight);}/*初始结点 = A结点 = B 边的权值 = 50结点 = E 边的权值 = 40结点 = D 边的权值 = 50结点 = F 边的权值 = 30结点 = G 边的权值 = 42结点 = C 边的权值 = 45Press any key to continue*/

最短路径

在一个图中，若从一个结点到另一个结点存在着路径，定义路径长度为一条路径上所经过的边的数目。图中从一个结点到另一个结点可能存在着多条路径，则路径长度最短的那条称作最短路径，其长度称作最短路径长度或最短距离。

狄克斯特拉算法

Dijkstra.h

/*函数共有四个参数，其中两个为输入参数，分别为带权图G和源点序号v0;另外两个为输出参数，分别为distance[]和path[]。distance[]用来存放得到的从源点v0到其余各结点的最短距离数值，path[]用来存放得到的从源点v0到其余各结点的最短路径上到达目标结点的前一结点下标  */void Dijkstra(AdjMGraph G,int v0,int distance[],int path[]){    int n = G.Vertices.size;    int * s = (int*)malloc(sizeof(int)*n);    int minDis,i,j,u;    //初始化    for(i = 0;i < n; i++)    {        distance[i] = G.edge[v0][i];        s[i] = 0;        if(i != v0 && distance[i] < MaxWeight) path[i] = v0;        else            path[i] = -1;    }    s[v0] = 1;  //标记结点v0已从集合T加入到集合S中    //在当前还未找到最短 路径的结点集中选取具有最短距离的结点u    for(i = 1;i < n; i++)    {        minDis = MaxWeight;        for(j = 0;j < n; j++)            if(s[j] == 0 && distance[j] < minDis)            {                u = j;                minDis = distance[j];            }        //当已不再存在路径时算法结束；此语句对非连通图是必须的        if(minDis == MaxWeight) return ;        s[u] = 1;//标记结点u已从集合T加入到集合S中        for(j = 0;j < n; j++)            if(s[j] == 0 && G.edge[j][j] < MaxWeight && distance[u] + G.edge[u][j] < distance[j])            {                //结点v0经结点u到其他结点的最短距离和最短路径                distance[j] = distance[u] + G.edge[u][j];                path[j] = u;            }    }}

main.c

#include <stdio.h>#include <stdlib.h>#include <malloc.h>typedef char DataType;#define MaxSize 100#define MaxVertices 10#define MaxWeight 10000#include "AdjMGraph.h"#include "AdjMGraphCreate.h"#include "Dijkstra.h"void main(void){    AdjMGraph g;    char a[] = {        'A','B','C','D','E','F'    };    RowColWeight rcw[] = {        {0,2,5},        {0,3,30},        {1,0,2},        {1,4,8},        {2,1,15},        {2,5,7},        {4,3,4},        {5,3,10},        {5,4,18},    };    int i,n = 6,e = 9;    int distance[6],path[6];    CreatGraph(&g,a,n,rcw,e);    Dijkstra(g,0,distance,path);    printf("从结点%c到其他各结点的最短距离为：\n",g.Vertices.list[0]);    for(i = 0;i < n; i++)        printf("到结点%c的最短距离为：%d\n",g.Vertices.list[i],distance[i]);    printf("从结点%c到其他各结点最短路径的前一结点为：\n",g.Vertices.list[0]);    for(i = 0;i < n; i++)        if(path[i] != -1)            printf("到结点%c的前一结点为：%c\n",g.Vertices.list[i],g.Vertices.list[path[i]]);}/*从结点A到其他各结点的最短距离为：到结点A的最短距离为：0到结点B的最短距离为：20到结点C的最短距离为：5到结点D的最短距离为：22到结点E的最短距离为：28到结点F的最短距离为：12从结点A到其他各结点最短路径的前一结点为：到结点B的前一结点为：C到结点C的前一结点为：A到结点D的前一结点为：F到结点E的前一结点为：B到结点F的前一结点为：CPress any key to continue*/

弗洛伊德算法

void Floyd(int cost[][N],int n,int weight[][N],int path[][N]){    int i,j,k;    //初始化    for(i = 0;i < n; i++)        for(j = 0;j < n; j++)        {            weight[i][j] = cost[i][j];            path[i][j] = -1;        }    //n次递推    for(k = 0;k < n; k++)    {        for(i = 0;i < n; i++)        {            for(j = 0;j < n; j++)            {                if(weight[i][j] > weight[i][k] + weight[k][j])                {                    //得到新的最短路径长度数值                    weight[i][j] = weight[i][k] + weight[k][j];                    //得到该最短路径经过的结点序号                    path[i][j] = k;                }            }        }    }}