c-以编程方式访问GLPK中的MathProg解决方案

Question

我在MathProg中表达了一个问题，使用C API似乎无法描述。具体来说，我在变量之间有约束。我已经生成了一个MathProg文件并将其传递给GLPK。它找到了正确的解决方案，但我不知道如何以编程方式访问这个解决方案。返回的glp_prob结构没有行或列。我可以解析求解器打印的解决方案，但我希望有更好的方法。

或者，如果可以使用C API表达变量之间的约束，我怀疑这也会解决我的问题。下面是我的MathProg代码。

param T := 200;
set b, dimen 3;
set C, dimen 2;
set S, dimen 2;
set Q := setof{(i,j,c) in b : j == 1} i;
set I := setof{(i,s) in S} i;
set E := setof{(i,j) in Q cross I} (i, j);
var x{(i,j) in E}, >=0, <=1, binary;
var y{I}, >=0, <=1, binary;
maximize obj :
  sum{(i,j,c) in b} x[i,j] * c;
s.t. q1c:
  sum{(i,s) in S} x[1,i] <= 1;
s.t. q2c:
  sum{(i,s) in S} x[2,i] <= 1;
s.t. size :
  sum{(i,c) in C} c * y[i] <= T;
s.t. c111avail :
  x[1,1] <= y[1];
s.t. c122avail :
  x[1,2] <= y[2];
s.t. c131avail :
  x[1,3] <= y[1];
s.t. c132avail :
  x[1,3] <= y[2];
s.t. c243avail :
  x[2,4] <= y[3];
solve;
printf "set:";
printf {(i,s) in S: y[i] == 1} " %i", i;
printf "nnot set:";
printf {(i,s) in S: y[i] == 0} " %i", i;
printf "n";
data;
set C :=
  1 100
  2 100
  3 100
  4 100
;
set S :=
  1 100
  2 100
  3 200
  4 100
;
set b :=
  1 1 30
  1 2 30
  1 3 75
  1 4 0
  2 1 0
  2 2 0
  2 3 0
  2 4 35
;
end;

Answer 1

在解决问题后，您可以使用一些额外的GLPK API例程。像

glp_get_row_prim(lp, i)
glp_get_col_prim(lp, i)

它将返回行/列的原始值作为浮点数。

为了生成两个变量之间的不等式的矩阵表示（例如使用C API例程），您必须重新公式化方程，以便所有变量都在方程的左手边。例如

x[1, 1] <= y[1];

将转变为

y[1] - x[1, 1] >= 0;