lingo程序哪里错？

a147596

2000年b题
! h* H7 l# l) X) P  u% q

B题
% [' p9 L$ z- K( U' y+ }钢管订购和运输

: [& A/ t$ J7 ?: R. y; M3 ?要铺设一条 的输送天然气的主管道, 如图一所示(见下页)。经筛选后可以生产这种主管道钢管的钢厂有 。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。
为方便计，1km主管道钢管称为1单位钢管。
, t; l# U. Q* q: e+ H* o& z8 W一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂 在指定期限内能生产该钢管的最大数量为 个单位，钢管出厂销价1单位钢管为 万元，如下表：

	1 2 H! w% |7 _' k- m3 o6 R	2 5 n7 N: d" b8 s% T- F; d	3 ; B" x5 }9 q/ q& L0 [% Z* _	4 : ?) ?$ c3 f" p* l6 r	5	6 4 a" u. }- x0 N	7
3 G* E7 C# Y" N& }	800 7 I" p$ V3 x7 M5 O3 I	800 5 K5 g5 c  D7 E8 A" T$ k. Z	1000 4 ]9 x2 n0 L$ K* H) O( r, w	2000 / c2 O5 F2 L. o. }/ w: Q	2000	2000	3000 8 [! }4 K& ?- y; i3 c
	160 " G/ l- G8 r1 T6 x% D1 b	155	155 + w' C: }) X* ?* Y" e' t	160 % V  a; D: w+ q$ Y* z( A	155	150	160

6 z2 k& x+ X6 L2 {9 V0 n; L1单位钢管的铁路运价如下表：
% d2 H* r2 F" H. T/ u. l" n

里程(km) 1 ^3 x4 }5 p3 k	≤300	301～350	351～400 : j1 K0 i; z: H- D! Y0 f0 Z/ s  |	401～450 * y: n6 c8 [$ z, `; `/ z; [! g	451～500 ) `3 m: C# V6 R$ b* U
运价(万元) 3 Y: u+ k8 J  n6 [* m( [4 B	20	23 ' p3 `0 x1 Z/ y$ D) @5 b	26	29 ' j' [, s0 }4 J7 X7 q( U  `	32 5 ^3 B2 f' S1 j$ m0 Y6 w+ H

. B( R$ U* `# b$ D+ e

里程(km) ( q6 a2 F/ q/ C! h( s( T2 [	501～600	601～700   Z/ u2 J3 l1 F) c$ x	701～800 ( U/ w2 C2 O# ]& m( _	801～900	901～1000 + j  V5 t. P$ O+ D
运价(万元)	37 - s/ X9 B5 |& h- ^* p	44	50	55 ' l1 D* ^- @  [" w- x	60 4 B9 D2 y  Z( V' R; a

1000km以上每增加1至100km运价增加5万元。?
+ \/ e2 j3 O! [+ ^, y公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。
钢管可由铁路、公路运往铺设地点（不只是运到点 ，而是管道全线）。
（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。
3 |# F# J  e+ T( ?" Z（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。
# q' L. f  O& j! j3 v' q3 _（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果
model：
sets:
6 F& @/ g9 M/ Eone/1..7/:p,s;
8 P: Z; k4 a: x6 _& ftwo/1..14/:R,Y;
three(one,two):c,k;
p=160    155    155    160    155    150    160;
s=800    800    1000   2000   2000   2000   3000;
& @. O, F7 h+ F9 Dc=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;
% o% \2 I! X& t' S  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;
! _' U4 d; @( U5 i8 \  64.2   164.2  48.2   134.2  79.2   84.2   79.2;
" p0 n5 y' Q  _7 p  87     192    82     62     57     62     67;
  96     151    86     51     33     51     38;
  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;
' A, j0 g4 \; i9 M6 R  152    192    132    97     87     28     2;
2 H3 _% J4 v1 |7 q/ cR=104    405    1155   1761   1955   2160   2361  3041   3521   3821   4041   4251   4671   5171;
7 ]' ~' n% _+ G5 D7 n7 `5 w1 vY=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;
@for(one(i):k(i,j)<=s(i));
@sum(two(j):a(i))=5171;
@sum(two(j):b(j))=5171;
@for(two(J):a(j)<R(j));
@for(two(J):a(j)<Y(j));
end