gams:model_permutations

How can I model in GAMS the following restriction? A vector of integer variables x(i) can assume values 1,2,…,card(i), and each x(i) should be different, i.e.

x(i) <> x(j) for all i <> j

This looks like we need the help of a permutation matrix P=pij, where there is exactly one 1 in each row and column, and the other elements are zero. The identity matrix is a trivial example of a permutation matrix. Other ones just have some rows and columns of the identity matrix interchanged. To model this in a MIP use binary variables P and add the constraints:

alias (i,j); binary variables p(i,j);

equations rowsum(i) colsum(j) ;

rowsum(i).. sum(j, p(i,j)) =e= 1; colsum(j).. sum(i, p(i,j)) =e= 1;

The different values for our vector x can now be calculated as:

x(i) =e= sum(j, ord(j)*p(i,j));

Because x is automatically integer if p is integer, we can declare x as a normal continuous variable.

gams/model_permutations.txt · Last modified: 2007/08/10 09:04 (external edit)