图中的K阶邻居 - Python networkx

时间:2013-08-23 02:45:59

标签: python networkx adjacency-list

我有一个有向图,我想有效地找到一个节点的所有K阶邻居的列表。 K阶邻居被定义为可以从完全 K跳跃中的相关节点到达的所有节点。

我查看了networkx,唯一相关的功能是neighbors。但是,这只返回1个邻居的顺序。对于更高阶,我们需要迭代以确定完整集。我认为应该有一种更有效的方式来访问networkx中的K阶邻居。

是否有一个函数可以有效地返回第K个邻居,而无需逐步构建集合?

编辑:如果Python中存在其他可能有用的图形库,请提及。

5 个答案:

答案 0 :(得分:15)

您可以使用: nx.single_source_shortest_path_length(G, node, cutoff=K)

其中G是您的图表对象。

答案 1 :(得分:3)

对于NetworkX,最好的方法可能是在每个k上构建一组邻居。你没有发布你的代码,但似乎你可能已经这样做了:

import networkx as nx

def knbrs(G, start, k):
    nbrs = set([start])
    for l in range(k):
        nbrs = set((nbr for n in nbrs for nbr in G[n]))
    return nbrs

if __name__ == '__main__':
    G = nx.gnp_random_graph(50,0.1,directed=True)
    print(knbrs(G, 0, 3))

答案 2 :(得分:0)

使用修改后的BFS算法解决问题。当您将节点存储在队列中时,也存储它的级别(与根的距离)。处理完节点(所有访问的邻居 - 节点标记为黑色)后,您可以将其添加到其级别的节点列表中。以下是基于this simple implementation的示例:

#!/usr/bin/python
# -*- coding: utf-8 -*-

from collections import defaultdict 
from collections import deque

kth_step = defaultdict(list)

class BFS:
    def __init__(self, node,edges, source):
        self.node = node
        self.edges = edges
        self.source = source
        self.color=['W' for i in range(0,node)] # W for White
        self.graph =color=[[False for i in range(0,node)] for j in range(0,node)]
        self.queue = deque()

        # Start BFS algorithm
        self.construct_graph()
        self.bfs_traversal()

    def construct_graph(self):
        for u,v in self.edges:
            self.graph[u][v], self.graph[v][u] = True, True

    def bfs_traversal(self):
        self.queue.append((self.source, 1))
        self.color[self.source] = 'B' # B for Black
        kth_step[0].append(self.source)

        while len(self.queue):
            u, level =  self.queue.popleft()
            if level > 5: # limit searching there
                return
            for v in range(0, self.node):
                if self.graph[u][v] == True and self.color[v]=='W':
                    self.color[v]='B'
                    kth_step[level].append(v)
                    self.queue.append((v, level+1))

'''
0 -- 1---7
|    |
|    |
2----3---5---6
|
|
4

'''


node = 8 # 8 nodes from 0 to 7
edges =[(0,1),(1,7),(0,2),(1,3),(2,3),(3,5),(5,6),(2,4)] # bi-directional edge
source = 0 # set fist node (0) as source

bfs = BFS(node, edges, source)


for key, value in kth_step.items():
    print key, value

输出:

$ python test.py
0 [0]
1 [1, 2]
2 [3, 7, 4]
3 [5]
4 [6]

我不知道networkx,我也没有准备好在图形工具中使用算法。我相信这样的问题不足以拥有自己的功能。此外,我认为在图形实例中存储任何节点的第k个邻居列表会过于复杂,效率低下且冗余,因此无论如何这样的函数可能必须遍历节点。

答案 3 :(得分:0)

我有一个类似的问题,除了我有一个有向图,我需要维护边属性字典。如果需要,这种相互递归解决方案保留边缘属性字典

def neighbors_n(G, root, n):
    E = nx.DiGraph()

    def n_tree(tree, n_remain):
        neighbors_dict = G[tree]

        for neighbor, relations in neighbors_dict.iteritems():
          E.add_edge(tree, neighbor, rel=relations['rel'])

        #you can use this map if you want to retain functional purity
        #map(lambda neigh_rel: E.add_edge(tree, neigh_rel[0], rel=neigh_rel[1]['rel']), neighbors_dict.iteritems() )

        neighbors = list(neighbors_dict.iterkeys())
        n_forest(neighbors, n_remain= (n_remain - 1))

    def n_forest(forest, n_remain):
        if n_remain <= 0:
            return
        else:
            map(lambda tree: n_tree(tree, n_remain=n_remain), forest)

    n_forest( [root] , n)

    return E

答案 4 :(得分:0)

如前所述,以下解决方案为您提供所有次要邻居(邻居的邻居)并一次性列出所有邻居(该解决方案基于 BFS):

{n: path for n, path in nx.single_source_shortest_path(G, 'a', cutoff=2).items() if len(path)==3}

另一种稍快的解决方案(6.68 µs ± 191 ns vs. 13.3 µs ± 32.1 ns,使用 timeit 测量)包括在无向图中,邻居的邻居可以再次成为源:

def k_neighbors(G, source, cutoff):
    neighbors = {}
    neighbors[0] = {source}
    for k in range(1, cutoff+1):
        neighbors[k] = set()
        for node in level[k-1]:
            neighbors[k].update(set(G.neighbors(node)))
    return neighbors

k_neighbors(B, 'a', 2) #dict keyed with level until `cutoff`, in this case 2

两种解决方案都将源本身作为 0 阶邻居。

所以这取决于你喜欢哪一个。