R -ROI LPSolve插件通过二进制约束获得错误的结果

Question

我试图通过ROI接口(ROI.plugin.lpsolve(使用LP_SOLVE求解器，并在MIP问题上获取随机结果。

这是我期望0/1二进制解决方案的示例：

library(slam)
library(ROI)
library(ROI.plugin.lpsolve)
library(lpSolve)
solver = "lpsolve"
volume <- c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 3L, 3L, 
            3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 10L, 10L, 10L, 
            10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L)
profit <- c(20L, 15L, 12L, 10L, 7L, 5L, 5L, 5L, 4L, 3L, 2L, 1L, 0L, 40L, 
            38L, 36L, 35L, 30L, 26L, 25L, 22L, 20L, 17L, 15L, 10L, 5L, 0L, 
            100L, 60L, 40L, 20L, 10L, 5L, 5L, 5L, 5L, 5L, 0L)
cap <- 65
lp <- OP(objective = profit,
         constraints = L_constraint(L = rbind(volume),
                                    dir = c("<="),
                                    rhs = cap),
         maximum = TRUE)
types(lp) <- rep("B", length(volume))
mod <- ROI_solve(lp, solver)
res <- mod$solution
print(res)
# 65  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  
#  0  0  0  0  0  0  0  0
#  0  0  0  0  0  0  0

正确的解决方案是：

# 1 1 1 1 1 1 1 1 1 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 0 0 0 0 0 0 0 0 0

似乎插件无法正确设置type vector。

Answer 1

谢谢您指出的！这是在新版本的ROI.plugin.lpsolve，版本0.3-2中修复的自2018-12-20以来的Cran上。

问题是，在lpSolveAPI软件包中，如果首先设置界限或类型，则会有所作为。在lpSolveAPI中，应该首先设置边界，然后设置类型。