0

我正在尝试通过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

插件似乎没有正确设置类型向量。

4

1 回答 1

1

感谢您指出这一点！ROI.plugin.lpsolve自 2018 年 12 月 20 日以来，此问题已在 CRAN的新版本 0.3-2 版本中得到修复。

问题在于，lpSolveAPI如果您首先设置边界或类型，则在包中会有所不同。lpSolveAPI首先应该设置边界，然后设置类型。

于 2019-03-01T13:05:36.280 回答