model:
sets:
city/1,2,3/:;
department/1..5/:;
link(city,department):c,x;
link1(department,department)|&2#gt#&1:a;
link2(city,city):b;
endsets
data:
a= 0 1000 1500 0
1400 1200 0
0 2000
700;
b=
100 130 90
130 50 140
90 140 50;
c=
0 0 0 0 0
10 15 10 20 5
10 20 15 15 15 ;
enddata
min=@sum(link1(m,n):a(m,n)*@sum(link2(p,q):b(p,q)*x(p,m)*x(q,n)))/10000-@sum(link:c*x);
@for(department(j):@sum(city(i):x(i,j))=1);
@for(city(i)|i#gt#1:@sum(department(j):x(i,j))<=3);
@for(link:@bin(x));
end