model:
sets:
pt/1..26/;
road(pt,pt):x,a;
endsets
data:
a=@file('a.txt');
enddata
min=@sum(road(i,j):a*x);
@for(pt(i)|i#ne#22#and#i#ne#16:@sum(pt(k):x(k,i))=@sum(pt(j):x(i,j)));
@sum(pt(j)|j#ne#22:x(22,j))=1;
@sum(pt(k)|k#ne#22:x(k,22))=0;
@sum(pt(j)|j#ne#16:x(16,j))=0;
@for(road(i,j):x(i,j)<=a(i,j));
@for(road(i,j):@bin(x(i,j)));
end