Set n / n1*n15 /; Set s / a1 /; yupp('n1','n1','a1')=0; Parameter fcost(m,n,s) / $ondelim n1 n2 a1 1334 n2 n1 a1 1724 n1 n3 a1 1962 n3 n1 a1 1145 n1 n4 a1 1961 n4 n1 a1 1942 n1 n5 a1 1391 n5 n1 a1 1153 n1 n6 a1 1421 n6 n1 a1 1895 n1 n7 a1 1771 n7 n1 a1 1912 n1 n8 a1 1035 n8 n1 a1 1811 n1 n9 a1 1673 n9 n1 a1 1711 n1 n10 a1 1547 n10 n1 a1 1757 n1 n11 a1 1723 n11 n1 a1 1778 n1 n12 a1 1190 n12 n1 a1 1106 n1 n13 a1 1264 n13 n1 a1 1805 n1 n14 a1 1370 n14 n1 a1 1101 n1 n15 a1 1629 n15 n1 a1 1954 n2 n3 a1 1966 n3 n2 a1 1308 n2 n4 a1 1626 n4 n2 a1 1538 n2 n5 a1 1929 n5 n2 a1 1115 n2 n6 a1 1704 n6 n2 a1 1306 n2 n7 a1 1021 n7 n2 a1 1072 n2 n8 a1 1777 n8 n2 a1 1512 n2 n9 a1 1161 n9 n2 a1 1767 n2 n10 a1 1031 n10 n2 a1 1150 n2 n11 a1 1966 n11 n2 a1 1191 n2 n12 a1 1457 n12 n2 a1 1383 n2 n13 a1 1209 n13 n2 a1 1588 n2 n14 a1 1506 n14 n2 a1 1168 n2 n15 a1 1762 n15 n2 a1 1359 n3 n4 a1 1548 n4 n3 a1 1041 n3 n5 a1 1291 n5 n3 a1 1020 n3 n6 a1 1348 n6 n3 a1 1484 n3 n7 a1 1053 n7 n3 a1 1938 n3 n8 a1 1127 n8 n3 a1 1893 n3 n9 a1 1807 n9 n3 a1 1617 n3 n10 a1 1309 n10 n3 a1 1451 n3 n11 a1 1519 n11 n3 a1 1303 n3 n12 a1 1844 n12 n3 a1 1702 n3 n13 a1 1093 n13 n3 a1 1587 n3 n14 a1 1448 n14 n3 a1 1618 n3 n15 a1 1798 n15 n3 a1 1798 n4 n5 a1 1472 n5 n4 a1 1292 n4 n6 a1 1190 n6 n4 a1 1191 n4 n7 a1 1156 n7 n4 a1 1634 n4 n8 a1 1328 n8 n4 a1 1886 n4 n9 a1 1869 n9 n4 a1 1416 n4 n10 a1 1322 n10 n4 a1 1699 n4 n11 a1 1892 n11 n4 a1 1712 n4 n12 a1 1003 n12 n4 a1 1688 n4 n13 a1 1255 n13 n4 a1 1585 n4 n14 a1 1088 n14 n4 a1 1757 n4 n15 a1 1169 n15 n4 a1 1189 n5 n6 a1 1368 n6 n5 a1 1555 n5 n7 a1 1441 n7 n5 a1 1060 n5 n8 a1 1139 n8 n5 a1 1996 n5 n9 a1 1549 n9 n5 a1 1949 n5 n10 a1 1297 n10 n5 a1 1105 n5 n11 a1 1455 n11 n5 a1 1701 n5 n12 a1 1786 n12 n5 a1 1355 n5 n13 a1 1912 n13 n5 a1 1945 n5 n14 a1 1321 n14 n5 a1 1982 n5 n15 a1 1196 n15 n5 a1 1161 n6 n7 a1 1173 n7 n6 a1 1659 n6 n8 a1 1253 n8 n6 a1 1510 n6 n9 a1 1186 n9 n6 a1 1022 n6 n10 a1 1787 n10 n6 a1 1391 n6 n11 a1 1477 n11 n6 a1 1824 n6 n12 a1 1372 n12 n6 a1 1070 n6 n13 a1 1518 n13 n6 a1 1270 n6 n14 a1 1192 n14 n6 a1 1480 n6 n15 a1 1802 n15 n6 a1 1625 n7 n8 a1 1023 n8 n7 a1 1181 n7 n9 a1 1505 n9 n7 a1 1031 n7 n10 a1 1222 n10 n7 a1 1900 n7 n11 a1 1413 n11 n7 a1 1170 n7 n12 a1 1711 n12 n7 a1 1667 n7 n13 a1 1140 n13 n7 a1 1624 n7 n14 a1 1576 n14 n7 a1 1302 n7 n15 a1 1678 n15 n7 a1 1484 n8 n9 a1 1119 n9 n8 a1 1087 n8 n10 a1 1010 n10 n8 a1 1315 n8 n11 a1 1043 n11 n8 a1 1109 n8 n12 a1 1565 n12 n8 a1 1474 n8 n13 a1 1629 n13 n8 a1 1520 n8 n14 a1 1123 n14 n8 a1 1261 n8 n15 a1 1264 n15 n8 a1 1116 n9 n10 a1 1011 n10 n9 a1 1547 n9 n11 a1 1790 n11 n9 a1 1763 n9 n12 a1 1662 n12 n9 a1 1713 n9 n13 a1 1365 n13 n9 a1 1200 n9 n14 a1 1303 n14 n9 a1 1324 n9 n15 a1 1113 n15 n9 a1 1850 n10 n11 a1 1993 n11 n10 a1 1540 n10 n12 a1 1835 n12 n10 a1 1350 n10 n13 a1 1556 n13 n10 a1 1357 n10 n14 a1 1337 n14 n10 a1 1361 n10 n15 a1 1617 n15 n10 a1 1696 n11 n12 a1 1423 n12 n11 a1 1565 n11 n13 a1 1296 n13 n11 a1 1962 n11 n14 a1 1654 n14 n11 a1 1369 n11 n15 a1 1607 n15 n11 a1 1635 n12 n13 a1 1675 n13 n12 a1 1142 n12 n14 a1 1741 n14 n12 a1 1825 n12 n15 a1 1626 n15 n12 a1 1783 n13 n14 a1 1279 n14 n13 a1 1734 n13 n15 a1 1556 n15 n13 a1 1705 n14 n15 a1 1881 n15 n14 a1 1641 $offdelim /; Parameter vcost(m,n,s) / $ondelim n1 n2 a1 9 n2 n1 a1 9 n1 n3 a1 1 n3 n1 a1 9 n1 n4 a1 7 n4 n1 a1 8 n1 n5 a1 1 n5 n1 a1 6 n1 n6 a1 8 n6 n1 a1 9 n1 n7 a1 1 n7 n1 a1 7 n1 n8 a1 2 n8 n1 a1 3 n1 n9 a1 4 n9 n1 a1 4 n1 n10 a1 2 n10 n1 a1 2 n1 n11 a1 8 n11 n1 a1 8 n1 n12 a1 3 n12 n1 a1 1 n1 n13 a1 6 n13 n1 a1 6 n1 n14 a1 7 n14 n1 a1 4 n1 n15 a1 3 n15 n1 a1 1 n2 n3 a1 6 n3 n2 a1 9 n2 n4 a1 4 n4 n2 a1 3 n2 n5 a1 4 n5 n2 a1 1 n2 n6 a1 4 n6 n2 a1 1 n2 n7 a1 4 n7 n2 a1 6 n2 n8 a1 7 n8 n2 a1 4 n2 n9 a1 7 n9 n2 a1 9 n2 n10 a1 5 n10 n2 a1 9 n2 n11 a1 8 n11 n2 a1 4 n2 n12 a1 7 n12 n2 a1 7 n2 n13 a1 9 n13 n2 a1 3 n2 n14 a1 8 n14 n2 a1 7 n2 n15 a1 3 n15 n2 a1 5 n3 n4 a1 9 n4 n3 a1 2 n3 n5 a1 6 n5 n3 a1 6 n3 n6 a1 1 n6 n3 a1 4 n3 n7 a1 1 n7 n3 a1 4 n3 n8 a1 9 n8 n3 a1 3 n3 n9 a1 2 n9 n3 a1 1 n3 n10 a1 2 n10 n3 a1 9 n3 n11 a1 3 n11 n3 a1 2 n3 n12 a1 2 n12 n3 a1 5 n3 n13 a1 2 n13 n3 a1 5 n3 n14 a1 9 n14 n3 a1 4 n3 n15 a1 6 n15 n3 a1 6 n4 n5 a1 5 n5 n4 a1 8 n4 n6 a1 6 n6 n4 a1 3 n4 n7 a1 2 n7 n4 a1 3 n4 n8 a1 4 n8 n4 a1 2 n4 n9 a1 3 n9 n4 a1 4 n4 n10 a1 4 n10 n4 a1 5 n4 n11 a1 1 n11 n4 a1 9 n4 n12 a1 9 n12 n4 a1 6 n4 n13 a1 7 n13 n4 a1 8 n4 n14 a1 8 n14 n4 a1 7 n4 n15 a1 9 n15 n4 a1 9 n5 n6 a1 1 n6 n5 a1 6 n5 n7 a1 8 n7 n5 a1 5 n5 n8 a1 1 n8 n5 a1 8 n5 n9 a1 2 n9 n5 a1 4 n5 n10 a1 3 n10 n5 a1 9 n5 n11 a1 2 n11 n5 a1 7 n5 n12 a1 8 n12 n5 a1 3 n5 n13 a1 2 n13 n5 a1 6 n5 n14 a1 1 n14 n5 a1 4 n5 n15 a1 5 n15 n5 a1 2 n6 n7 a1 3 n7 n6 a1 3 n6 n8 a1 7 n8 n6 a1 1 n6 n9 a1 4 n9 n6 a1 2 n6 n10 a1 6 n10 n6 a1 4 n6 n11 a1 8 n11 n6 a1 9 n6 n12 a1 9 n12 n6 a1 8 n6 n13 a1 2 n13 n6 a1 8 n6 n14 a1 5 n14 n6 a1 7 n6 n15 a1 5 n15 n6 a1 9 n7 n8 a1 8 n8 n7 a1 3 n7 n9 a1 1 n9 n7 a1 1 n7 n10 a1 8 n10 n7 a1 6 n7 n11 a1 9 n11 n7 a1 6 n7 n12 a1 9 n12 n7 a1 6 n7 n13 a1 5 n13 n7 a1 5 n7 n14 a1 2 n14 n7 a1 6 n7 n15 a1 3 n15 n7 a1 2 n8 n9 a1 7 n9 n8 a1 2 n8 n10 a1 1 n10 n8 a1 5 n8 n11 a1 6 n11 n8 a1 1 n8 n12 a1 5 n12 n8 a1 4 n8 n13 a1 5 n13 n8 a1 8 n8 n14 a1 5 n14 n8 a1 3 n8 n15 a1 9 n15 n8 a1 2 n9 n10 a1 5 n10 n9 a1 4 n9 n11 a1 2 n11 n9 a1 3 n9 n12 a1 8 n12 n9 a1 4 n9 n13 a1 2 n13 n9 a1 3 n9 n14 a1 5 n14 n9 a1 8 n9 n15 a1 6 n15 n9 a1 2 n10 n11 a1 1 n11 n10 a1 2 n10 n12 a1 4 n12 n10 a1 7 n10 n13 a1 5 n13 n10 a1 7 n10 n14 a1 2 n14 n10 a1 2 n10 n15 a1 9 n15 n10 a1 1 n11 n12 a1 3 n12 n11 a1 2 n11 n13 a1 5 n13 n11 a1 3 n11 n14 a1 8 n14 n11 a1 2 n11 n15 a1 3 n15 n11 a1 3 n12 n13 a1 7 n13 n12 a1 1 n12 n14 a1 5 n14 n12 a1 4 n12 n15 a1 6 n15 n12 a1 5 n13 n14 a1 2 n14 n13 a1 8 n13 n15 a1 3 n15 n13 a1 7 n14 n15 a1 1 n15 n14 a1 5 $offdelim /; Parameter demand(n) / n1 5 n2 14 n3 7 n4 6 n5 11 n6 3 n7 -100 n8 7 n9 7 n10 4 n11 9 n12 12 n13 5 n14 7 n15 3 /; u = -sum(n$(demand(n)<0), demand(n)); xupp(m,n,s)$(fcost(m,n,s) or vcost(m,n,s)) = u; loop((m,n,s)$(fcost(m,n,s)<0), fcost(m,n,s)=0; vcost(m,n,s)=0; );