R で lpSolveAPI を使用して、ソース ロケーションで製品製造をセットアップするコスト負担を含む施設の最適化問題を解決したいと考えています。lprec オブジェクトでこれをモデル化する方法についてのアドバイスを探しています。ここで 80 ページの問題例を使用します。
Minimize {1300000 x11 + 1040000 x12 + 780000 x13 +
780000 x21 + 1300000 x22 + 1040000 23 +
1040000 x31 + 780000 x32 + 1300000 x33 +
1300000 x41 + 780000 x42 + 780000 x43 +
500000 y1 + 500000 y2 + 500000 y3}
Subject to:
x 11 + x 12 + x 13 = l
x 21 + x 22 + x 23 = 1
x 31 + x 32 + x 33 = 1
x 41 + x 42 + x 43 = 1
lprec オブジェクトを次のように設定しました。
library(lpSolveAPI)
lprec <- make.lp(0, 15)
set.objfn(lprec, c(1300000, 1040000, 780000,
780000, 1300000, 1040000,
1040000, 780000, 1300000,
1300000, 780000, 780000,
500000, 500000, 500000
))
add.constraint(lprec, c(1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), "=", 1)
add.constraint(lprec, c(0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0), "=", 1)
add.constraint(lprec, c(0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0), "=", 1)
add.constraint(lprec, c(0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0), "=", 1)
set.type(lprec, c(1:15), "binary")
ただし、lpSolveAPI を使用してこれらの制約をモデル化する方法がわかりません。
1000 x11 + 1000 x21 +500 x31 +500 x41 <= 1500 y1
1000 x12 + 1000 x22 + 500 x32 + 500 x42 <= 1500 y2
1000 x13 + 1000 x23 + 500 x33 + 500 x43 <= 1500 y3
どんなアドバイスでも大歓迎です!