<P>model:
sets:
aa/1..15/;
ss/1..7/:s;
cost(ss,aa):c,x,vaule;
pp/1..7/:p;
yy/1..14/:y;
ll/1..14/:l;
tt:t;
endsets
data:
c=
264.3 160.3 140.2 98.6 38 20.5 3.1 21.2 64.2 92 96 105 121.2 128 142
309.3 205.3 100.2 171.6 111 95.5 86 71.2 114.2 142 146 156 171.2 178 192
324.3 220.3 200.2 181.6 121 105.5 96 86.2 48.2 82 86 96 111.2 118 132
354.3 250.3 235.2 216.6 156 140.5 131 116.2 84.2 62 51 61 76.2 83 97
349.3 245.3 225.2 206.6 146 130.5 121 111.2 79.2 57 33 51 71.2 73 87
359.3 255.3 235.2 216.6 156 140.5 131 121.2 84.2 62 51 45 26.2 11 28
369.3 265.3 245.2 226.6 166 150.5 141 131.2 99.2 77 66 56 38.2 26 2;
p=160 155 155 160 155 150 160;
l=104 301 750 606 194 205 201 680 480 300 220 210 420 500;
s=800 800 1000 2000 2000 2000 3000;
enddata
min=@sum(cost(i,j):p(i)*x(i,j)+c(i,j)*x(i,j))+@sum(yy(i):((y(i)^2+y(i))+(l(i)-y(i))*(l(i)-y(i)+1))*0.05);
@sum(cost(i,j)|2#le#j:x(i,j))=y(j)+l(j-1)-y(j-1);
@sum(ss(i):x(i,15))=l(14)-y(14);
@for(yy(i):@bnd(0,y(i),l(i)));
@for(aa(i):@for(ss(j):x(i,j)>=0));
y(1)=0;
@for(tt:@bin(t));
@for(ss(i):@bnd(500*t(i),@sum(aa(j):x(i,j)),s(i)*t(i)));
end
</P>
自己改正了 可是无最优解
<P>我自己改过来了 可是约束成整数都没有最优解改后代码如下:</P><P>model:
sets:
aa/1..15/;
ss/1..7/:s;
cost(ss,aa):c,x;
pp/1..7/:p;
yy/1..14/:y;
ll/1..14/:l;
tt/1..7/:t;
endsets
data:
c=
264.3 160.3 140.2 98.6 38 20.5 3.1 21.2 64.2 92 96 105 121.2 128 142
309.3 205.3 100.2 171.6 111 95.5 86 71.2 114.2 142 146 156 171.2 178 192
324.3 220.3 200.2 181.6 121 105.5 96 86.2 48.2 82 86 96 111.2 118 132
354.3 250.3 235.2 216.6 156 140.5 131 116.2 84.2 62 51 61 76.2 83 97
349.3 245.3 225.2 206.6 146 130.5 121 111.2 79.2 57 33 51 71.2 73 87
359.3 255.3 235.2 216.6 156 140.5 131 121.2 84.2 62 51 45 26.2 11 28
369.3 265.3 245.2 226.6 166 150.5 141 131.2 99.2 77 66 56 38.2 26 2;
p=160 155 155 160 155 150 160;
l=104 301 750 606 194 205 201 680 480 300 220 210 420 500;
s=800 800 1000 2000 2000 2000 3000;
enddata
min=@sum(cost(i,j):(p(i)+c(i,j))*x(i,j))+@sum(yy(i):((y(i)^2+y(i))+(l(i)-y(i))*(l(i)-y(i)+1))*0.05);
@for(ss(i):@for(aa(j):@gin(x(i,j))));
@for(ss(i):@for(yy(j)|2#le#j:@sum(cost(i,j):x(i,j))=y(j)+l(j-1)-y(j-1)));
@sum(ss(i):x(i,15))=l(14)-y(14);
@for(yy(i):y(i)<=l(i));
@for(yy(i):y(i)>=0);
@for(ss(i):@for(aa(j):x(i,j)>=0));
y(1)=0;
@for(tt:@bin(t));
@for(ss(i):500*t(i)<=@sum(aa(j):x(i,j)));
@for(ss(i):@sum(aa(j):x(i,j))<=s(i)*t(i));
end
</P> 我打不开啊 !
回第3楼的!
model:sets:
tou/1..4/:s,r,p,u,x,y,q;
endsets
data:
r=0.28 0.21 0.23 0.25;
p=0.01 0.02 0.045 0.065;
u=103 198 52 40;
q=0.025 0.015 0.055 0.026;
t=?;
v=?;
enddata
max <a href="mailt=@sum" target="_blank" >=@sum</A>( tou(i):s(i)*r(i)-x(i)*(p(i)*u(i)+(s(i)-u(i))*p(i)*y(i)))*(1-v)-v*@max(tou(i):s(i)*q(i));
@sum(tou(i):s(i)+x(i)*(p(i)*u(i)+(s(i)-u(i))*p(i)*y(i)))<=t;
@for(tou(i):2*x(i)=@sign(s(i))-@sign(-s(i)));
@for(tou(i):2*y(i)=@sign(s(i)-u(i))-@sign(u(i)-s(i)));
@for(tou(i):@bin(x(i)));
@for(tou(i):@bin(y(i)));
end
<P>可否麻烦楼主也往我的信箱发一份教程,谢谢了!</P> 我在作《内部网信息组织规划问题》的题目,除了层次分析法,还有没有其他方法,请赐教! <P>不好意思,这个主题主要是LINGO的。你可到数学竞赛上问问,我对此也不太熟。</P> <P> 可否麻烦楼主也往我的信箱发一份教程,谢谢了!<a href="mailthjen_118@163.com" target="_blank" >hjen_118@163.com</A></P> 我的邮箱<a href="mailtxz---zc@163.com" target="_blank" >xz---zc@163.com</A> 多谢了 终于可以看到了楼主太棒了