どうすればPythonでグラフをきれいに表現できますか?(最初から、つまりライブラリがない場合!)どのデータ構造(たとえば、dicts / tuples / dict(tuples))が高速で、メモリ効率も良いでしょうか?さまざまなグラフ操作を実行できる必要があります。
指摘したように、さまざまなグラフ表現が役立つ場合があります。それらをPythonに実装するにはどうすればよいですか?ライブラリに関しては、この質問はかなり良い答えを持っています。
どうすればPythonでグラフをきれいに表現できますか?(最初から、つまりライブラリがない場合!)どのデータ構造(たとえば、dicts / tuples / dict(tuples))が高速で、メモリ効率も良いでしょうか?さまざまなグラフ操作を実行できる必要があります。
指摘したように、さまざまなグラフ表現が役立つ場合があります。それらをPythonに実装するにはどうすればよいですか?ライブラリに関しては、この質問はかなり良い答えを持っています。
回答:
これはやや古い質問ですが、これに遭遇した人には実用的な答えを出そうと思いました。
接続の入力データを次のようなタプルのリストとして取得するとします。
[('A', 'B'), ('B', 'C'), ('B', 'D'), ('C', 'D'), ('E', 'F'), ('F', 'C')]
Pythonのグラフで最も便利で効率的であることがわかったデータ構造は、セットの辞書です。これは、Graph
クラスの基本的な構造になります。また、これらの接続が円弧(有向、一方向に接続)であるか、エッジ(無向、両方向に接続)であるかを知る必要もあります。メソッドにdirected
パラメーターを追加することでそれを処理しますGraph.__init__
。その他の便利なメソッドも追加します。
import pprint
from collections import defaultdict
class Graph(object):
""" Graph data structure, undirected by default. """
def __init__(self, connections, directed=False):
self._graph = defaultdict(set)
self._directed = directed
self.add_connections(connections)
def add_connections(self, connections):
""" Add connections (list of tuple pairs) to graph """
for node1, node2 in connections:
self.add(node1, node2)
def add(self, node1, node2):
""" Add connection between node1 and node2 """
self._graph[node1].add(node2)
if not self._directed:
self._graph[node2].add(node1)
def remove(self, node):
""" Remove all references to node """
for n, cxns in self._graph.items(): # python3: items(); python2: iteritems()
try:
cxns.remove(node)
except KeyError:
pass
try:
del self._graph[node]
except KeyError:
pass
def is_connected(self, node1, node2):
""" Is node1 directly connected to node2 """
return node1 in self._graph and node2 in self._graph[node1]
def find_path(self, node1, node2, path=[]):
""" Find any path between node1 and node2 (may not be shortest) """
path = path + [node1]
if node1 == node2:
return path
if node1 not in self._graph:
return None
for node in self._graph[node1]:
if node not in path:
new_path = self.find_path(node, node2, path)
if new_path:
return new_path
return None
def __str__(self):
return '{}({})'.format(self.__class__.__name__, dict(self._graph))
find_shortest_path
その他のメソッドを作成する「読者のための演習」として残しておきます。
これを実際に見てみましょう...
>>> connections = [('A', 'B'), ('B', 'C'), ('B', 'D'),
('C', 'D'), ('E', 'F'), ('F', 'C')]
>>> g = Graph(connections, directed=True)
>>> pretty_print = pprint.PrettyPrinter()
>>> pretty_print.pprint(g._graph)
{'A': {'B'},
'B': {'D', 'C'},
'C': {'D'},
'E': {'F'},
'F': {'C'}}
>>> g = Graph(connections) # undirected
>>> pretty_print = pprint.PrettyPrinter()
>>> pretty_print.pprint(g._graph)
{'A': {'B'},
'B': {'D', 'A', 'C'},
'C': {'D', 'F', 'B'},
'D': {'C', 'B'},
'E': {'F'},
'F': {'E', 'C'}}
>>> g.add('E', 'D')
>>> pretty_print.pprint(g._graph)
{'A': {'B'},
'B': {'D', 'A', 'C'},
'C': {'D', 'F', 'B'},
'D': {'C', 'E', 'B'},
'E': {'D', 'F'},
'F': {'E', 'C'}}
>>> g.remove('A')
>>> pretty_print.pprint(g._graph)
{'B': {'D', 'C'},
'C': {'D', 'F', 'B'},
'D': {'C', 'E', 'B'},
'E': {'D', 'F'},
'F': {'E', 'C'}}
>>> g.add('G', 'B')
>>> pretty_print.pprint(g._graph)
{'B': {'D', 'G', 'C'},
'C': {'D', 'F', 'B'},
'D': {'C', 'E', 'B'},
'E': {'D', 'F'},
'F': {'E', 'C'},
'G': {'B'}}
>>> g.find_path('G', 'E')
['G', 'B', 'D', 'C', 'F', 'E']
heapq
libではなく、セットではなくタプルのリストをヒープ化することです。たとえば、グラフは次のようなヒープのディクテーションになります_graph = {'A': heapify([(0.3, 'D'), (0.5, 'B'), (0.75, 'A'), (0.9, 'C')])}
(注:heapify
このように実際に使用するのではなく、libのヘルプを読んでください)。次に、heapq
関数を使用して重み付きエッジを挿入および取得できます。
log
時間アクセスを意味します。しかし、nodeIDとweightの両方をマップするために使用した辞書を拡張する方法は?
NetworkXは素晴らしいPythonグラフライブラリです。あなたはそれがまだしていないあなたが必要とする何かを見つけるのに苦労するでしょう。
そしてそれはオープンソースなので、彼らがアルゴリズムをどのように実装したかを見ることができます。アルゴリズムを追加することもできます。
https://github.com/networkx/networkx/tree/master/networkx/algorithms
graph.py --> class Graph
。そして、私が見たいのは、彼らの使い方__iter__
です。
まず、古典的なリストと行列表現のどちらを選択するかは、目的(表現をどのように処理するか)によって異なります。よく知られている問題とアルゴリズムは、選択に関連しています。抽象表現の種類の選択によって、実装方法が決まります。
第二に、問題は、頂点とエッジが存在の観点からのみ表現されるべきか、それともそれらが何らかの追加情報を運ぶかどうかです。
Pythonの組み込みデータ型の視点から、他の場所に含まれる値はすべて、ターゲットオブジェクトへの(非表示の)参照として表されます。変数(つまり、名前付き参照)の場合、名前と参照は常に(内部)辞書に格納されます。名前が不要な場合は、参照を独自のコンテナに保存できます。ここでは、リストとして抽象化として常にPythonリストが使用されます。
Pythonリストは参照の動的配列として実装され、Pythonタプルは一定の内容を持つ参照の静的配列として実装されます(参照の値は変更できません)。そのため、簡単にインデックスを付けることができます。このように、リストは行列の実装にも使用できます。
行列を表現するもう1つの方法は、標準モジュールによって実装される配列array
です。これは、保存された型、同種の値に関してより制約されます。要素は値を直接格納します。(リストには、代わりに値オブジェクトへの参照が格納されます)。これにより、メモリ効率が向上し、値へのアクセスも高速になります。
ときどき、のようなさらに限定された表現が役立つことがありますbytearray
。