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

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

Parameter fcost(m,n,s) /
$ondelim
   n1      n2      a1             63
   n2      n1      a1             77
   n1      n3      a1             62
   n3      n1      a1             67
   n1      n4      a1             75
   n4      n1      a1             83
   n1      n5      a1             66
   n5      n1      a1             61
   n1      n6      a1             81
   n6      n1      a1             86
   n1      n7      a1             94
   n7      n1      a1             67
   n1      n8      a1             61
   n8      n1      a1             88
   n1      n9      a1             67
   n9      n1      a1             53
   n1      n10     a1             53
   n10     n1      a1             53
   n1      n11     a1             72
   n11     n1      a1             95
   n1      n12     a1             51
   n12     n1      a1             59
   n1      n13     a1             79
   n13     n1      a1             90
   n1      n14     a1             64
   n14     n1      a1             75
   n1      n15     a1             93
   n15     n1      a1             86
   n1      n16     a1             71
   n16     n1      a1             61
   n1      n17     a1             79
   n17     n1      a1             93
   n1      n18     a1             88
   n18     n1      a1             71
   n1      n19     a1             82
   n19     n1      a1             53
   n1      n20     a1             51
   n20     n1      a1             84
   n2      n3      a1             98
   n3      n2      a1             99
   n2      n4      a1             66
   n4      n2      a1             60
   n2      n5      a1             76
   n5      n2      a1             82
   n2      n6      a1             52
   n6      n2      a1             80
   n2      n7      a1             85
   n7      n2      a1             65
   n2      n8      a1             76
   n8      n2      a1             97
   n2      n9      a1             73
   n9      n2      a1             56
   n2      n10     a1             63
   n10     n2      a1             67
   n2      n11     a1             76
   n11     n2      a1             66
   n2      n12     a1             78
   n12     n2      a1             87
   n2      n13     a1             52
   n13     n2      a1             83
   n2      n14     a1             78
   n14     n2      a1             59
   n2      n15     a1             75
   n15     n2      a1             85
   n2      n16     a1             63
   n16     n2      a1             76
   n2      n17     a1             53
   n17     n2      a1             79
   n2      n18     a1             57
   n18     n2      a1             57
   n2      n19     a1             81
   n19     n2      a1             74
   n2      n20     a1             97
   n20     n2      a1             65
   n3      n4      a1             92
   n4      n3      a1             63
   n3      n5      a1             82
   n5      n3      a1             74
   n3      n6      a1             89
   n6      n3      a1             89
   n3      n7      a1             64
   n7      n3      a1             85
   n3      n8      a1             50
   n8      n3      a1             68
   n3      n9      a1             81
   n9      n3      a1             86
   n3      n10     a1             61
   n10     n3      a1             90
   n3      n11     a1             58
   n11     n3      a1             51
   n3      n12     a1             64
   n12     n3      a1             57
   n3      n13     a1             55
   n13     n3      a1             75
   n3      n14     a1             88
   n14     n3      a1             75
   n3      n15     a1             94
   n15     n3      a1             63
   n3      n16     a1             61
   n16     n3      a1             69
   n3      n17     a1             89
   n17     n3      a1             64
   n3      n18     a1             85
   n18     n3      a1             56
   n3      n19     a1             88
   n19     n3      a1             57
   n3      n20     a1             90
   n20     n3      a1             59
   n4      n5      a1             96
   n5      n4      a1             63
   n4      n6      a1             77
   n6      n4      a1             90
   n4      n7      a1             79
   n7      n4      a1             64
   n4      n8      a1             62
   n8      n4      a1             93
   n4      n9      a1             83
   n9      n4      a1             92
   n4      n10     a1             64
   n10     n4      a1             76
   n4      n11     a1             78
   n11     n4      a1             65
   n4      n12     a1             81
   n12     n4      a1             99
   n4      n13     a1             50
   n13     n4      a1             83
   n4      n14     a1             74
   n14     n4      a1             92
   n4      n15     a1             85
   n15     n4      a1             84
   n4      n16     a1             84
   n16     n4      a1             89
   n4      n17     a1             53
   n17     n4      a1             63
   n4      n18     a1             85
   n18     n4      a1             55
   n4      n19     a1             53
   n19     n4      a1             72
   n4      n20     a1             94
   n20     n4      a1             55
   n5      n6      a1             58
   n6      n5      a1             65
   n5      n7      a1             80
   n7      n5      a1             87
   n5      n8      a1             81
   n8      n5      a1             79
   n5      n9      a1             85
   n9      n5      a1             99
   n5      n10     a1             76
   n10     n5      a1             80
   n5      n11     a1             69
   n11     n5      a1             99
   n5      n12     a1             88
   n12     n5      a1             76
   n5      n13     a1             97
   n13     n5      a1             96
   n5      n14     a1             71
   n14     n5      a1             89
   n5      n15     a1             76
   n15     n5      a1             52
   n5      n16     a1             53
   n16     n5      a1             82
   n5      n17     a1             68
   n17     n5      a1             62
   n5      n18     a1             56
   n18     n5      a1             65
   n5      n19     a1             93
   n19     n5      a1             59
   n5      n20     a1             94
   n20     n5      a1             70
   n6      n7      a1             63
   n7      n6      a1             70
   n6      n8      a1             93
   n8      n6      a1             90
   n6      n9      a1             59
   n9      n6      a1             97
   n6      n10     a1             97
   n10     n6      a1             62
   n6      n11     a1             52
   n11     n6      a1             98
   n6      n12     a1             96
   n12     n6      a1             96
   n6      n13     a1             68
   n13     n6      a1             61
   n6      n14     a1             54
   n14     n6      a1             97
   n6      n15     a1             88
   n15     n6      a1             91
   n6      n16     a1             77
   n16     n6      a1             67
   n6      n17     a1             63
   n17     n6      a1             65
   n6      n18     a1             96
   n18     n6      a1             61
   n6      n19     a1             79
   n19     n6      a1             88
   n6      n20     a1             66
   n20     n6      a1             52
   n7      n8      a1             76
   n8      n7      a1             82
   n7      n9      a1             55
   n9      n7      a1             55
   n7      n10     a1             80
   n10     n7      a1             72
   n7      n11     a1             55
   n11     n7      a1             96
   n7      n12     a1             54
   n12     n7      a1             59
   n7      n13     a1             59
   n13     n7      a1             97
   n7      n14     a1             67
   n14     n7      a1             97
   n7      n15     a1             61
   n15     n7      a1             67
   n7      n16     a1             91
   n16     n7      a1             68
   n7      n17     a1             70
   n17     n7      a1             95
   n7      n18     a1             59
   n18     n7      a1             69
   n7      n19     a1             70
   n19     n7      a1             77
   n7      n20     a1             90
   n20     n7      a1             51
   n8      n9      a1             74
   n9      n8      a1             98
   n8      n10     a1             61
   n10     n8      a1             62
   n8      n11     a1             83
   n11     n8      a1             84
   n8      n12     a1             77
   n12     n8      a1             61
   n8      n13     a1             84
   n13     n8      a1             69
   n8      n14     a1             99
   n14     n8      a1             84
   n8      n15     a1             66
   n15     n8      a1             60
   n8      n16     a1             68
   n16     n8      a1             73
   n8      n17     a1             73
   n17     n8      a1             66
   n8      n18     a1             91
   n18     n8      a1             78
   n8      n19     a1             68
   n19     n8      a1             59
   n8      n20     a1             55
   n20     n8      a1             59
   n9      n10     a1             94
   n10     n9      a1             89
   n9      n11     a1             81
   n11     n9      a1             97
   n9      n12     a1             52
   n12     n9      a1             79
   n9      n13     a1             52
   n13     n9      a1             60
   n9      n14     a1             95
   n14     n9      a1             85
   n9      n15     a1             58
   n15     n9      a1             53
   n9      n16     a1             68
   n16     n9      a1             79
   n9      n17     a1             58
   n17     n9      a1             62
   n9      n18     a1             99
   n18     n9      a1             74
   n9      n19     a1             63
   n19     n9      a1             57
   n9      n20     a1             55
   n20     n9      a1             60
   n10     n11     a1             54
   n11     n10     a1             82
   n10     n12     a1             68
   n12     n10     a1             71
   n10     n13     a1             55
   n13     n10     a1             98
   n10     n14     a1             61
   n14     n10     a1             52
   n10     n15     a1             99
   n15     n10     a1             59
   n10     n16     a1             50
   n16     n10     a1             79
   n10     n17     a1             85
   n17     n10     a1             76
   n10     n18     a1             72
   n18     n10     a1             50
   n10     n19     a1             63
   n19     n10     a1             61
   n10     n20     a1             66
   n20     n10     a1             99
   n11     n12     a1             71
   n12     n11     a1             90
   n11     n13     a1             66
   n13     n11     a1             65
   n11     n14     a1             50
   n14     n11     a1             62
   n11     n15     a1             60
   n15     n11     a1             94
   n11     n16     a1             89
   n16     n11     a1             70
   n11     n17     a1             88
   n17     n11     a1             88
   n11     n18     a1             74
   n18     n11     a1             55
   n11     n19     a1             72
   n19     n11     a1             97
   n11     n20     a1             97
   n20     n11     a1             55
   n12     n13     a1             67
   n13     n12     a1             98
   n12     n14     a1             52
   n14     n12     a1             76
   n12     n15     a1             70
   n15     n12     a1             84
   n12     n16     a1             87
   n16     n12     a1             74
   n12     n17     a1             73
   n17     n12     a1             70
   n12     n18     a1             97
   n18     n12     a1             53
   n12     n19     a1             68
   n19     n12     a1             66
   n12     n20     a1             57
   n20     n12     a1             76
   n13     n14     a1             65
   n14     n13     a1             64
   n13     n15     a1             78
   n15     n13     a1             93
   n13     n16     a1             77
   n16     n13     a1             58
   n13     n17     a1             61
   n17     n13     a1             90
   n13     n18     a1             73
   n18     n13     a1             75
   n13     n19     a1             92
   n19     n13     a1             87
   n13     n20     a1             62
   n20     n13     a1             66
   n14     n15     a1             58
   n15     n14     a1             60
   n14     n16     a1             70
   n16     n14     a1             95
   n14     n17     a1             98
   n17     n14     a1             54
   n14     n18     a1             74
   n18     n14     a1             92
   n14     n19     a1             79
   n19     n14     a1             64
   n14     n20     a1             63
   n20     n14     a1             52
   n15     n16     a1             61
   n16     n15     a1             60
   n15     n17     a1             67
   n17     n15     a1             67
   n15     n18     a1             93
   n18     n15     a1             77
   n15     n19     a1             85
   n19     n15     a1             54
   n15     n20     a1             53
   n20     n15     a1             81
   n16     n17     a1             89
   n17     n16     a1             96
   n16     n18     a1             88
   n18     n16     a1             98
   n16     n19     a1             51
   n19     n16     a1             69
   n16     n20     a1             80
   n20     n16     a1             81
   n17     n18     a1             56
   n18     n17     a1             80
   n17     n19     a1             85
   n19     n17     a1             98
   n17     n20     a1             58
   n20     n17     a1             62
   n18     n19     a1             91
   n19     n18     a1             68
   n18     n20     a1             73
   n20     n18     a1             58
   n19     n20     a1             52
   n20     n19     a1             80
