lingo程序哪里错？

a147596

2000年b题
& @" S+ c0 [' P) M. I! F3 O

B题
钢管订购和运输

; ?4 i" d0 ]1 ]; P0 n- t6 H5 S) \要铺设一条 的输送天然气的主管道, 如图一所示(见下页)。经筛选后可以生产这种主管道钢管的钢厂有 。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。
为方便计，1km主管道钢管称为1单位钢管。
5 e* r! W6 k9 E0 s一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂 在指定期限内能生产该钢管的最大数量为 个单位，钢管出厂销价1单位钢管为 万元，如下表：

) S6 }- y8 v/ U9 ]	1	2 - ~, l! p$ _3 H1 U! R) r9 |; p	3	4 - O& ]; \) K& a/ U& u	5 3 j7 a' G7 O2 a" f3 G	6	7
	800	800 2 C5 T; f" {+ D9 x3 M	1000 + [0 b$ ~  a9 D	2000 8 O% @4 X: o! Z9 ?$ g	2000	2000	3000
" O# R7 k1 Q4 U( j	160   \1 U6 i8 L$ }! x" Z7 v	155 0 q2 l/ U0 a; v	155	160	155	150 % w# \: v1 X% H1 C( r/ d	160 # |; n- e+ k- u

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

+ J8 X. ^3 `1 X5 l- Z- `

里程(km) ; S( @* e- t+ q( A" z+ N	≤300 ' `4 Z3 m1 f' m6 K. {# |	301～350	351～400	401～450 & z8 F* ~( y( X5 b	451～500
运价(万元)	20	23	26   v; E* ?' e* y- H+ h	29 ) a* q* O# p3 w2 j% }	32

里程(km)	501～600 0 H# W- e2 C5 G/ D8 x	601～700	701～800	801～900 , s6 K0 |6 Z+ b6 W* d	901～1000
运价(万元) $ e( m1 {2 a7 a1 a2 ?	37	44	50 # [, r  y9 r3 P9 I* @	55 / d3 |4 k  b3 A3 {	60

1000km以上每增加1至100km运价增加5万元。?
公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。
4 }" q, R5 M; H7 @" `& a1 H2 A' j钢管可由铁路、公路运往铺设地点（不只是运到点 ，而是管道全线）。
; i8 S/ \' u. u( j% u1 h/ D8 d（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。
" C5 _+ d4 q& W4 Q* h9 E（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。
) X+ q: d+ o/ c8 }（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果
model：
sets:
& r% q* {& T3 Z2 d, Wone/1..7/:p,s;
! b- b/ ?& r+ C/ ~two/1..14/:R,Y;
! E$ D5 \9 Q8 e% F! x* c4 A; P" Qthree(one,two):c,k;
p=160    155    155    160    155    150    160;
s=800    800    1000   2000   2000   2000   3000;
. Y. U* M7 K% h0 Z% n0 ~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;
: W( i- T! t9 ]  135    235    200    230    225    245    240;
" [+ B# K1 ?" N, @: j3 R' T1 b  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;
. C+ ~, n" B3 H) N  F( l; Y  64.2   164.2  48.2   134.2  79.2   84.2   79.2;
  87     192    82     62     57     62     67;
  96     151    86     51     33     51     38;
  106    156    96     61     51     45     56;
- |) D, f1 @9 m. z& D  121.2  171.2  111.2  76.2   71.2   26.2   38.2;
0 [' H3 m1 f; h- A* X2 F7 ~  128    173    118    83     73     7      26;
# A1 Z' h0 e: {. C( r  G  152    192    132    97     87     28     2;
$ J' @) C: }& V7 \R=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;
  [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]));
@for(two(j)sum(ones(i):k(i,j))<a(i)+b(j));
@sum(two(j):k(i,j))>=500;
$ i- i+ S) L. |/ O3 C, _  K. X@for(one(i):k(i,j)<=s(i));
@sum(two(j):a(i))=5171;
5 J% [+ R9 J) D/ p. O@sum(two(j):b(j))=5171;
@for(two(J):a(j)<R(j));
@for(two(J):a(j)<Y(j));
, n3 `0 m+ _8 s* L) Z& [$ yend