model:
sets:
set1/0,1,2,3,4,5/:;
set2/1/:;
set3/1..22/:;
set4/1..6/:h,m,n,om,on;
set5(set2,set3):c;
set7(set2,set3,set4):p,l,v,s,u,x,x1,s1,j1,z1,q1;
set9(set4,set1):g,d,r;
set16(set2,set3,set4,set1),q,a;
endsets
data:
c=15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15;
p=0,0,2,0,0,0,
0,0,0,0,2,0,
0,2,2,0,0,0,
0,0,0,1,0,1,
0,0,0,0,0,0,
0,0,0,0,2,0,
2,0,0,0,0,0,
0,0,0,1,1,2,
2,2,0,0,0,0,
0,0,0,0,3,0,
1,0,0,0,0,1,
0,0,0,1,0,0,
0,0,3,0,0,4,
2,2,2,0,0,0,
0,0,2,1,0,0,
0,0,0,2,3,1,
0,3,2,1,0,0,
0,0,0,1,2,1,
0,1,0,0,0,0,
0,0,0,2,2,0,
1,2,3,0,0,0,
0,0,0,0,0,2;
l=0,0,1,0,0,0,
0,1,0,0,0,0,
0,0,0,0,1,0,
0,0,1,0,0,0,
0,0,0,0,0,0,
0,0,2,0,0,0,
2,1,0,0,0,0,
0,0,0,0,0,1,
0,0,0,0,0,0,
3,0,0,0,0,0,
0,2,1,0,0,0,
1,1,0,0,0,0,
1,0,0,0,0,0,
2,0,0,0,0,0,
0,0,0,0,2,2,
0,0,1,1,0,0,
3,0,0,0,0,0,
0,0,0,0,0,1,
2,1,1,0,0,0,
0,0,0,0,1,1,
0,0,0,1,0,0,
1,0,0,0,0,0;
g=6,5,5,5,5,3,
3,2,5,5,6,0,
3,3,1,3,0,0,
5,2,2,0,0,0,
6,1,0,0,0,0,
6,0,0,0,0,0;
d=3,6,2,4,8,3,
6,1,5,4,5,0,
11,2,4,2,0,0,
4,6,3,0,0,0,
5,6,0,0,0,0,
12,0,0,0,0,0;
r=3,2,2,2,2,0,
0,1,2,2,3,0,
1,1,1,1,0,0,
2,1,2,0,0,0,
3,1,0,0,0,0,
3,0,0,0,0,0;
h=3,2,2,1,3,2;
m=4,3,3,0,2,1;
n=0,1,1,0,1,2;
enddata
@for(set4(t)sum(set5(i,j):u(i,j,t))=m(t););
@for(set9(t,k)sum(set5(i,j)(i,j,t,k))=d(t,k););
@for(set4(t)sum(set5(i,j):s(i,j,t))=h(t););
@for(set9(t,k)sum(set5(i,j):q(i,j,t,k))=g(t,k););
@for(set4(t)sum(set5(i,j):x(i,j,t))=n(t););
@for(set9(t,k)sum(set5(i,j):a(i,j,t,k))=r(t,k););
min=@sum(set4(t)m(t)-on(t));
@for(set4(t)m(t)=@smax(x1(1,1,t)+s1(1,1,t)+j1(1,1,t)+z1(1,1,t),x1(1,2,t)+s1(1,2,t)+j1(1,2,t)+z1(1,2,t),x1(1,3,t)+s1(1,3,t)+j1(1,3,t)+z1(1,3,t),
x1(1,4,t)+s1(1,4,t)+j1(1,4,t)+z1(1,4,t),x1(1,5,t)+s1(1,5,t)+j1(1,5,t)+z1(1,5,t),x1(1,6,t)+s1(1,6,t)+j1(1,6,t)+z1(1,6,t),
x1(1,7,t)+s1(1,7,t)+j1(1,7,t)+z1(1,7,t),x1(1,8,t)+s1(1,8,t)+j1(1,8,t)+z1(1,8,t),x1(1,9,t)+s1(1,9,t)+j1(1,9,t)+z1(1,9,t),
x1(1,10,t)+s1(1,10,t)+j1(1,10,t)+z1(1,10,t),x1(1,11,t)+s1(1,11,t)+j1(1,11,t)+z1(1,11,t),x1(1,12,t)+s1(1,12,t)+j1(1,12,t)+z1(1,12,t),
x1(1,13,t)+s1(1,13,t)+j1(1,13,t)+z1(1,13,t),x1(1,14,t)+s1(1,14,t)+j1(1,14,t)+z1(1,14,t),x1(1,15,t)+s1(1,15,t)+j1(1,15,t)+z1(1,15,t),
x1(1,16,t)+s1(1,16,t)+j1(1,16,t)+z1(1,16,t),x1(1,17,t)+s1(1,17,t)+j1(1,17,t)+z1(1,17,t),x1(1,18,t)+s1(1,18,t)+j1(1,18,t)+z1(1,18,t),
x1(1,19,t)+s1(1,19,t)+j1(1,19,t)+z1(1,19,t),x1(1,20,t)+s1(1,20,t)+j1(1,20,t)+z1(1,20,t),x1(1,21,t)+s1(1,21,t)+j1(1,21,t)+z1(1,21,t),
x1(1,22,t)+s1(1,22,t)+j1(1,22,t)+z1(1,22,t)));
@for(set4(t)n(t)=@smin(x1(1,1,t)+s1(1,1,t)+j1(1,1,t)+z1(1,1,t),x1(1,2,t)+s1(1,2,t)+j1(1,2,t)+z1(1,2,t),x1(1,3,t)+s1(1,3,t)+j1(1,3,t)+z1(1,3,t),
x1(1,4,t)+s1(1,4,t)+j1(1,4,t)+z1(1,4,t),x1(1,5,t)+s1(1,5,t)+j1(1,5,t)+z1(1,5,t),x1(1,6,t)+s1(1,6,t)+j1(1,6,t)+z1(1,6,t),
x1(1,7,t)+s1(1,7,t)+j1(1,7,t)+z1(1,7,t),x1(1,8,t)+s1(1,8,t)+j1(1,8,t)+z1(1,8,t),x1(1,9,t)+s1(1,9,t)+j1(1,9,t)+z1(1,9,t),
x1(1,10,t)+s1(1,10,t)+j1(1,10,t)+z1(1,10,t),x1(1,11,t)+s1(1,11,t)+j1(1,11,t)+z1(1,11,t),x1(1,12,t)+s1(1,12,t)+j1(1,12,t)+z1(1,12,t),
x1(1,13,t)+s1(1,13,t)+j1(1,13,t)+z1(1,13,t),x1(1,14,t)+s1(1,14,t)+j1(1,14,t)+z1(1,14,t),x1(1,15,t)+s1(1,15,t)+j1(1,15,t)+z1(1,15,t),
x1(1,16,t)+s1(1,16,t)+j1(1,16,t)+z1(1,16,t),x1(1,17,t)+s1(1,17,t)+j1(1,17,t)+z1(1,17,t),x1(1,18,t)+s1(1,18,t)+j1(1,18,t)+z1(1,18,t),
x1(1,19,t)+s1(1,19,t)+j1(1,19,t)+z1(1,19,t),x1(1,20,t)+s1(1,20,t)+j1(1,20,t)+z1(1,20,t),x1(1,21,t)+s1(1,21,t)+j1(1,21,t)+z1(1,21,t),
x1(1,22,t)+s1(1,22,t)+j1(1,22,t)+z1(1,22,t)));
@for(set2(i)for(set3(j)for(set4(t):x1(i,j,t)=@sum(set1(k)|k#le#(6-t):q(i,j,t,k))+s(i,j,t))));
@for(set2(i)for(set3(j)for(set4(t):j1(i,j,t)=@sum(set1(k)|k#le#(6-t)(i,j,t,k)+a(i,j,t,k))+u(i,j,t)+x(i,j,t))));
@for(set2(i)for(set3(j)for(set4(t):z1(i,j,t)=p(i,j,t)+@sum(set1(k)|k#le#(t-1)(i,j,t-k,k)))));
@for(set2(i)for(set3(j)for(set4(t):s1(i,j,t)=l(i,j,t)+@sum(set1(k)|k#le#(t-1):q(i,j,t-k,k)+a(i,j,t-k,k)))));
@for(set2(i)for(set3(j)for(set4(t):v(i,j,t)<=15)));
@for(set2(i)for(set3(j)for(set4(t)|t#ge#t+1:v(i,j,t+1)=v(i,j,t)+x1(i,j,t)+j1(i,j,t)-z1(i,j,t)-s1(i,j,t))));
v(1,1,1)=3;v(1,2,1)=3;v(1,3,1)=6;v(1,4,1)=4;v(1,5,1)=0;v(1,6,1)=4;v(1,7,1)=5;v(1,8,1)=5;v(1,9,1)=4;v(1,10,1)=6;v(1,11,1)=5;
v(1,12,1)=3;v(1,13,1)=8;v(1,14,1)=8;v(1,15,1)=7;v(1,16,1)=8;v(1,17,1)=9;v(1,18,1)=5;v(1,19,1)=5;v(1,20,1)=6;v(1,21,1)=7;v(1,22,1)=3;
@for(set16(i,j,t,k)gin(o(i,j,t,k)));
@for(set16(i,j,t,k)gin(q(i,j,t,k)));
@for(set16(i,j,t,k)gin(a(i,j,t,k)));
@for(set7(i,j,t)gin(s(i,j,t)));
@for(set7(i,j,t)gin(u(i,j,t)));
@for(set7(i,j,t)gin(x(i,j,t)));
@for(set7(i,j,t)gin(x1(i,j,t)));
@for(set7(i,j,t)gin(s1(i,j,t)));
@for(set7(i,j,t)gin(j1(i,j,t)));
@for(set7(i,j,t)gin(z1(i,j,t)));
@for(set7(i,j,t)gin(v(i,j,t)));
@for(set16(i,j,t,k)(i,j,t,k)<=15);
@for(set16(i,j,t,k):q(i,j,t,k)<=15);
@for(set16(i,j,t,k):a(i,j,t,k)<=15);
end
|