Model:
sets:
first/1..3/:L;
second/1..8/:x,z;
third/1..5/:c;
mode(first,second):y;
endsets
data:
c=106.333,129.923,149.81,165.7,179.52;
z=5,7,5,4,6,5,5,3;
x=25,36,32,15,31,28,22,12;
enddata
min=@sum(mode(i,j):y(i,j)*z(j))+@sum(first(i):0.66*L(i)*Q^0.51);
@sum(second(j):y(1,j)*x(j))>=c(1);
@sum(second(j):y(2,j)*x(j))+@sum(second(j):y(1,j)*x(j))*0.9>=c(2);
@sum(second(j):y(3,j)*x(j))+@sum(second(j):y(2,j)*x(j)*0.9+@sum(second(j):y(1,j)*x(j))*(0.9^2)>=c(3);
Q+@sum(second(j):y(4,j)*x(j))+@sum(second(j):y(3,j)*x(j))*0.9+@sum(second(j):y(2,j)*x(j))*(0.9^2)+@sum(second(j):y(1,j)*x(j))*(0.9^3)>=c(4);
Q+@sum(second(j):y(5,j)*x(j))+@sum(second(j):y(4,j)*x(j))*0.9+@sum(second(j):y(3,j)*x(j))*(0.9^2)+@sum(second(j):y(2,j)*x(j))*(0.9^3)+@sum(second(j):y(1,j)*x(j))*(0.9^4)>=c(5);
@for(first(i):@sum(second(j):y(i,j)*z(j))+(0.66*L(i)*Q^0.51)=0&y(i,j)=c(4);这一行，y的下标不对，y是一个三行八列的矩阵，不会有y(4,j) |
