Question

我的问题与此问题非常相似：Find a subset from a set of integer whose sum is closest to a value

它只讨论了算法，但我想用R解决它。我对R很新，并试图找出解决方案，但我想知道是否有更有效的方法。

以下是我的例子：

# Define a vector, to findout a subset whose sum is closest to the reference number 20. 
A <- c(2,5,6,3,7)

# display all the possible combinations
y1 <- combn(A,1)
y2 <- combn(A,2)
y3 <- combn(A,3)
y4 <- combn(A,4)
y5 <- combn(A,5)
Y <- list(y1,y2,y3,y4,y5)

# calculate the distance to the reference number of each combination
s1 <- abs(apply(y1,2,sum)-20)
s2 <- abs(apply(y2,2,sum)-20)
s3 <- abs(apply(y3,2,sum)-20)
s4 <- abs(apply(y4,2,sum)-20)
s5 <- abs(apply(y5,2,sum)-20)
S <- list(s1,s2,s3,s4,s5)

# find the minimum difference
M <- sapply(S,FUN=function(x) list(which.min(x),min(x)))
Mm <- which.min(as.numeric(M[2,]))

# return the right combination
data.frame(Y[Mm])[as.numeric(M[,Mm[1]])]

所以答案是2,5,6,7。

如何改进此计划？特别是五个combn（）和五个apply（）s，有没有一种方法可以同时使用它们？我希望当A中有更多项目时，我可以使用长度（A）来覆盖它。

Answer 1

这是另一种方法，

l1 <- sapply(seq_along(A), function(i) combn(A, i))
l2 <- sapply(l1, function(i) abs(colSums(i) - 20))

Filter(length, Map(function(x, y)x[,y], l1, sapply(l2, function(i) i == Reduce(min, l2))))
#[[1]]
#[1] 2 5 6 7

最后一行使用Map根据通过查找列表l1中的最小值创建的逻辑列表来索引l2。

Answer 2

combiter库有isubsetv迭代器，它遍历向量的所有子集。结合foreach简化了代码。

library(combiter)
library(foreach)
A <- c(2,5,6,3,7)

res <- foreach(x = isubsetv(A), .combine = c) %do% sum(x)
absdif <- abs(res-20)
ind <- which(absdif==min(absdif))
as.list(isubsetv(A))[ind]

如何找出总和最接近给定数字的给定向量的最佳组合

2 个答案: