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

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

Parameter fcost(m,n,s) /
$ondelim
   n9      n11     a1             33
   n11     n9      a1             23
   n9      n14     a1             36
   n14     n9      a1             45
   n13     n21     a1             34
   n21     n13     a1             44
   n10     n32     a1             42
   n32     n10     a1             41
   n33     n30     a1             21
   n30     n33     a1             37
   n15     n14     a1             25
   n14     n15     a1             41
   n34     n25     a1             46
   n25     n34     a1             20
   n7      n37     a1             29
   n37     n7      a1             38
   n18     n3      a1             38
   n3      n18     a1             44
   n7      n22     a1             21
   n22     n7      a1             24
   n29     n20     a1             28
   n20     n29     a1             21
   n21     n17     a1             37
   n17     n21     a1             20
   n12     n43     a1             34
   n43     n12     a1             24
   n35     n28     a1             44
   n28     n35     a1             41
   n35     n36     a1             34
   n36     n35     a1             21
   n2      n38     a1             43
   n38     n2      a1             47
   n21     n39     a1             43
   n39     n21     a1             38
   n30     n1      a1             41
   n1      n30     a1             23
   n31     n33     a1             48
   n33     n31     a1             25
   n32     n6      a1             23
   n6      n32     a1             20
   n12     n35     a1             47
   n35     n12     a1             40
   n12     n36     a1             31
   n36     n12     a1             49
   n20     n10     a1             29
   n10     n20     a1             48
   n7      n2      a1             48
   n2      n7      a1             32
   n24     n4      a1             38
   n4      n24     a1             25
   n9      n19     a1             49
   n19     n9      a1             24
   n10     n13     a1             28
   n13     n10     a1             37
   n14     n36     a1             47
   n36     n14     a1             40
   n26     n44     a1             38
   n44     n26     a1             47
   n35     n47     a1             45
   n47     n35     a1             30
   n20     n40     a1             43
   n40     n20     a1             27
   n8      n46     a1             25
   n46     n8      a1             46
   n35     n25     a1             32
   n25     n35     a1             37
   n25     n45     a1             30
   n45     n25     a1             26
   n20     n21     a1             48
   n21     n20     a1             47
   n13     n2      a1             38
   n2      n13     a1             42
   n1      n50     a1             39
   n50     n1      a1             21
   n29     n42     a1             22
   n42     n29     a1             27
   n32     n30     a1             29
   n30     n32     a1             43
   n7      n3      a1             45
   n3      n7      a1             47
   n17     n29     a1             30
   n29     n17     a1             23
   n29     n11     a1             32
   n11     n29     a1             36
   n3      n49     a1             20
   n49     n3      a1             20
   n15     n44     a1             31
   n44     n15     a1             21
   n16     n9      a1             34
   n9      n16     a1             25
   n28     n33     a1             26
   n33     n28     a1             35
   n5      n17     a1             45
   n17     n5      a1             22
   n20     n48     a1             27
   n48     n20     a1             47
   n15     n38     a1             25
   n38     n15     a1             46
   n7      n31     a1             48
   n31     n7      a1             21
   n20     n32     a1             37
   n32     n20     a1             46
   n19     n30     a1             44
   n30     n19     a1             22
   n8      n44     a1             48
   n44     n8      a1             25
   n34     n19     a1             23
   n19     n34     a1             46
   n26     n19     a1             43
   n19     n26     a1             48
   n3      n36     a1             25
   n36     n3      a1             21
   n19     n3      a1             25
   n3      n19     a1             32
   n12     n47     a1             45
   n47     n12     a1             35
   n5      n43     a1             41
   n43     n5      a1             49
   n16     n17     a1             47
   n17     n16     a1             27
   n34     n21     a1             24
   n21     n34     a1             48
   n2      n5      a1             27
   n5      n2      a1             26
   n33     n36     a1             22
   n36     n33     a1             32
   n16     n22     a1             33
   n22     n16     a1             45
   n27     n34     a1             24
   n34     n27     a1             29
   n16     n45     a1             23
   n45     n16     a1             22
   n13     n7      a1             40
   n7      n13     a1             21
   n9      n21     a1             24
   n21     n9      a1             26
   n16     n27     a1             44
   n27     n16     a1             44
   n18     n12     a1             24
   n12     n18     a1             21
   n20     n3      a1             26
   n3      n20     a1             42
   n24     n25     a1             25
   n25     n24     a1             28
   n22     n13     a1             25
   n13     n22     a1             20
   n33     n10     a1             35
   n10     n33     a1             25
   n18     n22     a1             45
   n22     n18     a1             21
   n17     n48     a1             36
   n48     n17     a1             42
   n26     n37     a1             43
   n37     n26     a1             22
   n20     n38     a1             30
   n38     n20     a1             20
   n34     n28     a1             26
   n28     n34     a1             44
   n25     n41     a1             35
   n41     n25     a1             27
   n20     n31     a1             49
   n31     n20     a1             45
   n21     n38     a1             34
   n38     n21     a1             37
   n11     n25     a1             44
   n25     n11     a1             22
   n11     n18     a1             27
   n18     n11     a1             37
   n17     n42     a1             44
   n42     n17     a1             38
   n24     n29     a1             43
   n29     n24     a1             49
   n27     n46     a1             43
   n46     n27     a1             21
   n17     n36     a1             33
   n36     n17     a1             43
   n22     n38     a1             43
   n38     n22     a1             28
   n19     n20     a1             32
   n20     n19     a1             45
   n6      n3      a1             27
   n3      n6      a1             26
   n11     n41     a1             26
   n41     n11     a1             43
   n14     n35     a1             46
   n35     n14     a1             49
   n9      n31     a1             42
   n31     n9      a1             32
   n12     n27     a1             40
   n27     n12     a1             30
   n8      n9      a1             40
   n9      n8      a1             24
   n35     n3      a1             35
   n3      n35     a1             22
   n33     n8      a1             31
   n8      n33     a1             46
   n20     n47     a1             28
   n47     n20     a1             20
   n33     n24     a1             23
   n24     n33     a1             32
   n2      n10     a1             43
   n10     n2      a1             43
   n4      n41     a1             42
   n41     n4      a1             32
   n13     n48     a1             45
   n48     n13     a1             39
   n2      n28     a1             36
   n28     n2      a1             29
   n5      n9      a1             42
   n9      n5      a1             33
   n17     n44     a1             39
   n44     n17     a1             44
   n3      n41     a1             32
   n41     n3      a1             24
   n2      n25     a1             33
   n25     n2      a1             24
   n32     n34     a1             35
   n34     n32     a1             42
   n30     n36     a1             41
   n36     n30     a1             37
   n21     n3      a1             42
   n3      n21     a1             48
   n27     n26     a1             36
   n26     n27     a1             41
   n29     n8      a1             25
   n8      n29     a1             29
   n7      n36     a1             44
   n36     n7      a1             44
   n12     n32     a1             35
   n32     n12     a1             37
   n14     n12     a1             42
   n12     n14     a1             25
   n21     n10     a1             35
   n10     n21     a1             35
   n9      n26     a1             22
   n26     n9      a1             27
   n32     n37     a1             49
   n37     n32     a1             35
   n28     n21     a1             45
   n21     n28     a1             22
   n12     n39     a1             48
   n39     n12     a1             46
   n16     n18     a1             49
   n18     n16     a1             38
   n17     n18     a1             33
   n18     n17     a1             31
   n8      n36     a1             20
   n36     n8      a1             25
   n28     n17     a1             45
   n17     n28     a1             30
   n25     n42     a1             41
   n42     n25     a1             45
   n33     n11     a1             45
   n11     n33     a1             44
   n11     n46     a1             29
   n46     n11     a1             26
   n23     n21     a1             34
   n21     n23     a1             35
   n3      n39     a1             41
   n39     n3      a1             25
   n7      n43     a1             37
   n43     n7      a1             27
   n16     n28     a1             26
   n28     n16     a1             20
   n19     n12     a1             22
   n12     n19     a1             25
   n2      n31     a1             44
   n31     n2      a1             41
   n34     n47     a1             29
   n47     n34     a1             39
   n29     n15     a1             44
   n15     n29     a1             27
   n4      n37     a1             46
   n37     n4      a1             45
   n13     n19     a1             34
   n19     n13     a1             27
   n26     n23     a1             36
   n23     n26     a1             25
   n6      n43     a1             30
   n43     n6      a1             39
   n12     n28     a1             40
   n28     n12     a1             33
   n31     n13     a1             46
   n13     n31     a1             27
   n19     n25     a1             23
   n25     n19     a1             30
   n12     n50     a1             41
   n50     n12     a1             23
   n17     n40     a1             40
   n40     n17     a1             21
   n30     n25     a1             43
   n25     n30     a1             39
   n6      n13     a1             31
   n13     n6      a1             25
   n7      n21     a1             40
   n21     n7      a1             24
   n23     n41     a1             35
   n41     n23     a1             46
   n16     n42     a1             28
   n42     n16     a1             28
   n31     n8      a1             21
   n8      n31     a1             22
   n34     n20     a1             23
   n20     n34     a1             49
   n11     n4      a1             23
   n4      n11     a1             47
   n32     n43     a1             28
   n43     n32     a1             34
   n15     n26     a1             34
   n26     n15     a1             47
   n7      n33     a1             36
   n33     n7      a1             48
   n28     n47     a1             26
   n47     n28     a1             30
   n21     n8      a1             22
   n8      n21     a1             27
   n25     n36     a1             40
   n36     n25     a1             34
   n32     n5      a1             27
   n5      n32     a1             30
   n4      n22     a1             44
   n22     n4      a1             47
   n35     n32     a1             23
   n32     n35     a1             45
   n16     n7      a1             22
   n7      n16     a1             29
   n7      n38     a1             33
   n38     n7      a1             21
   n20     n30     a1             30
   n30     n20     a1             24
   n18     n19     a1             28
   n19     n18     a1             37
   n3      n50     a1             21
   n50     n3      a1             27
   n18     n15     a1             39
   n15     n18     a1             45
   n16     n34     a1             27
   n34     n16     a1             28
   n8      n35     a1             33
   n35     n8      a1             23
   n4      n45     a1             27
   n45     n4      a1             34
   n25     n20     a1             45
   n20     n25     a1             40
   n26     n39     a1             21
   n39     n26     a1             34
   n21     n5      a1             44
   n5      n21     a1             45
   n35     n23     a1             32
   n23     n35     a1             37
   n13     n50     a1             30
   n50     n13     a1             38
   n10     n50     a1             22
   n50     n10     a1             25
   n6      n16     a1             36
   n16     n6      a1             20
   n15     n1      a1             22
   n1      n15     a1             23
   n24     n17     a1             22
   n17     n24     a1             28
