最小生成树-普利姆和克鲁斯卡尔算法
目录
- 最小生成树
- 普利姆算法
- 算法介绍
- 代码
- 克鲁斯卡尔算法
- 算法介绍
- 步骤解析
- 回路
- 代码实现
最小生成树主要是用于解决修路问题等类似问题,要将所有顶点连通,并且权值之和最小。
最小生成树
- 给定一个带权的无向连通图,如何选取一棵生成树,使树上所有边上权的总和为最小,这叫最小生成树
- N个顶点,一定有N-1条边
包含全部顶点 - N-1条边都在图中
- 求最小生成树的算法主要是普里姆算法和克鲁斯卡尔算法
普利姆算法
算法介绍
普利姆(Prim)算法求最小生成树,也就是在包含n个顶点的连通图中,找出只有(n-1)条边包含所有n个顶点的连通子图,也就是所谓的极小连通子图
普利姆的算法如下:
- 设G=(V,E)是连通网,T=(U,D)是最小生成树,V,U是顶点集合,E,D是边的集合
- 若从顶点u开始构造最小生成树,则从集合V中取出顶点u放入集合U中,标记顶点v的visited[u]=1
- 若集合U中顶点ui与集合V-U中的顶点vj之间存在边,则寻找这些边中权值最小的边,但不能构成回路,将顶点vj加入集合U中,将边(ui,vj)加入集合D中,标记visited[vj]=1
- 重复步骤②,直到U与V相等,即所有顶点都被标记为访问过,此时D中有n-1条边
其实,按照我个人的理解,可以简单概括为:对所有已经加入(已经访问过)的顶点,在它们的所有未访问过的边中找到权值weight最小的边,加入结果集,并标记边的另一个顶点为已访问过,一直循环这歌过程,直到加入了n-1条边。
代码
首先要有一个创建和存储图信息的类
class Graph:"""图的构造类"""def __init__(self, vertex_num):""":param vertex_num: 顶点的数量"""self.vertex_num = vertex_numself.vertexs = Noneself.weights = Nonedef create(self, vertexs, weights):""":param vertexs: 所有顶点集合:param weights: 边的邻接矩阵,二维数组:return:"""self.vertexs = vertexsself.weights = weights
然后,开始我们的普利姆算法
class MinimumSpanningTree:"""通过构建最小生成树,解决修路:要将所有顶点连通,并且使用的费用最低即权值之和最小"""def import_graph(self, graph: Graph):self.graph = graphdef prim(self, start):"""利用普利姆算法构建最小生成树:param start::return:"""# 用于判断每个顶点是否已经访问的listvisited = [False] * self.graph.vertex_numvisited[start] = Truefor k in range(1, self.graph.vertex_num):min_weight = 10000i1 = -1i2 = -1# 对所有已经加入(已经访问过)的顶点,在它们的所有未访问过的边中找到权值weight最小的边,加入结果集,并标记为已访问过for i in range(self.graph.vertex_num):if not visited[i]:continuefor j in range(self.graph.vertex_num):if (self.graph.weights[i][j] < min_weight) & (not visited[j]):min_weight = self.graph.weights[i][j]i1 = ii2 = jvisited[i2] = Trueprint('%s -> %s : %d' % (self.graph.vertexs[i1], self.graph.vertexs[i2], min_weight))
测试图例
if __name__ == '__main__':# 普利姆算法测试graph = Graph(7)vertexts = ['A', 'B', 'C', 'D', 'E', 'F', 'G']weights = [[10000, 5, 7, 10000, 10000, 10000, 2],[5, 10000, 10000, 9, 10000, 10000, 3],[7, 10000, 10000, 10000, 8, 10000, 10000],[10000, 9, 10000, 10000, 10000, 4, 10000],[10000, 10000, 8, 10000, 10000, 5, 4],[10000, 10000, 10000, 4, 5, 10000, 6],[2, 3, 10000, 10000, 4, 6, 10000]]graph.create(vertexts, weights)tree = MinimumSpanningTree()tree.import_graph(graph)tree.prim(0)
克鲁斯卡尔算法
算法介绍
克鲁斯卡尔(Kruskal)算法,是用来求加权连通图的最小生成树的算法。
基本思想:按照权值从小到大的顺序选择n-1条边,并保证这n-1条边不构成回路
具体做法:首先构造一个只含n个顶点的森林,然后依权值从小到大从连通网中选择边加入到森林中,并使森林中不产生回路,直至森林变成一棵树为止。
这里其实最关键的就是两点:
- 按照权值从小到大的顺序选择n-1条边
- 选择的边不构成回路
步骤解析
第1步:将边<E,F>加入R中。
边<E,F>的权值最小,因此将它加入到最小生成树结果R中。
第2步:将边<C,D>加入R中。
上一步操作之后,边<C,D>的权值最小,因此将它加入到最小生成树结果R中。
第3步:将边<D,E>加入R中。
上一步操作之后,边<D,E>的权值最小,因此将它加入到最小生成树结果R中。
第4步:将边<B,F>加入R中。
上一步操作之后,边<C,E>的权值最小,但<C,E>会和已有的边构成回路;因此,跳过边<C,E>。同理,跳过边<C,F>。将边<B,F>加入到最小生成树结果R中。
第5步:将边<E,G>加入R中。
上一步操作之后,边<E,G>的权值最小,因此将它加入到最小生成树结果R中。
第6步:将边<A,B>加入R中。
上一步操作之后,边<F,G>的权值最小,但<F,G>会和已有的边构成回路;因此,跳过边<F,G>。同理,跳过边<B,C>。将边<A,B>加入到最小生成树结果R中。
此时,最小生成树构造完成!它包括的边依次是:
回路
我们加入的边的两个顶点不能都指向同一个终点,否则将构成回路。
在将<E,F> <C,D> <D,E>加入到最小生成树R中之后,这几条边的顶点就都有了终点:
(01) C的终点是F。
(02) D的终点是F。
(03) E的终点是F。
(04) F的终点是F。
这里需要注意的是,只有当边加入连通了了之后,才能有所谓的终点
例如仅仅只有<E,F> <C,D> <D,E>加入后,像A、B这样的顶点是还没有终点,也可以认为终点是它们自己。
代码实现
def kruskal(self):"""利用克鲁斯卡尔算法构建最小生成树:return:"""ends = [0] * self.graph.vertex_num # 用于存储每个顶点的终点edges = [] # 用于存储排序后的边# 开始进行排序# 遍历邻接矩阵,将边提取出来的同时插入edges,并保证有序visited = [False] * self.graph.vertex_numfor i in range(self.graph.vertex_num):for j in range(self.graph.vertex_num):if visited[j]:continueif self.graph.weights[i][j] != 10000:index = len(edges)for n in range(len(edges)):if self.graph.weights[i][j] < weights[edges[n][0]][edges[n][1]]:index = nbreakedges.insert(index, (i, j))visited[i] = True# 开始克鲁斯卡尔算法# 选取权值最小的n-1条不构成回路的边,即为最小生成树selects = []for edge in edges:if len(selects) > self.graph.vertex_num - 1:break# 首先判断是否构成回路end1 = self.get_end(edge[0], ends)end2 = self.get_end(edge[1], ends)if end1 == end2:continueelif end1 > end2:ends[end2] = end1 # 更新终点数组selects.append((self.graph.vertexs[edge[0]], self.graph.vertexs[edge[1]]))else:ends[end1] = end2selects.append((self.graph.vertexs[edge[0]], self.graph.vertexs[edge[1]]))print(selects)def get_end(self, v, ends):"""计算输入顶点的终点:param v::param ends: 存储每个顶点的终点:return:"""while ends[v] != 0:v = ends[v]return v
测试代码
graph = Graph(7)
vertexts = ['A', 'B', 'C', 'D', 'E', 'F', 'G']
weights = [[10000, 12, 10000, 10000, 10000, 16, 14],
[12, 10000, 110000, 10000, 10000, 7, 10000],
[10000, 110000, 10000, 3, 5, 6, 10000],
[10000, 10000, 3, 10000, 4, 10000, 10000],
[10000, 10000, 5, 4, 10000, 2, 8],
[16, 7, 6, 10000, 2, 10000, 9],
[14, 10000, 10000, 10000, 8, 9, 10000]]
graph.create(vertexts, weights)
tree = MinimumSpanningTree()
tree.import_graph(graph)
tree.kruskal()
欢迎关注同名公众号:“我就算饿死也不做程序员”。
交个朋友,一起交流,一起学习,一起进步。
最小生成树-普利姆和克鲁斯卡尔算法相关推荐
- 最小生成树————普利姆和克鲁斯卡尔
最小生成树定义: 一个无向图,任意两个顶点都是联通的,并且是一个树,这棵树就叫生成树.如果边上有权值,使边的权值和最小的生成树叫做最小生成树. 求解最小生成树有两个算法,克鲁斯卡尔(Kruskal)算 ...
- 普利姆(prim)算法与迪杰斯特拉(dijestella)算法
一:普利姆(prim)算法 对于临界矩阵表示的图,时间复杂度为O(|V|2) 对于邻接表表示的图,时间复杂度为O(|V|+|E|) 代码: int dist[];//用来记录从原点到各个节点的距离 b ...
- 普利姆(Prim)算法与克鲁斯卡尔(Kruskal)算法
普利姆算法就是求出最小生成树的算法之一 最小生成树:给定一个大权的无向连接图如何选取一颗生成树,使树上的所有边上权的总和为最小,这就叫做最小生成树 例子: 用Prim算法从p城出发 求出生成树: 结果 ...
- 【算法基础12】最小生成树的两种解法(普里姆、克鲁斯卡尔)
一.稠密图:朴素版prim算法 主要思想:每次从树外的结点中找到一个距离树最近的点加入树,将这段路径长度计入最小生成树的路径长度中,然后依据新的树更新树外结点距离树的距离,再次找到最近点加入,直到所有 ...
- [算法]最小生成树-普利姆算法
2019独角兽企业重金招聘Python工程师标准>>> 目前正在看<大话数据结构>,其中介绍了普利姆算法,自己对算法理解能力太差,能够手写求出最小生成树,但是写出算法代码 ...
- 【数据结构】最小生成树(Prim算法,普里姆算法,普利姆)、最短路径(Dijkstra算法,迪杰斯特拉算法,单源最短路径)
文章目录 前置问题 问题解答 一.基础概念:最小生成树的定义和性质 (1)最小生成树(Minimal Spanning Tree)的定义 (2)最小生成树(MST)的性质 二.如何利用MST性质寻找最 ...
- 四、最小生成树——普里姆(Prim)算法
一个连通图的生成树是一个极小的连通子图,它含有图中全部的顶点,但只有足以构成一棵树的n-1条边,那么我们把构造连通图网的最小代价生成树称为最小生成树(就是n个顶点,用n-1条边全部连接起来,并且使得权 ...
- 最小生成树算法普利姆算法和克鲁斯卡尔算法实现
最小生成树算法: 普里姆算法:顶点集合N,辅助顶点集合S,初始化中,将出发点vi加入S,并从N中删除 1.从顶点集合N中找到一条到集合S最近的边(vi,vj),存储该边,并将vj从N移到S中 2.重复 ...
- 【数据结构】图的应用(普利姆算法、克鲁斯卡尔算法、迪杰斯特拉算法、弗洛伊德算法、拓扑排序)
最小生成树 什么是最小生成树 是一棵树 - 无回路 - |V|个顶点一定有|V|-1条边 是生成树 - 包含全部顶点 - |V|-1条边全在图里 贪心算法 什么是"贪":每一步都要 ...
最新文章
- 瑞幸咖啡生,互联网咖啡生;数据死,互联网咖啡死
- php框架 css文件引用,yii框架中怎么引入css文件
- android okhttpclient设置编码,Android之okhttp实现socket通讯(非原创)
- 改变世界的十位算法大师
- jenkins 脚本编写_使用Jenkins 2编写以代码形式编写管道和实现DevOps的简介
- 软硬交互代码示例_matlab交互式程序设计示例:[6]GUI界面编程1
- x11获得窗口名的代码
- Sql Server 2000 挂起怎么办
- Windows操作系统原理笔记
- ShipConstructor.2006v1.00.rar
- Node.js连接MySQL数据库
- Mac读写NTFS驱动程序神器Tuxera NTFS2022
- QQ空间官方账号被黑产利用漏洞分析
- js数组操作方法 concat
- delphi报表_顶级4 Delphi报表工具
- 用sync toy做增量备份
- 假定在使用CSMA/CD协议的10Mbit/s以太网中某个站在发送数据时检测到碰撞,执行退避算法时选择了随机数r=100.试问这个站需要等待多长时间后才能再次发送数据?如果是100Mbit/s的以太网
- 小众软件:录屏局部放大神器 ZoomIt
- NC23053月月查华华的手机
- 给企业微信成员不同人发不同的文件