[9] | 1 | /* TRICK, A Transportation Design Problem */ |
---|
| 2 | |
---|
| 3 | /* Translated from the Mosel modeling language to GNU MathProg by |
---|
| 4 | Andrew Makhorin <mao@gnu.org> */ |
---|
| 5 | |
---|
| 6 | /* This example model is described in the article "Formulations and |
---|
| 7 | Reformulations in Integer Programming" by Michael Trick (it is |
---|
| 8 | publicly available at http://mat.gsia.cmu.edu/trick/formul04.pdf). |
---|
| 9 | |
---|
| 10 | This model demonstrates an amazing effect when including in the |
---|
| 11 | formulation an additional constraint, which is redundant even for |
---|
| 12 | LP relaxation, makes the model easy for solving with the B&B. */ |
---|
| 13 | |
---|
| 14 | set TRUCKS := 1..10; |
---|
| 15 | |
---|
| 16 | set PACKAGES := 1..20; |
---|
| 17 | |
---|
| 18 | param capacity{TRUCKS}; |
---|
| 19 | |
---|
| 20 | param size{PACKAGES}; |
---|
| 21 | |
---|
| 22 | param cost{TRUCKS}; |
---|
| 23 | |
---|
| 24 | param can_use{PACKAGES, TRUCKS}; |
---|
| 25 | |
---|
| 26 | var x{PACKAGES, TRUCKS}, binary; |
---|
| 27 | |
---|
| 28 | var y{TRUCKS}, binary; |
---|
| 29 | |
---|
| 30 | minimize total: sum{i in TRUCKS} cost[i] * y[i]; |
---|
| 31 | |
---|
| 32 | f1{i in TRUCKS}: |
---|
| 33 | sum{j in PACKAGES} size[j] * x[j,i] <= capacity[i] * y[i]; |
---|
| 34 | |
---|
| 35 | f2{i in TRUCKS, j in PACKAGES}: |
---|
| 36 | x[j,i] <= y[i]; |
---|
| 37 | |
---|
| 38 | f3{j in PACKAGES}: |
---|
| 39 | sum{i in TRUCKS} can_use[j,i] * x[j,i] = 1; |
---|
| 40 | |
---|
| 41 | redundant_constraint: |
---|
| 42 | sum{i in TRUCKS} capacity[i] * y[i] >= sum{j in PACKAGES} size[j]; |
---|
| 43 | |
---|
| 44 | data; |
---|
| 45 | |
---|
| 46 | param capacity := |
---|
| 47 | [1] 100 [2] 200 [3] 100 [4] 200 [5] 100 [6] 200 [7] 100 [8] 200 |
---|
| 48 | [9] 100 [10] 200; |
---|
| 49 | |
---|
| 50 | param size := |
---|
| 51 | [1] 17 [2] 21 [3] 54 [4] 45 [5] 87 [6] 34 [7] 23 [8] 45 [9] 12 |
---|
| 52 | [10] 43 [11] 54 [12] 39 [13] 31 [14] 26 [15] 75 [16] 48 [17] 16 |
---|
| 53 | [18] 32 [19] 45 [20] 55; |
---|
| 54 | |
---|
| 55 | param cost := |
---|
| 56 | [1] 1 [2] 1.8 [3] 1 [4] 1.8 [5] 1 [6] 1.8 [7] 1 [8] 1.8 [9] 1 |
---|
| 57 | [10] 1.8; |
---|
| 58 | |
---|
| 59 | param can_use (tr): |
---|
| 60 | 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 := |
---|
| 61 | 1 1 1 1 1 1 1 0 0 0 0 1 1 1 1 0 0 0 0 0 0 |
---|
| 62 | 2 1 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 0 0 0 |
---|
| 63 | 3 0 1 1 1 1 0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 |
---|
| 64 | 4 0 0 1 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 0 0 |
---|
| 65 | 5 0 0 1 1 1 1 0 0 0 0 0 0 0 1 1 1 1 1 1 0 |
---|
| 66 | 6 0 0 0 1 1 1 1 0 0 0 0 0 0 1 1 1 0 0 0 0 |
---|
| 67 | 7 0 0 0 0 1 1 1 1 1 0 0 0 0 0 1 1 1 1 0 0 |
---|
| 68 | 8 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 1 1 1 1 1 |
---|
| 69 | 9 0 0 0 0 0 1 1 1 1 0 0 0 0 0 0 0 1 1 1 1 |
---|
| 70 | 10 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 1 1; |
---|
| 71 | |
---|
| 72 | end; |
---|