Set n / n1*n121 /;
Set s / a1 /;

yupp('n1','n1','a1')=0;

Parameter fcost(m,n,s) /
$ondelim
   n1      n2      a1              0
   n16     n1      a1             23
   n16     n17     a1              0
   n31     n16     a1             37
   n31     n32     a1              0
   n46     n31     a1             32
   n46     n47     a1              0
   n61     n46     a1             37
   n61     n62     a1              0
   n76     n61     a1             25
   n76     n77     a1              0
   n91     n76     a1             23
   n91     n92     a1              0
   n106    n91     a1             46
   n106    n107    a1              0
   n121    n106    a1             41
   n2      n3      a1              0
   n17     n2      a1             41
   n17     n18     a1              0
   n32     n17     a1             46
   n32     n33     a1              0
   n47     n32     a1             34
   n47     n48     a1              0
   n62     n47     a1             37
   n62     n63     a1              0
   n77     n62     a1             41
   n77     n78     a1              0
   n92     n77     a1             28
   n92     n93     a1              0
   n107    n92     a1             27
   n107    n108    a1              0
   n121    n107    a1             43
   n3      n4      a1              0
   n18     n3      a1             23
   n18     n19     a1              0
   n33     n18     a1             43
   n33     n34     a1              0
   n48     n33     a1             32
   n48     n49     a1              0
   n63     n48     a1             45
   n63     n64     a1              0
   n78     n63     a1             21
   n78     n79     a1              0
   n93     n78     a1             29
   n93     n94     a1              0
   n108    n93     a1             49
   n108    n109    a1              0
   n121    n108    a1             30
   n4      n5      a1              0
   n19     n4      a1             24
   n19     n20     a1              0
   n34     n19     a1             45
   n34     n35     a1              0
   n49     n34     a1             43
   n49     n50     a1              0
   n64     n49     a1             26
   n64     n65     a1              0
   n79     n64     a1             21
   n79     n80     a1              0
   n94     n79     a1             31
   n94     n95     a1              0
   n109    n94     a1             29
   n109    n110    a1              0
   n121    n109    a1             33
   n5      n6      a1              0
   n20     n5      a1             28
   n20     n21     a1              0
   n35     n20     a1             41
   n35     n36     a1              0
   n50     n35     a1             32
   n50     n51     a1              0
   n65     n50     a1             43
   n65     n66     a1              0
   n80     n65     a1             21
   n80     n81     a1              0
   n95     n80     a1             24
   n95     n96     a1              0
   n110    n95     a1             48
   n110    n111    a1              0
   n121    n110    a1             29
   n6      n7      a1              0
   n21     n6      a1             46
   n21     n22     a1              0
   n36     n21     a1             20
   n36     n37     a1              0
   n51     n36     a1             26
   n51     n52     a1              0
   n66     n51     a1             22
   n66     n67     a1              0
   n81     n66     a1             32
   n81     n82     a1              0
   n96     n81     a1             30
   n96     n97     a1              0
   n111    n96     a1             45
   n111    n112    a1              0
   n121    n111    a1             35
   n7      n8      a1              0
   n22     n7      a1             26
   n22     n23     a1              0
   n37     n22     a1             37
   n37     n38     a1              0
   n52     n37     a1             23
   n52     n53     a1              0
   n67     n52     a1             36
   n67     n68     a1              0
   n82     n67     a1             43
   n82     n83     a1              0
   n97     n82     a1             47
   n97     n98     a1              0
   n112    n97     a1             36
   n112    n113    a1              0
   n121    n112    a1             36
   n8      n9      a1              0
   n23     n8      a1             38
   n23     n24     a1              0
   n38     n23     a1             37
   n38     n39     a1              0
   n53     n38     a1             22
   n53     n54     a1              0
   n68     n53     a1             43
   n68     n69     a1              0
   n83     n68     a1             38
   n83     n84     a1              0
   n98     n83     a1             49
   n98     n99     a1              0
   n113    n98     a1             25
   n113    n114    a1              0
   n121    n113    a1             45
   n9      n10     a1              0
   n24     n9      a1             33
   n24     n25     a1              0
   n39     n24     a1             26
   n39     n40     a1              0
   n54     n39     a1             23
   n54     n55     a1              0
   n69     n54     a1             29
   n69     n70     a1              0
   n84     n69     a1             47
   n84     n85     a1              0
   n99     n84     a1             27
   n99     n100    a1              0
   n114    n99     a1             31
   n114    n115    a1              0
   n121    n114    a1             34
   n10     n11     a1              0
   n25     n10     a1             27
   n25     n26     a1              0
   n40     n25     a1             35
   n40     n41     a1              0
   n55     n40     a1             22
   n55     n56     a1              0
   n70     n55     a1             23
   n70     n71     a1              0
   n85     n70     a1             22
   n85     n86     a1              0
   n100    n85     a1             24
   n100    n101    a1              0
   n115    n100    a1             39
   n115    n116    a1              0
   n121    n115    a1             29
   n11     n12     a1              0
   n26     n11     a1             44
   n26     n27     a1              0
   n41     n26     a1             45
   n41     n42     a1              0
   n56     n41     a1             30
   n56     n57     a1              0
   n71     n56     a1             28
   n71     n72     a1              0
   n86     n71     a1             31
   n86     n87     a1              0
   n101    n86     a1             46
   n101    n102    a1              0
   n116    n101    a1             21
   n116    n117    a1              0
   n121    n116    a1             30
   n12     n13     a1              0
   n27     n12     a1             48
   n27     n28     a1              0
   n42     n27     a1             21
   n42     n43     a1              0
   n57     n42     a1             44
   n57     n58     a1              0
   n72     n57     a1             37
   n72     n73     a1              0
   n87     n72     a1             25
   n87     n88     a1              0
   n102    n87     a1             35
   n102    n103    a1              0
   n117    n102    a1             48
   n117    n118    a1              0
   n121    n117    a1             25
   n13     n14     a1              0
   n28     n13     a1             34
   n28     n29     a1              0
   n43     n28     a1             43
   n43     n44     a1              0
   n58     n43     a1             41
   n58     n59     a1              0
   n73     n58     a1             39
   n73     n74     a1              0
   n88     n73     a1             49
   n88     n89     a1              0
   n103    n88     a1             34
   n103    n104    a1              0
   n118    n103    a1             35
   n118    n119    a1              0
   n121    n118    a1             46
   n14     n15     a1              0
   n29     n14     a1             48
   n29     n30     a1              0
   n44     n29     a1             27
   n44     n45     a1              0
   n59     n44     a1             20
   n59     n60     a1              0
   n74     n59     a1             29
   n74     n75     a1              0
   n89     n74     a1             36
   n89     n90     a1              0
   n104    n89     a1             23
   n104    n105    a1              0
   n119    n104    a1             37
   n119    n120    a1              0
   n121    n119    a1             30
   n30     n15     a1             35
   n45     n30     a1             40
   n60     n45     a1             24
   n75     n60     a1             23
   n90     n75     a1             29
   n105    n90     a1             23
   n120    n105    a1             45
   n121    n120    a1             24
