Set n / n1*n150 /; Set s / a1*a3 /; yupp('n1','n1','a1')=0; Parameter fcost(m,n,s) / $ondelim n2 n7 a1 5977 n2 n7 a2 7985 n7 n6 a1 5639 n7 n6 a2 5501 n6 n8 a1 5861 n8 n5 a1 7661 n8 n5 a2 7478 n6 n9 a1 4826 n9 n5 a1 7028 n2 n10 a1 5749 n2 n10 a2 7427 n10 n11 a1 5409 n11 n5 a1 6775 n10 n5 a1 5645 n2 n13 a1 7423 n13 n12 a1 7061 n12 n14 a1 6820 n14 n5 a1 4772 n12 n5 a1 5483 n2 n16 a1 6876 n16 n15 a1 4421 n2 n17 a1 7706 n17 n15 a1 4819 n2 n15 a1 6405 n2 n18 a1 7388 n18 n15 a1 7657 n2 n15 a2 5892 n2 n15 a3 5977 n15 n19 a1 7589 n15 n19 a2 4430 n19 n20 a1 5039 n20 n5 a1 7323 n19 n5 a1 5605 n15 n22 a1 4984 n22 n21 a1 5300 n21 n5 a1 7601 n21 n5 a2 5131 n5 n24 a1 6711 n5 n25 a1 7056 n25 n24 a1 5971 n24 n23 a1 5764 n24 n23 a2 5579 n5 n26 a1 4120 n26 n23 a1 5833 n5 n27 a1 4799 n27 n23 a1 5248 n5 n29 a1 4898 n29 n28 a1 5843 n29 n28 a2 4655 n5 n30 a1 5279 n30 n28 a1 4236 n28 n23 a1 6857 n28 n23 a2 4564 n28 n31 a1 4578 n31 n23 a1 5603 n23 n34 a1 7906 n23 n34 a2 5140 n34 n35 a1 6195 n34 n35 a2 7203 n35 n33 a1 6336 n33 n32 a1 4452 n33 n32 a2 7583 n33 n36 a1 7669 n36 n32 a1 7242 n33 n32 a3 6237 n32 n37 a1 6480 n37 n4 a1 4948 n37 n4 a2 5307 n32 n38 a1 6866 n38 n4 a1 6503 n32 n39 a1 4007 n39 n40 a1 4191 n40 n4 a1 5692 n39 n4 a1 4595 n4 n43 a1 6983 n4 n43 a2 4572 n43 n44 a1 5415 n43 n44 a2 4459 n44 n42 a1 5401 n4 n45 a1 5990 n45 n42 a1 7063 n4 n46 a1 6075 n46 n42 a1 4097 n42 n47 a1 5588 n42 n48 a1 6248 n48 n47 a1 7465 n47 n41 a1 4450 n47 n41 a2 6912 n42 n49 a1 7294 n49 n41 a1 4790 n42 n50 a1 4566 n50 n41 a1 5855 n42 n41 a1 5875 n4 n52 a1 6941 n4 n53 a1 4030 n53 n52 a1 6909 n4 n52 a2 4791 n4 n52 a3 6716 n52 n54 a2 6855 n52 n54 a1 4812 n54 n55 a1 5497 n55 n51 a1 4734 n54 n51 a1 6244 n51 n57 a1 4752 n57 n56 a1 4010 n51 n56 a1 4198 n51 n56 a2 4815 n56 n58 a1 4222 n58 n41 a1 5838 n56 n59 a1 4630 n59 n41 a1 5874 n56 n41 a1 5673 n41 n62 a1 4525 n41 n62 a2 4717 n62 n61 a1 5398 n62 n63 a1 7070 n63 n61 a1 6468 n61 n60 a1 6372 n61 n60 a2 4335 n61 n64 a1 7501 n64 n60 a1 6609 n64 n60 a2 6155 n60 n66 a1 6849 n60 n66 a2 6942 n66 n67 a1 5634 n67 n65 a1 6170 n66 n65 a1 4130 n65 n3 a1 5307 n65 n68 a1 5313 n65 n68 a2 6560 n68 n3 a1 5413 n41 n70 a1 7927 n41 n70 a2 4586 n70 n71 a1 7726 n70 n71 a2 7802 n71 n69 a1 6546 n41 n72 a1 5045 n72 n69 a1 4410 n41 n73 a1 6544 n73 n69 a1 7003 n69 n75 a1 4409 n69 n75 a2 7741 n75 n76 a1 5383 n76 n74 a1 4690 n76 n74 a2 6648 n74 n3 a1 5385 n74 n3 a2 4596 n74 n3 a3 6928 n74 n77 a1 7348 n77 n3 a1 7935 n3 n79 a1 6356 n79 n78 a1 7378 n78 n80 a1 6975 n78 n80 a2 7627 n80 n1 a1 7037 n3 n81 a1 7628 n3 n81 a2 6634 n81 n82 a1 5561 n82 n1 a1 4097 n82 n1 a2 7050 n3 n85 a1 6970 n3 n85 a2 5504 n85 n84 a1 7147 n84 n83 a1 4090 n84 n83 a2 4818 n83 n86 a1 6424 n83 n86 a2 4841 n86 n1 a1 6089 n83 n1 a1 4246 n3 n88 a1 7882 n3 n88 a2 5221 n88 n87 a1 6759 n88 n87 a2 4595 n87 n89 a1 4391 n89 n1 a1 5122 n89 n1 a2 4465 n87 n90 a1 6163 n90 n1 a1 6886 n3 n92 a1 6275 n3 n92 a2 4413 n92 n93 a1 4133 n93 n91 a1 7411 n92 n91 a1 5735 n91 n94 a1 5872 n94 n1 a1 7404 n91 n1 a1 5167 n91 n95 a1 6247 n95 n1 a1 5569 n3 n98 a1 4903 n98 n97 a1 6028 n98 n97 a2 6187 n97 n96 a1 6369 n97 n96 a2 4056 n3 n100 a1 5667 n100 n99 a1 6699 n99 n96 a1 6367 n99 n96 a2 4001 n3 n102 a1 4400 n102 n101 a1 7649 n101 n96 a1 5523 n101 n96 a2 4188 n3 n104 a1 7888 n3 n104 a2 6024 n104 n103 a1 4958 n103 n105 a1 6986 n105 n96 a1 5810 n105 n96 a2 4407 n96 n107 a1 4045 n96 n107 a2 6926 n96 n108 a1 6649 n108 n107 a1 7144 n107 n109 a1 6402 n107 n109 a2 5198 n109 n106 a1 4685 n107 n110 a1 7824 n107 n110 a2 7938 n110 n106 a1 6034 n106 n112 a1 5362 n112 n111 a1 5717 n112 n111 a2 6051 n106 n111 a1 5594 n111 n1 a1 5112 n111 n113 a1 6779 n113 n1 a1 7128 n111 n114 a1 5991 n114 n1 a1 7794 n3 n117 a1 6841 n117 n116 a1 7445 n116 n115 a1 5071 n116 n118 a1 5850 n118 n115 a1 4635 n3 n120 a1 4614 n120 n119 a1 6572 n120 n119 a2 4325 n119 n115 a1 6253 n119 n115 a2 6698 n3 n121 a1 5297 n3 n121 a2 7850 n121 n122 a1 4681 n122 n115 a1 4103 n3 n123 a1 7919 n3 n124 a1 4729 n124 n123 a1 4033 n123 n115 a1 7396 n123 n115 a2 4037 n115 n125 a1 4266 n115 n125 a2 4092 n115 n126 a1 4152 n126 n125 a1 4127 n125 n128 a1 5425 n128 n127 a1 6813 n128 n127 a2 5979 n127 n1 a1 6986 n127 n1 a2 4977 n115 n129 a1 4378 n115 n130 a1 6810 n130 n129 a1 5319 n129 n1 a1 7250 n115 n131 a1 6469 n131 n1 a1 7637 n131 n1 a2 4086 n115 n132 a1 7596 n132 n1 a1 6337 n3 n134 a1 4754 n3 n134 a2 4221 n3 n135 a1 7722 n135 n134 a1 7336 n134 n136 a1 4697 n136 n133 a1 6370 n136 n133 a2 6988 n134 n137 a1 5513 n137 n133 a1 6938 n3 n138 a1 7248 n138 n133 a1 7810 n3 n139 a1 5903 n139 n133 a1 5520 n3 n140 a1 7275 n140 n133 a1 7288 n3 n141 a1 4142 n3 n141 a2 6723 n141 n142 a1 6072 n142 n133 a1 5192 n133 n144 a1 7022 n144 n143 a1 5885 n133 n145 a1 4480 n133 n145 a2 4413 n145 n143 a1 7614 n143 n146 a1 4526 n146 n1 a1 5708 n143 n1 a1 6029 n143 n1 a2 4602 n133 n149 a1 6881 n149 n148 a1 6232 n133 n148 a1 7083 n148 n147 a1 5416 n148 n147 a2 4961 n147 n1 a1 5089 n147 n1 a2 7225 n147 n150 a1 6825 n150 n1 a1 7300 $offdelim /; Parameter vcost(m,n,s) / $ondelim n2 n7 a1 5 n2 n7 a2 5 n7 n6 a1 7 n7 n6 a2 7 n6 n8 a1 1 n8 n5 a1 9 n8 n5 a2 9 n6 n9 a1 8 n9 n5 a1 5 n2 n10 a1 2 n2 n10 a2 3 n10 n11 a1 3 n11 n5 a1 10 n10 n5 a1 10 n2 n13 a1 10 n13 n12 a1 4 n12 n14 a1 10 n14 n5 a1 10 n12 n5 a1 7 n2 n16 a1 5 n16 n15 a1 9 n2 n17 a1 3 n17 n15 a1 8 n2 n15 a1 7 n2 n18 a1 5 n18 n15 a1 6 n2 n15 a2 9 n2 n15 a3 3 n15 n19 a1 3 n15 n19 a2 9 n19 n20 a1 6 n20 n5 a1 7 n19 n5 a1 7 n15 n22 a1 4 n22 n21 a1 9 n21 n5 a1 3 n21 n5 a2 1 n5 n24 a1 7 n5 n25 a1 1 n25 n24 a1 4 n24 n23 a1 6 n24 n23 a2 10 n5 n26 a1 1 n26 n23 a1 6 n5 n27 a1 9 n27 n23 a1 2 n5 n29 a1 8 n29 n28 a1 5 n29 n28 a2 2 n5 n30 a1 5 n30 n28 a1 2 n28 n23 a1 3 n28 n23 a2 9 n28 n31 a1 8 n31 n23 a1 6 n23 n34 a1 9 n23 n34 a2 6 n34 n35 a1 10 n34 n35 a2 3 n35 n33 a1 10 n33 n32 a1 3 n33 n32 a2 4 n33 n36 a1 1 n36 n32 a1 10 n33 n32 a3 4 n32 n37 a1 10 n37 n4 a1 5 n37 n4 a2 9 n32 n38 a1 7 n38 n4 a1 1 n32 n39 a1 8 n39 n40 a1 9 n40 n4 a1 7 n39 n4 a1 8 n4 n43 a1 4 n4 n43 a2 7 n43 n44 a1 10 n43 n44 a2 6 n44 n42 a1 9 n4 n45 a1 6 n45 n42 a1 2 n4 n46 a1 3 n46 n42 a1 6 n42 n47 a1 5 n42 n48 a1 1 n48 n47 a1 3 n47 n41 a1 6 n47 n41 a2 10 n42 n49 a1 9 n49 n41 a1 3 n42 n50 a1 4 n50 n41 a1 2 n42 n41 a1 5 n4 n52 a1 4 n4 n53 a1 5 n53 n52 a1 4 n4 n52 a2 5 n4 n52 a3 5 n52 n54 a2 7 n52 n54 a1 7 n54 n55 a1 10 n55 n51 a1 8 n54 n51 a1 8 n51 n57 a1 8 n57 n56 a1 9 n51 n56 a1 6 n51 n56 a2 10 n56 n58 a1 5 n58 n41 a1 6 n56 n59 a1 2 n59 n41 a1 1 n56 n41 a1 7 n41 n62 a1 3 n41 n62 a2 9 n62 n61 a1 6 n62 n63 a1 9 n63 n61 a1 8 n61 n60 a1 10 n61 n60 a2 9 n61 n64 a1 9 n64 n60 a1 1 n64 n60 a2 7 n60 n66 a1 8 n60 n66 a2 5 n66 n67 a1 2 n67 n65 a1 1 n66 n65 a1 5 n65 n3 a1 8 n65 n68 a1 6 n65 n68 a2 9 n68 n3 a1 2 n41 n70 a1 8 n41 n70 a2 5 n70 n71 a1 10 n70 n71 a2 7 n71 n69 a1 3 n41 n72 a1 4 n72 n69 a1 8 n41 n73 a1 1 n73 n69 a1 4 n69 n75 a1 10 n69 n75 a2 4 n75 n76 a1 2 n76 n74 a1 7 n76 n74 a2 9 n74 n3 a1 1 n74 n3 a2 2 n74 n3 a3 5 n74 n77 a1 10 n77 n3 a1 10 n3 n79 a1 7 n79 n78 a1 2 n78 n80 a1 2 n78 n80 a2 6 n80 n1 a1 7 n3 n81 a1 9 n3 n81 a2 7 n81 n82 a1 4 n82 n1 a1 8 n82 n1 a2 3 n3 n85 a1 1 n3 n85 a2 6 n85 n84 a1 5 n84 n83 a1 4 n84 n83 a2 8 n83 n86 a1 5 n83 n86 a2 1 n86 n1 a1 6 n83 n1 a1 8 n3 n88 a1 10 n3 n88 a2 5 n88 n87 a1 1 n88 n87 a2 2 n87 n89 a1 4 n89 n1 a1 10 n89 n1 a2 4 n87 n90 a1 9 n90 n1 a1 8 n3 n92 a1 10 n3 n92 a2 7 n92 n93 a1 6 n93 n91 a1 10 n92 n91 a1 2 n91 n94 a1 7 n94 n1 a1 4 n91 n1 a1 5 n91 n95 a1 5 n95 n1 a1 8 n3 n98 a1 8 n98 n97 a1 6 n98 n97 a2 5 n97 n96 a1 4 n97 n96 a2 2 n3 n100 a1 9 n100 n99 a1 3 n99 n96 a1 1 n99 n96 a2 10 n3 n102 a1 8 n102 n101 a1 8 n101 n96 a1 9 n101 n96 a2 10 n3 n104 a1 6 n3 n104 a2 7 n104 n103 a1 7 n103 n105 a1 1 n105 n96 a1 10 n105 n96 a2 4 n96 n107 a1 6 n96 n107 a2 2 n96 n108 a1 10 n108 n107 a1 2 n107 n109 a1 4 n107 n109 a2 1 n109 n106 a1 4 n107 n110 a1 5 n107 n110 a2 10 n110 n106 a1 1 n106 n112 a1 2 n112 n111 a1 5 n112 n111 a2 2 n106 n111 a1 3 n111 n1 a1 9 n111 n113 a1 2 n113 n1 a1 10 n111 n114 a1 9 n114 n1 a1 6 n3 n117 a1 1 n117 n116 a1 1 n116 n115 a1 2 n116 n118 a1 1 n118 n115 a1 5 n3 n120 a1 2 n120 n119 a1 1 n120 n119 a2 5 n119 n115 a1 6 n119 n115 a2 9 n3 n121 a1 3 n3 n121 a2 7 n121 n122 a1 2 n122 n115 a1 9 n3 n123 a1 3 n3 n124 a1 8 n124 n123 a1 3 n123 n115 a1 9 n123 n115 a2 1 n115 n125 a1 1 n115 n125 a2 10 n115 n126 a1 5 n126 n125 a1 2 n125 n128 a1 6 n128 n127 a1 3 n128 n127 a2 2 n127 n1 a1 7 n127 n1 a2 1 n115 n129 a1 10 n115 n130 a1 5 n130 n129 a1 7 n129 n1 a1 5 n115 n131 a1 7 n131 n1 a1 7 n131 n1 a2 2 n115 n132 a1 2 n132 n1 a1 10 n3 n134 a1 4 n3 n134 a2 7 n3 n135 a1 7 n135 n134 a1 5 n134 n136 a1 4 n136 n133 a1 1 n136 n133 a2 5 n134 n137 a1 3 n137 n133 a1 8 n3 n138 a1 8 n138 n133 a1 8 n3 n139 a1 6 n139 n133 a1 8 n3 n140 a1 8 n140 n133 a1 7 n3 n141 a1 9 n3 n141 a2 6 n141 n142 a1 8 n142 n133 a1 5 n133 n144 a1 10 n144 n143 a1 6 n133 n145 a1 2 n133 n145 a2 6 n145 n143 a1 6 n143 n146 a1 5 n146 n1 a1 5 n143 n1 a1 3 n143 n1 a2 7 n133 n149 a1 9 n149 n148 a1 7 n133 n148 a1 7 n148 n147 a1 1 n148 n147 a2 4 n147 n1 a1 3 n147 n1 a2 2 n147 n150 a1 4 n150 n1 a1 8 $offdelim /; Parameter demand(n) / n1 2903 n3 491 n4 161 n5 41 n6 -1524 n7 -519 n8 16 n9 12 n22 -1 n46 -11 n49 -1 n63 -1 n66 -3 n69 -2 n74 1 n84 -3 n91 53 n92 -122 n93 4 n94 2 n95 2 n96 133 n99 -1049 n106 38 n107 6 n110 1 n113 -118 n115 90 n116 -11 n121 -421 n122 18 n125 9 n126 6 n133 51 n137 -2 n140 -56 n141 -206 n142 15 n143 8 n144 2 n147 12 n148 1 n149 -27 n150 1 /; * Upper bound for flow on any arc u = -sum(n$(demand(n)<0), demand(n)); xupp(m,n,s)$(fcost(m,n,s) or vcost(m,n,s)) = u;