一辈子都报答不尽的一个问题，各位大侠来帮忙解决下哈

denniceyang

sets:
Order/1..3/:;
Location/1..3/:Cltl,ml;
Time/1,2/:;
Route/1..3/:Ctl;
AA(Location,Route):m;
BB(Route,Time):x;
CC(Order,Location,Time):TL,LTL,d,e,a;
endsets
max=@sum(Time(t):@sum(Route(r):x(r,t)*Ctl(r)))+@sum(Order(i):@sum(Time(t):@sum(Location(j):LTL(i,j,t) *Cltl(j))));
@for(Order(i):@for(Location(j):@for(Time(t):TL(i,j,t)+LTL(i,j,t)=1)));
@for(Time(t):x(1,t)*PTL+@sum(Order(i):@sum(Location(j)|j#LE#1:LTL(i,j,t)*PLTL>=@sum(Order(i):@sum(Location(j)|j#LE#1:d(i,j,t))))));
@for(Time(t):x(2,t)*PTL+@sum(Order(i):@sum(Location(j)|j#LE#2:LTL(i,j,t)*PLTL>=@sum(Order(i):@sum(Location(j)|j#LE#2:d(i,j,t))))));
@for(Time(t):x(3,t)*PTL+@sum(Order(i):@sum(Location(j)|j#LE#3:LTL(i,j,t)*PLTL>=@sum(Order(i):@sum(Location(j)|j#LE#3:d(i,j,t))))));
@for(Time(t):x(1,t)<=@sum(Order (i):@sum(Location(j)|j#LE#1:TL(i,j,t)*M)));
@for(Time(t):x(2,t)<=@sum(Order (i):@sum(Location(j)|j#LE#2:TL(i,j,t)*M)));
@for(Time(t):x(3,t)<=@sum(Order (i):@sum(Location(j)|j#LE#3:TL(i,j,t)*M)));
@for(Route(r):@for(Time(t):@for(Order(i):@for (Location(j):e(i,j,t)*TL(i,j,t)>=a(i,j,t)+m(j,r));
@for(Time(t):@for(Order(i):@for(Location(j):e(i,j,t) *LTL(i,j,t)>=a(i,j,t)+ml(j));
@for(Time(t):@for(Order(i):@for(Location(j):@bin(TL (i,j,t)));
@for(Time(t):@for(Order(i):@for(Location(j):@bin(LTL (i,j,t)));
@for(Time(t):@for(Route(r):@gin(x(r,t)));
data:
PTL=300;
PLTL=150;
Ctl=2 5 4;
Cltl=1 2 3;
d=60 90 160 200 50 100 100 230 260
   80 100 150 20 60 100 200 230 90;
e=100 150 200 50 100 150 150 200 250
   110 160 210 40 90 140 130 180 230;
a=60 120 160 30 90 120 110 150 200
   50 110 170 10 50 100 90 150 200;
m=15 20 25 20 25 30;
ml=25 30 35;
M=999999999;  
enddata