shortest_path#
- scipy.sparse.csgraph.shortest_path(csgraph, method='auto', directed=True, return_predecessors=False, unweighted=False, overwrite=False, indices=None)#
在正向有向图或无向图上执行最短路径图搜索。
在版本 0.11.0 中添加。
- 参数:
- csgraph数组、矩阵或稀疏矩阵,2 维
表示输入图的 N x N 距离数组。
- method字符串 ['auto'|'FW'|'D'],可选
用于最短路径的算法。选项是
- ‘auto’ – (默认) 在 ‘FW’、‘D’、‘BF’ 或 ‘J’ 中选择最佳选项
基于输入数据。
- ‘FW’ – Floyd-Warshall 算法。
计算成本大约为
O[N^3]
。输入 csgraph 将转换为密集表示。- ‘D’ – 带有斐波那契堆的 Dijkstra 算法。
计算成本大约为
O[N(N*k + N*log(N))]
,其中k
是每个节点的平均连接边数。输入 csgraph 将转换为 csr 表示。- ‘BF’ – Bellman-Ford 算法。
当权重为负数时,可以使用此算法。如果遇到负环,将引发错误。计算成本大约为
O[N(N^2 k)]
,其中k
是每个节点的平均连接边数。输入 csgraph 将转换为 csr 表示。- ‘J’ – Johnson 算法。
与 Bellman-Ford 算法类似,Johnson 算法适用于权重为负数的情况。它将 Bellman-Ford 算法与 Dijkstra 算法结合起来,以实现更快的计算。
- directed布尔值,可选
如果为 True (默认),则在有向图上找到最短路径:仅沿路径 csgraph[i, j] 从点 i 移动到点 j。如果为 False,则在无向图上找到最短路径:该算法可以沿 csgraph[i, j] 或 csgraph[j, i] 从点 i 移动到 j
- return_predecessors布尔值,可选
如果为 True,则返回大小为 (N, N) 的前驱矩阵。
- unweighted布尔值,可选
如果为 True,则查找非加权距离。也就是说,不是找到每个点之间的路径,使权重之和最小化,而是找到使边数最小化的路径。
- overwrite布尔值,可选
如果为 True,则用结果覆盖 csgraph。这仅适用于 method == 'FW' 且 csgraph 是一个密集的、c 顺序的数组,其 dtype=float64。
- indices类数组或整数,可选
如果指定,则仅计算给定索引处的点的路径。与 method == 'FW' 不兼容。
- 返回:
- dist_matrixndarray
图节点之间距离的 N x N 矩阵。dist_matrix[i,j] 给出了沿图从点 i 到点 j 的最短距离。
- predecessorsndarray
仅当 return_predecessors == True 时返回。N x N 的前驱矩阵,可用于重建最短路径。前驱矩阵的第 i 行包含有关从点 i 出发的最短路径的信息:每个条目 predecessors[i, j] 给出了从点 i 到点 j 的路径中前一个节点的索引。如果点 i 和 j 之间不存在路径,则 predecessors[i, j] = -9999
- 引发:
- NegativeCycleError
如果图中存在负环
注释
正如目前所实现的那样,当 directed == False 时,Dijkstra 算法和 Johnson 算法不适用于具有方向依赖距离的图。即,如果 csgraph[i,j] 和 csgraph[j,i] 是不相等的边,则 method='D' 可能会产生错误的结果。
如果可能有多个有效的解决方案,输出可能会因 SciPy 和 Python 版本而异。
示例
>>> from scipy.sparse import csr_matrix >>> from scipy.sparse.csgraph import shortest_path
>>> graph = [ ... [0, 1, 2, 0], ... [0, 0, 0, 1], ... [2, 0, 0, 3], ... [0, 0, 0, 0] ... ] >>> graph = csr_matrix(graph) >>> print(graph) (np.int32(0), np.int32(1)) 1 (np.int32(0), np.int32(2)) 2 (np.int32(1), np.int32(3)) 1 (np.int32(2), np.int32(0)) 2 (np.int32(2), np.int32(3)) 3
>>> dist_matrix, predecessors = shortest_path(csgraph=graph, directed=False, indices=0, return_predecessors=True) >>> dist_matrix array([0., 1., 2., 2.]) >>> predecessors array([-9999, 0, 0, 1], dtype=int32)