附上我的代码：
obj <- c(110,55,80,65,400,220,35,10,20,0,0,0,0,0,0,0,0,0)

mat <- matrix(c(80,65,250,225,500,300,55,250,300,0,0,0,0,0,0,0,0,0,
                40,20,5,8,70,40,0,75,65,0,0,0,0,0,0,0,0,0,
                1,1,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,1,1,1,1,1,1,1,1,1,
                1,0,0,0,0,0,0,0,0,-30,0,0,0,0,0,0,0,0,
                0,1,0,0,0,0,0,0,0,0,-30,0,0,0,0,0,0,0,
                0,0,1,0,0,0,0,0,0,0,0,-24,0,0,0,0,0,0,
                0,0,0,1,0,0,0,0,0,0,0,0,-24,0,0,0,0,0,
                0,0,0,0,1,0,0,0,0,0,0,0,0,-10,0,0,0,0,
                0,0,0,0,0,1,0,0,0,0,0,0,0,0,-10,0,0,0,
                0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,-15,0,0,
                0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,-12,0,
                0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,-12), nrow = 13, byrow = TRUE)

dir <- c("<=",">=","<=","<=","<=","<=","<=","<=","<=","<=","<=","<=","<=")

rhs<- c(20000,4000,100,6,0,0,0,0,0,0,0,0,0)

type <- c("I","I","I","I","I","I","I","I","I",
          "B","B","B","B","B","B","B","B","B")

res <- Rglpk_solve_LP(obj, mat, dir, rhs, types = type, max = TRUE)

sol <- matrix(res$solution, nrow = 2, byrow = TRUE)

sol
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9]
[1,]   30   30    0    0   10   10    0    8   12
[2,]    1    1    0    0    1    1    0    1    1


手动输入矩阵元素是真的累…………