最小生成树(Prim、Kruskal)算法,秒懂!

简介: 在数据结构与算法的图论中,(生成)最小生成树算法是一种常用并且和生活贴切比较近的一种算法。但是可能很多人对概念不是很清楚,什么是最小生成树?

前言



在数据结构与算法的图论中,(生成)最小生成树算法是一种常用并且和生活贴切比较近的一种算法。但是可能很多人对概念不是很清楚,什么是最小生成树?


一个有 n 个结点的连通图的生成树是原图的极小连通子图,且包含原图中的所有 n 个结点,并且有保持图连通的最少的边。 最小生成树可以用kruskal(克鲁斯卡尔)算法或prim(普里姆)算法求出。


通俗易懂的讲就是最小生成树包含原图的所有节点而只用最少的边和最小的权值距离。因为n个节点最少需要n-1个边联通,而距离就需要采取某种策略选择恰当的边。


学习最小生成树实现算法之前我们要先高清最小生成树的结构和意义所在。咱么首先根据一些图更好的祝你理解。


一个故事


在城市道路规划中,是一门很需要科学的研究(只是假设学习不必当真)。在公路时代城市联通的主要矛盾是时间慢,而造价相比运输时间是次要矛盾。所以在公路时代我们尽量使得城市能够直接联通,缩短城市联系时间,而稍微考虑建路成本!随着科技发展、高级铁路、信息传输相比公路运输快非常非常多,从而事件的主要矛盾从运输时间转变为造价成本,故有时会关注联通所有点的路程(最短),这就用到最小生成树算法。


城市道路铺设可能经历以下几个阶段。


初始,各个城市没有高速公路(铁路)。

政府打算各个城市铺设公路(铁路),每个城市都想成为交通枢纽,快速到达其他城市,但每个城市都有这种想法,如果实现下去造价太昂贵。并且造成巨大浪费。

最终国家选择一些主要城市进行联通,有个别城市只能稍微绕道而行,而绕道太远的、人流量多的国家考虑新建公路(铁路),适当提高效率。



不过上面铁路规划上由于庞大的人口可能不能够满足与"有铁路"这个需求,人们对速度、距离、直达等条件一直在追求中……


但是你可以想象这个场景:有些东西造假非常非常昂贵,使用效率非常高,我这里假设成黄金镶钻电缆 铺设,所以各个城市只要求不给自己落下,能通上就行(没人敢跳了吧)。


要从有环图中选取代价和最小的路线一方面代价最小(总距离最小最省黄金)另一方面联通所有城市。


然而根据上图我们可以得到以下最小生成树,但是最么生成这个最小生成树,就是下面要讲的了。


20191007130031340.png


而类似的还有局部区域岛屿联通修桥,海底通道这些高成本的都多多少少会运用。


Kruskal算法



上面介绍了最小生成树是什么,现在需要掌握和理解最小生成树如何形成。给你一个图,用一个规则生成一个最小生成树。而在实现最小生成树方面有prim和kruskal算法,这两种算法的策略有所区别,但是时间复杂度一致。


Kruskal算法,和前面讲到的并查集关系很大,它的主要思想为:


先构造一个只含 n 个顶点、而边集为空的子图,把子图中各个顶点看成各棵树上的根结点,之后,从网的边集 E 中选取一条权值最小的边,若该条边的两个顶点分属不同的树,则将其加入子图,即把两棵树合成一棵树,反之,若该条边的两个顶点已落在同一棵树上,则不可取,而应该取下一条权值最小的边再试之。依次类推,直到森林中只有一棵树,也即子图中含有 n-1 条边为止。


简而言之,Kruskal算法进行调度的单位是边,它的信仰为:所有边能小则小,算法的实现方面要用到并查集判断两点是否在同一集合。


而算法的具体步骤为:


将图中所有边对象(边长、两端点)依次加入集合(优先队列)q1中。初始所有点相互独立。


取出集合(优先队列)q1最小边,判断边的两点是否联通。


