lingo程序哪里错？

a147596

2000年b题

B题
9 D; [- e! J% |7 H1 P6 i钢管订购和运输

+ v* G6 N8 r8 ?1 c5 i
- m, k( t8 h( ~: p5 X# W要铺设一条 的输送天然气的主管道, 如图一所示(见下页)。经筛选后可以生产这种主管道钢管的钢厂有 。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。
  o/ q: ~& X/ g; T7 c3 ?) v! L为方便计，1km主管道钢管称为1单位钢管。
! N& c" ~: d7 {$ T7 a一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂 在指定期限内能生产该钢管的最大数量为 个单位，钢管出厂销价1单位钢管为 万元，如下表：

	1	2 8 h- O3 t% w: H7 n( z	3	4	5	6	7 / x9 ?, u; O. }4 D! p6 U  n1 J
5 }) t7 s+ B9 S% a) C	800	800	1000 ; A8 G% X- o( W	2000 " ~$ M/ Q  f; B% L! U/ z9 W2 H	2000	2000 4 X2 P8 k4 ~9 C	3000 . K5 [3 D( l- |9 v4 \6 h7 f
	160 5 Q% {( I! ^9 K5 k2 {" h. r% ^	155 9 V# c0 i# r( O# i$ z% ?- p5 H1 ]	155	160	155	150 ) r& H4 P) ~: ~+ v% r6 s3 w# M	160

1单位钢管的铁路运价如下表：

里程(km) 0 ^! R" v1 Q1 m. e8 e  p	≤300	301～350	351～400 8 h0 B7 ^1 q' x; r	401～450	451～500 % t. |* }) s" O0 P8 [% K
运价(万元)	20	23	26 # b8 R* q$ _* r6 }	29	32

里程(km)	501～600 , Z8 M" i9 P5 M: S8 y  o2 F- {	601～700 ' j5 r3 S2 b2 Z. O  d! b9 K' G	701～800	801～900	901～1000 5 F. n9 G" D# j5 g1 ^7 K7 u7 q
运价(万元) - u7 r; x( H# M( o& x9 b( t	37 8 N) W( I/ A1 a% |& y	44 ; ~2 @+ ?8 N- L6 K' {	50	55 # b% z5 H( h4 i8 w# d6 O. n% A	60 ' C; F! H" F7 j9 G6 e! z

; n, K1 z! _+ I7 C1 t# t7 l1000km以上每增加1至100km运价增加5万元。?
公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。
钢管可由铁路、公路运往铺设地点（不只是运到点 ，而是管道全线）。
5 K: v& ]) b0 b- O& ^（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。
（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。
（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果
model：
sets:
one/1..7/:p,s;
# `* c) Z, X1 s! o, \! Q- U5 Ftwo/1..14/:R,Y;
) X. M4 H  ?( I' Nthree(one,two):c,k;
2 J/ u1 ]0 d" f1 l& k1 }p=160    155    155    160    155    150    160;
s=800    800    1000   2000   2000   2000   3000;
c=110.3  205.3  220.3  255.3  245.3  255.3  260.3;
  90.2   190.2  200.2  235.2  225.2  235.2  240.2;
  135    235    200    230    225    245    240;
  38     161    121    156    146    156    160.1;
  20.5   145.5  105.5  150.5  130.5  140.5  145.5;
  3.1    138.1  96     131    121    131    136;
  21.2   121.2  86.2   126.5  111.2  121.2  126.2;
- Z1 q$ o% @6 E4 h1 H. _2 h5 {4 K+ a  64.2   164.2  48.2   134.2  79.2   84.2   79.2;
: K: r8 \% J5 T$ d% J# M  U  87     192    82     62     57     62     67;
  96     151    86     51     33     51     38;
% ^' N- T! F5 [  G! w1 r$ L  106    156    96     61     51     45     56;
  121.2  171.2  111.2  76.2   71.2   26.2   38.2;
  128    173    118    83     73     7      26;
  152    192    132    97     87     28     2;
. S# k; t7 Z2 U- h( aR=104    405    1155   1761   1955   2160   2361  3041   3521   3821   4041   4251   4671   5171;
Y=5067   4766   4016   3410   3216   3216   2810  2130   1650   1350   1130   920    500    0;
4 _' G& L& Q% `7 T: L  [email=min=@sum(three(i,j):p(i)*k(i,j))+@sum(three(i,j):c(i,j)*k(i,j))+@sum(two(j):h*(n))+@sum(two(j):h*g(m]min=@sum(three(i,j):p(i)*k(i,j))+@sum(three(i,j):c(i,j)*k(i,j))+@sum(two(j):h*(n))+@sum(two(j):h*g(m[/email]));
; l# i: j! ?& I* w: Z@for(two(j)sum(ones(i):k(i,j))<a(i)+b(j));
! x: h7 F9 }; A$ T@sum(two(j):k(i,j))>=500;
@for(one(i):k(i,j)<=s(i));
3 ^: s+ |! ~7 L0 Z@sum(two(j):a(i))=5171;
@sum(two(j):b(j))=5171;
@for(two(J):a(j)<R(j));
1 H# c7 C6 X7 w" B" Z@for(two(J):a(j)<Y(j));
end