我正在使用 networkx 来处理图表。我有相当大的图表(其中有近 200 个节点),我尝试找到两个节点之间的所有可能路径。但是,据我了解,networkx 只能找到最短路径。我怎样才能不仅获得最短路径,而且获得所有可能的路径?
UPD:路径只能包含每个节点一次。
UPD2:我需要类似 find_all_paths() 函数,在此处描述:python.org/doc/essays/graphs.html 但此函数不适用于大量节点和边缘 =(
igraph是 Python 的另一个图形模块,可以计算给定节点对之间的所有最短路径。计算所有路径没有意义,因为您有无数这样的路径。
计算从顶点 0 开始的所有最短路径的示例:
>>> from igraph import Graph
>>> g = Graph.Lattice([10, 10], circular=False)
>>> g.get_all_shortest_paths(0)
[...a list of 3669 shortest paths starting from vertex 0...]
如果您有 igraph 0.6 或更高版本(这是撰写本文时的开发版本),您也可以将结果限制为get_all_shortest_paths
给定的结束顶点:
>>> g.get_all_shortest_paths(0, 15)
[[0, 1, 2, 3, 4, 14, 15],
[0, 1, 2, 12, 13, 14, 15],
[0, 10, 11, 12, 13, 14, 15],
[0, 1, 11, 12, 13, 14, 15],
[0, 1, 2, 3, 13, 14, 15],
[0, 1, 2, 3, 4, 5, 15]]
当然,你必须小心;例如,假设您有一个 100 x 100 的网格图(可以通过Graph.Lattice([100, 100], circular=False)
in igraph 轻松生成)。从左上角节点到右下角节点的最短路径的数量等于从 200 个元素中选择 100 个元素的可能性的数量(证明:最短路径的长度有 200 条边,其中 100 条将“水平”在网格中,其中 100 个将“垂直”移动)。这可能不适合您的记忆,因此即使计算这两个节点之间的所有最短路径在这里也不可行。
如果您真的需要两个节点之间的所有路径,您可以使用 igraph 重写您提到的网页上给出的函数,这可能会比纯 Python 解决方案更快,因为 igraph 的核心是用 C 实现的:
def find_all_paths(graph, start, end, path=[]):
path = path + [start]
if start == end:
return [path]
paths = []
for node in set(graph.neighbors(start)) - set(path):
paths.extend(find_all_paths(graph, node, end, path))
return paths
它可以通过首先将图形转换为邻接表表示来进行更多优化,因为它可以避免重复调用graph.neighbors
:
def find_all_paths(graph, start, end):
def find_all_paths_aux(adjlist, start, end, path):
path = path + [start]
if start == end:
return [path]
paths = []
for node in adjlist[start] - set(path):
paths.extend(find_all_paths_aux(adjlist, node, end, path))
return paths
adjlist = [set(graph.neighbors(node)) for node in xrange(graph.vcount())]
return find_all_paths_aux(adjlist, start, end, [])
编辑:修复了第一个示例以在 igraph 0.5.3 中工作,而不仅仅是在 igraph 0.6 中。
这个实际上适用于networkx,它是非递归的,这对于大型图可能很好。
def find_all_paths(graph, start, end):
path = []
paths = []
queue = [(start, end, path)]
while queue:
start, end, path = queue.pop()
print 'PATH', path
path = path + [start]
if start == end:
paths.append(path)
for node in set(graph[start]).difference(path):
queue.append((node, end, path))
return paths
Dijkstra 的算法将以类似于广度优先搜索的方式找到最短路径(它将按深度加权的优先级队列替换为 BFS 的朴素队列)。如果您需要一些替代方案,您可以相当简单地扩展它以产生“N”条最短路径,但如果您需要路径大不相同(例如安排安全货车的路线),您可能需要更聪明地选择彼此显着不同的路径。