如果联通说明两个点已经有其它边将两点联通了,跳过,如果不连通,则使用union(并查集合并)将两个顶点合并,这条边被使用(可以储存或者计算数值)。

重复2,3操作直到集合(优先队列)q1为空。此时被选择的边构成最小生成树。


c90f303bd8fb435faaa2cfa7a1925ea2.png


Prim算法



除了Kruskal算法以外,普里姆算法(Prim算法)也是常用的最小生成树算法。虽然在效率上差不多。但是贪心的方式和Kruskal完全不同。


prim算法的核心信仰是:从已知扩散寻找最小。它的实现方式和Dijkstra算法相似但稍微有所区别,Dijkstra是求单源最短路径,而每计算一个点需要对这个点重新更新距离,而prim不用更新距离。直接找已知点的邻边最小加入即可!prim和kruskal算法都是从边入手处理。


对于具体算法具体步骤,大致为:


寻找图中任意点,以它为起点,它的所有边V加入集合(优先队列)q1,设置一个boolean数组bool[]标记该位置(边有两个点,每次加入没有被标记那个点的所有边)。


从集合q1找到距离最小的那个边v1并 判断边是否存在未被标记的一点p ,如果p不存在说明已经确定过那么跳过当前边处理,如果未被标(访问)记那么标记该点p,并且与p相连的未知点(未被标记)构成的边加入集合q1, 边v1(可以进行计算距离之类,该边构成最小生成树) .


重复1,2直到q1为空,构成最小生成树 !

大体步骤图解为:


5d1226d9c1664b759c06a1a8e20fd493.pngedc7bf34d4ce444abd0f6427afdcce8a.png


因为prim从开始到结束一直是一个整体在扩散,所以不需要考虑两棵树合并的问题,在这一点实现上稍微方便了一点。


当然,要注意的是最小生成树并不唯一,甚至同一种算法生成的最小生成树都可能有所不同,但是相同的是无论生成怎样的最小生成树:


能够保证所有节点连通(能够满足要求和条件)

能够保证所有路径之和最小(结果和目的相同)

最小生成树不唯一,可能多样的



代码实现



上面分析了逻辑实现。下面我们用代码简单实现上述的算法。


prim


package 图论;
import java.util.ArrayList;
import java.util.Arrays;
import java.util.Comparator;
import java.util.List;
import java.util.PriorityQueue;
import java.util.Queue;
public class prim {
  public static void main(String[] args) {
    int minlength=0;//最小生成树的最短路径长度
    int max=66666;
    String cityname[]= {"北京","武汉","南京","上海","杭州","广州","深圳"};
    int city[][]= {
        { max, 8, 7, max, max, max, max }, //北京和武汉南京联通
        { 8, max,6, max,9, 8,max }, //武汉——北京、南京、杭州、广州
        { 7, 6, max, 3,4, max,max }, //南京——北京、武汉、上海、杭州
        { max, max,3, max,2, max,max }, //上海——南京、杭州
        { max, 9,4, 2,max, max,10 }, //杭州——武汉、南京、上海、深圳
        { max, 8,max, max,max, max,2 }, //广州——武汉、深圳
        { max, max,max, max,10,2,max }//深圳——杭州、广州
    };// 地图
    boolean istrue[]=new boolean[7];
    //南京
    Queue<side>q1=new PriorityQueue<side>(new Comparator<side>() {
      public int compare(side o1, side o2) {
        // TODO Auto-generated method stub
        return o1.lenth-o2.lenth;
      }
    });
    for(int i=0;i<7;i++)
    {
      if(city[2][i]!=max)
      {
        istrue[2]=true;
        q1.add(new side(city[2][i], 2, i));
      }
    }   
    while(!q1.isEmpty())
    {
      side newside=q1.poll();//抛出
      if(istrue[newside.point1]&&istrue[newside.point2])
      {
        continue;
      }
      else {
        if(!istrue[newside.point1])
        {
          istrue[newside.point1]=true;
          minlength+=city[newside.point1][newside.point2];
          System.out.println(cityname[newside.point1]+" "+cityname[newside.point2]+" 联通");
          for(int i=0;i<7;i++)
          {
            if(!istrue[i])
            {
              q1.add(new side(city[newside.point1][i],newside.point1,i));
            }
          }
        }
        else {
          istrue[newside.point2]=true;
          minlength+=city[newside.point1][newside.point2];
          System.out.println(cityname[newside.point2]+" "+cityname[newside.point1]+" 联通");
          for(int i=0;i<7;i++)
          {
            if(!istrue[i])
            {
              q1.add(new side(city[newside.point2][i],newside.point2,i));
            }
          }
        }
      }
    }
    System.out.println(minlength);    
  }
  static class side//边
  {
    int lenth;
    int point1;
    int point2;
    public side(int lenth,int p1,int p2) {
      this.lenth=lenth;
      this.point1=p1;
      this.point2=p2;
    }
  }
}