$offdelim
/;

Parameter vcost(m,n,s) /
$ondelim
   n1      n2      a1              8
   n2      n1      a1              2
   n1      n3      a1              4
   n3      n1      a1              5
   n1      n4      a1              1
   n4      n1      a1              7
   n1      n5      a1              7
   n5      n1      a1              3
   n1      n6      a1              4
   n6      n1      a1              2
   n1      n7      a1              3
   n7      n1      a1              3
   n1      n8      a1              8
   n8      n1      a1              3
   n1      n9      a1              4
   n9      n1      a1              1
   n1      n10     a1              4
   n10     n1      a1              5
   n1      n11     a1              5
   n11     n1      a1              9
   n1      n12     a1              5
   n12     n1      a1              5
   n1      n13     a1              9
   n13     n1      a1              5
   n1      n14     a1              7
   n14     n1      a1              8
   n1      n15     a1              4
   n15     n1      a1              2
   n1      n16     a1              8
   n16     n1      a1              7
   n1      n17     a1              3
   n17     n1      a1              8
   n1      n18     a1              1
   n18     n1      a1              1
   n1      n19     a1              4
   n19     n1      a1              6
   n1      n20     a1              2
   n20     n1      a1              3
   n2      n3      a1              8
   n3      n2      a1              8
   n2      n4      a1              6
   n4      n2      a1              9
   n2      n5      a1              9
   n5      n2      a1              2
   n2      n6      a1              6
   n6      n2      a1              6
   n2      n7      a1              6
   n7      n2      a1              3
   n2      n8      a1              5
   n8      n2      a1              2
   n2      n9      a1              3
   n9      n2      a1              3
   n2      n10     a1              4
   n10     n2      a1              9
   n2      n11     a1              5
   n11     n2      a1              3
   n2      n12     a1              5
   n12     n2      a1              8
   n2      n13     a1              5
   n13     n2      a1              6
   n2      n14     a1              9
   n14     n2      a1              9
   n2      n15     a1              6
   n15     n2      a1              4
   n2      n16     a1              5
   n16     n2      a1              6
   n2      n17     a1              5
   n17     n2      a1              2
   n2      n18     a1              1
   n18     n2      a1              3
   n2      n19     a1              5
   n19     n2      a1              2
   n2      n20     a1              4
   n20     n2      a1              8
   n3      n4      a1              8
   n4      n3      a1              1
   n3      n5      a1              8
   n5      n3      a1              3
   n3      n6      a1              3
   n6      n3      a1              4
   n3      n7      a1              9
   n7      n3      a1              1
   n3      n8      a1              4
   n8      n3      a1              9
   n3      n9      a1              2
   n9      n3      a1              2
   n3      n10     a1              9
   n10     n3      a1              9
   n3      n11     a1              6
   n11     n3      a1              4
   n3      n12     a1              7
   n12     n3      a1              9
   n3      n13     a1              8
   n13     n3      a1              6
   n3      n14     a1              2
   n14     n3      a1              6
   n3      n15     a1              1
   n15     n3      a1              8
   n3      n16     a1              5
   n16     n3      a1              1
   n3      n17     a1              7
   n17     n3      a1              8
   n3      n18     a1              5
   n18     n3      a1              8
   n3      n19     a1              2
   n19     n3      a1              5
   n3      n20     a1              2
   n20     n3      a1              7
   n4      n5      a1              3
   n5      n4      a1              7
   n4      n6      a1              3
   n6      n4      a1              9
   n4      n7      a1              1
   n7      n4      a1              9
   n4      n8      a1              4
   n8      n4      a1              1
   n4      n9      a1              2
   n9      n4      a1              5
   n4      n10     a1              3
   n10     n4      a1              8
   n4      n11     a1              2
   n11     n4      a1              7
   n4      n12     a1              9
   n12     n4      a1              1
   n4      n13     a1              9
   n13     n4      a1              7
   n4      n14     a1              3
   n14     n4      a1              8
   n4      n15     a1              2
   n15     n4      a1              2
   n4      n16     a1              7
   n16     n4      a1              7
   n4      n17     a1              8
   n17     n4      a1              8
   n4      n18     a1              9
   n18     n4      a1              3
   n4      n19     a1              1
   n19     n4      a1              8
   n4      n20     a1              1
   n20     n4      a1              6
   n5      n6      a1              5
   n6      n5      a1              6
   n5      n7      a1              9
   n7      n5      a1              2
   n5      n8      a1              6
   n8      n5      a1              9
   n5      n9      a1              6
   n9      n5      a1              9
   n5      n10     a1              6
   n10     n5      a1              7
   n5      n11     a1              7
   n11     n5      a1              4
   n5      n12     a1              9
   n12     n5      a1              8
   n5      n13     a1              2
   n13     n5      a1              4
   n5      n14     a1              8
   n14     n5      a1              7
   n5      n15     a1              9
   n15     n5      a1              2
   n5      n16     a1              9
   n16     n5      a1              9
   n5      n17     a1              5
   n17     n5      a1              3
   n5      n18     a1              5
   n18     n5      a1              7
   n5      n19     a1              9
   n19     n5      a1              9
   n5      n20     a1              2
   n20     n5      a1              7
   n6      n7      a1              1
   n7      n6      a1              4
   n6      n8      a1              6
   n8      n6      a1              2
   n6      n9      a1              8
   n9      n6      a1              1
   n6      n10     a1              8
   n10     n6      a1              8
   n6      n11     a1              9
   n11     n6      a1              7
   n6      n12     a1              9
   n12     n6      a1              8
   n6      n13     a1              6
   n13     n6      a1              4
   n6      n14     a1              9
   n14     n6      a1              6
   n6      n15     a1              2
   n15     n6      a1              3
   n6      n16     a1              7
   n16     n6      a1              3
   n6      n17     a1              9
   n17     n6      a1              2
   n6      n18     a1              7
   n18     n6      a1              5
   n6      n19     a1              6
   n19     n6      a1              6
   n6      n20     a1              5
   n20     n6      a1              3
   n7      n8      a1              6
   n8      n7      a1              4
   n7      n9      a1              1
   n9      n7      a1              1
   n7      n10     a1              8
   n10     n7      a1              3
   n7      n11     a1              1
   n11     n7      a1              7
   n7      n12     a1              4
   n12     n7      a1              9
   n7      n13     a1              5
   n13     n7      a1              2
   n7      n14     a1              3
   n14     n7      a1              6
   n7      n15     a1              3
   n15     n7      a1              3
   n7      n16     a1              5
   n16     n7      a1              7
   n7      n17     a1              3
   n17     n7      a1              9
   n7      n18     a1              8
   n18     n7      a1              1
   n7      n19     a1              2
   n19     n7      a1              8
   n7      n20     a1              1
   n20     n7      a1              9
   n8      n9      a1              5
   n9      n8      a1              7
   n8      n10     a1              8
   n10     n8      a1              8
   n8      n11     a1              4
   n11     n8      a1              4
   n8      n12     a1              8
   n12     n8      a1              6
   n8      n13     a1              5
   n13     n8      a1              5
   n8      n14     a1              4
   n14     n8      a1              9
   n8      n15     a1              6
   n15     n8      a1              1
   n8      n16     a1              8
   n16     n8      a1              7
   n8      n17     a1              5
   n17     n8      a1              5
   n8      n18     a1              2
   n18     n8      a1              5
   n8      n19     a1              9
   n19     n8      a1              9
   n8      n20     a1              3
   n20     n8      a1              7
   n9      n10     a1              6
   n10     n9      a1              2
   n9      n11     a1              6
   n11     n9      a1              2
   n9      n12     a1              5
   n12     n9      a1              3
   n9      n13     a1              4
   n13     n9      a1              6
   n9      n14     a1              7
   n14     n9      a1              6
   n9      n15     a1              2
   n15     n9      a1              2
   n9      n16     a1              8
   n16     n9      a1              6
   n9      n17     a1              2
   n17     n9      a1              2
   n9      n18     a1              8
   n18     n9      a1              6
   n9      n19     a1              6
   n19     n9      a1              6
   n9      n20     a1              1
   n20     n9      a1              5
   n10     n11     a1              9
   n11     n10     a1              1
   n10     n12     a1              7
   n12     n10     a1              2
   n10     n13     a1              7
   n13     n10     a1              6
   n10     n14     a1              4
   n14     n10     a1              6
   n10     n15     a1              2
   n15     n10     a1              9
   n10     n16     a1              1
   n16     n10     a1              8
   n10     n17     a1              9
   n17     n10     a1              3
   n10     n18     a1              9
   n18     n10     a1              1
   n10     n19     a1              9
   n19     n10     a1              2
   n10     n20     a1              1
   n20     n10     a1              9
   n11     n12     a1              2
   n12     n11     a1              6
   n11     n13     a1              8
   n13     n11     a1              6
   n11     n14     a1              4
   n14     n11     a1              3
   n11     n15     a1              1
   n15     n11     a1              5
   n11     n16     a1              1
   n16     n11     a1              4
   n11     n17     a1              1
   n17     n11     a1              2
   n11     n18     a1              3
   n18     n11     a1              2
   n11     n19     a1              9
   n19     n11     a1              3
   n11     n20     a1              6
   n20     n11     a1              3
   n12     n13     a1              4
   n13     n12     a1              9
   n12     n14     a1              2
   n14     n12     a1              4
   n12     n15     a1              1
   n15     n12     a1              2
   n12     n16     a1              4
   n16     n12     a1              2
   n12     n17     a1              9
   n17     n12     a1              8
   n12     n18     a1              4
   n18     n12     a1              8
   n12     n19     a1              8
   n19     n12     a1              4
   n12     n20     a1              1
   n20     n12     a1              7
   n13     n14     a1              9
   n14     n13     a1              6
   n13     n15     a1              4
   n15     n13     a1              5
   n13     n16     a1              4
   n16     n13     a1              6
   n13     n17     a1              2
   n17     n13     a1              5
   n13     n18     a1              7
   n18     n13     a1              4
   n13     n19     a1              3
   n19     n13     a1              8
   n13     n20     a1              1
   n20     n13     a1              8
   n14     n15     a1              7
   n15     n14     a1              8
   n14     n16     a1              2
   n16     n14     a1              6
   n14     n17     a1              9
   n17     n14     a1              6
   n14     n18     a1              8
   n18     n14     a1              2
   n14     n19     a1              2
   n19     n14     a1              5
   n14     n20     a1              7
   n20     n14     a1              3
   n15     n16     a1              3
   n16     n15     a1              3
   n15     n17     a1              2
   n17     n15     a1              2
   n15     n18     a1              7
   n18     n15     a1              9
   n15     n19     a1              5
   n19     n15     a1              9
   n15     n20     a1              1
   n20     n15     a1              5
   n16     n17     a1              6
   n17     n16     a1              8
   n16     n18     a1              1
   n18     n16     a1              2
   n16     n19     a1              4
   n19     n16     a1              7
   n16     n20     a1              9
   n20     n16     a1              4
   n17     n18     a1              9
   n18     n17     a1              8
   n17     n19     a1              4
   n19     n17     a1              9
   n17     n20     a1              9
   n20     n17     a1              8
   n18     n19     a1              3
   n19     n18     a1              6
   n18     n20     a1              6
   n20     n18     a1              9
   n19     n20     a1              1
   n20     n19     a1              2
$offdelim
/;

Parameter demand(n) /
   n1             63
   n3            -36
   n4           -118
   n5            -47
   n6            -67
   n7            -49
   n8            -42
   n9            -50
   n10            77
   n11           -18
   n14           -57
   n16            34
   n17           -16
   n18            85
   n19           110
   n20           131
/;

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