Prim算法

Question

我正在尝试实现近似算法来解决旅行商问题（TSP），当三角不等式适用于边权重时，可以使用该问题。如Cormen等人所述，算法导论（第3版），伪代码为：

这是一个例子：

我正在努力的是如何在Prim算法生成的树上实现预订树行走。由于这不是二叉搜索树，因此https://en.wikipedia.org/wiki/Tree_traversal#Pre-order_2处给出的伪代码似乎不适用。

相反，节点不具有left和right属性，而是具有key和parent属性。以下是我在Prim算法实现中生成它们的方法（使用小测试用例）：

import math
import copy
import pytest
import pandas as pd
from cached_property import cached_property


class Node(object):
    def __init__(self, key=math.inf, parent=None):
        self.key = key
        self.parent = parent

    def __lt__(self, other):
        return self.key < other.key


class Graph(object):
    def __init__(self, edges):
        self.edges = edges

    @cached_property
    def nodes(self):
        _nodes = set()
        for edge in self.edges:
            _nodes.add(edge[0])
            _nodes.add(edge[1])
        return {node: Node() for node in list(_nodes)}

    @cached_property
    def adj(self):
        A = {node: [] for node in self.nodes}
        for edge in self.edges:
            u, v, _ = edge
            A[u].append(v)
            A[v].append(u)
        return A

    @cached_property
    def w(self):
        N = len(self.nodes)
        none_array = [[None for _ in range(N)] for _ in range(N)]
        df = pd.DataFrame(none_array, index=sorted(self.nodes), columns=sorted(self.nodes))
        for edge in self.edges:
            u, v, weight = edge
            df.set_value(u, v, weight)
            df.set_value(v, u, weight)
        return df

    def mst_prim(self, root):
        r = self.nodes[root]
        r.key = 0
        Q = copy.copy(self.nodes)
        while Q:
            u = min(Q, key=Q.get)
            u_node = Q.pop(u)
            for v in self.adj[u]:
                if v in Q and self.w[u][v] < self.nodes[v].key:
                    self.nodes[v].parent = u
                    self.nodes[v].key = self.w[u][v]


@pytest.fixture
def edges_simple():
    return [('a', 'b', 4),
            ('a', 'h', 8),
            ('b', 'h', 11),
            ('h', 'i', 7),
            ('b', 'c', 8),
            ('h', 'g', 1),
            ('i', 'c', 2),
            ('i', 'g', 6),
            ('c', 'd', 7),
            ('g', 'f', 2),
            ('c', 'f', 4),
            ('d', 'f', 14),
            ('d', 'e', 9),
            ('f', 'e', 10)
            ]

def test_mst_prim(edges_simple):
    graph = Graph(edges_simple)
    graph.mst_prim(root='a')
    # print("\n")
    # for u, node in graph.nodes.items():
    #   print(u, node.__dict__)
    assert graph.nodes['a'].parent is None
    assert graph.nodes['i'].parent == 'c'
    assert graph.nodes['d'].parent == 'c'



if __name__ == "__main__":
    # pytest.main([__file__+"::test_mst_prim", "-s"])
    pytest.main([__file__, "-s"])

如何在此图表上执行预订树遍历？ （请注意，这个问题类似于pre-order traversal of a Minimum spanning tree，但我发现那里的答案相当高级。）

Answer 1

我建议您在名为Node的{​​{1}}班级中添加新列表。

在children算法后，您可以浏览获得的节点并将其添加到其父Prim's。复杂性是children，所以没什么大不了的。之后O(n)遍历将很容易。

但同样，正如您提到的post一样，您必须为您的孩子选择DFS遍历的订单。在您只引用preorder的情况下，例如，无法知道parent孩子是什么。

Answer 2

我有点疑惑为什么Cormen的书中的算法提到了一个前序遍历，因为在最小生成树MST中节点的“子节点”之间没有顺序。

据我了解，您只需在MST上执行深度优先搜索（DFS），如here和here所述。这意味着，如果在某个节点u，您无需特定顺序访问其中一个邻居或“子女”。

您可以使用图表的邻接列表表示来实现DFS，该图表在您的班级中表示为adj。