创建clique图以确定独立集

时间:2014-03-14 01:01:21

标签: python graph-theory

尝试将s-clique转换为s独立集合。下面是代码和最底层的函数" independent_set_decision(H,s)'是我正在努力的方向。我知道我需要创建图的补充,然后检查该图是否是一个集团。但是,我写的功能并不是按照预期创建图形。任何人对该代码有什么问题都有任何建议吗?

# Returns a list of all the subsets of a list of size k
def k_subsets(lst, k):
    if len(lst) < k:
        return []
    if len(lst) == k:
        return [lst]
    if k == 1:
        return [[i] for i in lst]
    return k_subsets(lst[1:],k) + map(lambda x: x + [lst[0]], k_subsets(lst[1:], k-1))

# Checks if the given list of nodes forms a clique in the given graph.
def is_clique(G, nodes):
    for pair in k_subsets(nodes, 2):
        if pair[1] not in G[pair[0]]:
            return False
    return True

# Determines if there is clique of size k or greater in the given graph.
def k_clique_decision(G, k):
    nodes = G.keys()
    for i in range(k, len(nodes) + 1):
        for subset in k_subsets(nodes, i):
            if is_clique(G, subset):
                return True
    return False

def make_link(G, node1, node2):
    if node1 not in G:
        G[node1] = {}
    (G[node1])[node2] = 1
    if node2 not in G:
        G[node2] = {}
    (G[node2])[node1] = 1
    return G

def break_link(G, node1, node2):
    if node1 not in G:
        print "error: breaking link in a non-existent node"
        return
    if node2 not in G:
        print "error: breaking link in a non-existent node"
        return
    if node2 not in G[node1]:
        print "error: breaking non-existent link"
        return
    if node1 not in G[node2]:
        print "error: breaking non-existent link"
        return
    del G[node1][node2]
    del G[node2][node1]
    return G

# This function should use the k_clique_decision function
# to solve the independent set decision problem
def independent_set_decision(H, s):
    nodes = H.keys()
    I = {}
    for node1 in nodes:
        for node2 in H[node1]:
            if (H[node1])[node2] != 1:
                make_link(I,node1,node2)

    return k_clique_decision(I, s)

1 个答案:

答案 0 :(得分:1)

            if (H[node1])[node2] != 1:

您的图表表示不代表缺少非1值的链接。它表示根本不包括相关的dict条目而缺少链接。迭代所有节点而不是仅包含链接的节点,并检查node2中的H[node1]是否为密钥:

for node1 in H:
    for node2 in H:
        if node2 not in H[node1]:
            make_link(I, node1, node2)