图算法--最小生成树算法的实现与分析

时间  2019-12-10
标签 算法 最小 生成树 实现 分析 繁體版
原文   原文链接

图是一种灵活的数据结构，它多用于描述对象之间的关系和链接模型。html

关于图的算法：最小生成树、最短路径、旅行商问题以及许多其余算法大都会使用到广度优先搜索和深度优先搜索，由于它们提供了一套系统地访问图数据结构的方法。算法

带权图，是指图的每条边上带有一个值或权，这些权用一个小的数字标记在边上。不少条件因素均可以做为权值，但一般它表示遍历这条边所产生的代价。数据结构

最小生成树简述

咱们作一个简单的模型，在一块木板上钉上一些钉子并用一些细绳链接起来，假设每个钉子都经由一根或多根细绳链接起来。如今咱们一根一根拿走细绳，直到用最少的细绳将全部钉子链接起来。这个模型的背后思想就是最小生成树。函数

正式的表述法是，给定一个无方向的带权图G=（V，E），最小生成树为集合T，T是以最小代价链接V中全部顶点所用边E的最小集合。集合T中的边能造成一颗树，这是由于每一个顶点都能向上找到它的一个父节点（根结点除外）。spa

Prim算法

Prim算法是一种产生最小生成树的方法。3d

Prim算法从任意一个顶点开始，每次选择一个与当前顶点最近的一个顶点，并将两个顶点之间的边加入到树中。code

从根本上讲，Prim算法就是不断的选择顶点，并计算边的权值，同时判断是否还有更有效的链接方式。该算法相似广度优先搜索算法，由于在往图中更深的顶点探索以前，它首先要遍历与此顶点相关的全部顶点。每一个阶段都要决定选择哪一个顶点，因此须要维护顶点的颜色和键值。htm

开始，将全部顶点的色值设置为白色，键值设置为∞（它表明一个足够大的数，大于图中全部边的权值）。同时，将起始顶点的键值设置为0。随着算法的不断演进，在最小生成树中为每一个顶点（除起始顶点外）指派一个父节点。只有当顶点的色值变为黑色时，此顶点才是最小生成树的一部分。对象

Prim算法的运行过程以下：blog

首先，在图中全部的白色顶点中，选择键值最小的顶点u。开始，键值被设置为0的那一顶点将做为起始顶点。

当选择此顶点以后，将其标记为黑色。

接下来，对于每一个与u相邻的顶点v，设置v的键值为边（u，v）的权值，同时将u设置为v的父结点。

重复这个过程，直到全部的顶点都标记为黑色。

随着最小生成树的增加，该树会包含图中的全部的边（链接全部顶点最少数量边），且每条边的两端都有一个黑色的顶点。

下图展现了最小生成树的产生过程。在图中，键值和父结点都显示在每一个顶点的旁边，用斜线分开。键值显示在斜线的左边，父结点显示在斜线的右边。浅灰色的边是最小生成树增加过程当中的边。图中最小生成树总的权值为17。

 

最小生成树的接口定义

mst

int mst(Graph *graph, const MstVertex *start, List *span, int (*match)(const void *key1, const void *key2));

返回值：若是计算最小生成树成功，返回0，不然返回-1。

描述：为一个无方向的带权图graph计算最小生成树。

最小生成树从顶点start开始。

此操做会改变graph，因此若是有必要，在调用此操做以前先对图进行备份。

graph中的每一个顶点必须包含MstVertex类型的数据。经过设置MstVertex结构体中的成员weight的值来指定每一个边的权值，weight的值由传入graph_ins_edge的参数data2决定。用MstVertex结构体的成员data来保存与顶点相关的数据。

graph的match函数（此函数在用graph_init对图进行初始化时调用）用来比较MstVertex结构体中的data成员。此函数与传入mst中的参数match相同。

一旦计算完成，最小生成树的相关数据将会返回到span。span是存储MstVertex结构体的列表。在span中，父结点为NULL的顶点为最小生成树的根结点。其余每一个顶点的parent成员都指向span中位于该顶点以前的那个顶点。

