Lingo二维矩阵

model:

sets:

factory /1..6/:a;
plant /1..8/:d;
coo(factory,plant):c,x;

endsets

data:

a=60,55,51,43,41,52;
d=35,37,22,32,41,32,43,38;
c=6,2,6,7,4,2,5,8
  4,9,5,3,8,5,8,2
  5,2,1,9,7,4,3,3
  7,6,7,3,9,2,7,1
  2,3,9,5,7,2,6,5
  5,5,2,2,8,1,4,3;
enddata


min=@sum(coo(i,j):c(i,j)*x(i,j));
@for(factory(i):@sum(plant(j):x(i,j))<=a(i));
@for(plant(j):@sum(factory(i):x(i,j))=d(j));

end