model:
sets:
num1/1..8/:a,b;
num2/1..3/:l;
aa(num2,num1):x;
endsets
data:
a=5,7,5,4,6,5,5,3;
b=25,36,32,15,31,28,22,12;
enddata
min=0.66*20*Q^0.51+@sum(aa(i,j):x(i,j)*a(j));
m1=@sum(num1(j):x(1,j)*b(j));
m2=0.9*m1+@sum(num1(j):x(2,j)*b(j));
m3=0.9*m2+@sum(num1(j):x(3,j)*b(j));
46.5+m1>=150;
40.7+m2>=160;
34.7+m3>=170;
28.4+0.9*m3+Q>=180;
22.4+0.81*m3+Q>=190;
@for(num2(i):@bnd(0,l(i),8.68);
@bnd(100,Q,167.8);
@for(aa(i,j):@bin(x(i,j)));
@for(num1(j):@sum(num2(i):x(i,j))<=1);
@for(num2(i):@gin(0.66*l(i)*Q^0.51));
@sum(num2(i):l(i))=20;
end

error code 1017
主要看看min=0.66*20*Q^0.51+@sum(aa(i,j):x(i,j)*a(j));
m1=@sum(num1(j):x(1,j)*b(j));
m2=0.9*m1+@sum(num1(j):x(2,j)*b(j));
m3=0.9*m2+@sum(num1(j):x(3,j)*b(j));
有没有错，我想xij*aj，i=1..3,j=1..8
和x1j*bj,j=1..8