在OMPR优化中,是否可以使用产品表达式作为约束?

我知道 sum_expr 中的函数 ompr 包作为创建带有动态和的约束的方法。但是,我想知道是否有一种方法可以创建一个使用乘积而不是求和的约束。或者这在 线性的 优化?

例如:

library(dplyr)
library(ROI)
library(ROI.plugin.glpk)
library(ompr)
library(ompr.roi)

n <- 20
score <- round(runif(n, 0, 25))
penalties <- round(runif(n, 0, 25))

model <- MIPModel() %>%
  add_variable(x[i], i = 1:n, type = "binary") %>%
  set_objective(sum_expr(score[i] * x[i], i = 1:n), "max") %>%
  add_constraint(sum_expr(penalties[i] * x[i], i = 1:n) <= 100)

result <- solve_model(model, with_ROI(solver = "glpk", verbose = TRUE))
result$solution

而不是 add_constraint(sum_expr()) 有什么办法吗 add_constraint(product_expr()) ?

如果线性优化无法实现,我应该在哪里寻找替代方案?