model:
sets:
input/1..3/;
inputt/1..3/;
qiye/1..10/;
link1(input,qiye):x;
link2(inputt,qiye):z;
link3(input,inputt):g;
endsets
data:
enddata
submodel mylinprog:
min=@sum(qiye(k):a0+@sum(input:a*x(n,k))+c1*y(k)+r1*b(k)+(1/2)*@sum(link3:g(n,m)*
x(n,k)*z(m,k))+(1/2)*c2*y(k)^2+(1/2)*r2*b(k)^2+@sum(input:d*x(n,k)*b(k))+@sum(input:
e*x(n,k)*y(k))+f*y(k)*b(k));
@for(qiye(k):a0+@sum(input:a*x(n,k))+c1*y(k)+r1*b(k)+(1/2)*@sum(link3:g(n,m)*
x(n,k)*z(m,k))+(1/2)*c2*y(k)^2+(1/2)*r2*b(k)^2+@sum(input:d*x(n,k)*b(k))+@sum(input:
e*x(n,k)*y(k))+f*y(k)*b(k)>0);
@for(qiye(k):a0+@sum(input:a*x(n,k))+c1*y(k)+(1/2)*@sum(link3:g(n,m)*x(n,k)*z(m,k))
+(1/2)*c2*y(k)^2+@sum(input:e*x(n,k)*y(k))<0);
@for(qiye(k):@for(input:a+@sum(input:g(n,m)*x(n,k))+d*b+e*y>0));
@for(qiye(k):c1+c2*y(k)+@sum(input:e*x(n,k))+f*b(k)<0);
@for(qiye(k):r1+r2*b(k)+@sum(input:d*x(n,k))+f*y(k)>0);
@for(input:d=e;
c1-r1=-1;
c2=r2;
r2=f;
g(1,2)=g(2,1);
g(1,3)=g(3,1);
g(2,3)=g(3,2));
endsubmodel