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

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

Parameter fcost(m,n,s) /
$ondelim
   n1      n7      a1             88
   n7      n1      a1             86
   n1      n13     a1             88
   n13     n1      a1             90
   n1      n40     a1             51
   n40     n1      a1             54
   n1      n3      a1             58
   n3      n1      a1             79
   n1      n6      a1             97
   n6      n1      a1             67
   n1      n43     a1             71
   n43     n1      a1             53
   n1      n36     a1             56
   n36     n1      a1             53
   n1      n5      a1             61
   n5      n1      a1             65
   n1      n17     a1             64
   n17     n1      a1             77
   n1      n4      a1             62
   n4      n1      a1             72
   n1      n12     a1             89
   n12     n1      a1             97
   n1      n45     a1             80
   n45     n1      a1             92
   n1      n22     a1             72
   n22     n1      a1             69
   n1      n25     a1             77
   n25     n1      a1             89
   n1      n26     a1             51
   n26     n1      a1             58
   n1      n30     a1             69
   n30     n1      a1             51
   n1      n8      a1             74
   n8      n1      a1             56
   n1      n34     a1             68
   n34     n1      a1             87
   n1      n9      a1             70
   n9      n1      a1             76
   n1      n24     a1             93
   n24     n1      a1             64
   n1      n10     a1             72
   n10     n1      a1             93
   n2      n1      a1             64
   n1      n2      a1             94
   n2      n10     a1             71
   n10     n2      a1            100
   n2      n33     a1             77
   n33     n2      a1             87
   n2      n9      a1             89
   n9      n2      a1             99
   n2      n8      a1             50
   n8      n2      a1             91
   n2      n4      a1             79
   n4      n2      a1             84
   n2      n3      a1             85
   n3      n2      a1             55
   n2      n21     a1             65
   n21     n2      a1             63
   n2      n42     a1             68
   n42     n2      a1            100
   n2      n38     a1             86
   n38     n2      a1             86
   n2      n7      a1             92
   n7      n2      a1             69
   n3      n4      a1             61
   n4      n3      a1             62
   n3      n35     a1             61
   n35     n3      a1             97
   n3      n16     a1             55
   n16     n3      a1             57
   n3      n27     a1             96
   n27     n3      a1             80
   n3      n11     a1             52
   n11     n3      a1             64
   n3      n5      a1             68
   n5      n3      a1             55
   n3      n6      a1             73
   n6      n3      a1             52
   n3      n40     a1             63
   n40     n3      a1             92
   n3      n13     a1             97
   n13     n3      a1             50
   n3      n20     a1             84
   n20     n3      a1             58
   n3      n15     a1             74
   n15     n3      a1             76
   n3      n37     a1             57
   n37     n3      a1             73
   n3      n7      a1             55
   n7      n3      a1             79
   n3      n21     a1             66
   n21     n3      a1             87
   n4      n8      a1             92
   n8      n4      a1             99
   n4      n14     a1             56
   n14     n4      a1             97
   n4      n12     a1             58
   n12     n4      a1             51
   n4      n17     a1             79
   n17     n4      a1             66
   n4      n32     a1             63
   n32     n4      a1             87
   n4      n5      a1             77
   n5      n4      a1             56
   n4      n11     a1             96
   n11     n4      a1             81
   n4      n16     a1             71
   n16     n4      a1             85
   n4      n35     a1             53
   n35     n4      a1             61
   n5      n11     a1             97
   n11     n5      a1             92
   n5      n32     a1             76
   n32     n5      a1             86
   n5      n17     a1             76
   n17     n5      a1             97
   n5      n36     a1             70
   n36     n5      a1             74
   n5      n6      a1             73
   n6      n5      a1             85
   n6      n36     a1             57
   n36     n6      a1             76
   n6      n44     a1             92
   n44     n6      a1             78
   n6      n43     a1            100
   n43     n6      a1             78
   n7      n38     a1             51
   n38     n7      a1             61
   n7      n21     a1             60
   n21     n7      a1             66
   n7      n37     a1             82
   n37     n7      a1             76
   n7      n15     a1             60
   n15     n7      a1             65
   n7      n19     a1             76
   n19     n7      a1             67
   n7      n41     a1             65
   n41     n7      a1             78
   n7      n39     a1             92
   n39     n7      a1             70
   n7      n23     a1             94
   n23     n7      a1             76
   n7      n18     a1             78
   n18     n7      a1             66
   n7      n13     a1             55
   n13     n7      a1             81
   n8      n30     a1             93
   n30     n8      a1             69
   n8      n26     a1             52
   n26     n8      a1             80
   n8      n25     a1             58
   n25     n8      a1             92
   n8      n46     a1             78
   n46     n8      a1             64
   n8      n22     a1             72
   n22     n8      a1             85
   n8      n12     a1             53
   n12     n8      a1             95
   n8      n49     a1             57
   n49     n8      a1            100
   n8      n14     a1             62
   n14     n8      a1            100
   n8      n9      a1             67
   n9      n8      a1             75
   n8      n34     a1             51
   n34     n8      a1             82
   n9      n34     a1             54
   n34     n9      a1             72
   n9      n33     a1             61
   n33     n9      a1             55
   n9      n10     a1             88
   n10     n9      a1             74
   n9      n24     a1             54
   n24     n9      a1             96
   n10     n33     a1             87
   n33     n10     a1             75
   n10     n24     a1             89
   n24     n10     a1             76
   n11     n27     a1             91
   n27     n11     a1             96
   n11     n16     a1             58
   n16     n11     a1             92
   n12     n22     a1             93
   n22     n12     a1             59
   n12     n45     a1             91
   n45     n12     a1             94
   n12     n14     a1             71
   n14     n12     a1             90
   n12     n49     a1            100
   n49     n12     a1             80
   n13     n40     a1             73
   n40     n13     a1             56
   n13     n18     a1             78
   n18     n13     a1             93
   n13     n28     a1             50
   n28     n13     a1             68
   n13     n48     a1             79
   n48     n13     a1             89
   n13     n47     a1             54
   n47     n13     a1             55
   n13     n29     a1             89
   n29     n13     a1             53
   n13     n15     a1             92
   n15     n13     a1             56
   n13     n20     a1             86
   n20     n13     a1             86
   n14     n49     a1             56
   n49     n14     a1             98
   n15     n37     a1             56
   n37     n15     a1             80
   n15     n20     a1             50
   n20     n15     a1             99
   n15     n29     a1             58
   n29     n15     a1             58
   n15     n28     a1             64
   n28     n15     a1             58
   n15     n18     a1             56
   n18     n15     a1             83
   n15     n31     a1             92
   n31     n15     a1             64
   n15     n19     a1             82
   n19     n15     a1             82
   n16     n35     a1             85
   n35     n16     a1             56
   n16     n27     a1             99
   n27     n16     a1             52
   n17     n32     a1             76
   n32     n17     a1             62
   n18     n23     a1             58
   n23     n18     a1             63
   n18     n50     a1             62
   n50     n18     a1             60
   n18     n19     a1             77
   n19     n18     a1             54
   n18     n31     a1             67
   n31     n18     a1             71
   n18     n28     a1             93
   n28     n18     a1             76
   n19     n31     a1             98
   n31     n19     a1             52
   n19     n50     a1             79
   n50     n19     a1             92
   n19     n23     a1             98
   n23     n19     a1             87
   n19     n39     a1             97
   n39     n19     a1             54
   n19     n41     a1             99
   n41     n19     a1             92
   n21     n38     a1             69
   n38     n21     a1             85
   n21     n42     a1             50
   n42     n21     a1             75
   n22     n46     a1             71
   n46     n22     a1             72
   n22     n25     a1             58
   n25     n22     a1             59
   n22     n45     a1             90
   n45     n22     a1             88
   n23     n39     a1             55
   n39     n23     a1             54
   n23     n50     a1             78
   n50     n23     a1             88
   n25     n46     a1             56
   n46     n25     a1             82
   n25     n26     a1             64
   n26     n25     a1             65
   n26     n30     a1             98
   n30     n26     a1             60
   n28     n29     a1            100
   n29     n28     a1             53
   n28     n47     a1             83
   n47     n28     a1             57
   n28     n48     a1             56
   n48     n28     a1             88
   n29     n47     a1             83
   n47     n29     a1             90
   n36     n43     a1             89
   n43     n36     a1             73
   n36     n44     a1             83
   n44     n36     a1             60
   n38     n42     a1             56
   n42     n38     a1             80
   n39     n41     a1             60
   n41     n39     a1             96
   n43     n44     a1             70
   n44     n43     a1             85
   n47     n48     a1             81
   n48     n47     a1             85
$offdelim
/;

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

Parameter demand(n) /
   n1             -4
   n4             39
   n5             40
   n6             44
   n8             27
   n9            -21
   n11            30
   n12            33
   n13            32
   n15           -64
   n18           -59
   n24           -75
   n25            32
   n26           -47
   n27            32
   n31            24
   n33           -45
   n36            67
   n39           -73
   n40           -27
   n41            33
   n42            21
   n48           -85
   n49            20
   n50            26
/;

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