求问Lingo错误29

1255684595

以下是旅游线路规划模型求解代码
sets:
jingdian/1..16/:c,t;
其中：1,2,...,16表示西安火车站、秦始皇陵兵马俑、华清池、半坡遗址、骊山风景区、陕西省历史
博物馆、大唐芙蓉园、曲江海洋馆、西安古城墙、大雁塔广场、钟鼓楼、曲江寒窑遗址公园、大唐不夜
城、回民街、青龙寺和翠华山这16个地方
c，t分别表示游客在各景点的吃住消费和逗留时间；
links(jingdian,jingdian):r,cc,tt;
!其中：r为0-1变量（0表示两景点不相连，1表示两景点相通）；cc为两景点之间的交通费用；tt为两
景点之间的交通时间;
endsets
data:
！十六个地点之间的出租车时间表
tt =
0 1.5 1.2 0.5 1.5 0.4 0.5 0.3 0.1 0.3 0.2 0.4 0.3 0.2 0.3 1.4
1.5 0 0.5 1 0.5 1.5 2 1.5 1.5 1.5 1.5 2 1.5 1.3 1.5 2.5
1.2 0.5 0 1 1.5 1.5 1.5 1.4 1.2 1 1 1.1 1.2 1 1 2.5
0.5 1 1 0 0.5 0.5 0.4 0.5 0.5 0.6 0.6 0.8 0.5 0.4 0.4 1.5
1.5 0.5 1.5 0.5 0 1.6 1.7 1.8 1.5 1.5 1.4 1.6 1.5 1.2 1.5 2.6
0.4 1.5 1.5 0.6 1.6 0 0.2 0.4 0.4 0.2 0.6 0.2 0.4 0.3 0.5 1.5
0.5 2 1.5 0.4 1.7 0.2 0 0.1 0.5 0.2 0.4 0.2 0.2 0.4 0.3 1.2
0.4 1.5 1.4 0.5 1.8 0.4 0.1 0 0.4 0.2 0.4 0.2 0.2 0.4 0.2 0.9
0.1 1.6 1 0.5 1.6 0.5 0.5 0.4 0 0.3 0.1 0.5 0.2 0.3 0.4 1.5
0.4 1.5 1 0.6 1.5 0.2 0.2 0.2 0.4 0 0.4 0.3 0.1 0.4 0.3 0.8
0.2 1.5 1 0.6 1.4 0.6 0.4 0.4 0.2 0.4 0 0.4 0.4 0.1 0.4 0.9
0.4 2 1.1 0.8 1.6 0.2 0.2 0.2 0.4 0.3 0.4 0 0.3 0.5 0.3 0.6
0.3 1.5 1.2 0.5 1.5 0.4 0.2 0.2 0.3 0.1 0.4 0.3 0 0.4 0.2 0.8
0.2 1.3 1 0.4 1.2 0.3 0.4 0.4 0.2 0.4 0.1 0.5 0.4 0 0.4 0.9
0.3 1.5 1 0.4 1.5 0.5 0.3 0.2 0.3 0.3 0.4 0.3 0.2 0.4 0 0.7
1.4 2.5 2.5 1.5 2.6 1.5 1.2 0.9 1.4 0.8 0.9 0.6 0.8 0.9 0.7 0
！十六个地点之间的公交车费用表
cc =
0 20 10 3 5 4 4 6 1 5 3 4 5 2 4 16
20 0 2 15 2 15 18 20 20 18 25 23 20 24 23 30
10 2 0 12 25 30 16 18 10 18 16 12 12 15 14 30
3 15 2 0 15 6 8 9 3 6 2 8 10 6 5 18
5 2 25 15 0 15 15 16 5 6 14 15 14 13 15 30
4 15 30 6 15 0 2 3 4 2 2 4 4 6 4 10
4 18 16 8 15 2 0 2 4 5 6 4 2 4 3 10
6 20 18 9 16 3 2 0 6 4 5 5 3 4 6 15
1 19 10 2 6 4 5 5 1 5 5 4 5 3 5 15
5 18 18 6 6 2 4 4 5 0 4 4 4 5 8 18
2 25 16 2 14 2 6 5 2 3 0 4 5 0 3 15
4 23 12 8 15 4 4 5 4 4 4 0 3 4 6 15
5 20 12 10 14 4 2 3 5 4 5 3 0 5 7 14
2 24 15 6 13 2 4 4 2 5 1 4 5 0 4 15
4 23 14 5 15 4 3 6 4 8 3 6 7 4 0 17
16 30 30 18 30 10 10 15 16 18 15 15 14 15 17 0
n=?;
!其中：n表示计划游玩的景点数目;
enddata
min=@sum(jingdian(j):@sum(jingdian(i):r(i,j)*(cc(i,j)+0.5*(c(i)+c(j)))));
min=@sum(jingdian(j):@sum(jingdian(i):100*w(j)*r(i,j)*(cc(i,j)+0.5*90*(c(i)
+c(j)))));
min=@sum(jingdian(j):@sum(jingdian(i):r(i,j)*(cc(i,j)+0.5*(c(i)+c(j)))));
@for(jingdian(i):r(i,i)=0);
@for(jingdian(i)|i#ge#2:@for(jingdian(j)|j#ge#2:r(i,j)+r(j,i)<1));
a=@sum(jingdian(j):@sum(jingdian(i):r(i,j)*(tt(i,j)+0.5*(t(i)+t(j))));
@sum(jingdian(j):@sum(jingdian(i):r(i,j)*(tt(i,j)+0.5*(t(i)+t(j)))))<120;
@for(jingdian(i):@sum(jingdian(j):r(i,j))=@sum(jingdian(j):r(j,i)));
@for(jingdian(i)|i#eq#1:@sum(jingdian(j):r(i,j))=1);
@for(jingdian(i)|i#ne#1:@sum(jingdian(j):r(i,j))<1);
@for(links:@bin(r));
@sum(jingdian(j):@sum(jingdian(i):r(i,j)))=n;
@for(jingdian(i):@for(jingdian(j)|j#gt#1#and#j#ne#i:l(j)>=l(i)+r(i,j)-(n-2)
*(1-r(i,j))+(n-3)*r(j,i)));
@for(jingdian(i)|i#gt#1:l(i)<n-1-(n-2)*r(1,i);l(i)>1+(n-2)*r(i,1));
min=@sum(jingdian(j):@sum(jingdian(i):w(j)*r(i,j)*(cc(i,j)+0.95*c(j)*(1-rrv
(j)*(1-z)))))+@sum(jingdian(j):@sum(jingdian(i):w2(j)*r2(i,j)*(cc(i,j)+0.95
*c(j)*(1-rrv(j)*(1-z)))));
feiyong=@sum(jingdian(j):@sum(jingdian(i):r(i,j)*(cc(i,j)+0.95*c(j)*(1-rrv(
j)*(1-z)))));
feiyong2=@sum(jingdian(j):@sum(jingdian(i):r2(i,j)*(cc(i,j)+0.95*c(j)*(1-rr
v(j)*(1-z)))));
@sum(jingdian(j):@sum(jingdian(i):w(j)*r(i,j)*(cc(i,j)+0.95*c(j)*(1-rrv(j)*
(1-z)))))+@sum(jingdian(j):@sum(jingdian(i):w2(j)*r2(i,j)*(cc(i,j)+0.95*c(j)
*(1-rrv(j)*(1-z)))))>1493;
@for(jingdian(i):r(i,i)=0);@for(jingdian(i):r2(i,i)=0);
@for(jingdian(i)|i#ge#2:@for(jingdian(j)|j#ge#2:r(i,j)+r(j,i)<1));
@for(jingdian(i)|i#ge#2:@for(jingdian(j)|j#ge#2:r2(i,j)+r2(j,i)<1));
@sum(jingdian(j):@sum(jingdian(i):r(i,j)*(tt(i,j)+t(j))))<120;
@sum(jingdian(j):@sum(jingdian(i):r2(i,j)*(tt(i,j)+t(j))))<120;
@for(jingdian(i):@sum(jingdian(j):r(i,j))=@sum(jingdian(j):r(j,i)));
@for(jingdian(i):@sum(jingdian(j):r2(i,j))=@sum(jingdian(j):r2(j,i));
@for(jingdian(i)|i#eq#1:@sum(jingdian(j):r(i,j))=1);
@for(jingdian(i)|i#eq#1:@sum(jingdian(j):r2(i,j))=1);
@for(jingdian(i)|i#ne#1:@sum(jingdian(j):r(i,j))<1);
@for(jingdian(i)|i#ne#1:@sum(jingdian(j):r2(i,j))<1);
@for(links:@bin(r));
@for(links:@bin(r2));
@for(jingdian:@bin(rrv));
@for(links:@bin(rrr));
@sum(jingdian(j):@sum(jingdian(i):r(i,j)))=n;
@sum(jingdian(j):@sum(jingdian(i):r2(i,j)))=n2;
!@sum(jingdian(j):@sum(jingdian(i):r(i,j)))>n-0.5;
!@sum(jingdian(j):@sum(jingdian(i):r(i,j)))<n+0.5;
!@sum(jingdian(j):@sum(jingdian(i):r2(i,j)))<n2+0.5;
!@sum(jingdian(j):@sum(jingdian(i):r2(i,j)))>n2-0.5;
@for(jingdian(i):@for(jingdian(j)|j#gt#1#and#j#ne#i:l(j)>=l(i)+r(i,j)-(n-2)
*(1-r(i,j))+(n-3)*r(j,i)));
@for(jingdian(i):@for(jingdian(j)|j#gt#1#and#j#ne#i:l2(j)>=l2(i)+r2(i,j)-(n
2-2)*(1-r2(i,j))+(n2-3)*r2(j,i)));
@for(jingdian(i)|i#gt#1:l(i)<n-1-(n-2)*r(1,i);l(i)>1+(n-2)*r(i,1));
@for(jingdian(i)|i#gt#1:l2(i)<n2-1-(n2-2)*r2(1,i);l2(i)>1+(n2-2)*r2(i,1));
@for(jingdian(i)|1#eq#i:v(i)=1;b(i)=0);
!@for(jingdian(k)|1#lt#k#and#k#le#n:@for(jingdian(i):@for(jingdian(j):x(k,i,
j)=@if(r(i,j)#eq#1#and#v(k-1)#eq#i,j,0)));!v(k)=@sum(jingdian(j):@sum(jingd
ian(i):x(k,i,j))));
@for(jingdian(k)|1#lt#k#and#k#le#n:@for(jingdian(i):@for(jingdian(j):x(k,i,
j)=@if(0.5#le#r(i,j)#and#r(i,j)#le#1.5#and#(i-0.5)#le#v(k-1)#and#v(k-1)#le#
(i+0.5),j,0)));v(k)=@sum(jingdian(j):@sum(jingdian(i):x(k,i,j))));
@for(jingdian(i)|1#eq#i:v2(i)=1;b2(i)=0);
!@for(jingdian(k)|1#lt#k#and#k#le#n2:@for(jingdian(i):@for(jingdian(j):x2(k,
i,j)=@if(r(i,j)#eq#1#and#v2(k-1)#eq#i,j,0)));!v2(k)=@sum(jingdian(j):@sum(j
ingdian(i):x2(k,i,j))));
@for(jingdian(k)|1#lt#k#and#k#le#n2:@for(jingdian(i):@for(jingdian(j):x2(k,
i,j)=@if(0.5#le#r2(i,j)#and#r2(i,j)#le#1.5#and#(i-0.5)#le#v2(k-1)#and#v2(k-
1)#le#(i+0.5),j,0)));v2(k)=@sum(jingdian(j):@sum(jingdian(i):x2(k,i,j))));
@for(jingdian(i)|2#le#i#and#i#le#(n):@for(jingdian(k):@for(jingdian(j):bb(k,
i,j)=@if(v(i)#eq#j#and#v(i-1)#eq#k,t(k)+tt(k,j),0)));b(i)=@sum(jingdian(j):
@sum(jingdian(k):bb(k,i,j))));
@for(jingdian(k)|k#le#n:tv(k)=@sum(jingdian(i)|i#le#(k):b(i))/12);
@for(jingdian(i)|2#le#i#and#i#le#(n2):@for(jingdian(k):@for(jingdian(j):bb2
(k,i,j)=@if(v2(i)#eq#j#and#v2(i-1)#eq#k,t(k)+tt(k,j),0)));b2(i)=@sum(jingdi
an(j):@sum(jingdian(k):bb2(k,i,j))));
@for(jingdian(k)|k#ne#1#and#k#le#n2:tv2(k)=@sum(jingdian(i)|i#le#(k):b2(i))
/12+4);
@for(jingdian(k)|k#eq#1#or#k#gt#n2:tv2(k)=0);
@for(jingdian(k)|k#le#n2:@for(jingdian(j)|j#le#n:rrr(j,k)=@if(v(j)#eq#v2(k)
#and#(tv2(k)-tv(j))#lt#1#and#(tv(j)-tv2(k))#lt#1,1,0)));
@for(links(j,k)|(j#gt#n#and#k#le#n2)#or#(j#le#n#and#k#gt#n)#or#(j#gt#n#and#
k#gt#n):rrr(j,k)=0);
@for(jingdian(j)|j#ge#2#and#j#le#n2:rrv(j)=@sum(jingdian(k):rrr(j,k)));
@for(jingdian(k)|k#eq#1#or#k#gt#n2:rrv(k)=0);
End



tt和cc矩阵我有数过，都是16*16矩阵，但还是一直报这个错：-------------------------------------------------------------------------------- [Error Code:   29]

  Invalid number of initialization values in a DATA/INIT/CALC
  statement.  The required number of values is:         256.

     25]  cc =
          ^

--------------------------------------------------------------------------------