我开始编写此 OPL 代码,但一直收到 5002 错误。我试图努力解决它,但对我来说没有任何效果。除了我找不到改造它的方法。代码中一定有一些错误,现在我无法弄清楚。有没有人可以帮助我解决这个问题。我需要在星期一提交我的项目。
range operation=1..8;
range machine=1..2;
range tool=1..10;
dvar boolean y[operation,machine];
dvar boolean x[operation,tool];
dvar boolean w1;
dvar boolean w2;
dvar boolean w3;
float coef1[operation,machine]=[[133,133],
[71,71],
[42,42],
[83,83],
[91,91],
[63,63],
[108,108],
[125,125]];
float coef2[operation,machine]=[[60,60],
[252,252],
[302,302],
[79,79],
[500,500],
[440,440],
[445,445],
[151,151]];
float coef3[operation,machine]=[[4,4],
[4,4],
[4,4],
[6,6],
[6,6],
[6,6],
[4,4],
[6,6]];
maximize sum(i in operation,j in machine) y[i,j]*coef1[i,j];
subject to
{
sum(i in operation,j in machine) y[i,1]*coef2[i,1]<=384;
sum(i in operation,j in machine) y[i,2]*coef2[i,2]<=492;
y[1,1]+y[1,2]==x[1,1]*x[1,7];
y[2,1]+y[2,2]==x[2,1]*x[2,4];
y[3,1]+y[3,2]==x[3,5]*x[3,8];
y[4,1]+y[4,2]==w1; //x[4,1]*x[4,2]*x[4,6]
y[5,1]+y[5,2]==x[5,3]*x[5,8];
y[6,1]+y[6,2]==w2; //x[6,3]*x[6,5]*x[6,8]
y[7,1]+y[7,2]==x[7,2]*x[7,4];
y[8,1]+y[8,2]==w3; //x[8,7]*x[8,8]*x[8,9]
w1<=x[4,1];
w1<=x[4,2];
w1<=x[4,6];
w1>=x[4,1]+x[4,2]+x[4,6]-2;
w2<=x[6,3];
w2<=x[6,5];
w2<=x[6,8];
w2>=x[6,3]+x[6,5]+x[6,8]-2;
w3<=x[8,7];
w3<=x[8,8];
w3<=x[8,9];
w3>=x[8,7]+x[8,8]+x[8,9]-2;
x[1,1]+x[2,1]+x[4,1]<=2;
x[4,2]+x[7,2]<=1;
x[6,3]<=1;
x[2,4]+x[7,4]<=2;
x[3,5]+x[5,5]+x[6,5]<=2;
x[4,6]<=2;
x[1,7]+x[8,7]<=1;
x[3,8]+x[6,8]+x[8,8]<=2;
x[8,9]<=1;
x[5,10]<=1;
y[1,1]+y[1,2]<=1;
y[2,1]+y[2,2]<=1;
y[3,1]+y[3,2]<=1;
y[4,1]+y[4,2]<=1;
y[5,1]+y[5,2]<=1;
y[6,1]+y[6,2]<=1;
y[7,1]+y[7,2]<=1;
y[8,1]+y[8,2]<=1;
sum (i in operation,j in machine) y[i,1]*coef3[i,j]<=6;
sum (i in operation,j in machine) y[i,2]*coef3[i,j]<=6;
}