Sets nn / n1*n89 / ss / a1*a2 /; yupp('n1','n1','a1')=0; Parameter fcost(nn,nn,ss) / $ondelim n3 n5 a1 4100 n5 n3 a1 4100 n5 n15 a1 15900 n6 n7 a1 6050 n7 n6 a1 6050 n7 n8 a1 10000 n8 n7 a1 10000 n8 n9 a1 3500 n9 n8 a1 3500 n9 n10 a1 4000 n9 n25 a1 4500 n10 n9 a1 4000 n10 n11 a1 2400 n11 n10 a1 2400 n11 n12 a1 2100 n12 n11 a1 2100 n12 n13 a1 6500 n13 n12 a1 6500 n13 n14 a1 11910 n15 n4 a1 6840 n15 n5 a1 15900 n15 n17 a1 2000 n16 n17 a1 19450 n16 n18 a1 8200 n17 n15 a1 2000 n17 n16 a1 19450 n17 n24 a1 8600 n18 n16 a1 8200 n18 n19 a1 3650 n19 n18 a1 3650 n19 n20 a1 2350 n20 n19 a1 2350 n20 n21 a1 6900 n21 n20 a1 6900 n21 n22 a1 12100 n21 n37 a1 7420 n22 n21 a1 12100 n22 n23 a1 1000 n22 n26 a1 3250 n22 n28 a1 15200 n23 n22 a1 1000 n23 n24 a1 2280 n24 n17 a1 8600 n24 n23 a1 2280 n25 n9 a1 4500 n25 n26 a1 4950 n26 n22 a1 3250 n26 n25 a1 4950 n27 n28 a1 2400 n27 n30 a1 3720 n28 n22 a1 15200 n28 n27 a1 2400 n28 n29 a1 7925 n29 n28 a1 7925 n29 n68 a1 9150 n31 n32 a1 1800 n32 n31 a2 1800 n32 n33 a1 1000 n33 n32 a2 1000 n33 n34 a1 1200 n34 n33 a2 1200 n34 n35 a1 3840 n34 n38 a1 8495 n35 n34 a1 3840 n35 n36 a1 5760 n36 n35 a1 5760 n36 n37 a1 6425 n36 n44 a1 8280 n37 n21 a1 7420 n37 n36 a1 6425 n38 n34 a1 8495 n38 n39 a1 3000 n39 n38 a1 3000 n39 n40 a1 1100 n40 n39 a1 1100 n40 n41 a1 2720 n41 n40 a1 2720 n41 n42 a1 1890 n42 n41 a1 1890 n42 n43 a1 5720 n43 n42 a1 5720 n43 n44 a1 3120 n44 n36 a1 8280 n44 n43 a1 3120 n44 n45 a1 4350 n45 n44 a1 4350 n45 n46 a1 2605 n46 n45 a1 2605 n46 n47 a1 6530 n46 n64 a1 9440 n47 n46 a1 6530 n47 n48 a1 2400 n48 n47 a1 2400 n48 n49 a1 2880 n49 n48 a1 2880 n49 n50 a1 1740 n50 n49 a1 1740 n50 n51 a1 1200 n51 n52 a2 3000 n52 n53 a2 2500 n52 n63 a2 1500 n53 n52 a2 2500 n53 n54 a2 3000 n54 n53 a2 3000 n54 n55 a2 1800 n55 n54 a2 1800 n55 n56 a2 2000 n56 n55 a2 2000 n56 n57 a1 6600 n57 n56 a1 6600 n57 n58 a1 7150 n58 n57 a1 7150 n58 n59 a1 8000 n62 n60 a2 1800 n63 n62 a2 1200 n64 n46 a1 9440 n64 n65 a1 7715 n65 n64 a1 7715 n65 n66 a1 4130 n66 n65 a1 4130 n66 n67 a1 2135 n67 n66 a1 2135 n67 n68 a1 1485 n68 n29 a1 9150 n68 n67 a1 1485 n68 n69 a1 10200 n70 n69 a1 6700 n70 n71 a1 8790 n71 n70 a1 8790 n71 n72 a1 2605 n72 n71 a1 2605 n72 n85 a1 1500 n73 n74 a1 3170 n74 n73 a1 3170 n74 n75 a1 3290 n74 n75 a2 6810 n74 n83 a1 1230 n75 n74 a1 3290 n75 n74 a2 6810 n75 n76 a1 9555 n76 n75 a1 9555 n76 n77 a1 2040 n77 n76 a1 2040 n77 n78 a1 10800 n78 n77 a1 10800 n78 n79 a1 9400 n82 n81 a1 5600 n82 n83 a1 2900 n83 n74 a1 1230 n83 n82 a1 2900 n85 n72 a1 1500 n85 n84 a1 6500 $offdelim / Parameter vcost(nn,nn,ss) / $ondelim n2 n1 a1 2.5 n2 n4 a1 0.9 n3 n2 a1 2.3 n3 n5 a1 1.4 n3 n6 a1 1.1 n4 n2 a1 0.9 n4 n31 a1 2.8 n5 n3 a1 1.4 n5 n15 a1 1.9 n6 n3 a1 1.1 n6 n7 a1 1.2 n7 n6 a1 1.0 n7 n8 a1 1.8 n8 n7 a1 1.8 n8 n9 a1 0.8 n9 n8 a1 0.6 n9 n10 a1 0.9 n9 n25 a1 1.4 n10 n9 a1 0.9 n10 n11 a1 1.4 n11 n10 a1 1.4 n11 n12 a1 0.7 n12 n11 a1 0.7 n12 n13 a1 0.9 n13 n12 a1 0.9 n13 n14 a1 1.1 n14 n1 a1 3.1 n14 n30 a1 1.8 n15 n4 a1 0.7 n15 n5 a1 1.9 n15 n17 a1 0.9 n16 n17 a1 2.0 n16 n18 a1 2.3 n17 n15 a1 0.7 n17 n16 a1 2.0 n17 n24 a1 2.4 n18 n16 a1 1.4 n18 n19 a1 1.3 n19 n18 a1 1.3 n19 n20 a1 0.8 n20 n19 a1 0.8 n20 n21 a1 2.0 n21 n20 a1 2.0 n21 n22 a1 2.4 n21 n37 a1 2.2 n22 n21 a1 2.4 n22 n23 a1 0.4 n22 n26 a1 1.1 n22 n28 a1 2.5 n23 n22 a1 0.4 n23 n24 a1 0.7 n24 n17 a1 2.1 n24 n23 a1 0.7 n25 n9 a1 1.4 n25 n26 a1 1.6 n26 n22 a1 1.1 n26 n25 a1 1.6 n27 n28 a1 1.0 n27 n30 a1 4.2 n28 n22 a1 2.5 n28 n27 a1 0.9 n28 n29 a1 2.2 n29 n28 a1 2.2 n29 n68 a1 2.5 n30 n87 a1 1.2 n31 n4 a1 2.8 n31 n32 a1 0.4 n31 n59 a1 1.2 n32 n31 a1 0.5 n32 n31 a2 0.4 n32 n33 a1 0.7 n33 n32 a1 0.9 n33 n32 a2 0.7 n33 n34 a1 0.5 n34 n33 a1 0.7 n34 n33 a2 0.5 n34 n35 a1 1.1 n34 n38 a1 1.4 n35 n34 a1 1.1 n35 n36 a1 1.3 n36 n35 a1 1.3 n36 n37 a1 1.6 n36 n44 a1 3.0 n37 n21 a1 2.2 n37 n36 a1 1.6 n38 n34 a1 1.4 n38 n39 a1 0.8 n39 n38 a1 0.8 n39 n40 a1 0.4 n40 n39 a1 0.4 n40 n41 a1 0.7 n41 n40 a1 0.7 n41 n42 a1 1.1 n42 n41 a1 1.1 n42 n43 a1 1.5 n43 n42 a1 0.9 n43 n44 a1 0.9 n44 n36 a1 3.0 n44 n43 a1 0.9 n44 n45 a1 1.5 n45 n44 a1 1.5 n45 n46 a1 1.3 n46 n45 a1 1.3 n46 n47 a1 1.6 n46 n64 a1 2.6 n47 n46 a1 1.6 n47 n48 a1 0.8 n48 n47 a1 0.8 n48 n49 a1 0.8 n49 n48 a1 0.8 n49 n50 a1 1.0 n50 n49 a1 1.0 n50 n51 a1 1.0 n51 n52 a1 2.4 n51 n52 a2 1.8 n52 n53 a1 2.3 n52 n53 a2 1.6 n52 n63 a1 1.4 n52 n63 a2 1.0 n53 n52 a1 2.3 n53 n52 a2 1.6 n53 n54 a1 3.2 n53 n54 a2 2.6 n54 n53 a1 3.2 n54 n53 a2 2.6 n54 n55 a1 0.9 n54 n55 a2 0.6 n55 n54 a1 0.9 n55 n54 a2 0.6 n55 n56 a1 1.8 n55 n56 a2 1.3 n56 n55 a1 1.8 n56 n55 a2 1.3 n56 n57 a1 1.1 n57 n56 a1 1.3 n57 n58 a1 1.1 n58 n57 a1 1.1 n58 n59 a1 1.8 n59 n60 a1 1.2 n60 n61 a1 0.8 n62 n60 a1 1.7 n62 n60 a2 1.3 n62 n79 a1 2.8 n63 n62 a1 1.3 n63 n62 a2 0.8 n64 n46 a1 2.6 n64 n65 a1 1.8 n65 n64 a1 1.8 n65 n66 a1 0.9 n66 n65 a1 0.9 n66 n67 a1 0.8 n67 n66 a1 0.8 n67 n68 a1 0.5 n68 n29 a1 2.5 n68 n67 a1 0.5 n68 n69 a1 1.0 n69 n86 a1 2.5 n70 n69 a1 0.7 n70 n71 a1 3.2 n71 n70 a1 3.2 n71 n72 a1 0.7 n72 n71 a1 0.7 n72 n85 a1 1.1 n73 n74 a1 0.9 n74 n73 a1 0.9 n74 n75 a1 0.9 n74 n75 a2 1.6 n74 n83 a1 0.7 n75 n74 a1 0.9 n75 n74 a2 1.6 n75 n76 a1 2.3 n76 n75 a1 2.3 n76 n77 a1 1.6 n77 n76 a1 1.6 n77 n78 a1 1.2 n78 n77 a1 1.2 n78 n79 a1 1.1 n79 n62 a1 2.8 n79 n80 a1 1.6 n80 n88 a1 3.3 n81 n80 a1 0.9 n81 n84 a1 1.0 n82 n81 a1 1.2 n82 n83 a1 1.0 n83 n74 a1 0.7 n83 n82 a1 1.0 n84 n81 a1 1.0 n84 n86 a1 1.3 n85 n72 a1 1.1 n85 n84 a1 1.1 n86 n87 a1 1.6 n87 n88 a1 0.1 $offdelim / Parameter demand(nn) / n3 -66 n5 -100 n6 -84 n7 -88 n8 -56 n9 -110 n10 -70 n11 -256 n12 -2136 n13 -127 n15 -70 n16 -65 n17 -225 n18 -138 n19 -228 n20 -680 n21 -227 n22 -73 n23 -263 n24 -321 n25 -3 n26 -131 n27 -86 n28 -11 n29 -40 n32 -46 n33 -832 n34 -29 n35 -63 n36 -270 n37 -119 n38 -89 n39 -552 n40 -94 n41 -85 n42 -139 n43 -480 n44 -490 n45 -219 n46 -175 n47 -374 n48 -265 n49 -51 n50 -374 n51 -259 n52 -203 n53 -306 n54 -480 n55 -431 n56 -315 n57 -200 n58 -24 n62 -83 n63 -205 n64 -100 n65 -124 n66 -90 n68 -40 n70 -72 n71 -149 n72 -31 n73 -145 n74 -24 n75 -156 n76 -325 n77 -100 n78 -106 n82 -38 n83 -49 n85 -50 n89 14505 $offdelim /; Parameter xupp(nn,nn,ss) / $ondelim n1 n89 a1 10754 n2 n1 a1 10754 n2 n4 a1 10754 n3 n2 a1 10754 n3 n5 a1 10754 n3 n6 a1 10754 n4 n2 a1 10754 n4 n31 a1 10754 n5 n3 a1 10754 n5 n15 a1 10754 n6 n3 a1 10754 n6 n7 a1 10754 n7 n6 a1 10754 n7 n8 a1 10754 n8 n7 a1 56 n8 n9 a1 10754 n9 n8 a1 10754 n9 n10 a1 10754 n9 n25 a1 10754 n10 n9 a1 10754 n10 n11 a1 10754 n11 n10 a1 10754 n11 n12 a1 10754 n12 n11 a1 10754 n12 n13 a1 10754 n13 n12 a1 10754 n13 n14 a1 10754 n14 n1 a1 10754 n14 n30 a1 10754 n15 n4 a1 10754 n15 n5 a1 10754 n15 n17 a1 10754 n16 n17 a1 65 n16 n18 a1 10754 n17 n15 a1 10754 n17 n16 a1 10754 n17 n24 a1 10754 n18 n16 a1 10754 n18 n19 a1 10754 n19 n18 a1 10754 n19 n20 a1 10754 n20 n19 a1 10754 n20 n21 a1 10754 n21 n20 a1 10754 n21 n22 a1 10754 n21 n37 a1 10754 n22 n21 a1 10754 n22 n23 a1 10754 n22 n26 a1 10754 n22 n28 a1 10754 n23 n22 a1 10754 n23 n24 a1 10754 n24 n17 a1 10754 n24 n23 a1 10754 n25 n9 a1 10754 n25 n26 a1 10754 n26 n22 a1 10754 n26 n25 a1 10754 n27 n28 a1 10754 n27 n30 a1 10754 n28 n22 a1 10754 n28 n27 a1 10754 n28 n29 a1 10754 n29 n28 a1 10754 n29 n68 a1 10754 n30 n87 a1 10754 n31 n4 a1 10754 n31 n32 a1 10754 n31 n59 a1 10754 n32 n31 a1 18000 n32 n31 a2 10754 n32 n33 a1 10754 n33 n32 a1 5000 n33 n32 a2 10754 n33 n34 a1 10754 n34 n33 a1 6000 n34 n33 a2 10754 n34 n35 a1 10754 n34 n38 a1 10754 n35 n34 a1 10754 n35 n36 a1 10754 n36 n35 a1 10754 n36 n37 a1 10754 n36 n44 a1 10754 n37 n21 a1 10754 n37 n36 a1 10754 n38 n34 a1 10754 n38 n39 a1 10754 n39 n38 a1 10754 n39 n40 a1 10754 n40 n39 a1 10754 n40 n41 a1 10754 n41 n40 a1 10754 n41 n42 a1 10754 n42 n41 a1 10754 n42 n43 a1 10754 n43 n42 a1 10754 n43 n44 a1 10754 n44 n36 a1 10754 n44 n43 a1 10754 n44 n45 a1 10754 n45 n44 a1 10754 n45 n46 a1 10754 n46 n45 a1 10754 n46 n47 a1 10754 n46 n64 a1 10754 n47 n46 a1 10754 n47 n48 a1 10754 n48 n47 a1 10754 n48 n49 a1 10754 n49 n48 a1 10754 n49 n50 a1 10754 n50 n49 a1 10754 n50 n51 a1 10754 n51 n52 a1 5000 n51 n52 a2 11013 n52 n53 a1 3571 n52 n53 a2 12972 n52 n63 a1 3750 n52 n63 a2 12972 n53 n52 a1 3571 n53 n52 a2 12972 n53 n54 a1 5000 n53 n54 a2 12972 n54 n53 a1 5000 n54 n53 a2 12972 n54 n55 a1 6000 n54 n55 a2 12972 n55 n54 a1 6000 n55 n54 a2 12972 n55 n56 a1 4000 n55 n56 a2 12972 n56 n55 a1 4000 n56 n55 a2 12972 n56 n57 a1 746 n57 n56 a1 12972 n57 n58 a1 12972 n58 n57 a1 12972 n58 n59 a1 12972 n59 n60 a1 12972 n60 n61 a1 14203 n61 n89 a1 14203 n62 n60 a1 4500 n62 n60 a2 14203 n62 n79 a1 14203 n63 n62 a1 2400 n63 n62 a2 13177 n64 n46 a1 10754 n64 n65 a1 10754 n65 n64 a1 10754 n65 n66 a1 10754 n66 n65 a1 10754 n66 n67 a1 10754 n67 n66 a1 10754 n67 n68 a1 10754 n68 n29 a1 254 n68 n67 a1 10754 n68 n69 a1 394 n69 n86 a1 11056 n70 n69 a1 302 n70 n71 a1 302 n71 n70 a1 302 n71 n72 a1 302 n72 n71 a1 302 n72 n85 a1 302 n73 n74 a1 943 n74 n73 a1 943 n74 n75 a1 943 n74 n75 a2 943 n74 n83 a1 943 n75 n74 a1 943 n75 n74 a2 943 n75 n76 a1 943 n76 n75 a1 943 n76 n77 a1 943 n77 n76 a1 943 n77 n78 a1 943 n78 n77 a1 943 n78 n79 a1 943 n79 n62 a1 14203 n79 n80 a1 14203 n80 n88 a1 14505 n81 n80 a1 1245 n81 n84 a1 1245 n82 n81 a1 943 n82 n83 a1 943 n83 n74 a1 943 n83 n82 a1 943 n84 n81 a1 1245 n84 n86 a1 1245 n85 n72 a1 302 n85 n84 a1 302 n86 n87 a1 11999 n87 n88 a1 11999 n88 n89 a1 14505 $offdelim /; ** This is how the xupp values for double arcs were actually calculated * arc(m,n,s)$(fcost(m,n,s) or vcost(m,n,s) or xupp(m,n,s)) = yes; * loop(arc(m,n,'a2'), * xupp(m,n,'a1') = (fcost(m,n,'a2')-fcost(m,n,'a1'))/(vcost(m,n,'a1')-vcost(m,n,'a2')); * ); u = -sum(n$(demand(n)<0), demand(n)); usetree = 1;