R中的Bron-Kerbosch算法

时间:2017-04-02 21:20:16

标签: r algorithm graph igraph

我想手动在igraph中实现Bron-Kerbosch算法。我想用Wikipedia page中的图表复制示例。这是我到目前为止提出的解决方案。

library(igraph)
g <- graph.formula(1-2,1-5,5-2,2-3,3-4,4-5,4-6)
P <- as.vector(V(g))
bron_kerbosch <- function(R,P,X) {

    if (length(P) == 0 && length(X) == 0) {
        print (R)
    } else {
        for (v in P) {
            neis <- as.vector(neighbors(g,v))
            bron_kerbosch(union(R,v),intersect(P,neis),intersect(X,neis))
            P <- c(P[-which(P==v)])
            X <- c(X,v)
        }
    }
}
> bron_kerbosch(c(),P,c())
[1] 1 2 3
[1] 2 4
[1] 3 5
[1] 4 5
[1] 5 6

答案应该是:  [[1,2,5],[2,3],[3,4],[4,5],[4,6]] 我从1开始索引,而不是从0索引。 知道我做错了吗?

1 个答案:

答案 0 :(得分:2)

我刚刚找到了解决方案。我使用了g graph的索引,而不是名字:/。 所以我换了两行:

 P <- V(g)$name
 neis <- neighbors(g,v)$name

现在可以使用字符而不是数字。结果:

 [1] "1" "2" "5"
 [1] "2" "3"
 [1] "5" "4"
 [1] "3" "4"
 [1] "4" "6"

所以顶点不是有序的(5,4)应该是(4,5),但我认为它并不重要。但是将它写成整数而不是字符应该是好的。