节点最短路径算法改进

2744557306

实现了一种基于邻接矩阵的图算法，具体是 Floyd-Warshall 算法的一个变种，用于计算每对节点之间的最短路径。下面对代码进行逐行分析并解释算法思想。

Floyd-Warshall 算法的基本思想Floyd-Warshall 算法是一种用于求解加权图中任意两点之间最短路径的动态规划算法。它的思想是通过中间节点逐步更新路径长度，最终得到所有节点对之间的最短路径。
2 Z( M( r4 M% t+ }/ h# M( g+ ~代码解析

[color=rgba(6, 8, 31, 0.88)]函数定义:

1. function a = dij2_m(a)

复制代码

- `a` 为输入的邻接矩阵，表示图中节点之间的边的权重。
- 输出结果也是更新后的邻接矩阵 `a`，其中的值代表任意两节点之间的最短路径距离。

0 d& o5 O. P3 B: K- p* A1 `$ T[color=rgba(6, 8, 31, 0.88)]将邻接矩阵变为对称矩阵[color=rgba(6, 8, 31, 0.88)]:

1. n = length(a);  + k2 h- v0 F- o1 ~! E+ c1 a3 U
   
2. for i = 2:n  4 {- O# Z7 s; i4 j- a6 f
   
3.     for j = 1:(i-1)  5 l+ r; N; y* w3 u, U( a& a
   
4.         a(i,j) = a(j,i);  , s/ f0 ]; c: F4 [
   
5.     end  
   ( [. ?% q1 W6 d: M; l9 @8 Z3 z8 V- j
6. end

复制代码

- 这部分代码确保输入的邻接矩阵是对称的，即 \(a(i, j) = a(j, i)\)。这是因为对于无向图，节点间的距离应该是相同的。
0 G! o( x+ F8 l! T0 a6 _
- \* {. l8 d6 l7 K[color=rgba(6, 8, 31, 0.88)]主程序[color=rgba(6, 8, 31, 0.88)]:

1. for k = 1:(n-1)  
   3 A5 f& o* F( @) k
2.     b = [1:(k-1),(k+1):n]; % 初始化剩余节点  8 x1 @) B% c% D4 \+ q
   
3.     kk = length(b);  
   * k9 L9 D3 N6 v; X# l  N2 q1 W
4.     a_id = k;              % 当前节点  
   / O& H! j+ q4 |& l  z; F
5.     b1 = [(k+1):n];       % 节点的一部分  
   - H% U& m: a+ t. @5 U
6.     kk1 = length(b1);

复制代码

- 这段代码是在外层循环中执行的。变量 `k` 表示当前考虑的中间节点。
- `b` 是未考虑的节点集合，初始为包含所有节点（除了当前节点 `k`）的数组。
- `b1` 是一组要更新的节点，将从 `k+1` 开始划分，以找到与节点 `k` 的路径。

3 L* s% P7 I" R5 ^4 r- z+ r# e[color=rgba(6, 8, 31, 0.88)]更新最短路径[color=rgba(6, 8, 31, 0.88)]:

1. while kk > 0  
   4 j4 d, r# [% _1 \- f# h! m: t
2.     for j = 1:kk1  7 w. v' ?' Z  w- g! X  `4 `# G
   
3.         te = a(k, a_id) + a(a_id, b1(j));  4 s8 t4 e3 q: p6 {
   
4.         if te < a(k, b1(j))  8 l) _2 ?7 W4 |4 [* b
   
5.             a(k, b1(j)) = te;  
   + j( k4 Q, f; {6 e! S4 Z- \& w
6.         end  5 v( K2 {( H\" ?8 p; U/ B, J
   
7.     end

复制代码

- 在 `while` 循环中，对于每个未访问节点，从当前节点 `k` 到所有在 `b1` 中节点的路径进行更新。
- 计算通过当前节点 `a_id` 访问 `b1` 中的每个节点的路径 `te`，如果这个新路径比已知的更短，就更新邻接矩阵。
) U8 e6 y$ ]' _+ x! Z4 U/ t
[color=rgba(6, 8, 31, 0.88)]选择下一个节点[color=rgba(6, 8, 31, 0.88)]:

1. miid = 1;  5 m/ s/ H  Z3 f/ {- C
   
2. for j = 2:kk  \" P) Z: s( {7 n/ a4 G% J
   
3.     if a(k, b(j)) < a(k, b(miid))  / [% }( I5 l/ r8 h; N( O4 L
   
4.         miid = j;  9 X9 E' W, E2 H+ o$ j1 J% t7 ^
   
5.     end  \" _( y% _& ^& \: z2 j; M* r. u
   
6. end  
   0 M3 i1 b% J$ t0 V
7. a_id = b(miid);          % 选取出最短路径的节点  
   6 t  n4 H6 |$ r
8. b = [b(1:(miid-1)), b((miid+1):kk)]; % 从该集合中移除选择的节点  $ O7 N$ m2 t' p5 x
   
9. kk = length(b);

复制代码

- 找到当前未访问节点中与 `k` 相连的最短距离节点，将其标记为 `a_id`，并更新集合 `b`，剔除已经访问的节点。

[color=rgba(6, 8, 31, 0.88)]更新所有节点间的距离[color=rgba(6, 8, 31, 0.88)]:

1. for j = (k+1):n  + w2 M. q' q$ N/ G( a9 @% G7 V
   
2.     a(j, k) = a(k, j);  
   ) B% a' @5 o# F  w* O
3. end

复制代码

- 更新 `a` 中关于节点 `k` 的距离，保证后续节点间的对称性。

总结整体来看，这段代码通过逐步选择中间节点，并更新已处理节点之间的最短路径距离，从而计算出无向图的任意两点之间的最短路径。其核心思想是动态规划和“逐步逼近”，和传统的 Floyd-Warshall 算法有相似之处，但实现上有所变化。
+ T  |! V' h6 v7 l! x4 ]
此算法的时间复杂度为 \(O(V^3)\)，其中 \(V\) 是图中节点的数量，适合用于小型图的最短路径计算。


6 E+ S$ t. J# j
/ F! z+ p+ V7 |/ i