41546156101511111111(2)()()2()ijijikijkkmjkmjjmjkmijikjkmjznbxncyngsdhns约束条件:ijijnxa41511ijkjkijyev61ijijkkxy61011ijkjkmimys61jkmjmknsh41011ijjmimnxh415461561015111111112180ijijikijkkmjkmijikjkmjbxcygsmodel:sets:si/1..4/:;sj/1..15/:d,e;sk/1..6/:v;sm/1..10/:;lij(si,sj):a,b,x;lik(si,sk):c;ljk(sj,sk):;lijk(si,sj,sk):y;lkm(sk,sm):g;ljm(sj,sm):h;ljkm(sj,sk,sm):s;endsetsdata:enddatamin=n*sum+2*@sum(ljkm(j,k,m):d(j)*(h(j,m)-n*s(j,k,m)));sum=@sum(lij(i,j):b(i,j)*x(i,j))+20000+@sum(lik(i,k):c(i,k)*@sum(sj(j):y(i,j,k)))+@sum(lkm(k,m):g(k,m)*@sum(sj(j):s(j,k,m)));@for(lij:n*x=a);@for(sk(k):@sum(lij(i,j):y(i,j,k)*e(j))=v(k));@for(lij(i,j):x(i,j)=@sum(sk(k):y(i,j,k)));@for(ljk(j,k):@sum(si(i):y(i,j,k))=@sum(sm(m):s(j,k,m)));@for(ljm(j,m):n*@sum(sk(k):s(j,k,m))=h(j,m));@for(sj(j):n*@sum(si(i):x(i,j))=@sum(sm(m):h(j,m)));sum=1800000;n=3;end