span中的顶点指向graph中的实际顶点，因此只要可以访问span，函数调用者就必须保证graph中内存空间有效。一旦再也不使用span，就调用list_destroy销毁span。

复杂度：O（E V2），其中V是图中顶点的个数，E是边的条数。

最小生成树的实现与分析

为了计算一个无方向的带权图的最小生成树，首先，咱们要用表示图的基本抽象数据类型来表示带权图。同时，Prim算法还须要一种追踪顶点和边信息的方法。这就用到了MstVertex结构体：它用来为图中的顶点计算最小生成树。此结构包含5个成员：data是与顶点相关的数据；weight是到达该顶点的边的权值；color是顶点的颜色；key是顶点的键值；parent是最小生成树中顶点的父结点。

创建一个包含MstVertex结构体的图的过程几乎与创建一个包含其余类型的图的过程同样：要将一个顶点插入图中，调用graph_ins_vertex，并将MstVertex结构体传入data。相似地，要将一条边插入图中，调用函数graph_ins_edge，并将MstVertex结构体传入data1和data2。当插入一个顶点时，只设置MstVertex结构体的data成员。当插入一条边时，设置data1的data成员，data2的data和weight成员。在data2中，weight是边的权值，此边是data1中的顶点到data2中顶点的链接线。在实际中，权值一般用浮点数进行存储和计算。因为键值是由权值计算来的，所以键值也用浮点数表示。

mst操做首先初始化邻接表结构链表中的每一个顶点。将每一个顶点的键值初始化为DBL_MAX（除超始顶点外，超始顶点初始值为0.0）。在图的抽象数据类型中，图由一个邻接表结构链表来表示。每一个邻接表包含一个顶点和一个相邻顶点的集合。用存储在邻接表结构中的顶点来维护顶点的色值、键值、和父结点。维护邻接表结构链表中信息的关键是能将这些信息存储起来，而不是仅仅列出与本身相邻的顶点。鉴于一个顶点可能会出如今众多的邻接表中，因此每一个顶点只能在邻接表结构链表中出现一次。

Prim算法的核心是用一个单循环为图中的每一个结点迭代一次。在每次迭代的过程当中：

首先，在全部的白色顶点中选择键值最小的顶点。同时，在邻接表结构链表把此顶点涂黑。

接下来，遍历与所选顶点相邻的顶点。在遍历每一个顶点时，检查它在邻接表结构链表中的颜色和键值。一旦获取了这些信息，就将它与所选顶点的颜色和键值进行比较。若是相邻顶点是白色，且其键值比所选顶点的小，就将所选顶点与相邻顶点之间边的权值设置为相邻顶点的键值；同时，将相邻顶点的父结点设置为所选顶点。

而后，更新存储在邻接表结构链表中相邻顶点的信息。

重复这个过程，直到全部顶点都涂黑。

一旦Prim算法中的主循环结束，最小生成树也就计算完成了。此时，将邻接表结构链表中的每一个黑色MstVertex结构体插入到链表span中。在span中，父结点为NULL的顶点就是最小生成树的根结点。其余每一个顶点的parent成员都指向span中位于该顶点以前的那个顶点。每一个MstVertex结构体的成员weight并不常用，由于它只有在存储到邻接表中时才用的到。

下图显示了上面的示例图中计算最小生成树所返回的MstVertex结构体链表：

 

示例：图算法的头文件（含最小生成树、最短路径、旅行商问题三种实现所需函数定义的头文件）

/*graphalg.h*/
#ifndef GRAPHALG_H
#define GRAPHALG_H

#include "graph.h"
#include "list.h"

/*定义最小生成树中结点的数据结构*/
typedef struct MstVertex_
{
    void    *data;
    double  weight;
    
    VertexColor color;
    double      kdy;
    
    struct MstVertex_ *parent;
}MstVertex; 

/*定义最短路径中结点的数据结构*/
typedef struct PathVertex_
{
    void    *data;
    double  weight;
    
    VertexColor  color;
    double       d;
    
    struct PathVertex_  *parent;
}PathVertex;