$offdelim
/;

Parameter vcost(m,n,s) /
$ondelim
   n1      n2      a1              9
   n16     n1      a1              8
   n16     n17     a1              2
   n31     n16     a1              2
   n31     n32     a1              7
   n46     n31     a1              4
   n46     n47     a1              1
   n61     n46     a1              5
   n61     n62     a1              4
   n76     n61     a1              1
   n76     n77     a1              3
   n91     n76     a1              7
   n91     n92     a1              1
   n106    n91     a1              7
   n106    n107    a1              2
   n121    n106    a1              3
   n2      n3      a1              1
   n17     n2      a1              4
   n17     n18     a1              6
   n32     n17     a1              2
   n32     n33     a1              6
   n47     n32     a1              3
   n47     n48     a1              8
   n62     n47     a1              3
   n62     n63     a1              2
   n77     n62     a1              8
   n77     n78     a1              7
   n92     n77     a1              3
   n92     n93     a1              1
   n107    n92     a1              4
   n107    n108    a1              6
   n121    n107    a1              1
   n3      n4      a1              3
   n18     n3      a1              4
   n18     n19     a1              3
   n33     n18     a1              5
   n33     n34     a1              7
   n48     n33     a1              5
   n48     n49     a1              8
   n63     n48     a1              9
   n63     n64     a1              8
   n78     n63     a1              5
   n78     n79     a1              7
   n93     n78     a1              5
   n93     n94     a1              4
   n108    n93     a1              9
   n108    n109    a1              5
   n121    n108    a1              5
   n4      n5      a1              3
   n19     n4      a1              7
   n19     n20     a1              7
   n34     n19     a1              8
   n34     n35     a1              5
   n49     n34     a1              4
   n49     n50     a1              4
   n64     n49     a1              2
   n64     n65     a1              5
   n79     n64     a1              8
   n79     n80     a1              5
   n94     n79     a1              7
   n94     n95     a1              3
   n109    n94     a1              3
   n109    n110    a1              4
   n121    n109    a1              8
   n5      n6      a1              9
   n20     n5      a1              1
   n20     n21     a1              8
   n35     n20     a1              1
   n35     n36     a1              7
   n50     n35     a1              4
   n50     n51     a1              3
   n65     n50     a1              6
   n65     n66     a1              9
   n80     n65     a1              2
   n80     n81     a1              4
   n95     n80     a1              3
   n95     n96     a1              4
   n110    n95     a1              8
   n110    n111    a1              5
   n121    n110    a1              8
   n6      n7      a1              6
   n21     n6      a1              6
   n21     n22     a1              4
   n36     n21     a1              9
   n36     n37     a1              8
   n51     n36     a1              9
   n51     n52     a1              2
   n66     n51     a1              2
   n66     n67     a1              5
   n81     n66     a1              6
   n81     n82     a1              6
   n96     n81     a1              6
   n96     n97     a1              7
   n111    n96     a1              6
   n111    n112    a1              2
   n121    n111    a1              3
   n7      n8      a1              2
   n22     n7      a1              5
   n22     n23     a1              9
   n37     n22     a1              2
   n37     n38     a1              4
   n52     n37     a1              3
   n52     n53     a1              5
   n67     n52     a1              3
   n67     n68     a1              2
   n82     n67     a1              4
   n82     n83     a1              2
   n97     n82     a1              9
   n97     n98     a1              4
   n112    n97     a1              5
   n112    n113    a1              5
   n121    n112    a1              3
   n8      n9      a1              6
   n23     n8      a1              5
   n23     n24     a1              1
   n38     n23     a1              8
   n38     n39     a1              2
   n53     n38     a1              5
   n53     n54     a1              5
   n68     n53     a1              6
   n68     n69     a1              6
   n83     n68     a1              9
   n83     n84     a1              7
   n98     n83     a1              9
   n98     n99     a1              8
   n113    n98     a1              6
   n113    n114    a1              4
   n121    n113    a1              4
   n9      n10     a1              5
   n24     n9      a1              5
   n24     n25     a1              4
   n39     n24     a1              6
   n39     n40     a1              6
   n54     n39     a1              5
   n54     n55     a1              7
   n69     n54     a1              2
   n69     n70     a1              5
   n84     n69     a1              1
   n84     n85     a1              4
   n99     n84     a1              3
   n99     n100    a1              1
   n114    n99     a1              5
   n114    n115    a1              6
   n121    n114    a1              2
   n10     n11     a1              5
   n25     n10     a1              4
   n25     n26     a1              9
   n40     n25     a1              8
   n40     n41     a1              8
   n55     n40     a1              8
   n55     n56     a1              4
   n70     n55     a1              1
   n70     n71     a1              9
   n85     n70     a1              8
   n85     n86     a1              2
   n100    n85     a1              3
   n100    n101    a1              4
   n115    n100    a1              3
   n115    n116    a1              7
   n121    n115    a1              4
   n11     n12     a1              3
   n26     n11     a1              9
   n26     n27     a1              4
   n41     n26     a1              1
   n41     n42     a1              7
   n56     n41     a1              4
   n56     n57     a1              7
   n71     n56     a1              9
   n71     n72     a1              3
   n86     n71     a1              2
   n86     n87     a1              8
   n101    n86     a1              2
   n101    n102    a1              6
   n116    n101    a1              9
   n116    n117    a1              7
   n121    n116    a1              9
   n12     n13     a1              6
   n27     n12     a1              6
   n27     n28     a1              2
   n42     n27     a1              4
   n42     n43     a1              7
   n57     n42     a1              7
   n57     n58     a1              9
   n72     n57     a1              9
   n72     n73     a1              9
   n87     n72     a1              8
   n87     n88     a1              6
   n102    n87     a1              6
   n102    n103    a1              3
   n117    n102    a1              2
   n117    n118    a1              4
   n121    n117    a1              6
   n13     n14     a1              4
   n28     n13     a1              1
   n28     n29     a1              1
   n43     n28     a1              8
   n43     n44     a1              5
   n58     n43     a1              5
   n58     n59     a1              8
   n73     n58     a1              1
   n73     n74     a1              6
   n88     n73     a1              7
   n88     n89     a1              3
   n103    n88     a1              8
   n103    n104    a1              3
   n118    n103    a1              5
   n118    n119    a1              3
   n121    n118    a1              8
   n14     n15     a1              6
   n29     n14     a1              2
   n29     n30     a1              7
   n44     n29     a1              5
   n44     n45     a1              6
   n59     n44     a1              2
   n59     n60     a1              1
   n74     n59     a1              7
   n74     n75     a1              3
   n89     n74     a1              3
   n89     n90     a1              6
   n104    n89     a1              7
   n104    n105    a1              8
   n119    n104    a1              3
   n119    n120    a1              7
   n121    n119    a1              9
   n30     n15     a1              3
   n45     n30     a1              9
   n60     n45     a1              9
   n75     n60     a1              2
   n90     n75     a1              5
   n105    n90     a1              1
   n120    n105    a1              5
   n121    n120    a1              3
$offdelim
/;

Parameter demand(n) /
   n1            322
   n2            142
   n3            395
   n4            112
   n5            364
   n6            410
   n7            485
   n8            126
   n9            488
   n10           185
   n11           478
   n12            50
   n13           476
   n14           365
   n15           126
   n121        -4524
/;

fcost(m,n,s) = 30 * fcost(m,n,s);
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;
);