Set n / n1*n16 /; Set s / a1 /; yupp('n1','n1','a1')=0; Parameter fcost(m,n,s) / $ondelim n1 n2 a1 834 n2 n1 a1 724 n1 n3 a1 962 n3 n1 a1 645 n1 n4 a1 961 n4 n1 a1 942 n1 n5 a1 891 n5 n1 a1 653 n1 n6 a1 921 n6 n1 a1 895 n1 n7 a1 771 n7 n1 a1 912 n1 n8 a1 535 n8 n1 a1 811 n1 n9 a1 673 n9 n1 a1 711 n1 n10 a1 547 n10 n1 a1 757 n1 n11 a1 723 n11 n1 a1 778 n1 n12 a1 690 n12 n1 a1 606 n1 n13 a1 764 n13 n1 a1 805 n1 n14 a1 870 n14 n1 a1 601 n1 n15 a1 629 n15 n1 a1 954 n1 n16 a1 966 n16 n1 a1 808 n2 n3 a1 626 n3 n2 a1 538 n2 n4 a1 929 n4 n2 a1 615 n2 n5 a1 704 n5 n2 a1 806 n2 n6 a1 521 n6 n2 a1 572 n2 n7 a1 777 n7 n2 a1 512 n2 n8 a1 661 n8 n2 a1 767 n2 n9 a1 531 n9 n2 a1 650 n2 n10 a1 966 n10 n2 a1 691 n2 n11 a1 957 n11 n2 a1 883 n2 n12 a1 709 n12 n2 a1 588 n2 n13 a1 506 n13 n2 a1 668 n2 n14 a1 762 n14 n2 a1 859 n2 n15 a1 548 n15 n2 a1 541 n2 n16 a1 791 n16 n2 a1 520 n3 n4 a1 848 n4 n3 a1 984 n3 n5 a1 553 n5 n3 a1 938 n3 n6 a1 627 n6 n3 a1 893 n3 n7 a1 807 n7 n3 a1 617 n3 n8 a1 809 n8 n3 a1 951 n3 n9 a1 519 n9 n3 a1 803 n3 n10 a1 844 n10 n3 a1 702 n3 n11 a1 593 n11 n3 a1 587 n3 n12 a1 948 n12 n3 a1 618 n3 n13 a1 798 n13 n3 a1 798 n3 n14 a1 972 n14 n3 a1 792 n3 n15 a1 690 n15 n3 a1 691 n3 n16 a1 656 n16 n3 a1 634 n4 n5 a1 828 n5 n4 a1 886 n4 n6 a1 869 n6 n4 a1 916 n4 n7 a1 822 n7 n4 a1 699 n4 n8 a1 892 n8 n4 a1 712 n4 n9 a1 503 n9 n4 a1 688 n4 n10 a1 755 n10 n4 a1 585 n4 n11 a1 588 n11 n4 a1 757 n4 n12 a1 669 n12 n4 a1 689 n4 n13 a1 868 n13 n4 a1 555 n4 n14 a1 941 n14 n4 a1 560 n4 n15 a1 639 n15 n4 a1 996 n4 n16 a1 549 n16 n4 a1 949 n5 n6 a1 797 n6 n5 a1 605 n5 n7 a1 955 n7 n5 a1 701 n5 n8 a1 786 n8 n5 a1 855 n5 n9 a1 912 n9 n5 a1 945 n5 n10 a1 821 n10 n5 a1 982 n5 n11 a1 696 n11 n5 a1 661 n5 n12 a1 673 n12 n5 a1 659 n5 n13 a1 753 n13 n5 a1 510 n5 n14 a1 686 n14 n5 a1 522 n5 n15 a1 787 n15 n5 a1 891 n5 n16 a1 977 n16 n5 a1 824 n6 n7 a1 872 n7 n6 a1 570 n6 n8 a1 518 n8 n6 a1 770 n6 n9 a1 692 n9 n6 a1 980 n6 n10 a1 802 n10 n6 a1 625 n6 n11 a1 523 n11 n6 a1 681 n6 n12 a1 505 n12 n6 a1 531 n6 n13 a1 722 n13 n6 a1 900 n6 n14 a1 913 n14 n6 a1 670 n6 n15 a1 711 n15 n6 a1 667 n6 n16 a1 640 n16 n6 a1 624 n7 n8 a1 576 n8 n7 a1 802 n7 n9 a1 678 n9 n7 a1 984 n7 n10 a1 619 n10 n7 a1 587 n7 n11 a1 510 n11 n7 a1 815 n7 n12 a1 543 n12 n7 a1 609 n7 n13 a1 565 n13 n7 a1 974 n7 n14 a1 629 n14 n7 a1 520 n7 n15 a1 623 n15 n7 a1 761 n7 n16 a1 764 n16 n7 a1 616 n8 n9 a1 511 n9 n8 a1 547 n8 n10 a1 790 n10 n8 a1 763 n8 n11 a1 662 n11 n8 a1 713 n8 n12 a1 865 n12 n8 a1 700 n8 n13 a1 803 n13 n8 a1 824 n8 n14 a1 613 n14 n8 a1 850 n8 n15 a1 993 n15 n8 a1 540 n8 n16 a1 835 n16 n8 a1 850 n9 n10 a1 556 n10 n9 a1 857 n9 n11 a1 837 n11 n9 a1 861 n9 n12 a1 617 n12 n9 a1 696 n9 n13 a1 923 n13 n9 a1 565 n9 n14 a1 796 n14 n9 a1 962 n9 n15 a1 654 n15 n9 a1 869 n9 n16 a1 607 n16 n9 a1 635 n10 n11 a1 675 n11 n10 a1 642 n10 n12 a1 741 n12 n10 a1 825 n10 n13 a1 626 n13 n10 a1 783 n10 n14 a1 779 n14 n10 a1 734 n10 n15 a1 556 n15 n10 a1 705 n10 n16 a1 881 n16 n10 a1 641 n11 n12 a1 642 n12 n11 a1 877 n11 n13 a1 752 n13 n11 a1 673 n11 n14 a1 875 n14 n11 a1 610 n11 n15 a1 612 n15 n11 a1 831 n11 n16 a1 685 n16 n11 a1 589 n12 n13 a1 853 n13 n12 a1 740 n12 n14 a1 835 n14 n12 a1 605 n12 n15 a1 503 n15 n12 a1 608 n12 n16 a1 549 n16 n12 a1 627 n13 n14 a1 629 n14 n13 a1 674 n13 n15 a1 715 n15 n13 a1 992 n13 n16 a1 892 n16 n13 a1 928 n14 n15 a1 575 n15 n14 a1 970 n14 n16 a1 604 n16 n14 a1 663 n15 n16 a1 510 n16 n15 a1 740 $offdelim /; Parameter vcost(m,n,s) / $ondelim n1 n2 a1 9 n2 n1 a1 9 n1 n3 a1 1 n3 n1 a1 9 n1 n4 a1 7 n4 n1 a1 8 n1 n5 a1 1 n5 n1 a1 6 n1 n6 a1 8 n6 n1 a1 9 n1 n7 a1 1 n7 n1 a1 7 n1 n8 a1 2 n8 n1 a1 3 n1 n9 a1 4 n9 n1 a1 4 n1 n10 a1 2 n10 n1 a1 2 n1 n11 a1 8 n11 n1 a1 8 n1 n12 a1 3 n12 n1 a1 1 n1 n13 a1 6 n13 n1 a1 6 n1 n14 a1 7 n14 n1 a1 4 n1 n15 a1 3 n15 n1 a1 1 n1 n16 a1 6 n16 n1 a1 9 n2 n3 a1 4 n3 n2 a1 3 n2 n4 a1 4 n4 n2 a1 1 n2 n5 a1 4 n5 n2 a1 1 n2 n6 a1 4 n6 n2 a1 6 n2 n7 a1 7 n7 n2 a1 4 n2 n8 a1 7 n8 n2 a1 9 n2 n9 a1 5 n9 n2 a1 9 n2 n10 a1 8 n10 n2 a1 4 n2 n11 a1 7 n11 n2 a1 7 n2 n12 a1 9 n12 n2 a1 3 n2 n13 a1 8 n13 n2 a1 7 n2 n14 a1 3 n14 n2 a1 5 n2 n15 a1 9 n15 n2 a1 2 n2 n16 a1 6 n16 n2 a1 6 n3 n4 a1 1 n4 n3 a1 4 n3 n5 a1 1 n5 n3 a1 4 n3 n6 a1 9 n6 n3 a1 3 n3 n7 a1 2 n7 n3 a1 1 n3 n8 a1 2 n8 n3 a1 9 n3 n9 a1 3 n9 n3 a1 2 n3 n10 a1 2 n10 n3 a1 5 n3 n11 a1 2 n11 n3 a1 5 n3 n12 a1 9 n12 n3 a1 4 n3 n13 a1 6 n13 n3 a1 6 n3 n14 a1 5 n14 n3 a1 8 n3 n15 a1 6 n15 n3 a1 3 n3 n16 a1 2 n16 n3 a1 3 n4 n5 a1 4 n5 n4 a1 2 n4 n6 a1 3 n6 n4 a1 4 n4 n7 a1 4 n7 n4 a1 5 n4 n8 a1 1 n8 n4 a1 9 n4 n9 a1 9 n9 n4 a1 6 n4 n10 a1 7 n10 n4 a1 8 n4 n11 a1 8 n11 n4 a1 7 n4 n12 a1 9 n12 n4 a1 9 n4 n13 a1 1 n13 n4 a1 6 n4 n14 a1 8 n14 n4 a1 5 n4 n15 a1 1 n15 n4 a1 8 n4 n16 a1 2 n16 n4 a1 4 n5 n6 a1 3 n6 n5 a1 9 n5 n7 a1 2 n7 n5 a1 7 n5 n8 a1 8 n8 n5 a1 3 n5 n9 a1 2 n9 n5 a1 6 n5 n10 a1 1 n10 n5 a1 4 n5 n11 a1 5 n11 n5 a1 2 n5 n12 a1 3 n12 n5 a1 3 n5 n13 a1 7 n13 n5 a1 1 n5 n14 a1 4 n14 n5 a1 2 n5 n15 a1 6 n15 n5 a1 4 n5 n16 a1 8 n16 n5 a1 9 n6 n7 a1 9 n7 n6 a1 8 n6 n8 a1 2 n8 n6 a1 8 n6 n9 a1 5 n9 n6 a1 7 n6 n10 a1 5 n10 n6 a1 9 n6 n11 a1 8 n11 n6 a1 3 n6 n12 a1 1 n12 n6 a1 1 n6 n13 a1 8 n13 n6 a1 6 n6 n14 a1 9 n14 n6 a1 6 n6 n15 a1 9 n15 n6 a1 6 n6 n16 a1 5 n16 n6 a1 5 n7 n8 a1 2 n8 n7 a1 6 n7 n9 a1 3 n9 n7 a1 2 n7 n10 a1 7 n10 n7 a1 2 n7 n11 a1 1 n11 n7 a1 5 n7 n12 a1 6 n12 n7 a1 1 n7 n13 a1 5 n13 n7 a1 4 n7 n14 a1 5 n14 n7 a1 8 n7 n15 a1 5 n15 n7 a1 3 n7 n16 a1 9 n16 n7 a1 2 n8 n9 a1 5 n9 n8 a1 4 n8 n10 a1 2 n10 n8 a1 3 n8 n11 a1 8 n11 n8 a1 4 n8 n12 a1 2 n12 n8 a1 3 n8 n13 a1 5 n13 n8 a1 8 n8 n14 a1 6 n14 n8 a1 2 n8 n15 a1 1 n15 n8 a1 2 n8 n16 a1 4 n16 n8 a1 7 n9 n10 a1 5 n10 n9 a1 7 n9 n11 a1 2 n11 n9 a1 2 n9 n12 a1 9 n12 n9 a1 1 n9 n13 a1 3 n13 n9 a1 2 n9 n14 a1 5 n14 n9 a1 3 n9 n15 a1 8 n15 n9 a1 2 n9 n16 a1 3 n16 n9 a1 3 n10 n11 a1 7 n11 n10 a1 1 n10 n12 a1 5 n12 n10 a1 4 n10 n13 a1 6 n13 n10 a1 5 n10 n14 a1 2 n14 n10 a1 8 n10 n15 a1 3 n15 n10 a1 7 n10 n16 a1 1 n16 n10 a1 5 n11 n12 a1 7 n12 n11 a1 9 n11 n13 a1 9 n13 n11 a1 1 n11 n14 a1 8 n14 n11 a1 3 n11 n15 a1 2 n15 n11 a1 3 n11 n16 a1 2 n16 n11 a1 4 n12 n13 a1 8 n13 n12 a1 4 n12 n14 a1 9 n14 n12 a1 6 n12 n15 a1 5 n15 n12 a1 7 n12 n16 a1 4 n16 n12 a1 2 n13 n14 a1 4 n14 n13 a1 4 n13 n15 a1 3 n15 n13 a1 6 n13 n16 a1 4 n16 n13 a1 7 n14 n15 a1 8 n15 n14 a1 1 n14 n16 a1 7 n16 n14 a1 9 n15 n16 a1 6 n16 n15 a1 4 $offdelim /; Parameter demand(n) / n1 -402 n2 24 n3 104 n4 138 n5 -15 n6 209 n7 169 n8 -77 n9 127 n10 66 n11 -4 n12 163 n13 -101 n14 -190 n15 -86 n16 -125 /; 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; );