lingo程序哪里错？

a147596

2000年b题
3 d9 t+ A5 Y- x' r/ R

B题
钢管订购和运输

6 y  ~: u2 t9 W0 G3 t4 Y4 j, d2 A& D
  u# I0 f( H9 z: |* t, {+ v要铺设一条 的输送天然气的主管道, 如图一所示(见下页)。经筛选后可以生产这种主管道钢管的钢厂有 。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。
6 m" g/ n) ~8 s$ j% g为方便计，1km主管道钢管称为1单位钢管。
一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂 在指定期限内能生产该钢管的最大数量为 个单位，钢管出厂销价1单位钢管为 万元，如下表：

	1	2	3	4 , E' G7 f# e% E7 G( T. V	5 1 Q4 s- B7 M. `  E6 }7 O	6	7 ' b& C/ J- A; o# h3 Q
	800 + {6 r3 @* H- o+ ?	800 3 d/ V, q% O; s# {	1000	2000 ! h; b: N, |  L8 D	2000 + S" H6 y' X" {$ r/ m8 Q	2000 * _: n) B" q" s5 U	3000
7 l; ~) p& q+ U  O6 H% b	160	155 4 D- u( O2 P, y4 S8 [	155 1 Y+ \3 w2 ^2 @9 E3 V	160 3 Y5 r3 R, t% x+ |: h, b	155	150	160 ! i- |) c) }, p# j

7 H0 K. P( A7 N5 L1单位钢管的铁路运价如下表：
! y5 |: C( f0 k1 D2 P* |8 e3 C
2 c9 c8 u- b  v# d0 n7 m' b( A0 o2 f

里程(km) 4 C6 t; S" F" }8 j; c% @	≤300 ; X  L. g5 A; o. g5 m	301～350   G7 `  u6 Q9 |2 x. ~" I	351～400	401～450	451～500
运价(万元) & O2 O! j0 y' u% r	20 , M2 r2 x0 X8 {$ n6 p	23	26 - u4 O( ~$ M3 z- J	29	32

4 W7 F+ O9 r9 A* j; D* r

里程(km) 1 v/ r+ W/ _: x9 O8 ~- j	501～600	601～700	701～800 5 [% ?# r' G' b7 _, j7 R$ m9 A	801～900	901～1000 ) \3 u! E4 t9 J- g9 M' ~5 V2 a+ S
运价(万元)	37	44 1 Z- `- F1 R" }) h' U$ k! v	50 3 X, Y+ g+ L3 N5 U) f8 q- X' E	55 % a- k6 y" ^+ l- g6 [& I) p	60

1000km以上每增加1至100km运价增加5万元。?
9 b0 @' W: H- C. j. I( e2 q3 n公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。
" \* V1 G/ ^' A3 Y! w( U钢管可由铁路、公路运往铺设地点（不只是运到点 ，而是管道全线）。
（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。
（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。
（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果
model：
sets:
one/1..7/:p,s;
' A8 I3 L2 m3 h' H' a+ S% b8 Ntwo/1..14/:R,Y;
three(one,two):c,k;
; c6 n, g) o/ U9 {  _9 Op=160    155    155    160    155    150    160;
s=800    800    1000   2000   2000   2000   3000;
. m3 C% K) t7 d9 P. d" ]c=110.3  205.3  220.3  255.3  245.3  255.3  260.3;
& T' y; l' Z# L  L9 |8 \  90.2   190.2  200.2  235.2  225.2  235.2  240.2;
( Q3 L3 [" l6 f$ a9 _/ Q- z  135    235    200    230    225    245    240;
3 f6 @* u2 n; h! E" t( J3 I) l- L0 d  38     161    121    156    146    156    160.1;
8 Z. F( n- W) H! C, ~; Q9 s  20.5   145.5  105.5  150.5  130.5  140.5  145.5;
  3.1    138.1  96     131    121    131    136;
: _% {% V( T+ n+ \  21.2   121.2  86.2   126.5  111.2  121.2  126.2;
  64.2   164.2  48.2   134.2  79.2   84.2   79.2;
! P4 c3 t( k/ Z  87     192    82     62     57     62     67;
% |2 Q2 i7 X: b, f# T  96     151    86     51     33     51     38;
1 Q4 ?% _" T4 t8 C1 ]$ c  G5 c: {  106    156    96     61     51     45     56;
' h! D$ c# i7 X9 W  121.2  171.2  111.2  76.2   71.2   26.2   38.2;
7 ]- o+ `% {- \$ ^  128    173    118    83     73     7      26;
  152    192    132    97     87     28     2;
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]));
6 S" T( V: X) P- U% V@for(two(j)sum(ones(i):k(i,j))<a(i)+b(j));
@sum(two(j):k(i,j))>=500;
, F5 W$ G# h# s$ N@for(one(i):k(i,j)<=s(i));
@sum(two(j):a(i))=5171;
@sum(two(j):b(j))=5171;
8 D$ r$ f; i* h$ B: L# b1 I@for(two(J):a(j)<R(j));
+ m$ H# {+ V* g+ o( G@for(two(J):a(j)<Y(j));
end