查找给定图形(Python)的所有完整子图的有效方法?

时间:2016-10-27 12:37:48

标签: python graph networkx

是否有一种有效的方法可以使用networkx查找给定(无向)图形的所有完全连接的组件(即完整的子图)?例如,我有以下邻接矩阵(没有自循环):

    |0 1 1 0 0|
    |1 0 1 0 0|
G = |1 1 0 1 0|
    |0 0 1 0 1|
    |0 0 0 1 0|

,对应于以下图表enter image description here 代码应该返回以下节点元组:

(0,1), (1,2), (0,2), (3,4), (2,3), (0,1,2)

我知道networkx有寻找周期,强连接组件等的例程,但是我找不到关于完全连接组件的任何内容。如果使用networkx无法实现,那么Numpy + Scipy也可以。非常感谢提前!

修改

这就是我所做的:

import networkx as nx
import itertools


def findsubsets(S, m):
    return set(itertools.combinations(S, m))



A = np.array([[0, 1, 1, 0, 0],
              [1, 0, 1, 0, 0],
              [1, 1, 0, 1, 0],
              [0, 0, 1, 0, 1],
              [0, 0, 0, 1, 0]])


G = nx.from_numpy_matrix(A)

M = np.sqrt(np.size(A))


for m in range(2, M+1):

    for a in findsubsets(range(0, M), m):

        if(nx.number_of_edges(G.subgraph(a)) == (m**2 - m)/2.):

            print nx.nodes(G.subgraph(a))

基本上找到给定的mXm子图的所有可能的mXm子图,然后检查它们是否具有最大(即(m ** 2 - m)/ 2)个连接数。但是我想知道是否有更有效的方法可以做到这一点,因为函数itertools.combinations的性能对于大图不是很好。

1 个答案:

答案 0 :(得分:2)

好的,我找到了。它只是list(nx.find_cliques(G)),因为我不知道在图论中,一个集团是一个完全连接的子图。

修改

更准确地说,list(nx.find_cliques(G))找到了最大的派系,因此它不是我需要的。我在this link找到了类似的帖子。

所以正确的答案是使用list(nx.enumerate_all_cliques(G))。但是,这个函数还返回大小为1的派系,我不喜欢,因为我的图中没有自循环。因此,最终的解决方案是使用以下代码行:

[s for s in nx.enumerate_all_cliques(G) if len(s) > 1]