$offdelim
/;

Parameter vcost(m,n,s) /
$ondelim
   n9      n11     a1              3
   n11     n9      a1              1
   n9      n14     a1              6
   n14     n9      a1              6
   n13     n21     a1              8
   n21     n13     a1              8
   n10     n32     a1              7
   n32     n10     a1              1
   n33     n30     a1              8
   n30     n33     a1              4
   n15     n14     a1              7
   n14     n15     a1              2
   n34     n25     a1              1
   n25     n34     a1              1
   n7      n37     a1              6
   n37     n7      a1              8
   n18     n3      a1              9
   n3      n18     a1              2
   n7      n22     a1              3
   n22     n7      a1              4
   n29     n20     a1              8
   n20     n29     a1              1
   n21     n17     a1              3
   n17     n21     a1              9
   n12     n43     a1              6
   n43     n12     a1              9
   n35     n28     a1              8
   n28     n35     a1              1
   n35     n36     a1              3
   n36     n35     a1              6
   n2      n38     a1              6
   n38     n2      a1              7
   n21     n39     a1              4
   n39     n21     a1              6
   n30     n1      a1              7
   n1      n30     a1              3
   n31     n33     a1              7
   n33     n31     a1              9
   n32     n6      a1              5
   n6      n32     a1              1
   n12     n35     a1              4
   n35     n12     a1              7
   n12     n36     a1              2
   n36     n12     a1              6
   n20     n10     a1              4
   n10     n20     a1              5
   n7      n2      a1              6
   n2      n7      a1              2
   n24     n4      a1              8
   n4      n24     a1              7
   n9      n19     a1              6
   n19     n9      a1              8
   n10     n13     a1              3
   n13     n10     a1              9
   n14     n36     a1              6
   n36     n14     a1              6
   n26     n44     a1              6
   n44     n26     a1              8
   n35     n47     a1              2
   n47     n35     a1              5
   n20     n40     a1              2
   n40     n20     a1              8
   n8      n46     a1              8
   n46     n8      a1              8
   n35     n25     a1              9
   n25     n35     a1              7
   n25     n45     a1              9
   n45     n25     a1              7
   n20     n21     a1              4
   n21     n20     a1              7
   n13     n2      a1              9
   n2      n13     a1              2
   n1      n50     a1              9
   n50     n1      a1              4
   n29     n42     a1              6
   n42     n29     a1              5
   n32     n30     a1              3
   n30     n32     a1              8
   n7      n3      a1              7
   n3      n7      a1              2
   n17     n29     a1              8
   n29     n17     a1              4
   n29     n11     a1              7
   n11     n29     a1              1
   n3      n49     a1              4
   n49     n3      a1              7
   n15     n44     a1              4
   n44     n15     a1              1
   n16     n9      a1              5
   n9      n16     a1              4
   n28     n33     a1              1
   n33     n28     a1              7
   n5      n17     a1              2
   n17     n5      a1              8
   n20     n48     a1              8
   n48     n20     a1              3
   n15     n38     a1              7
   n38     n15     a1              3
   n7      n31     a1              1
   n31     n7      a1              6
   n20     n32     a1              7
   n32     n20     a1              9
   n19     n30     a1              5
   n30     n19     a1              1
   n8      n44     a1              2
   n44     n8      a1              8
   n34     n19     a1              5
   n19     n34     a1              5
   n26     n19     a1              6
   n19     n26     a1              3
   n3      n36     a1              3
   n36     n3      a1              6
   n19     n3      a1              6
   n3      n19     a1              9
   n12     n47     a1              3
   n47     n12     a1              3
   n5      n43     a1              3
   n43     n5      a1              1
   n16     n17     a1              3
   n17     n16     a1              2
   n34     n21     a1              8
   n21     n34     a1              7
   n2      n5      a1              2
   n5      n2      a1              1
   n33     n36     a1              5
   n36     n33     a1              2
   n16     n22     a1              1
   n22     n16     a1              8
   n27     n34     a1              7
   n34     n27     a1              3
   n16     n45     a1              5
   n45     n16     a1              4
   n13     n7      a1              9
   n7      n13     a1              3
   n9      n21     a1              2
   n21     n9      a1              1
   n16     n27     a1              5
   n27     n16     a1              1
   n18     n12     a1              4
   n12     n18     a1              5
   n20     n3      a1              8
   n3      n20     a1              1
   n24     n25     a1              3
   n25     n24     a1              2
   n22     n13     a1              2
   n13     n22     a1              3
   n33     n10     a1              4
   n10     n33     a1              6
   n18     n22     a1              6
   n22     n18     a1              6
   n17     n48     a1              2
   n48     n17     a1              2
   n26     n37     a1              8
   n37     n26     a1              3
   n20     n38     a1              7
   n38     n20     a1              1
   n34     n28     a1              7
   n28     n34     a1              1
   n25     n41     a1              5
   n41     n25     a1              6
   n20     n31     a1              1
   n31     n20     a1              6
   n21     n38     a1              4
   n38     n21     a1              3
   n11     n25     a1              5
   n25     n11     a1              4
   n11     n18     a1              3
   n18     n11     a1              3
   n17     n42     a1              9
   n42     n17     a1              4
   n24     n29     a1              1
   n29     n24     a1              6
   n27     n46     a1              1
   n46     n27     a1              8
   n17     n36     a1              7
   n36     n17     a1              4
   n22     n38     a1              1
   n38     n22     a1              2
   n19     n20     a1              3
   n20     n19     a1              1
   n6      n3      a1              1
   n3      n6      a1              1
   n11     n41     a1              7
   n41     n11     a1              6
   n14     n35     a1              4
   n35     n14     a1              4
   n9      n31     a1              6
   n31     n9      a1              5
   n12     n27     a1              3
   n27     n12     a1              6
   n8      n9      a1              1
   n9      n8      a1              4
   n35     n3      a1              3
   n3      n35     a1              3
   n33     n8      a1              9
   n8      n33     a1              9
   n20     n47     a1              4
   n47     n20     a1              2
   n33     n24     a1              1
   n24     n33     a1              3
   n2      n10     a1              8
   n10     n2      a1              2
   n4      n41     a1              6
   n41     n4      a1              6
   n13     n48     a1              1
   n48     n13     a1              9
   n2      n28     a1              3
   n28     n2      a1              7
   n5      n9      a1              7
   n9      n5      a1              5
   n17     n44     a1              9
   n44     n17     a1              9
   n3      n41     a1              1
   n41     n3      a1              5
   n2      n25     a1              5
   n25     n2      a1              7
   n32     n34     a1              1
   n34     n32     a1              3
   n30     n36     a1              5
   n36     n30     a1              2
   n21     n3      a1              6
   n3      n21     a1              6
   n27     n26     a1              5
   n26     n27     a1              6
   n29     n8      a1              6
   n8      n29     a1              4
   n7      n36     a1              3
   n36     n7      a1              4
   n12     n32     a1              3
   n32     n12     a1              6
   n14     n12     a1              4
   n12     n14     a1              1
   n21     n10     a1              6
   n10     n21     a1              3
   n9      n26     a1              4
   n26     n9      a1              5
   n32     n37     a1              3
   n37     n32     a1              7
   n28     n21     a1              3
   n21     n28     a1              6
   n12     n39     a1              4
   n39     n12     a1              8
   n16     n18     a1              1
   n18     n16     a1              2
   n17     n18     a1              1
   n18     n17     a1              9
   n8      n36     a1              7
   n36     n8      a1              6
   n28     n17     a1              4
   n17     n28     a1              8
   n25     n42     a1              2
   n42     n25     a1              8
   n33     n11     a1              1
   n11     n33     a1              8
   n11     n46     a1              2
   n46     n11     a1              1
   n23     n21     a1              9
   n21     n23     a1              1
   n3      n39     a1              7
   n39     n3      a1              7
   n7      n43     a1              5
   n43     n7      a1              7
   n16     n28     a1              4
   n28     n16     a1              5
   n19     n12     a1              5
   n12     n19     a1              3
   n2      n31     a1              5
   n31     n2      a1              2
   n34     n47     a1              9
   n47     n34     a1              8
   n29     n15     a1              9
   n15     n29     a1              9
   n4      n37     a1              8
   n37     n4      a1              5
   n13     n19     a1              9
   n19     n13     a1              7
   n26     n23     a1              1
   n23     n26     a1              3
   n6      n43     a1              6
   n43     n6      a1              3
   n12     n28     a1              2
   n28     n12     a1              6
   n31     n13     a1              4
   n13     n31     a1              9
   n19     n25     a1              2
   n25     n19     a1              1
   n12     n50     a1              8
   n50     n12     a1              7
   n17     n40     a1              3
   n40     n17     a1              4
   n30     n25     a1              1
   n25     n30     a1              6
   n6      n13     a1              7
   n13     n6      a1              6
   n7      n21     a1              9
   n21     n7      a1              9
   n23     n41     a1              9
   n41     n23     a1              2
   n16     n42     a1              3
   n42     n16     a1              9
   n31     n8      a1              9
   n8      n31     a1              7
   n34     n20     a1              8
   n20     n34     a1              5
   n11     n4      a1              7
   n4      n11     a1              3
   n32     n43     a1              1
   n43     n32     a1              9
   n15     n26     a1              4
   n26     n15     a1              9
   n7      n33     a1              3
   n33     n7      a1              3
   n28     n47     a1              7
   n47     n28     a1              8
   n21     n8      a1              2
   n8      n21     a1              2
   n25     n36     a1              3
   n36     n25     a1              9
   n32     n5      a1              2
   n5      n32     a1              6
   n4      n22     a1              3
   n22     n4      a1              6
   n35     n32     a1              8
   n32     n35     a1              3
   n16     n7      a1              8
   n7      n16     a1              1
   n7      n38     a1              4
   n38     n7      a1              2
   n20     n30     a1              4
   n30     n20     a1              6
   n18     n19     a1              4
   n19     n18     a1              2
   n3      n50     a1              8
   n50     n3      a1              2
   n18     n15     a1              1
   n15     n18     a1              6
   n16     n34     a1              2
   n34     n16     a1              4
   n8      n35     a1              4
   n35     n8      a1              3
   n4      n45     a1              9
   n45     n4      a1              8
   n25     n20     a1              6
   n20     n25     a1              2
   n26     n39     a1              1
   n39     n26     a1              7
   n21     n5      a1              4
   n5      n21     a1              3
   n35     n23     a1              6
   n23     n35     a1              1
   n13     n50     a1              3
   n50     n13     a1              1
   n10     n50     a1              5
   n50     n10     a1              9
   n6      n16     a1              3
   n16     n6      a1              5
   n15     n1      a1              3
   n1      n15     a1              7
   n24     n17     a1              5
   n17     n24     a1              2
$offdelim
/;

Parameter demand(n) /
   n1             -1
   n2             14
   n3             -3
   n4            -16
   n5             -1
   n6              5
   n7             -4
   n8             -1
   n9             -2
   n12            -2
   n13             5
   n14             8
   n15            -3
   n16            -4
   n17            -7
   n18            -5
   n19            -1
   n20            -6
   n21            -2
   n23            -5
   n24            -4
   n25             4
   n26            -4
   n27            -2
   n28            -1
   n29            -1
   n30            -3
   n31            10
   n33            -1
   n34            -5
   n35            -2
   n36            -4
   n37             6
   n38             6
   n40             3
   n41             6
   n42            -2
   n43             5
   n44            -2
   n45            11
   n46             7
   n47            -1
   n48             4
   n49            -5
   n50             6
/;

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;
);