输出结果:


上海 南京 联通

杭州 上海 联通

武汉 南京 联通

北京 南京 联通

广州 武汉 联通

深圳 广州 联通

28


Kruskal:


package 图论;
import java.util.Comparator;
import java.util.PriorityQueue;
import java.util.Queue;
import 图论.prim.side;
/*
 * 作者:bigsai(公众号)
 */
public class kruskal {
  static int tree[]=new int[10];//bing查集
  public static void init() {
    for(int i=0;i<10;i++)//初始
    {
      tree[i]=-1;
    }
  }
  public static int search(int a)//返回头节点的数值
  {
    if(tree[a]>0)//说明是子节点
    {
      return tree[a]=search(tree[a]);//路径压缩
    }
    else
      return a;
  }
  public static void union(int a,int b)//表示 a,b所在的树合并小树合并大树(不重要)
  {
    int a1=search(a);//a根
    int b1=search(b);//b根
    if(a1==b1) {//System.out.println(a+"和"+b+"已经在一棵树上");
    }
    else {
    if(tree[a1]<tree[b1])//这个是负数,为了简单减少计算,不在调用value函数
    {
      tree[a1]+=tree[b1];//个数相加  注意是负数相加
      tree[b1]=a1;       //b树成为a的子树,直接指向a;
    }
    else
    {
      tree[b1]+=tree[a1];//个数相加  注意是负数相加
      tree[a1]=b1;       //b树成为a的子树,直接指向a;
    }
    }
  }
  public static void main(String[] args) {
    // TODO Auto-generated method stub
    init();
    int minlength=0;//最小生成树的最短路径长度
    int max=66666;
    String cityname[]= {"北京","武汉","南京","上海","杭州","广州","深圳"};
    boolean jud[][]=new boolean[7][7];//加入边需要防止重复 比如 ba和ab等价的
    int city[][]= {
        { max, 8, 7, max, max, max, max }, 
        { 8, max,6, max,9, 8,max }, 
        { 7, 6, max, 3,4, max,max }, 
        { max, max,3, max,2, max,max }, 
        { max, 9,4, 2,max, max,10 }, 
        { max, 8,max, max,max, max,2 }, 
        { max, max,max, max,10,2,max }
    };// 地图
    boolean istrue[]=new boolean[7];
    //南京
    Queue<side>q1=new PriorityQueue<side>(new Comparator<side>() {//优先队列存边+
      public int compare(side o1, side o2) {
        // TODO Auto-generated method stub
        return o1.lenth-o2.lenth;
      }
    });
    for(int i=0;i<7;i++)
    {
      for(int j=0;j<7;j++)
      {
        if(!jud[i][j]&&city[i][j]!=max)//是否加入队列
        {
          jud[i][j]=true;jud[j][i]=true;
          q1.add(new side(city[i][j], i, j));
        }
      }
    }
    while(!q1.isEmpty())//执行算法
    {
      side newside=q1.poll();
      int p1=newside.point1;
      int p2=newside.point2;
      if(search(p1)!=search(p2))
      {
        union(p1, p2);
        System.out.println(cityname[p1]+" "+cityname[p2]+" 联通");
        minlength+=newside.lenth;
      }
    }
    System.out.println(minlength);
  }
  static class side//边
  {
    int lenth;
    int point1;
    int point2;
    public side(int lenth,int p1,int p2) {
      this.lenth=lenth;
      this.point1=p1;
      this.point2=p2;
    }
  }
}


