我需要使用 GLPK 解决一个问题,实际上我需要的是计算两个数组中元素之间的差异但是,我在域外收到错误我希望你能修复它,因为我不知道这个错误这是代码
set cols1;
set cols2;
set w;
set range:=1..2;
#set Groups;
#set Items{Groups} ;
param p1{cols1}>=0;
param p2{cols2}>=0;
param uf{w}>=0;
var x>=0;
maximize regret: x;
s.t. a: sum { i in range } (p1[i] -p2[i]) * uf[i]>=x;
s.t. b: sum { i in range } p2[i] * uf[i] =1;
data;
set cols1:=40 110;
set cols2:=51 134;
set w:= 0.8 0.2;
#set Groups := {cols1, cols2, w};
end;
约束b
不包含变量,只包含参数:
s.t. b: sum { i in range } p2[i] * uf[i] =1;
也许这导致了glpk中的域外错误。最好将其重写为检查语句:
check: sum { i in range } p2[i] * uf[i] =1;
也没有参数p1
、p2
和uf
的数据。
外错误位于两个约束的sum{i in range}
语句中。
参数 p1{cols1}
的索引如下
p1[40]
p1[110]
在您的 sum 语句中,您将 p1 与既不40
也不110
的 range
索引。
在这个 gmpl 示例中还有其他一些内容。正确的索引如下所示:
set range:={1,2};
param p1{range}>=0;
param p2{range}>=0;
param uf{range}>=0;
var x>=0;
maximize regret: x;
s.t. a: sum { i in range } (p1[i] -p2[i]) * uf[i]>=x;
s.t. b: sum { i in range } p2[i] * uf[i] =1;
data;
param p1:=1 40 2 110;
param p2:=1 51 2 134;
param uf:= 1 0.8 2 0.2;
end;
但是,问题的技术问题仍然较少。就像 vitaut 已经提到的那样,您的 secound 约束看起来不对,并且会使您的问题仍然不可行。因为您的参数不满足equal 1
约束。
51 * 0.8 + 134 * 0.2 != 1