/*定义旅行商问题中结点的数据结构*/
typedef struct TspVertex_
{
    void    *data;
    
    double  x,y;
    
    VertexColor  color;
}TspVertex;

/*函数接口*/
int mst(Graph *graph, const MstVertex *start, List *span, int (*match)(const void *key1, const void *key2));

int shortest(Graph *graph, const PathVertex *start, List *paths, int(*match)(const void *key1, const void *key2));

int tsp(List *vertexs, const TspVertex *start, List *tour, int (match*)(const void *key1, const void *key2));

#endif // GRAPHALG_H

 

示例：计算最小生成树的实现

/*mst.c*/
#include <float.h>
#include <stdlib.h>
#include "graph.h"
#include "list.h"

/*mst 计算最小生成树函数*/
int mst(Graph *graph, const MstVertex *start, List *span,
        int (*match)(const void *key1, const void *key2))
{
    AdjList   *adjlist;
    MstVertex *mst_vertex,
              *adj_vertex;
    ListElmt  *element,
              *member;
    double    minmum;
    int       found,i;

    /*初始化图中的全部结点*/
    found=0;

    for(element=list_head(&graph_adjlists(graph)); element!=NULL; element = list_next(element))
    {
        mst_vertex = ((AdjList *)list_data(element))->vertex;

        if(match(mst_vertex,start))
        {
            /*匹配到起始顶点，并对其初始化*/
            mst_vertex->color = white;
            mst_vertex->key = 0;
            mst_vertex->parent = NULL;
            found = 1;
        }
        else
        {
            /*非起始顶点的初始化*/
            mst_vertex->color = white;
            mst_vertex->key = DBL_MAX;
            mst_vertex->parent = NULL;
        }
    }
    /*未找到起始顶点，函数返回*/
    if(!found )
        return -1;

    /*运用Prim算法计算最小生成树*/
    i=0;

    while(i<graph_vcount(graph))
    {
        /*选择拥有最小键值的白色顶点*/
        minimum = DBL_MAX;

        for(element = list_head(&graph_adjlists(graph)); element != NULL; element = list_next(element))
        {
            mst_vertex = ((AdjList *)list_data(element))->vertex;

            if(mst_vertex->color == white && mst_vertex->key < minmum)
            {
                minmum = mst_vertex->key;
                adjlist = list_data(element);
            }
        }
        /*将已选择的顶点涂成黑色*/
        ((MstVertex *)adjlist->vertex)->color = black;

        /*遍历被选中顶点的全部邻接顶点*/
        for(member = list_head(&adjlist->adjacent); member != NULL; member = list_next(member))
        {
            adj_vertex = list_data(member);

            for(element = list_head(&graph_adjlists(graph)); element != NULL; elemet = list_next(element))
            {
                mst_vertex = ((AdjList *)list_data(element))->vertex;

                if(match(mst_vertex,adj_vertex))
                {
                    /*决定是否改变该顶点的键值或父结点*/
                    if(mst_vertex->color == white && adj_vertex->weight < mst_vertex->key)
                    {
                        mst_vertex->key = adj_vertex->weight;
                        mst_vertex->parent = adjlist_vertex;
                    }
                    break;
                }
            }
        }
    /*准备选取下一个顶点*/
    i++;
    }


/*加载最小生成树到链表中*/
list_init(span,NULL);
/*从邻接表结构链表中加载每一个黑色结点*/
for(element = list_head(&graph_adjlists(graph)); elemet != NULL; element = list_next(element))
{
    /**/
    mst_vertex = ((AdjList *)list_data(element))->vertex;

    if(mst_vertex->color == black)
    {
        if(list_ins_next(span, list_tail(span),mst_vertex) != 0)
        {
            list_destroy(span);
            return -1;
        }
    }
}
    return 0;
}
相关文章
相关标签/搜索
每日一句
    每一个你不满意的现在，都有一个你没有努力的曾经。
本站公众号
   欢迎关注本站公众号,获取更多信息
联系我们 最近搜索 最新文章 沪ICP备13005482号-10 MyBatis教程 SQL 教程 MySQL教程 Java 教程 Thymeleaf 教程 Hibernate教程 Spring教程 Redis教程