sets:
x05dy/1..21/;
x05nl/1..7/:qz;
x05dw/1..7/:d;
x05link(dy,nl):C;
links(nl,dw):Z;
x05DD(dy,nl,dw):X;
x05link2(dy,dw):y;
endsets
data:
qz=0.350396277663283 0.237472604788531 0.158965769190678 0.105558337909803 0.069645448792050 0.046163178982660 0.031798382672994;
C=86 90 71 80 78 88 93
82 87 91 70 85 91 90
80 75 88 93 86 78 78
90 88 95 79 77 82 90
87 83 90 86 65 84 91
90 91 73 67 77 89 81
88 87 85 86 90 94 85
91 88 70 65 61 78 84
73 89 91 87 88 78 80
87 88 90 91 83 84 89
80 87 72 90 86 84 79
87 90 80 82 79 90 80
95 78 88 90 76 78 63
93 92 77 82 84 90 88
78 85 84 89 73 88 83
88 90 79 83 87 76 80
95 78 90 91 78 83 85
88 77 89 85 83 90 78
92 89 70 75 78 92 88
89 93 75 80 87 89 89
95 88 76 80 74 86 78;
enddata
max=@sum(nl(j):
x05x05qz(j)*@sum(dy(i):
x05x05x05x05x05C(i,j)*@sum(dw(k):
x05x05x05x05x05x05x05X(i,j,k))));
@for(DD:@bin(X));
@for(dw(k):@sum(link(i,j):X(i,j,k))=7);
@for(dy(i):@sum(dw(m):y(i,m))