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