model:
sets:
m/1..7/;
n/1..16/:c;
links(m,n):s,x;
endsets
data:
c=1 0.8454 0.4412 0.7787 0.6241 0.3899 0.5358 0.6101 0.6316 0.2059 0.1471 0.5219 0.0588 0.1546 0.3594 0.3300;
s=0.8328 0.7284 0.6503 0.7956 0.7225 0.6085 0.7607 0.7306 0.7688 0.5809 0.6099 0.7978 0.6503 0.6563 0.7607 0.7225
0.8700 0.8350 0.6853 0.8350 0.7978 0.7203 0.7978 0.8328 0.8328 0.6375 0.6085 0.8060 0.6563 0.6724 0.7978 0.7978
0.8120 0.7656 0.5681 0.8038 0.7284 0.7607 0.7284 0.7656 0.8328 0.6853 0.7257 0.7656 0.6375 0.7607 0.7284 0.7284
0.8410 0.7306 0.6449 0.8060 0.7574 0.7203 0.7688 0.7924 0.8328 0.6853 0.7257 0.8060 0.6853 0.7607 0.7688 0.7574;
enddata
max=@sum(links(i,j):c(j)*x(i,j)+s(i,j)*x(i,j));
@sum(links(i,j):x(i,j))=8;
@for(n(j):@sum(m(i):x(i,j))<=1);
@for(m(i):@sum(n(j):x(i,j))>=1);
@for(m(i):@sum(n(j):x(i,j))<=2);
@for(links:@bin(x));
end