数据结构与算法学习 (13)最小生成树

把构成连通⽹的最⼩代价的⽣成树称为最⼩⽣成树。
连通图:在无向图中,若任意两个顶点vi和vj之间都有路径想通,则称该无向图为连通图。
连通图的生成树:一个连通图的生成树是指一个连通子图,它包含有图中全部n个顶点,但只有足以构成一棵树的n-1条边。一颗有n个顶点的生成树有且仅有n-1条边,如果生成树中再添加一条边,必定成环。
最小生成树:在一给定的无向图G=(V,E)中,(u,v)代表连接顶点u与顶点v的边,而w(u,v)代表边的权重,若存在T为E的子集且为无循环图,使得W(T)最小,则此T为G的最小生成树。把构成连通网的最⼩代价的生成树称为最⼩生成树Minimum Spanning Tree,MST。

最小生成树性质:设G=(V,E)是一个连通网络,T1是顶点集V的一个非空真子集,T2是顶点集V-T1后剩余的子集。若(u,v)是G中一条一个端点在T1中(例如:u∈T1),另一个端点在T2中的边(例如:v∈T2),且(u,v)具有最小权值,则一定存在G的最小生成树包括此边(u,v)。

我们先把图转换成邻接矩阵,(https://juejin.im/post/5eae34e951882560d56d6588)

image

利用两个数组来辅助我们记录遍历顶点过程中的数据,两个数组lowcost和arjvex,

  • lowcost用于记录当前顶点与所有与它关联顶点间的权重值。下标:顶点下标,值:与当前顶点的权重。
  • arjvex用于记录当前顶点与哪个顶点相连接。下标:当前顶点下标,值:与当前顶点相连的前一个顶点。
  1. 第一次执行
    image

    <figcaption></figcaption>

  2. 第二次执行
    image

    <figcaption></figcaption>

  3. 第三次执行
    image

    <figcaption></figcaption>

  4. 第四次执行
    image

    <figcaption></figcaption>

  5. 第五次执行
    image

Prim算法
思路

定义两个数组lowcost和adjvex。其中:

lowcost保存相关顶点间边的权值,lowcost[i]=0代表顶点vi已经加入了最小生成子树中,lowcost[i]=INFINITYC表示i对应的顶点与当前的最小生成树中的顶点暂时不连通,在0<lowcost[i]<INFINITYC的情况中会生成一个与当前最小生成子树的各个顶点之间权值最小的顶点加入到最小生成子树当中。
adjvex保存相关顶点下标 adjvex[i]默认=-1.表示i对应的结点和最小生成子树暂时不连通,adjvex[i] >-1 表示i对应的结点和最小生成树相连的顶点为adjvex[i]

初始化上述数组,默认将v0加入到最小子树中,并更新数组中对应的值
以v0为最小子树,寻找剩余顶点中与v0路径最短的顶点加入v0当中,组成新的最小子树,并更新两个数组的值。
以新的最小子树开始,重复步骤3,直到找到最后一个顶点。

void MiniSpanTree_Prim_1(MGraph G) {
    int lowcost[MAXVEX];
    int adjvex[MAXVEX];
    
    //初始化两个数组,即默认将v0放入最小生成子树中
    adjvex[0] = 0;
    lowcost[0] = 0;
    int I;
    //将v0对应的连通顶点信息存入adjvex,将连通顶点的权值存入lowcost
    for(i = 1; i < G.numVertexes; i++) {
        if(G.arc[0][i] < INFINITYC) {
            adjvex[i] = 0;
        }else {
            adjvex[i] = -1;
        }
        lowcost[i] = G.arc[0][I];
    }
    
    int k = 0;
    int j = 0;
    int sum = 0;
    for(i = 1; i < G.numVertexes; i++) {
        int min = INFINITYC;
        //第一次遍历找到lowcost中除去0以外最小的顶点,作为下一个顶点加入最小生成树
        for(j = 1; j < G.numVertexes; j++) {
            if(lowcost[j] != 0 && lowcost[j] < min) {
                min = lowcost[j];
                k = j;
            }
        }
//        printf("lowcost为:\n");
//        for(int temp = 0; temp < G.numVertexes; temp++) {
//            printf("%d ",lowcost[temp]);
//        }
//        printf("\n");
//        printf("adjvex为:\n");
//        for(int temp = 0; temp < G.numVertexes; temp++) {
//            printf("%d ",adjvex[temp]);
//        }
//        printf("\n");
        
        /* 打印当前顶点边中权值最小的边 */
        printf("(V%d, V%d)=%d\n", adjvex[k], k ,G.arc[adjvex[k]][k]);
        sum += G.arc[adjvex[k]][k];
        lowcost[k] = 0;
        //将与新加入顶点连通的其他顶点加入lowcost中并更新对应的lowcost以及adjvex的值
        for(j = 0; j < G.numVertexes; j++) {
            if(lowcost[j] != 0 && G.arc[k][j] < lowcost[j]) {
                lowcost[j] = G.arc[k][j];
                adjvex[j] = k;
            }
        }
    }
    printf("sum = %d\n",sum);
}

克鲁斯卡尔算法
思路:

T是边的集合,初始状态为空
将图中的所有边以权值从小到大的顺序进行排序,预定排序后边的起点永远为比较小的顶点
取出权值最小的边,查看它与T中的边是否造成环路,
如果未构成环路,则加入T中;否则,抛弃改边
如果还有剩余边,则返回3中,否则程序结束

证明:
分两步来证明Kruskal算法得到的是最小生成树

Kruskal得到的是一棵生成树
Kruskal得到的生成树代价最小

假设改算法得到的不是生成树,那么有两种情况,第一种是得到的图是有环的,第二种情况是得到的图不连通。由于算法要求,每次加入的边都不能形成环,因此第一种情形不存在。第二种情况,首先我们知道原图是属于连通图,根据算法的要求,取出权值最小的边,如果未构成环路,则加入T中,因为原图是连通的,必然会存在可以连接所有顶点又不构成环路的路径,可以得到Kruskal得到的图肯定也是连通图,因此第二种情况也不存在。

接下来证明Kruskal得到的生成树代价最小。假设图有n个顶点,则生成树一定有n-1条边,最少存在一个最小生成树U。假设:

Kruskal得到的生成树为T。
U和T中不同的边数为k,其他n-1-k条边相同,相同的边构成集合E
在T中而不在U中的边按代价从小到大一次为a1<a2<..<ak
在U中而不在T中的边按代价从小到大一次为x1<x2<..<xk

我们通过把T中的边依次移动到U中来证明U和T是等价的。
首先将a1移动到U中,由于U是一棵生成树,加入一条新边后会行程环路,且这条回路一定会包括x1,x2...,xk中的边,否则a1和E中的边形成环路,则T不是生成树,产生矛盾。
在回路中删除属于x1,x2...,xk且代价最大的边xi构成一棵新的生成树V。
假设a1<xi,则w(V)<w(U),与最小生成树定义矛盾。
假设a1>xi, 则xi小于a1,a2...,ak,根据Kruskal算法,在考虑最小边的时候会优先考虑xi,而且xi和E也不会构成环路,但是T中并没有xi,说明a1也不会大于xi,即a1 = xi,新得到的w(V)=w(U)。
同理,我们将a2...,ak依次加入到U中,最终得到w(U)=w(T)。证明完毕。
判断是否生成环路

171f3fbbb3f1eeb8.png

上图中,边(v0,v1)在一条环路中,我们分别从v0和v1开始,因为是环路,我们通过某种规则,最终会达到同样的终点。
判断环路需要一个辅助数组parent[]={0},存储新加入顶点的路径关系,值为0代表没有新加入的顶点

新增加一个边的结构体来保存每一条边的信息
/* 对边集数组Edge结构的定义 */
typedef struct
{
    int begin;
    int end;
    int weight;
}Edge;

/* 查找连线顶点的尾部下标 */
//根据顶点f以及parent 数组,可以找到当前顶点的尾部下标; 帮助我们判断2点之间是否存在闭环问题;
int Find(int *parent, int f)
{
    while ( parent[f] > 0)
    {
        f = parent[f];
    }
    return f;
}

/* 交换权值以及头和尾 */
void Swapn(Edge *edges,int i, int j)
{
    int tempValue;
    
    //交换edges[i].begin 和 edges[j].begin 的值
    tempValue = edges[i].begin;
    edges[i].begin = edges[j].begin;
    edges[j].begin = tempValue;
    
    //交换edges[i].end 和 edges[j].end 的值
    tempValue = edges[i].end;
    edges[i].end = edges[j].end;
    edges[j].end = tempValue;
    
    //交换edges[i].weight 和 edges[j].weight 的值
    tempValue = edges[i].weight;
    edges[i].weight = edges[j].weight;
    edges[j].weight = tempValue;
}

/* 对权值进行排序 */
void sort(Edge edges[],MGraph *G)
{
    //对权值进行排序(从小到大)
    int i, j;
    for ( i = 0; i < G->numEdges; i++)
    {
        for ( j = i + 1; j < G->numEdges; j++)
        {
            if (edges[i].weight > edges[j].weight)
            {
                Swapn(edges, i, j);
            }
        }
    }
    
    printf("边集数组根据权值排序之后的为:\n");
    for (i = 0; i < G->numEdges; i++)
    {
        printf("(%d, %d) %d\n", edges[i].begin, edges[i].end, edges[i].weight);
    }
    
}

/* 生成最小生成树 */
void MiniSpanTree_Kruskal(MGraph G)
{
    int i, j, n, m;
    int sum = 0;
    int k = 0;
    /* 定义一数组用来判断边与边是否形成环路
     用来记录顶点间的连接关系. 通过它来防止最小生成树产生闭环;*/
    
    int parent[MAXVEX];
    /* 定义边集数组,edge的结构为begin,end,weight,均为整型 */
    Edge edges[MAXEDGE];
    
    /*1. 用来构建边集数组*/
    for ( i = 0; i < G.numVertexes-1; i++)
    {
        for (j = i + 1; j < G.numVertexes; j++)
        {
            //如果当前路径权值 != ∞
            if (G.arc[i][j]<INFINITYC)
            {
                //将路径对应的begin,end,weight 存储到edges 边集数组中.
                edges[k].begin = i;
                edges[k].end = j;
                edges[k].weight = G.arc[i][j];
                
                //边集数组计算器k++;
                k++;
            }
        }
    }
    
    //2. 对边集数组排序
    sort(edges, &G);
    
    
    //3.初始化parent 数组为0. 9个顶点;
    // for (i = 0; i < G.numVertexes; i++)
    for (i = 0; i < MAXVEX; i++)
        parent[i] = 0;
    
    //4. 计算最小生成树
    printf("打印最小生成树:\n");
    /* 循环每一条边 G.numEdges 有15条边*/
    for (i = 0; i < G.numEdges; i++)
    {
        //获取begin,end 在parent 数组中的信息;
        //如果n = m ,将begin 和 end 连接,就会产生闭合的环.
        n = Find(parent,edges[i].begin);
        m = Find(parent,edges[i].end);
        //printf("n = %d,m = %d\n",n,m);
        
        /* 假如n与m不等,说明此边没有与现有的生成树形成环路 */
        if (n != m)
        {
            /* 将此边的结尾顶点放入下标为起点的parent中。 */
            /* 表示此顶点已经在生成树集合中 */
            parent[n] = m;
            
            /*打印最小生成树路径*/
            printf("(%d, %d) %d\n", edges[i].begin, edges[i].end, edges[i].weight);
            sum += edges[i].weight;
        }
    }
    
    printf("sum = %d\n",sum);
}

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 206,126评论 6 481
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 88,254评论 2 382
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 152,445评论 0 341
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 55,185评论 1 278
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 64,178评论 5 371
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,970评论 1 284
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,276评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,927评论 0 259
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 43,400评论 1 300
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,883评论 2 323
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,997评论 1 333
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,646评论 4 322
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,213评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,204评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,423评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,423评论 2 352
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,722评论 2 345