Set n / n1*n20 /; Set s / a1 /; yupp('n1','n1','a1')=0; Parameter fcost(m,n,s) / $ondelim n1 n5 a1 944 n5 n1 a1 924 n7 n4 a1 978 n4 n7 a1 762 n12 n7 a1 868 n7 n12 a1 744 n8 n20 a1 881 n20 n8 a1 910 n3 n12 a1 766 n12 n3 a1 833 n6 n13 a1 863 n13 n6 a1 940 n4 n15 a1 609 n15 n4 a1 906 n6 n15 a1 910 n15 n6 a1 893 n10 n14 a1 752 n14 n10 a1 987 n11 n8 a1 861 n8 n11 a1 526 n2 n14 a1 911 n14 n2 a1 540 n1 n16 a1 841 n16 n1 a1 807 n6 n12 a1 907 n12 n6 a1 719 n7 n5 a1 790 n5 n7 a1 602 n1 n12 a1 837 n12 n1 a1 922 n12 n5 a1 645 n5 n12 a1 815 n2 n16 a1 995 n16 n2 a1 674 n8 n6 a1 889 n6 n8 a1 751 n1 n6 a1 952 n6 n1 a1 727 n10 n6 a1 988 n6 n10 a1 880 n4 n13 a1 757 n13 n4 a1 986 n7 n8 a1 769 n8 n7 a1 654 n12 n17 a1 703 n17 n12 a1 637 n11 n10 a1 803 n10 n11 a1 854 n4 n5 a1 677 n5 n4 a1 819 n7 n9 a1 759 n9 n7 a1 692 n4 n12 a1 650 n12 n4 a1 899 n3 n19 a1 621 n19 n3 a1 766 n4 n19 a1 838 n19 n4 a1 713 n9 n12 a1 788 n12 n9 a1 965 n11 n5 a1 731 n5 n11 a1 746 n1 n13 a1 591 n13 n1 a1 717 n9 n11 a1 529 n11 n9 a1 699 n1 n20 a1 807 n20 n1 a1 964 n2 n11 a1 616 n11 n2 a1 895 n11 n6 a1 645 n6 n11 a1 798 n3 n13 a1 839 n13 n3 a1 623 n10 n20 a1 706 n20 n10 a1 860 n5 n13 a1 845 n13 n5 a1 705 n3 n8 a1 887 n8 n3 a1 511 n7 n3 a1 507 n3 n7 a1 938 n7 n15 a1 574 n15 n7 a1 817 n9 n1 a1 706 n1 n9 a1 839 n9 n5 a1 924 n5 n9 a1 561 n8 n5 a1 939 n5 n8 a1 742 n6 n19 a1 879 n19 n6 a1 611 n9 n18 a1 728 n18 n9 a1 997 n6 n2 a1 846 n2 n6 a1 939 n8 n16 a1 773 n16 n8 a1 544 n4 n17 a1 941 n17 n4 a1 819 n2 n15 a1 708 n15 n2 a1 916 n3 n9 a1 719 n9 n3 a1 701 n9 n14 a1 801 n14 n9 a1 683 n6 n16 a1 920 n16 n6 a1 590 n10 n16 a1 975 n16 n10 a1 771 n4 n6 a1 858 n6 n4 a1 907 n10 n12 a1 995 n12 n10 a1 580 n1 n2 a1 893 n2 n1 a1 755 n12 n14 a1 848 n14 n12 a1 532 n2 n9 a1 979 n9 n2 a1 533 n7 n6 a1 998 n6 n7 a1 932 n4 n18 a1 789 n18 n4 a1 526 n7 n11 a1 750 n11 n7 a1 955 n7 n2 a1 609 n2 n7 a1 731 n2 n18 a1 984 n18 n2 a1 885 n1 n7 a1 855 n7 n1 a1 606 n10 n1 a1 513 n1 n10 a1 768 n12 n8 a1 808 n8 n12 a1 697 n1 n19 a1 999 n19 n1 a1 570 n6 n9 a1 522 n9 n6 a1 555 n4 n14 a1 569 n14 n4 a1 814 n12 n13 a1 733 n13 n12 a1 947 n8 n14 a1 661 n14 n8 a1 780 n8 n13 a1 987 n13 n8 a1 730 n12 n20 a1 585 n20 n12 a1 541 n2 n20 a1 819 n20 n2 a1 764 n1 n3 a1 993 n3 n1 a1 921 n7 n20 a1 584 n20 n7 a1 674 n6 n17 a1 875 n17 n6 a1 580 n4 n20 a1 838 n20 n4 a1 713 n10 n15 a1 713 n15 n10 a1 609 n3 n14 a1 888 n14 n3 a1 899 n10 n3 a1 658 n3 n10 a1 560 n11 n20 a1 829 n20 n11 a1 509 n8 n19 a1 716 n19 n8 a1 785 n4 n2 a1 754 n2 n4 a1 587 n8 n1 a1 626 n1 n8 a1 790 n9 n16 a1 693 n16 n9 a1 517 n8 n4 a1 891 n4 n8 a1 865 n12 n11 a1 742 n11 n12 a1 704 n7 n18 a1 859 n18 n7 a1 817 n12 n2 a1 909 n2 n12 a1 527 n1 n11 a1 750 n11 n1 a1 728 n2 n13 a1 762 n13 n2 a1 930 n8 n10 a1 898 n10 n8 a1 962 n4 n10 a1 735 n10 n4 a1 934 n5 n17 a1 568 n17 n5 a1 514 n12 n18 a1 977 n18 n12 a1 958 n10 n7 a1 969 n7 n10 a1 850 n3 n18 a1 947 n18 n3 a1 861 $offdelim /; Parameter vcost(m,n,s) / $ondelim n1 n5 a1 2 n5 n1 a1 5 n7 n4 a1 3 n4 n7 a1 5 n12 n7 a1 7 n7 n12 a1 7 n8 n20 a1 3 n20 n8 a1 9 n3 n12 a1 4 n12 n3 a1 6 n6 n13 a1 4 n13 n6 a1 6 n4 n15 a1 3 n15 n4 a1 2 n6 n15 a1 4 n15 n6 a1 3 n10 n14 a1 2 n14 n10 a1 7 n11 n8 a1 4 n8 n11 a1 3 n2 n14 a1 1 n14 n2 a1 9 n1 n16 a1 6 n16 n1 a1 4 n6 n12 a1 2 n12 n6 a1 8 n7 n5 a1 3 n5 n7 a1 9 n1 n12 a1 4 n12 n1 a1 9 n12 n5 a1 4 n5 n12 a1 9 n2 n16 a1 6 n16 n2 a1 3 n8 n6 a1 9 n6 n8 a1 1 n1 n6 a1 8 n6 n1 a1 5 n10 n6 a1 7 n6 n10 a1 3 n4 n13 a1 4 n13 n4 a1 1 n7 n8 a1 6 n8 n7 a1 3 n12 n17 a1 6 n17 n12 a1 3 n11 n10 a1 3 n10 n11 a1 1 n4 n5 a1 1 n5 n4 a1 9 n7 n9 a1 3 n9 n7 a1 4 n4 n12 a1 3 n12 n4 a1 1 n3 n19 a1 5 n19 n3 a1 4 n4 n19 a1 7 n19 n4 a1 3 n9 n12 a1 3 n12 n9 a1 5 n11 n5 a1 7 n5 n11 a1 5 n1 n13 a1 8 n13 n1 a1 1 n9 n11 a1 4 n11 n9 a1 7 n1 n20 a1 6 n20 n1 a1 6 n2 n11 a1 2 n11 n2 a1 2 n11 n6 a1 9 n6 n11 a1 2 n3 n13 a1 3 n13 n3 a1 7 n10 n20 a1 3 n20 n10 a1 7 n5 n13 a1 3 n13 n5 a1 4 n3 n8 a1 3 n8 n3 a1 1 n7 n3 a1 3 n3 n7 a1 1 n7 n15 a1 1 n15 n7 a1 8 n9 n1 a1 6 n1 n9 a1 6 n9 n5 a1 2 n5 n9 a1 2 n8 n5 a1 4 n5 n8 a1 1 n6 n19 a1 8 n19 n6 a1 6 n9 n18 a1 2 n18 n9 a1 3 n6 n2 a1 4 n2 n6 a1 5 n8 n16 a1 2 n16 n8 a1 8 n4 n17 a1 6 n17 n4 a1 5 n2 n15 a1 2 n15 n2 a1 6 n3 n9 a1 6 n9 n3 a1 1 n9 n14 a1 6 n14 n9 a1 1 n6 n16 a1 1 n16 n6 a1 4 n10 n16 a1 4 n16 n10 a1 5 n4 n6 a1 8 n6 n4 a1 7 n10 n12 a1 2 n12 n10 a1 1 n1 n2 a1 2 n2 n1 a1 4 n12 n14 a1 6 n14 n12 a1 4 n2 n9 a1 7 n9 n2 a1 5 n7 n6 a1 9 n6 n7 a1 3 n4 n18 a1 4 n18 n4 a1 8 n7 n11 a1 8 n11 n7 a1 2 n7 n2 a1 2 n2 n7 a1 7 n2 n18 a1 6 n18 n2 a1 4 n1 n7 a1 6 n7 n1 a1 7 n10 n1 a1 2 n1 n10 a1 9 n12 n8 a1 5 n8 n12 a1 4 n1 n19 a1 4 n19 n1 a1 5 n6 n9 a1 5 n9 n6 a1 8 n4 n14 a1 2 n14 n4 a1 8 n12 n13 a1 3 n13 n12 a1 5 n8 n14 a1 4 n14 n8 a1 3 n8 n13 a1 5 n13 n8 a1 9 n12 n20 a1 7 n20 n12 a1 6 n2 n20 a1 1 n20 n2 a1 9 n1 n3 a1 5 n3 n1 a1 2 n7 n20 a1 1 n20 n7 a1 3 n6 n17 a1 9 n17 n6 a1 8 n4 n20 a1 9 n20 n4 a1 2 n10 n15 a1 8 n15 n10 a1 2 n3 n14 a1 7 n14 n3 a1 9 n10 n3 a1 4 n3 n10 a1 7 n11 n20 a1 4 n20 n11 a1 5 n8 n19 a1 9 n19 n8 a1 7 n4 n2 a1 5 n2 n4 a1 1 n8 n1 a1 3 n1 n8 a1 5 n9 n16 a1 5 n16 n9 a1 6 n8 n4 a1 4 n4 n8 a1 9 n12 n11 a1 3 n11 n12 a1 7 n7 n18 a1 4 n18 n7 a1 2 n12 n2 a1 2 n2 n12 a1 1 n1 n11 a1 8 n11 n1 a1 8 n2 n13 a1 6 n13 n2 a1 1 n8 n10 a1 4 n10 n8 a1 2 n4 n10 a1 3 n10 n4 a1 5 n5 n17 a1 2 n17 n5 a1 9 n12 n18 a1 8 n18 n12 a1 8 n10 n7 a1 9 n7 n10 a1 6 n3 n18 a1 6 n18 n3 a1 8 $offdelim /; Parameter demand(n) / n1 -61 n2 30 n3 -189 n4 342 n5 39 n6 157 n7 -143 n8 -113 n11 -137 n12 -26 n13 125 n14 -199 n15 14 n17 200 n18 -132 n20 93 /; 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; );