输出结果


上海 杭州 联通

广州 深圳 联通

南京 上海 联通

武汉 南京 联通

北京 南京 联通

武汉 广州 联通

28


总结



最小生成树算法理解起来也相对简单,实现起来也不是很难。Kruskal和Prim主要是贪心算法的两种角度。一个从整体开始找最小边,遇到关联不断合并,另一个从局部开始扩散找身边的最小不断扩散直到生成最小生成树。在学习最小生成树之前最好学习一下dijkstra算法和并查集,这样在实现起来能够快一点,清晰一点。


力扣1584就是一个最小生成树的入门题,不过哪个有点区别的就是默认所有点是联通的,所以需要你剪枝优化。这里就不带大家一起看啦,有问题下面也可交流!


最后,如果你那天真的获得一大笔资金去修建这么一条昂贵的黄金路线,可以适当采取此方法,另外剩下的大批,苟富贵,勿相忘。。


目录
相关文章
|
5月前
|
算法 索引
class061 最小生成树【算法】
class061 最小生成树【算法】
53 0
|
4天前
|
算法 索引
数据结构与算法-最小生成树入门
数据结构与算法-最小生成树入门
10 0
|
15天前
|
存储 机器学习/深度学习 算法
上机实验三 图的最小生成树算法设计 西安石油大学数据结构
上机实验三 图的最小生成树算法设计 西安石油大学数据结构
19 1
|
3月前
|
算法 Java C++
第十四届蓝桥杯集训——练习解题阶段(无序阶段)-ALGO-6 算法训练 安慰奶牛 最小生成树
第十四届蓝桥杯集训——练习解题阶段(无序阶段)-ALGO-6 算法训练 安慰奶牛 最小生成树
23 0
|
5月前
|
算法 搜索推荐
Kruskal算法
Kruskal算法
|
6月前
|
算法 C++
用prim和kruskal算法求最小生成树问题
用prim和kruskal算法求最小生成树问题
47 0
|
10天前
|
机器学习/深度学习 人工智能 算法
基于DCT和扩频的音频水印嵌入提取算法matlab仿真
本文介绍了结合DCT和扩频技术的音频水印算法,用于在不降低音质的情况下嵌入版权信息。在matlab2022a中实现,算法利用DCT进行频域处理,通过扩频增强水印的隐蔽性和抗攻击性。核心程序展示了水印的嵌入与提取过程,包括DCT变换、水印扩频及反变换步骤。该方法有效且专业,未来研究将侧重于提高实用性和安全性。
|
25天前
|
机器学习/深度学习 算法
【MATLAB】GA_BP神经网络时序预测算法
【MATLAB】GA_BP神经网络时序预测算法
35 8
|
29天前
|
机器学习/深度学习 算法 Serverless
【MATLAB】PSO_BP神经网络回归预测算法(适用光伏发电回归预测等)
【MATLAB】PSO_BP神经网络回归预测算法(适用光伏发电回归预测等)
30 1
|
1天前
|
算法 TensorFlow 算法框架/工具
基于直方图的图像阈值计算和分割算法FPGA实现,包含tb测试文件和MATLAB辅助验证
这是一个关于图像处理的算法实现摘要,主要包括四部分:展示了四张算法运行的效果图;提到了使用的软件版本为VIVADO 2019.2和matlab 2022a;介绍了算法理论,即基于直方图的图像阈值分割,通过灰度直方图分布选取阈值来区分图像区域;并提供了部分Verilog代码,该代码读取图像数据,进行处理,并输出结果到&quot;result.txt&quot;以供MATLAB显示图像分割效果。