r - R lpSolve を使用して線形計画法に 1000 個の制約を追加する

Question

大きな線形モデルを実装しようとしていますが、変数/制約を追加するのに問題があります。たとえば、次の例のように、マトリックス (f.con) に他のゼロ (0) を追加せずに制約を追加したいとします。

library(linprog)
f.obj <- c(3,24,79,140,230,306,338,
           26,78,145,226,309,336,354,
           146,250,312,364,417,496,508,
           314,382,424,472,503,525,548)
f.con <- matrix (c(1,1,1,1,1,1,1, 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, 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, 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)             
                 ,nrow=4, byrow=TRUE)
f.rhs <- c(10,22,8,30)
f.dir <- c("<=", "<=","<=","<=")
lp ("max", f.obj, f.con, f.dir, f.rhs)

n<-lp ("max", f.obj, f.con, f.dir, f.rhs)$solution
n

ありがとうございました！