本文详细介绍了图的最短路径的概念,然后介绍了求最短路径的两种算法:dijkstra算法和floyd算法的原理,最后提供了基于邻接矩阵和邻接表的图对两种算法的java实现。
阅读本文需要一定的图的基础,如果对于图不是太明白的可以看看这篇文章:java数据结构之图的原理与实现。
在生活中,图形结构的应用是最广泛的。比如常见的交通路线选择,站点可以看作顶点,站点之间如果有路径,则算作两点之间的边或者弧,站点之间的通行时间,可以看作边或者弧的权值。
上图就是生活中出行路线的选择映射到图形结构的案例。顶点作为站点,站点之间能够到达则拥有边,站点的之间的通行时间则是边的权值。
对于出行路线的选择,不同的人有不同的选择。其中有一种很常见的选择是要求出发地和目的地之间的总通行时间最短,而不在乎中途到底有几站。毕竟通行时间对于很多人来说是宝贵的!
这样的问题转转换为数学模型,就是求带权图的最短路径,就是求带权图形两顶点之间的权值最小的路径。即如果从图中某一顶点(源点)到达另一顶点(终点)的路径可能不止一条,如何找到一条路径使得沿此路径上各边的权值总和(称为路径长度)达到最小。
实际上最短路径有两重含义,一个两顶点之间的路径数量最少,另一个是两顶点之间的路径距离最短,本次主要解决路径距离最短的问题,即最小权值和。常见的解决算法一般是两种,迪杰斯特拉(dijkstra)算法和弗洛伊德(floyd)算法。
迪杰斯特拉(dijkstra)算法是由荷兰计算机科学家狄克斯特拉于1959 年提出的,因此又叫狄克斯特拉算法。是从一个顶点到其余各顶点的最短路径算法,解决的是寻找给定的加权图中指定顶点间最短路径问题的算法
dijkstra算法并不是一下子就求出了起点到终点的最短路径,而是采用的是贪心算法策略,一步步求出它们之间顶点的最短路径,过程中都是基于已经求出的最短路径的基础上,求得更远顶点的最短路径,最终得到起点和终点的最短路径。
通用步骤如下:
1.指定两个集合s和u。s的作用是记录已求出最短路径的顶点,而u则是记录还未求出最短路径的顶点,以及这些顶点到起始顶点的权。
2.指定一个起始顶点a,存入集合s中,其他顶点以及到顶点a的权存入集合u中,从u中找出并移除路径最短的顶点b,并将其加入到s中,并且更新u中对应的路径权值(更新源点将新加入节点作为中间节点到达其它节点的距离);重复该操作直到遍历所有顶点,此时s中的集合就是起点a到其他各个顶点的最短路径。
迪杰斯特拉算法只支持非负权图,它计算的是单源最短路径,即单个源点到剩余节点的最短路径,时间复杂度为o(n²),对稀疏图运行更快。如果想知道所有顶点到所有顶点的最短路径,那么等于在原有算法的基础上,再来一次循环,此时整个算法的时间复杂度就成了o(n³)。
该案例对应着下面实现代码中的案例,设起始点为a,初始化a到其他点的路径数组{0, 99, 8, 2, 99, 3, 99},初始化标志位数组{true, fal, fal, fal, fal, fal, fal}。
开始第一轮循环,排除已找到的路径,即排除0,寻找到a的最短路径,找到了a-d,即索引为3的顶点,设置对应位置的最短路径标志位为true,更新未获取最短路径的顶点的最短路径,因为其他已找到的顶点的最短路径已经找到了,这里只需要更新新找到的d可达的顶点路径到a点的最短路径,如果经过d点的路径比数组中已存在的最短路径小,那么更新值。
这里d点可达c、b,d-c+a-d=7<8,因此更新a-c的最短路径为7;d-b+a-d=11<99,因此更新a-b的最短路径为11,第一轮循环结束,此时最短路径数组为{0, 11, 7, 2, 99, 3, 99},标志位数组为{true, fal, fal, true, fal, fal, fal}:
开始第二轮循环,排除已找到的路径,即排除0、2,寻找到a的最短路径,这里找到3,即索引为5的顶点,即a-f,设置对应位置的最短路径标志位为true,更新未获取最短路径的顶点的最短路径,因为其他已找到的顶点的最短路径已经找到了,这里只需要更新新找到的f可达的顶点路径到a点的最短路径,如果经过f点的路径比数组中已存在的最短路径小,那么更新值。
这里f点可达g,f-g+a-f=12<99,因此更新a-g的最短路径为12,第二轮循环结束,此时最短路径数组为{0, 11, 7, 2, 99, 3, 12},标志位数组为{true, fal, fal, true, fal, true, fal}。
开始第三轮循环,排除已找到的路径,即排除0、2、3,寻找到a的最短路径,这里找到7,即索引为2的顶点,即a-c,设置对应位置的最短路径标志位为true,更新未获取最短路径的顶点的最短路径,因为其他已找到的顶点的最短路径已经找到了,这里只需要更新新找到的c可达的顶点路径到a点的最短路径,如果经过c点的路径比数组中已存在的最短路径小,那么更新值。
这里c点可达b,c-b+a-c=11 = 11,因此不更新a-b的最短路径,第三轮循环结束,此时最短路径数组为{0, 11, 7, 2, 99, 3, 12},标志位数组为{true, fal, true, true, fal, true, fal}。
开始第四轮循环,排除已找到的路径,即排除0、2、3、7,寻找到a的最短路径,这里找到11,即索引为1的顶点,即a-b,设置对应位置的最短路径标志位为true,更新未获取最短路径的顶点的最短路径,因为其他已找到的顶点的最短路径已经找到了,这里只需要更新新找到的b可达的顶点路径到a点的最短路径,如果经过b点的路径比数组中已存在的最短路径小,那么更新值。
这里b点可达e,b-e+a-b=18 < 99,因此更新a-e的最短路径为18,第四轮循环结束,此时最短路径数组为{0, 11, 7, 2, 18, 3, 12},标志位数组为{true, true, true, true, fal, true, fal}。
开始第五轮循环,排除已找到的路径,即排除0、2、3、7、11,寻找到a的最短路径,这里找到12,即索引为6的顶点,即a-g,设置对应位置的最短路径标志位为true,更新未获取最短路径的顶点的最短路径,因为其他已找到的顶点的最短路径已经找到了,这里只需要更新新找到的g可达的顶点路径到a点的最短路径,如果经过g点的路径比数组中已存在的最短路径小,那么更新值。
排除已找到的顶点,这里g点可达e,g-e+a-g=18 = 18,因此不更新最短路径,第五轮循环结束,此时最短路径数组为{0, 11, 7, 2, 18, 3, 12},标志位数组为{true, true, true, true, fal, true, true}。
开始第六轮循环,排除已找到的路径,即排除0、2、3、7、11、12,寻找到a的最短路径,这里找到18,即索引为4的顶点,即a-e,设置对应位置的最短路径标志位为true,更新未获取最短路径的顶点的最短路径,因为其他已找到的顶点的最短路径已经找到了,这里只需要更新新找到的e可达的顶点路径到a点的最短路径,如果经过e点的路径比数组中已存在的最短路径小,那么更新值。
排除已找到的顶点,这里不更新最短路径,第四轮循环结束,此时最短路径数组为{0, 11, 7, 2, 18, 3, 12},标志位数组为{true, true, true, true, true, true, true}。
此时大循环结束,dijkstra算法结束,顶点a到各个顶点的最短路径已经找到,即a到{a,b,c,d,e,f,g}的最短路径为{0, 11, 7, 2, 18, 3, 12}。
弗洛伊德(floyd)算法又称插点法,是一种利用动态规划的思想寻找给定的加权图中多源点之间最短路径的算法。算出来的结果是所有的节点到其余各节点之间的最短距离。
通用步骤如下:
1.设图顶点数为n。首先需要准备一个长度为n的距离矩阵s,矩阵s中的元素a[i][j]=sum的表示顶点i到顶点j的最短路径为sum;
2.然后对s矩阵进行初始化,距离矩阵s中顶点a[i][j]的值为顶点i到顶点j的直接权值;
3.然后对s矩阵循环进行n次更新,在第k次更新时,如果s矩阵的a[i][j] > a[i][k]+a[k][j],那么更新a[i][j]=a[i][k]+a[k][j]。循环更新完毕,则算法完成,所有的节点到其余各节点之间的最短距离已经找到了。
相比于dijkstra 算法,floyd算法支持带有负权边的图,但是不能解决带有“负权回路”(或者叫“负权环”)的图,实际上如果一个图中带有“负权回路”那么这个图则没有最短路径。
floyd算法的时间复杂度同样是时间复杂度o(n³),空间复杂度是o(n²),代码非常简单,但是思想相却是非常的巧妙,将所有的可能都枚举出来一一对比,取最小值,这样最终会得到最小值。
该案例对应着下面实现代码中的案例:
首先初始化距离矩阵s如下:
然后就是三层嵌套循环,开始第一轮大循环,即当k=0,循环遍历s矩阵,判断是否小于 shortestpath[i][j],即所有的路径都经过a点中转,如果经过a中转后的路径shortestpath[i][0] + shortestpath[k][0]< shortestpath[i][j],自然更新路径:shortestpath[i][j]= shortestpath[i][0] + shortestpath[k][0]。一轮大循环之后的数组如下:
然后经过一共经过n次的大循环,表示经过所有的顶点,最终取得的矩阵如下:
这里的实现能够构造一个基于邻接矩阵实现无向加权图的类,并且提供深度优先遍历和广度优先遍历的方法,提供获取边集数组的方法,提供prim和kruskal两种求最小生成树的方法,提供dijkstra和floyd两种求最短路径的方法。
/*** 无向加权图邻接矩阵实现* {@link matrixdijkstraandfloyd#matrixdijkstraandfloyd(object[], edge[])} 构建无向加权图* {@link matrixdijkstraandfloyd#dfs()} 深度优先遍历无向加权图* {@link matrixdijkstraandfloyd#bfs()} 广度优先遍历无向加权图* {@link matrixdijkstraandfloyd#tostring()} 宾语从句的引导词输出无向加权图* {@link matrixdijkstraandfloyd#prim()} prim算法实现最小生成树* {@link matrixdijkstraandfloyd#kruskal()} kruskal算法实现最小生成树* {@link matrixdijkstraandfloyd#kruskalandprim()} kruskal算法结合prim算法实现最小生成树* {@link matrixdijkstraandfloyd#getedges()} 获取边集数组* {@link matrixdijkstraandfloyd#dijkstra(int)} ()} dijkstra算法获取指定顶点到所有顶点的最短路径* {@link matrixdijkstraandfloyd#dijkstra(int, int)} dijkstra算法获取指定顶点到指定顶点的最短路径* {@link matrixdijkstraandfloyd#floyd()} floyd获取所有顶点到所有顶点的最短路径** @author lx*/public class matrixdijkstraandfloyd<e> {/*** 顶点数组*/private object[] vertexs;/*** 邻接矩阵*/private int[][] matrix;/****/private edge<e>[] edges;/*** 由于是加权图,这里设置一个边的权值上限,任何边的最大权值不能大于等于该值,在实际应用中,该值应该根据实际情况确定*/private static final int no_edge = 99;/*** 边对象,具有权值,在构建加权无向图时使用*/private static class edge<e> {private e from;private e to;private int weight;public edge(e from, e to, int weight) {this.from = from;this.to = to;this.weight = weight;}@overridepublic string tostring() {return "edge{" +"from=" + from +", to=" + to +", weight=" + weight +'}';}}/*** 创建无向加权图** @param vertexs 顶点数组* @param edges 边对象数组*/public matrixdijkstraandfloyd(object[] vertexs, edge<e>[] edges) {//初始化边数组this.edges = edges;// 初始化顶点数组,并添加顶点this.vertexs = arrays.copyof(vertexs, vertexs.length);// 初始化边矩阵,并预先填充边信息this.matrix = new int[vertexs.length][vertexs.length];for (int i = 0; i < vertexs.length; i++) {for (int j = 0; j < vertexs.length; j++) {if (i == j) {this.matrix[i][j] = 0;} el {this.matrix[i][j] = no_edge;}}}for (edge<e> edge : edges) {// 读取一条边的起始顶点和结束顶点索引值int p1 = getposition(edge.from);int p2 = getposition(edge.to);//对称的两个点位都置为edge.weight,无向图可以看作相互可达的有向图this.matrix[p1][p2] = edge.weight;this.matrix[p2][p1] = edge.weight;}}/*** 获取某条边的某个顶点所在顶点数组的索引位置** @param e 顶点的值* @return 所在顶点数组的索引位置, 或者-1 - 表示不存在*/private int getposition(e e) {for (int i = 0; i < vertexs.length; i++) {if (vertexs[i] == e) {return i;}}return -1;}/*** 深度优先搜索遍历图,类似于树的前序遍历,*/public void dfs() {//新建顶点访问标记数组,对应每个索引对应相同索引的顶点数组中的顶点boolean[] visited = new boolean[vertexs.length];//初始化所有顶点都没有被访问for (int i = 0; i < vertexs.length; i++) {visited[i] = fal;}system.out.println("dfs: ");system.out.print("\t");for (int i = 0; i < vertexs.length; i++) {if (!visited[i]) {dfs(i, visited);}}system.out.println();}/*** 深度优先搜索遍历图的递归实现,类似于树的先序遍历* 因此模仿树的先序遍历,同样借用栈结构,这里使用的是方法的递归,隐式的借用栈** @param i 顶点索引* @param visited 访问标志数组*/private void dfs(int i, boolean[] visited) {visited[i] = true;system.out.print(vertexs[i] + " ");// 遍历该顶点的所有邻接点。若该邻接点是没有访问过,那么继续递归遍历领接点for (int w = firstvertex(i); w >= 0; w = nextvertex(i, w)) {if (!visited[w]) {dfs(w, visited);}}}/*** 广度优先搜索图,类似于树的层序遍历* 因此模仿树的层序遍历,同样借用队列结构*/public void bfs() {// 辅组队列queue<integer> indexlinkedlist = new linkedlist<>();//新建顶点访问标记数组,对应每个索引对应相同索引的顶点数组中的顶点boolean[] visited = new boolean[vertexs.length];for (int i = 0; i < vertexs.length; i++) {visited[i] = fal;}system.out.println("bfs: ");system.out.print("\t");for (int i = 0; i < vertexs.length; i++) {if (!visited[i]) {visited[i] = true;system.out.print(vertexs[i] + " ");indexlinkedlist.add(i);}if (!indexlinkedlist.impty()) {//j索引出队列integer j = indexlinkedlist.poll();//继续访问j的邻接点for (int k = firstvertex(j); k >= 0; k = nextvertex(j, k)) {if (!visited[k]) {visited[k] = true;system.out.print(vertexs[k] + " ");//继续入队列indexlinkedlist.add(k);}}}}system.out.println();}/*** 返回顶点v的第一个邻接顶点的索引,失败则返回-1** @param v 顶点v在数组中的索引* @return 返回顶点v的第一个邻接顶点的索引,失败则返回-1*/private int firstvertex(int v) {//如果索引超出范围,则返回-1if (v < 0 || v > (vertexs.length - 1)) {return -1;}/*根据邻接矩阵的规律:顶点索引v对应着边二维矩阵的matrix[v][i]一行记录* 从i=0开始*/for (int i = 0; i < vertexs.length; i++) {if (matrix[v][i] != 0 && matrix[v][i] != no_edge) {return i;}}return -1;}/*** 返回顶点v相对于w的下一个邻接顶点的索引,失败则返回-1** @param v 顶点索引* @param w 第一个邻接点索引* @return 返回顶点v相对于w的下一个邻接顶点的索引,失败则返回-1*/private int nextvertex(int v, int w) {//如果索引超出范围,则返回-1if (v < 0 || v > (vertexs.length - 1) || w < 0 || w > (vertexs.length - 1)) {return -1;}/*根据邻接矩阵的规律:顶点索引v对应着边二维矩阵的matrix[v][i]一行记录* 由于邻接点w的索引已经获取了,所以从i=w+1开始寻找*/for (int i = w + 1; i < vertexs.length; i++) {if (matrix[v][i] != 0 && matrix[v][i] != no_edge) {return i;}}return -1;}/*** 输出图** @return 输出图字符串*/@overridepublic string tostring() {stringbuilder stringbuilder = new stringbuilder();for (int i = 0; i < vertexs.length; i++) {for (int j = 0; j < vertexs.length; j++) {stringbuilder.append(matrix[i][j]).append("\t");}stringbuilder.append("\n");}return stringbuilder.tostring();}/*** prim算法求最小生成树*/public void prim() {system.out.println("prim: ");//对应节点应该被连接的前驱节点,用来输出//默认为0,即前驱结点为第一个节点int[] mid = new int[matrix.length];//如果某顶点作为末端顶点被连接,对应位置应该为true//第一个顶点默认被连接boolean[] connected = new boolean[matrix.length];connected[0] = true;//存储未连接顶点到已连接顶点的最短距离(最小权)int[] dis = new int[matrix.length];//首先将矩阵第一行即其他顶点到0索引顶点的权值拷贝进去system.arraycopy(matrix[0], 0, dis, 0, matrix.length);//存储路径长度int sum = 0;//最小权值int min;/*默认第一个顶点已经找到了,因此最多还要需要大循环n-1次*/for (int k = 1; k < matrix.length; k++) {min = no_edge;//最小权值的顶点的索引int minindex = 0;/*寻找权值最小的且未被连接的顶点索引*/for (int i = 1; i < matrix.length; i++) {//排除已连接的顶点,排除权值等于0的值,这里权值等于0表示已生成的最小生成树的顶点都未能与该顶点连接if (!connected[i] && dis[i] != 0 && dis[i] < min) {min = dis[i];minindex = i;}}//如果没找到,那么该图可能不是连通图,直接返回了,此时最小生成树没啥意义if (minindex == 0) {return;}//权值和增加sum += min;//该新连接顶点对应的索引值变成true,表示已被连接,后续判断时跳过该顶点connected[minindex] = true;//输出对应的前驱顶点到该最小顶点的权值system.out.println("\t" + vertexs[mid[minindex]] + " ---> " + vertexs[minindex] + " 权值:" + min);/*在新顶点minindex加入之前的其他所有顶点到连接顶点最小的权值已经计算过了因此只需要更新其他未连接顶点到新连接顶点minindex是否还有更短的权值,有的话就更新找到距离已连接的顶点权最小的顶点*/for (int i = 1; i < matrix.length; i++) {//如果该顶点未连接if (!connected[i]) {/*如果新顶点到未连接顶点i的权值不为0,并且比原始顶点到未连接顶点i的权值还要小,那么更新对应位置的最小权值*/if (matrix[minindex][i] != 0 && dis[i] > matrix[minindex][i]) {//更新最小权值dis[i] = matrix[minindex][i];//更新前驱节点索引为新加入节点索引mid[i] = minindex;}}}}system.out.println("\t" + "sum: " + sum);}/*** kruskal算法求最小生成树传统实现,要求知道边集数组,和顶点数组*/public void kruskal() {system.out.println("kruskal: ");//由于创建图的时候保存了边集数组,这里直接使用就行了//edge[] edges = getedges();//this.edges=edges;//对边集数组进行排序arrays.sort(this.edges, comparator.comparingint(o -> o.weight));// 用于保存已有最小生成树中每个顶点在该最小树中的最终终点的索引int[] vends = new int[this.edges.length];//能够知道终点索引范围是[0,this.edges.length-1],因此填充edges.length表示没有终点arrays.fill(vends, this.edges.length);int sum = 0;for (edge<e> edge : this.edges) {// 获取第i条边的起点索引fromint from = getposition(edge.from);// 获取第i条边的终点索引toint to = getposition(edge.to);// 获取顶点from在"已有的最小生成树"中的终点int m = getendindex(vends, from);// 获取顶点to在"已有的最小生成树"中的终点int n = getendindex(vends, to);// 如果m!=n,意味着没有形成环路,则可以添加,否则直接跳过,进行下一条边的判断if (m != n) {//添加设置原始终点索引m在已有的最小生成树中的终点为nvends[m] = n;system.out.println("\t" + vertexs[from] + " ---> " + vertexs[to] + " 权值:" + edge.weight);sum += edge.weight;}}system.out.println("\t" + "sum: " + sum);//system.out.println(arrays.tostring(this.edges));}/*** 获取顶点索引i的终点如果没有终点则返回顶点索引本身** @param vends 顶点在最小生成树中的终点* @param i 顶点索引* @return 顶点索引i的终点如果没有终点则返回顶点索引本身*/private int getendindex(int[] vends, int i) {//这里使用循环查找的逻辑,寻找的是最终的终点while (vends[i] != this.edges.length) {i = vends[i];}return i;}/*** 如果没有现成的边集数组,那么根据邻接矩阵结构获取图中的边集数组** @return 图的边集数组*/private edge[] getedges() {list<edge> edges = new arraylist<>();/*遍历矩阵数组 只需要遍历一半就行了*/for (int i = 0; i < vertexs.length; i++) {for (int j = i + 1; j < vertexs.length; j++) {//如果存在边if (matrix[i][j] != no_edge && matrix[i][j] != 0) {edges.add(new edge<>(vertexs[i], vertexs[j], matrix[i][j]));//edges[index++] = new edge(vertexs[i], vertexs[j], matrix[i][j]);}}}return edges.toarray(new edge[0]);}/*** kruskal结合prim算法.不需要知道边集,只需要矩阵数组,和顶点数组* 同样是求最小权值的边,但是有一个默认起点顶点,该起点可以是要求[0,顶点数量-1]之间的任意值,同时查找最小权的边。* 可能会有bug,目前未发现*/public void kruskalandprim() {system.out.println("kruskalandprim: ")极限挑战第三季第十期;//已经找到的边携带的顶点对应的索引将变为true,其余未找到边对应的顶点将是falboolean[] connected = new boolean[matrix.length];//这里选择第一个顶点为起点,表示以该顶点开始寻找包含该顶点的最小边connected[0] = true;int sum = 0, n1 = 0, n2 = 0;//最小权值int min;while (true) {min = no_edge;/*找出所有带有已找到顶点的边中,最小权值的边,只需要寻找对称矩阵的一半即可*///第一维for (int i = 0; i < matrix.length; i++) {//第二维for (int j = i + 1; j < matrix.length; j++) {//排除等于0的,排除两个顶点都找到了的,这里实际上已经隐含了排除环的逻辑,如果某条边的两个顶点都找到了,那么如果算上该条边,肯定会形成环//寻找剩下的最小的权值的边if (matrix[i][j] != 0 && connected[i] != connected[j] && matrix[i][j] < min) {min = matrix[i][j];n1 = i;n2 = j;}}}//如果没找到最小权值,该图可能不是连通图,或者已经寻找完毕,直接返回if (min == no_edge) {system.out.println("\t" + "sum:" + sum);return;}//已经找到的边对应的两个顶点都置为trueconnected[n1] = true;connected[n2] = true;//输出找到的边和最小权值system.out.println("\t" + vertexs[n1] + " ---> " + vertexs[n2] + " 权值:" + min);sum += min;}}/*** dijkstra算法求最短路径。** @param start 起始顶点索引。即计算"顶点vs"到其它顶点的最短路径。*/public void dijkstra(int start) {checkindex(start);int[] shortestpathance = getshortestdistance(start, vertexs.length);// 打印dijkstra最短路径的结果system.out.println("dijkstra(" + vertexs[start] + "):");for (int i = 0; i < vertexs.length; i++) {system.out.println("\t(" + vertexs[start] + " ---> " + vertexs[i] + ")最短路径:" + shortestpathance[i]);}}/*** dijkstra算法求最短路径** @param start 起始点* @param end 终点,如果end=vertexs.length说明是遍历查找所有的最短路径* @return 起始顶点到其他点或者指定点的最短权值*/private int[] getshortestdistance(int start, int end) {/*1、该数组存放起始顶点到其他点的权值*/int[] shortestpathance = new int[vertexs.length];//初始化数据//首先设置起始点到顶点i到的最短路径为起始点到顶点i的权。system.arraycopy(matrix[start], 0, shortestpathance, 0, matrix.length);/*2、标志位数组.某个位置如果为true表示对应位置的顶点到起始顶点的最短路径已成功获取。*/boolean[] shortest = new boolean[vertexs.length];//首先设置起始点到自己的的路径已经找到了,为0shortest[start] = true;/*3、最多遍历vertexs.length-1次;每次找出起始点到一个顶点的最短路径。*/int k;int min;for (int i = 1; i < vertexs.length; i++) {k = 0;// 寻找当前最小的路径;min = no_edge;for (int j = 0; j < vertexs.length; j++) {//排除已经找到的最短路径之后,找到离start最近的顶点(k)。if (!shortest[j] && shortestpathance[j] < min) {min = shortestpathance[j];k = j;}}//先设置起始点到新顶点k的最短路径已经找到shortest[k] = true;if (end != vertexs.length && k == end) {break;}//更新未获取最短路径的顶点的最短路径,因为其他已找到的顶点的最短路径已经找到了,这里指需要更新新加入的已找到的可达顶点的路径.for (int j = 0; j < vertexs.length; j++) {int tmp = matrix[k][j];//排除已经找到的最短路径,排除未连接的路径,排除等于0的路径(连接自己)之后//找到离start最如果新的最短路径比以前的最短路径还要短,则更新最短路径。if (!shortest[j] && tmp != no_edge && tmp != 0 && ((tmp = min + tmp) < shortestpathance[j])) {shortestpathance[j] = tmp;}}}return shortestpathance;}/*** 索引检查** @param index 多个索引*/private void checkindex(int... index) {for (int i : index) {if (i < 0 || i >= vertexs.length) {throw new arrayindexoutofboundxception("索引越界:" + i);}}}/*** dijkstra算法求最短路径。** @param start 起始顶点索引。* @param end 结束点索引*/public void dijkstra(int start, int end) {checkindex(start, end);int[] shortestpathance = getshortestdistance(start, end);// 打印dijkstra最短路径的结果system.out.println("dijkstra(" + vertexs[start] + " ---> " + vertexs[end] + ")最短路径:" + shortestpathance[end]);}/*** floyd算法获取所有顶点到所有顶点的最短路径,代码很简单,思想很巧妙*/public void floyd() {//路径矩阵(两顶点最短路径,即最小权值)int[][] shortestpath = new int[matrix.length][matrix.length];/*初始化数据*/for (int i = 0; i < matrix.length; i++) {system.arraycopy(matrix[i], 0, shortestpath[i], 0, vertexs.length);}// 计算最短路径for (int k = 0; k < matrix.length; k++) {for (int i = 0; i < matrix.length; i++) {for (int j = 0; j < matrix.length; j++) {//要求经过下标k顶点的两个路径都不能等于no_edge,否则就是没有路径,no_edge应该选取的足够的大,否则可能出错int tmp = (shortestpath[i][k] == no_edge || shortestpath[k][j] == no_edge) ? no_edge : (shortestpath[i][k] + shortestpath[k][j]);// 如果经过下标为k顶点路径比原两点间路径更短,则更新shortestpath[i][j]if (shortestpath[i][j] > tmp) {// i到j最短路径对应的值设为经过k的更小的一个shortestpath[i][j] = tmp;}}}}/*输出路径矩阵*/system.out.println("floyd: ");for (int i = 0; i < matrix.length; i++) {for (int j = 0; j < matrix.length; j++) {system.out.print("\t" + shortestpath[i][j]);}system.out.println();}}public static void main(string[] args) {//顶点数组character[] vexs = {'a', 'b', 'c', 'd', 'e', 'f', 'g'};//边数组,加权值edge[] edges = {new edge<>('a', 'c', 8),new edge<>('d', 'a', 2),new edge<>('a', 'f', 3),new edge<>('b', 'c', 4),new edge<>('c', 'd', 5),new edge<>('e', 'g', 6),new edge<>('e', 'b', 7),new edge<>('d', 'b', 9),new edge<>('f', 'g', 9)};//构建图matrixdijkstraandfloyd<character> matrixdijkstraandfloyd = new matrixdijkstraandfloyd<character>(vexs, edges);//输出图system.out.println(matrixdijkstraandfloyd);//深度优先遍历matrixdijkstraandfloyd.dfs();//广度优先遍历matrixdijkstraandfloyd.bfs();//prim算法输出最小生成树matrixdijkstraandfloyd.prim();//kruskal算法输出最小生成树matrixdijkstraandfloyd.kruskal();//kruskal算法结合prim算法输出最小生成树,可能会有bug,目前未发现matrixdijkstraandfloyd.kruskalandprim();// dijkstra算法获取某个索引的顶点到其它各个顶点的最短距离// 这里参数是索引,也可以是一个顶点,需要稍微修改代码获取顶点的索引,比较简单这里就不做了matrixdijkstraandfloyd.dijkstra(0);// dijkstra算法获取一个顶点到另一个顶点的最短距离matrixdijkstraandfloyd.dijkstra(2, 0);// floyd算法获取所有顶点到所有顶点的最短路径matrixdijkstraandfloyd.floyd();}}
这里的实现能够构造一个基于邻接表实现无向加权图的类;并且提供深度优先遍历和广度优先遍历的方法,提供获取边集数组的方法,提供prim和kruskal两种求最小生成树的方法,提供dijkstra和floyd两种求最短路径的方法。
/*** 无向加权图邻接表实现* {@link listdijkstraandfloyd#listdijkstraandfloyd(object[], edge[])} 构建无向加权图* {@link listdijkstraandfloyd#dfs()} 深度优先遍历无向加权图* {@link listdijkstraandfloyd#bfs()} 广度优先遍历无向加权图* {@link listdijkstraandfloyd#tostring()} 输出无向加权图* {@link listdijkstraandfloyd#prim()} prim算法实现最小生成树* {@link listdijkstraandfloyd#kruskal()} kruskal算法实现最小生成树* {@link listdijkstraandfloyd#getedges()} 获取边集数组* {@link listdijkstraandfloyd#dijkstra(int)} ()} 获取指定顶点到所有顶点的最短路径* {@link listdijkstraandfloyd#dijkstra(int, int)} 获取指定顶点到指定顶点的最短路径* {@link listdijkstraandfloyd#floyd()} floyd获取所有顶点到所有顶点的最短路径** @author lx*/public class listdijkstraandfloyd<e> {/*** 顶点类** @param <e>*/private class node<e> {/*** 顶点信息*/e data;/*** 指向第一条依附该顶点的边*/lnode firstlnode;public node(e data, lnode firstlnode) {this.data = data;this.firstlnode = firstlnode;}}/*** 边表节点类*/private class lnode {/*** 该边所指向的顶点的索引位置*/int vertex;/*** 该边的权值*/int weight;/*** 指向下一条边的指针*/lnode nextlnode;}/*** 边对象,具有权值,在构建加权无向图时使用*/private static class edge<e> {private e from;private e to;private int weight;public edge(e from, e to, int weight) {this.from = from;this.to = to;this.weight = weight;}@overridepublic string tostring() {return "edge{" +"from=" + from +", to=" + to +", weight=" + weight +'}';}}/*** 顶点数组*/private node<e>[] vertexs;/*** 边数组*/private edge<e>[] edges;/*** 由于是加权图,这里设置一个边的权值上限,任何边的最大权值不能大于等于该值,在实际应用中,该值应该根据实际情况确定*/private static final int no_edge = 99;/*** 创建无向加权图** @param vexs 顶点数组* @param edges 边二维数组*/public listdijkstraandfloyd(e[] vexs, edge<e>[] edges) {this.edges = edges;/*初始化顶点数组,并添加顶点*/vertexs = new node[vexs.length];for (int i = 0; i < vertexs.length; i++) {vertexs[i] = new node<>(vexs[i], null);}/*初始化边表,并添加边节点到边表尾部,即采用尾插法*/for (edge<e> edge : edges) {// 读取一条边的起始顶点和结束顶点索引值int p1 = getposition(edge.from);int p2 = getposition(edge.to);int weight = edge.weight;/*这里需要相互添加边节点,无向图可以看作相互可达的有向图*/// 初始化lnode1边节点lnode lnode1 = new lnode();lnode1.vertex = p2;lnode1.weight = weight;// 将lnode链接到"p1所在链表的末尾"if (vertexs[p1].firstlnode == null) {vertexs[p1].firstlnode = lnode1;} el {linklast(vertexs[p1].firstlnode, lnode1);}// 初始化lnode2边节点lnode lnode2 = new lnode();lnode2.vertex = p1;lnode2.weight = weight;// 将lnode2链接到"p2所在链表的末尾"if (vertexs[p2].firstlnode == null) {vertexs[p2].firstlnode = lnode2;} el {linklast(vertexs[p2].firstlnode, lnode2);}}}/*** 获取某条边的某个顶点所在顶点数组的索引位置** @param e 顶点的值* @return 所在顶点数组的索引位置, 或者-1 - 表示不存在*/private int getposition(e e) {for (int i = 0; i < vertexs.length; i++) {if (vertexs[i].data == e) {return i;}}return -1;}/*** 将lnode节点链接到边表的最后,采用尾插法** @param first 边表头结点* @param node 将要添加的节点*/private void linklast(lnode first, lnode node) {while (true) {if (first.vertex == node.vertex) {return;}if (first.nextlnode == null) {break;}first = first.nextlnode;}first.nextlnode = node;}@overridepublic string tostring() {stringbuilder stringbuilder = new stringbuilder();for (int i = 0; i < vertexs.length; i++) {stringbuilder.append(i).append("(").append(vertexs[i].data).append("): ");lnode node = vertexs[i].firstlnode;while (node != null) {stringbuilder.append(node.vertex).append("(").append(vertexs[node.vertex].data).append("-").append(node.weight).append(")");node = node.nextlnode;if (node != null) {stringbuilder.append("->");} el {break;}}stringbuilder.ap四年级语文教学总结pend("\n");}return stringbuilder.tostring();}/*** 深度优先搜索遍历图的递归实现,类似于树的先序遍历* 因此模仿树的先序遍历,同样借用栈结构,这里使用的是方法的递归,隐式的借用栈** @param i 顶点索引* @param visited 访问标志数组*/private void dfs(int i, boolean[] visited) {//索引索引标记为true ,表示已经访问了visited[i] = true;system.out.print(vertexs[i].data + " ");//获取该顶点的边表头结点lnode node = vertexs[i].firstlnode;//循环遍历该顶点的邻接点,采用同样的方式递归搜索while (node != null) {if (!visited[node.vertex]) {dfs(node.vertex, visited);}node = node.nextlnode;}}/*** 深度优先搜索遍历图,类似于树的前序遍历,*/public void dfs() {//新建顶点访问标记数组,对应每个索引对应相同索引的顶点数组中的顶点boolean[] visited = new boolean[vertexs.length];//初始化所有顶点都没有被访问for (int i = 0; i < vertexs.length; i++) {visited[i] = fal;}system.out.println("dfs: ");system.out.print("\t");/*循环搜索*/for (int i = 0; i < vertexs.length; i++) {//如果对应索引的顶点的访问标记为fal,则搜索该顶点if (!visited[i]) {dfs(i, visited);}}/*走到这一步,说明顶点访问标记数组全部为true,说明全部都访问到了,深度搜索结束*/system.out.println();}/*** 广度优先搜索图,类似于树的层序遍历* 因此模仿树的层序遍历,同样借用队列结构*/public void bfs() {// 辅组队列queue<integer> indexlinkedlist = new linkedlist<>();//新建顶点访问标记数组,对应每个索引对应相同索引的顶点数组中的顶点boolean[] visited = new boolean[vertexs.length];//初始化所有顶点都没有被访问for (int i = 0; i < vertexs.length; i++) {visited[i] = fal;}system.out.println("bfs: ");system.out.print("\t");for (int i = 0; i < vertexs.length; i++) {//如果访问方剂为fal,则设置为true,表示已经访问,然后开始访问if (!visited[i]) {visited[i] = true;system.out.print(vertexs[i].data + " ");indexlinkedlist.add(i);}//判断队列是否有值,有就开始遍历if (!indexlinkedlist.impty()) {//出队列integer j = indexlinkedlist.poll();lnode node = vertexs[j].firstlnode;while (node != null) {int k = node.vertex;if (!visited[k]) {visited[k] = true;system.out.print(vertexs[k].data + " ");//继续入队列indexlinkedlist.add(k);}node = node.nextlnode;}}}system.out.println();}/*** prim算法求最小生成树*/public void prim() {system.out.println("prim: ");//对应节点应该被连接的前驱节点,用来输出//默认为0,即前驱结点为第一个节点int[] mid = new int[vertexs.length];int start = 0;int min, tmp, sum = 0;int num = vertexs.length;//顶点间边的权值//存储未连接顶点到已连接顶点的最短距离(最小权)int[] dis = new int[num];// 初始化"顶点的权值数组",// 将每个顶点的权值初始化为"第start个顶点"到"该顶点"的权值。//首先将其他顶点到0索引顶点的权值存储进去for (int i = 0; i < num; i++) {dis[i] = getweight(start, i);}//如果某顶点作为末端顶点被连接,对应位置应该为true//第一个顶点默认被连接boolean[] connected = new boolean[vertexs.length];connected[0] = true;/*默认第一个顶点已经找到了,因此最多还要需要大循环n-1次*/for (int k = 1; k < num; k++) {min = no_edge;//最小权值的顶点的索引int minindex = 0;// 在未被加入到最小生成树的顶点中,找出权值最小的顶点。for (int i = 1; i < vertexs.length; i++) {//排除已连接的顶点,排除权值等于0的值,因为这里默认顶点指向自己的权值为0if (!connected[i] && dis[i] != 0 && dis[i] < min) {min = dis[i];minindex = i;}}//如果没找到,那么该图可能不是连通图,直接返回了,此时最小生成树没啥意义if (minindex == 0) {return;}//权值和增加sum += min;//该新连接顶点对应的索引值变成true,表示已被连接,后续判断时跳过该顶点connected[minindex] = true;//输出对应的前驱顶点到该最小顶点的权值system.out.println("\t" + vertexs[mid[minindex]].data + " ---> " + vertexs[minindex].data + " 权值:" + min);/*在新顶点minindex加入之前的其他所有顶点到连接顶点最小的权值已经计算过了因此只需要更新其他顶点到新连接顶点mininde基督教圣诞节节目x是否还有更短的权值,有的话就更新找到距离已连接的顶点权最小的顶点*/for (int i = 1; i < num; i++) {//如果该顶点未连接if (!connected[i]) {// 获取minindex顶点到未连接顶点i的权值tmp = getweight(minindex, i);/*如果新顶点到未连接顶点i的权值不为0,并且比原始顶点到未连接顶点i的权值还要小,那么更新对应位置的最小权值*/if (tmp != 0 && dis[i] > tmp) {dis[i] = tmp;//更新前驱节点索引为新加入节点索引mid[i] = minindex;}}}}system.out.println("\t" + "sum: " + sum);}/*** 尝试获取边起点start到边终点end的边的权值,当然可能获取不到** @param start 边起点* @param end 边终点* @return 返回权值; 如果起点和终点相同则返回0;如果边起点和边终点之间并没有边, 则返回no_edge*/private int getweight(int start, int end) {//如果start=end,则返回0if (start == end) {return 0;}//获取该顶点的边表的第一个值lnode node = vertexs[start].firstlnode;//循环查找边表,看能否找到对应的索引=end,找不到就返回no_edge,表示两个顶点未连接。while (node != null) {if (end == node.vertex) {return node.weight;}node = node.nextlnode;}return no_edge;}/*** kruskal算法求最小生成树,可以说邻接矩阵和邻接链表的实现方式是完全一致的*/public void kruskal() {system.out.println("kruskal: ");//由于创建图的时候保存了边集数组,这里直接使用就行了//edge[] edges = getedges();//this.edges=edges;//对边集数组进行排序arrays.sort(this.edges, comparator.comparingint(o -> o.weight));// 用于保存已有最小生成树中每个顶点在该最小树中的最终终点的索引int[] vends = new int[this.edges.length];//能够知道终点索引范围是[0,this.edges.length-1],因此填充edges.length表示没有终点arrays.fill(vends, this.edges.length);int sum = 0;for (edge<e> edge : this.edges) {// 获取第i条边的起点索引fromint from = getposition(edge.from);// 获取第i条边的终点索引toint to = getposition(edge.to);// 获取顶点from在"已有的最小生成树"中的终点int m = getendindex(vends, from);// 获取顶点to在"已有的最小生成树"中的终点int n = getendindex(vends, to);// 如果m!=n,意味着没有形成环路,则可以添加,否则直接跳过,进行下一条边的判断if (m != n) {//添加设置原始终点索引m在已有的最小生成树中的终点为nvends[m] = n;system.out.println("\t" + vertexs[from].data + " ---> " + vertexs[to].data + " 权值:" + edge.weight);sum += edge.weight;}}system.out.println("\t" + "sum: " + sum);//system.out.println(arrays.tostring(this.edges));}/*** 获取顶点索引i的终点如果没有终点则返回顶点索引本身** @param vends 顶点在最小生成树中的终点* @param i 顶点索引* @return 顶点索引i的终点如果没有终点则返回顶点索引本身*/private int getendindex(int[] vends, int i) {//这里使用循环查找的逻辑,寻找的是最终的终点while (vends[i] != this.edges.length) {i = vends[i];}return i;}/*** 如果没有现成的边集数组,那么根据邻接表结构获取图中的边集数组** @return 图的边集数组*/private edge[] getedges() {list<edge> edges = new arraylist<>();//遍历顶点数组for (int i = 0; i < vertexs.length; i++) {lnode node = vertexs[i].firstlnode;while (node != null) {//只需添加起点索引小于终点索引的边就行了if (node.vertex > i) {edges.add(new edge<>(vertexs[i].data, vertexs[node.vertex].data, node.weight));}node = node.nextlnode;}}return edges.toarray(new edge[0]);}/*** dijkstra算法求最短路径。** @param start 起始顶点索引。即计算"顶点vs"到其它顶点的最短路径。*/public void dijkstra(int start) {checkindex(start);int[] distance = getshortestdistance(start, vertexs.length);// 打印di鸡肝怎么做好吃jkstra最短路径的结果system.out.println("dijkstra(" + vertexs[start].data + "):");for (int i = 0; i < vertexs.length; i++) {system.out.println("\t(" + vertexs[start].data + " ---> " + vertexs[i].data + ")最短路径:" + distance[i]);}}/*** dijkstra算法求最短路径。** @param start 起始顶点索引。* @param end 结束点索引*/public void dijkstra(int start, int end) {checkindex(start, end);int[] shortestpathance = getshortestdistance(start, end);// 打印dijkstra最短路径的结果system.out.println("dijkstra(" + vertexs[start].data + " ---> " + vertexs[end].data + ")最短路径:" + shortestpathance[end]);}/*** dijkstra算法求最短路径** @param start 起始点* @param end 终点,如果end=vertexs.length说明是遍历查找所有的最短路径* @return 起始顶点到其他点或者指定点的最短权值*/private int[] getshortestdistance(int start, int end) {/*1、该数组存放起始顶点到其他点的权值*/int[] distance = new int[vertexs.length];//初始化数据for (int i = 0; i < vertexs.length; i++) {//首先设置起始点到顶点i到的最短路径为起始点到顶点i的权。distance[i] = getweight(start, i);}/*2、标志位数组.某个位置表示true表示,对应位置的顶点到起始顶点的最短路径已成功获取。*/boolean[] shortest = new boolean[vertexs.length];//首先设置起始点到自己的的路径已经找到了,为0shortest[start] = true;/*3、最多遍历vertexs.length-1次;每次找出起始点到一个顶点的最短路径。*/int k;int min;for (int i = 1; i < vertexs.length; i++) {k = 0;// 寻找当前最小的路径;min = no_edge;for (int j = 0; j < vertexs.length; j++) {//排除已经找到的最短路径之后,找到离start最近的顶点(k)。if (!shortest[j] && distance[j] < min) {min = distance[j];k = j;}}//先设置起始点到新顶点k的最短路径已经找到shortest[k] = true;if (end != vertexs.length && k == end) {break;}//更新未获取最短路径的顶点的最短路径,因为其他已找到的顶点的最短路径已经找到了,这里指需要更新新加入的已找到的可达顶点的路径.for (int j = 0; j < vertexs.length; j++) {int tmp = getweight(k, j);//排除已经找到的最短路径,排除未连接的路径,排除等于0的路径(连接自己)之后//找到离start最如果新的最短路径比以前的最短路径还要短,则更新最短路径。if (!shortest[j] && tmp != no_edge && tmp != 0 && ((tmp = min + tmp) < distance[j])) {distance[j] = tmp;}}}return distance;}/*** 索引检查** @param index 多个索引*/private void checkindex(int... index) {for (int i : index) {if (i < 0 || i >= vertexs.length) {throw new arrayindexoutofboundxception("索引越界:" + i);}}}/*** floyd算法获取所有顶点到所有顶点的最短路径,与邻接矩阵的实现基本一致*/public void floyd() {//路径矩阵(两顶点最短路径,即最小权值)int[][] shortestpath = new int[vertexs.length][vertexs.length];/*初始化数据*/for (int i = 0; i < vertexs.length; i++) {for (int j = 0; j < vertexs.length; j++) {//获取两点的直接权值//如果是频繁调用该方法,因此可以创建一个属于对象的权值矩阵用来保存权值,这里为了简单没做shortestpath[i][j] = getweight(i, j);}}// 计算最短路径for (int k = 0; k < vertexs.length; k++) {for (int i = 0; i < vertexs.length; i++) {for (int j = 0; j < vertexs.length; j++) {//要求经过下标k顶点的两个路径都不能等于no_edge,否则就是没有路径,no_edge应该选取的足够的大,否则可能出错int tmp = (shortestpath[i][k] == no_edge || shortestpath[k][j] == no_edge) ? no_edge : (shortestpath[i][k] + shortestpath[k][j]);// 如果经过下标为k顶点路径比原两点间路径更短,则更新shortestpath[i][j]if (shortestpath[i][j] > tmp) {// i到j最短路径对应的值设为经过k的更小的一个shortestpath[i][j] = tmp;}}}}/*输出路径矩阵*/system.out.println("floyd: ");for (int i = 0; i < vertexs.length; i++) {for (int j = 0; j < vertexs.length; j++) {system.out.print("\t" + shortestpath[i][j]);}system.out.println();}}public static void main(string[] args) {//顶点数组character[] vexs = {'a', 'b', 'c', 'd', 'e', 'f', 'g'};//边数组,加权值edge[] edges = {new edge<>('a', 'c', 8),new edge<>('d', 'a', 2),new edge<>('a', 'f', 3),new edge<>('b', 'c', 4),new edge<>('c', 'd', 5),new edge<>('e', 'g', 6),new edge<>('e', 'b', 7),new edge<>('d', 'b', 9),new edge<>('f', 'g', 9)};//构建图listdijkstraandfloyd<character> listdijkstraandfloyd = new listdijkstraandfloyd<character>(vexs, edges);//输出图system.out.println(listdijkstraandfloyd);//深度优先遍历//dfs://a c b e g f dlistdijkstraandfloyd.dfs();//广度优先遍历//bfs://a c d f b g elistdijkstraandfloyd.bfs();//prim算法求最小生成树listdijkstraandfloyd.prim();//kruskal算法求最小生成树listdijkstraandfloyd.kruskal();// dijkstra算法获取某个索引的顶点到其它各个顶点的最短距离// 这里参数是索引,也可以是一个顶点,需要稍微修改代码获取顶点的索引,比较简单这里就不做了listdijkstraandfloyd.dijkstra(0);// dijkstra算法获取一个顶点到另一个顶点的最短距离listdijkstraandfloyd.dijkstra(2, 0);// floyd算法获取所有顶点到所有顶点的最短路径listdijkstraandfloyd.floyd();}}
以上就是java利用dijkstra和floyd分别求取图的最短路径的详细内容,更多关于java求最短路径的资料请关注www.887551.com其它相关文章!
本文发布于:2023-04-04 18:49:57,感谢您对本站的认可!
本文链接:https://www.wtabcd.cn/fanwen/zuowen/ac149d090fcee02eb3e91a3fa9edefd9.html
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。
本文word下载地址:Java利用Dijkstra和Floyd分别求取图的最短路径.doc
本文 PDF 下载地址:Java利用Dijkstra和Floyd分别求取图的最短路径.pdf
留言与评论(共有 0 条评论) |