sets:
hangye/1..19/:y,b;
qi/1..3/;
touru/1..3/:v,s,l;
link(hangye,touru):x;
link2(touru,qi):h;
Endsets
data:
x=115.89         7.97          110.51
10.04         1.58          25.11
10.82         2.35          40.14
54.36         5.05          61.10
4.32          0.24          3.06
11.64         2.86          20.98
1.74          0.38          1.47
8.29          1.67          13.10
87.62         2.05          273.08
37.04          3.47          175.90
8.56          1.21          8.94
6.90          1.25          9.77
46.04         8.41          265.89
187.14         15.12         908.43
4.77          1.05          7.25
55.73         9.75          58.86
34.54         5.57          42.91
301.94         7.06          435.48
13.20         0.81          7.;
y=118.72
15.03
11.32
62.87
27.07
21.40
2.31
8.84
194.55
42.93
18.63
13.86
48.17
232.92
11.78
52.63
49.53
107.52
3.29 ;
b=1093.32
69.21
121.23
164.25
8.42
55.15
4.02
35.52
2552.54
351.91
24.53
25.64
713.80
2695.56
20.10
124.89
112.92
7033.80
20.99 ;
enddata
Min=@sum(hangye(i):a+@sum(touru(j):l(j)*@log(x(i,j)))+0.5*@sum(touru(j):@sum(link2(j,t):h(j,t)*@log(x(i,j))*@log(x(i,t))))+r*@log(y(i))+0.5*e*(@log(y(i)))^2+p*@log(b(i))+0.5*q*(@log(b(i)))^2+u*@log(y(i))*@log(b(i))+@sum(touru(j):v(j)*@log((x(i,j))*@log(y(i))))+@sum(touru(j):s(j)*@log(x(i,j))* @log(b(i))));
@sum(hangye(i):a+@sum(touru(j):l(j)*@log(x(i,j)))+0.5*@sum(touru(j):@sum(link2(j,t):h(j,t)*@log(x(i,j))*@log(x(i,t))))+r*@log(y(i))+0.5*e*(@log(y(i)))^2+p*@log(b(i))+0.5*q*(@log(b(i)))^2+u*@log(y(i))*@log(b(i))+@sum(touru(j):v(j)*@log((x(i,j))*@log(y(i))))+@sum(touru(j):s(j)*@log(x(i,j))* @log(b(i))))>=0;
@for(hangye(i):p+q*@log(b(i))+u*@log(y(i))+@sum(touru(j):s(j)*@log(x(i,j)))<=0);
@for(hangye(i):r+e*@log(y(i))+u*@log(b(i))+@sum(touru(j):v(j)*@log(x(i,j)))>=0);
r-p=-1;
e-u=0;
u-q=0;
@for(qi(j):v(j)-s(j)=0);
@for(link2(j,t) | j #ne# t:h(j,t)=h(j,t));
end