The problem considered is problem 8 of Heragu and Kusiak [3, p. 9]. The computer program listed below follows the mathematical formulation appearing on their page 4 [3] and on page 139 of Heragu [4]. The lengths of the thirty facilities are from their page 9 [3] and are used in line 561 through line 995 of the computer program below; the inter-facility flows are from Nugent et al. [5, p. 170] and are used in line 1238 through line 1347 of the computer program below.
The following computer program benefits from the computer programs of the present blog and the computer program of Conley [2, pp. 229-232]; the output is also presented here.
2 DEFINT I,J,K
5 DIM B(99),N(99),A(99),H(99),L(99),U(99),X(1111),D(111),P(522)
12 FOR JJJJ=-32000 TO 32000
14 RANDOMIZE JJJJ
16 M=-1D+17
25 B(1)=0:B(2)=0:B(3)=0:B(4)=0:B(5)=0:B(6)=0
26 B(7)=0:B(8)=0:B(9)=0:B(10)=0:B(11)=0
27 B(12)=0:B(13)=0:B(14)=0:B(15)=0
31 B(16)=0:B(17)=0:B(18)=0:B(19)=0:B(20)=0
32 B(21)=0:B(22)=0:B(23)=0:B(24)=0:B(25)=0
33 B(26)=0:B(27)=0:B(28)=0:B(29)=0:B(30)=0
37 N(1)=99200!:N(2)=99200!:N(3)=99200!:N(4)=99200!:N(5)=99200!:N(6)=99200!
38 N(7)=99200!:N(8)=99200!:N(9)=99200!:N(10)=99200!:N(11)=92900!
39 N(12)=99200!:N(13)=99200!:N(14)=99200!:N(15)=99200!
41 N(16)=99200!:N(17)=99200!:N(18)=99200!:N(19)=99200!:N(20)=92900!
42 N(21)=99200!:N(22)=99200!:N(23)=99200!:N(24)=99200!:N(25)=92900!
43 N(26)=99200!:N(27)=99200!:N(28)=99200!:N(29)=99200!:N(30)=92900!
48 GOTO 61
51 FOR IOCTT=1 TO 20
53 N(IOCTT)=999999!
57 NEXT IOCTT
61 FOR KLQ=1 TO 30
62 A(KLQ)=B(KLQ)+RND*(N(KLQ)-B(KLQ))
63 NEXT KLQ
71 FOR KLR=1 TO 30
72 H(KLR)=3
73 NEXT KLR
88 FOR JJJ=1 TO 1000 STEP 10
90 FOR INEW=1 TO 13
94 FOR J=INEW*2 TO 0 STEP -1
102 FOR JJ=0 TO 30
111 REM
128 FOR I=1 TO 2
129 FOR KKQQ=1 TO 30
130 X(KKQQ)=A(KKQQ)
131 NEXT KKQQ
132 FOR K=1 TO 30
133 IF RND<=.5 THEN 298 ELSE 230
230 IF A(K)-(N(K)-B(K))/H(K)^J
255 GOTO 265
260 L(K)=A(K)-(N(K)-B(K))/H(K)^J
265 IF A(K)+(N(K)-B(K))/H(K)^J>N(K) THEN 266 ELSE 268
266 U(K)=N(K)-L(K)
267 GOTO 272
268 U(K)=A(K)+(N(K)-B(K))/H(K)^J-L(K)
272 GOTO 300
298 X(K)=A(K)+2*RND*(2*RND-1)*(1/(1+RND*JJJ))*.05*A(K)
299 GOTO 302
300 X(K)=(L(K)+2*RND*RND*U(K))
302 NEXT K
304 X(JJ)=A(JJ)
305 REM
306 KLPS=FIX(1+RND*90)
307 FOR KLA1=1 TO KLPS
308 KLA2=FIX(1+30*RND)
309 X(KLA2)=A(KLA2)
310 NEXT KLA1
343 IF RND>.1 GOTO 461
352 IF RND<=.9 THEN 353 ELSE 355
353 X(JJ)=B(JJ)
354 GOTO 461
355 X(JJ)=N(JJ)
461 IF RND>.1 GOTO 561
471 IOCT1=1+FIX(RND*30)
474 IOCT2=1+FIX(RND*30)
477 X(IOCT1)=A(IOCT2)
480 X(IOCT2)=A(IOCT1)
561 P(11)=ABS(X(1)-X(2))-6
562 P(12)=ABS(X(1)-X(3))-3
563 P(13)=ABS(X(1)-X(4))-5
564 P(14)=ABS(X(1)-X(5))-3
565 P(15)=ABS(X(1)-X(6))-5
566 P(16)=ABS(X(1)-X(7))-4
567 P(17)=ABS(X(1)-X(8))-6
568 P(18)=ABS(X(1)-X(9))-4.5
569 P(19)=ABS(X(1)-X(10))-4
570 P(20)=ABS(X(1)-X(11))-3
571 P(21)=ABS(X(1)-X(12))-6
572 P(22)=ABS(X(1)-X(13))-3
573 P(23)=ABS(X(1)-X(14))-5
574 P(24)=ABS(X(1)-X(15))-3
575 P(25)=ABS(X(1)-X(16))-5
576 P(26)=ABS(X(1)-X(17))-4
577 P(27)=ABS(X(1)-X(18))-6
578 P(28)=ABS(X(1)-X(19))-4.5
579 P(29)=ABS(X(1)-X(20))-4
580 P(30)=ABS(X(2)-X(3))-6
581 P(31)=ABS(X(2)-X(4))-8
582 P(32)=ABS(X(2)-X(5))-6
583 P(33)=ABS(X(2)-X(6))-8
584 P(34)=ABS(X(2)-X(7))-7
585 P(35)=ABS(X(2)-X(8))-9
586 P(36)=ABS(X(2)-X(9))-7.5
587 P(37)=ABS(X(2)-X(10))-7
588 P(38)=ABS(X(2)-X(11))-6
589 P(39)=ABS(X(2)-X(12))-9
590 P(40)=ABS(X(2)-X(13))-6
591 P(41)=ABS(X(2)-X(14))-8
592 P(42)=ABS(X(2)-X(15))-6
593 P(43)=ABS(X(2)-X(16))-8
594 P(44)=ABS(X(2)-X(17))-7
595 P(45)=ABS(X(2)-X(18))-9
596 P(46)=ABS(X(2)-X(19))-7.5
597 P(47)=ABS(X(2)-X(20))-7
598 P(48)=ABS(X(3)-X(4))-5
599 P(49)=ABS(X(3)-X(5))-3
600 P(50)=ABS(X(3)-X(6))-5
601 P(51)=ABS(X(3)-X(7))-4
602 P(52)=ABS(X(3)-X(8))-6
603 P(53)=ABS(X(3)-X(9))-4.5
604 P(54)=ABS(X(3)-X(10))-4
605 P(55)=ABS(X(3)-X(11))-3
606 P(56)=ABS(X(3)-X(12))-6
607 P(57)=ABS(X(3)-X(13))-3
608 P(58)=ABS(X(3)-X(14))-5
609 P(59)=ABS(X(3)-X(15))-3
610 P(60)=ABS(X(3)-X(16))-5
611 P(61)=ABS(X(3)-X(17))-4
612 P(62)=ABS(X(3)-X(18))-6
613 P(63)=ABS(X(3)-X(19))-4.5
614 P(64)=ABS(X(3)-X(20))-4
615 P(65)=ABS(X(4)-X(5))-5
616 P(66)=ABS(X(4)-X(6))-7
617 P(67)=ABS(X(4)-X(7))-6
618 P(68)=ABS(X(4)-X(8))-8
619 P(69)=ABS(X(4)-X(9))-6.5
620 P(70)=ABS(X(4)-X(10))-6
621 P(71)=ABS(X(4)-X(11))-5
622 P(72)=ABS(X(4)-X(12))-8
623 P(73)=ABS(X(4)-X(13))-5
624 P(74)=ABS(X(4)-X(14))-7
625 P(75)=ABS(X(4)-X(15))-5
626 P(76)=ABS(X(4)-X(16))-7
627 P(77)=ABS(X(4)-X(17))-6
628 P(78)=ABS(X(4)-X(18))-8
629 P(79)=ABS(X(4)-X(19))-6.5
630 P(80)=ABS(X(4)-X(20))-6
631 P(81)=ABS(X(5)-X(6))-5
632 P(82)=ABS(X(5)-X(7))-4
633 P(83)=ABS(X(5)-X(8))-6
634 P(84)=ABS(X(5)-X(9))-4.5
635 P(85)=ABS(X(5)-X(10))-4
636 P(86)=ABS(X(5)-X(11))-3
637 P(87)=ABS(X(5)-X(12))-6
638 P(88)=ABS(X(5)-X(13))-3
639 P(89)=ABS(X(5)-X(14))-5
640 P(90)=ABS(X(5)-X(15))-3
641 P(91)=ABS(X(5)-X(16))-5
642 P(92)=ABS(X(5)-X(17))-4
643 P(93)=ABS(X(5)-X(18))-6
644 P(94)=ABS(X(5)-X(19))-4.5
645 P(95)=ABS(X(5)-X(20))-4
646 P(96)=ABS(X(6)-X(7))-6
647 P(97)=ABS(X(6)-X(8))-8
648 P(98)=ABS(X(6)-X(9))-6.5
649 P(99)=ABS(X(6)-X(10))-6
650 P(100)=ABS(X(6)-X(11))-5
651 P(101)=ABS(X(6)-X(12))-8
652 P(102)=ABS(X(6)-X(13))-5
653 P(103)=ABS(X(6)-X(14))-7
654 P(104)=ABS(X(6)-X(15))-5
655 P(105)=ABS(X(6)-X(16))-7
656 P(106)=ABS(X(6)-X(17))-6
657 P(107)=ABS(X(6)-X(18))-8
658 P(108)=ABS(X(6)-X(19))-6.5
659 P(109)=ABS(X(6)-X(20))-6
660 P(110)=ABS(X(7)-X(8))-7
661 P(111)=ABS(X(7)-X(9))-5.5
662 P(112)=ABS(X(7)-X(10))-5
663 P(113)=ABS(X(7)-X(11))-4
664 P(114)=ABS(X(7)-X(12))-7
665 P(115)=ABS(X(7)-X(13))-4
666 P(116)=ABS(X(7)-X(14))-6
667 P(117)=ABS(X(7)-X(15))-4
668 P(118)=ABS(X(7)-X(16))-6
669 P(119)=ABS(X(7)-X(17))-5
670 P(120)=ABS(X(7)-X(18))-7
671 P(121)=ABS(X(7)-X(19))-5.5
672 P(122)=ABS(X(7)-X(20))-5
673 P(123)=ABS(X(8)-X(9))-7.5
674 P(124)=ABS(X(8)-X(10))-7
675 P(125)=ABS(X(8)-X(11))-6
676 P(126)=ABS(X(8)-X(12))-9
677 P(127)=ABS(X(8)-X(13))-6
678 P(128)=ABS(X(8)-X(14))-8
679 P(129)=ABS(X(8)-X(15))-6
680 P(130)=ABS(X(8)-X(16))-8
681 P(131)=ABS(X(8)-X(17))-7
682 P(132)=ABS(X(8)-X(18))-9
683 P(133)=ABS(X(8)-X(19))-7.5
684 P(134)=ABS(X(8)-X(20))-7
685 P(135)=ABS(X(9)-X(10))-5.5
686 P(136)=ABS(X(9)-X(11))-4.5
687 P(137)=ABS(X(9)-X(12))-7.5
688 P(138)=ABS(X(9)-X(13))-4.5
689 P(139)=ABS(X(9)-X(14))-6.5
690 P(140)=ABS(X(9)-X(15))-4.5
691 P(141)=ABS(X(9)-X(16))-6.5
692 P(142)=ABS(X(9)-X(17))-5.5
693 P(143)=ABS(X(9)-X(18))-7.5
694 P(144)=ABS(X(9)-X(19))-6
695 P(145)=ABS(X(9)-X(20))-5.5
696 P(146)=ABS(X(10)-X(11))-4
697 P(147)=ABS(X(10)-X(12))-7
698 P(148)=ABS(X(10)-X(13))-4
699 P(149)=ABS(X(10)-X(14))-6
700 P(150)=ABS(X(10)-X(15))-4
701 P(151)=ABS(X(10)-X(16))-6
702 P(152)=ABS(X(10)-X(17))-5
703 P(153)=ABS(X(10)-X(18))-7
704 P(154)=ABS(X(10)-X(19))-5.5
705 P(155)=ABS(X(10)-X(20))-5
706 P(156)=ABS(X(11)-X(12))-6
707 P(157)=ABS(X(11)-X(13))-3
708 P(158)=ABS(X(11)-X(14))-5
709 P(159)=ABS(X(11)-X(15))-3
710 P(160)=ABS(X(11)-X(16))-5
711 P(161)=ABS(X(11)-X(17))-4
712 P(162)=ABS(X(11)-X(18))-6
713 P(163)=ABS(X(11)-X(19))-4.5
714 P(164)=ABS(X(11)-X(20))-4
715 P(165)=ABS(X(12)-X(13))-6
716 P(166)=ABS(X(12)-X(14))-8
717 P(167)=ABS(X(12)-X(15))-6
718 P(168)=ABS(X(12)-X(16))-8
719 P(169)=ABS(X(12)-X(17))-7
720 P(170)=ABS(X(12)-X(18))-9
721 P(171)=ABS(X(12)-X(19))-7.5
722 P(172)=ABS(X(12)-X(20))-7
723 P(173)=ABS(X(13)-X(14))-5
724 P(174)=ABS(X(13)-X(15))-3
725 P(175)=ABS(X(13)-X(16))-5
726 P(176)=ABS(X(13)-X(17))-4
727 P(177)=ABS(X(13)-X(18))-6
728 P(178)=ABS(X(13)-X(19))-4.5
729 P(179)=ABS(X(13)-X(20))-4
730 P(180)=ABS(X(14)-X(15))-5
731 P(181)=ABS(X(14)-X(16))-7
732 P(182)=ABS(X(14)-X(17))-6
733 P(183)=ABS(X(14)-X(18))-8
734 P(184)=ABS(X(14)-X(19))-6.5
735 P(185)=ABS(X(14)-X(20))-6
736 P(186)=ABS(X(15)-X(16))-5
737 P(187)=ABS(X(15)-X(17))-4
738 P(188)=ABS(X(15)-X(18))-6
739 P(189)=ABS(X(15)-X(19))-4.5
740 P(190)=ABS(X(15)-X(20))-4
741 P(191)=ABS(X(16)-X(17))-6
742 P(192)=ABS(X(16)-X(18))-8
743 P(193)=ABS(X(16)-X(19))-6.5
744 P(194)=ABS(X(16)-X(20))-6
745 P(195)=ABS(X(17)-X(18))-7
746 P(196)=ABS(X(17)-X(19))-5.5
747 P(197)=ABS(X(17)-X(20))-5
748 P(198)=ABS(X(18)-X(19))-7.5
749 P(199)=ABS(X(18)-X(20))-7
750 P(200)=ABS(X(19)-X(20))-5.5
751 P(201)=ABS(X(1)-X(21))-3
752 P(202)=ABS(X(1)-X(22))-6
753 P(203)=ABS(X(1)-X(23))-3
754 P(204)=ABS(X(1)-X(24))-5
755 P(205)=ABS(X(1)-X(25))-3
756 P(206)=ABS(X(1)-X(26))-5
757 P(207)=ABS(X(1)-X(27))-4
758 P(208)=ABS(X(1)-X(28))-6
759 P(209)=ABS(X(1)-X(29))-4.5
760 P(210)=ABS(X(1)-X(30))-4
761 P(211)=ABS(X(2)-X(21))-6
762 P(212)=ABS(X(2)-X(22))-9
763 P(213)=ABS(X(2)-X(23))-6
764 P(214)=ABS(X(2)-X(24))-8
765 P(215)=ABS(X(2)-X(25))-6
766 P(216)=ABS(X(2)-X(26))-8
767 P(217)=ABS(X(2)-X(27))-7
768 P(218)=ABS(X(2)-X(28))-9
769 P(219)=ABS(X(2)-X(29))-7.5
770 P(220)=ABS(X(2)-X(30))-7
771 P(221)=ABS(X(3)-X(21))-3
772 P(222)=ABS(X(3)-X(22))-6
773 P(223)=ABS(X(3)-X(23))-3
774 P(224)=ABS(X(3)-X(24))-5
775 P(225)=ABS(X(3)-X(25))-3
776 P(226)=ABS(X(3)-X(26))-5
777 P(227)=ABS(X(3)-X(27))-4
778 P(228)=ABS(X(3)-X(28))-6
779 P(229)=ABS(X(3)-X(29))-4.5
780 P(230)=ABS(X(3)-X(30))-4
781 P(231)=ABS(X(4)-X(21))-5
782 P(232)=ABS(X(4)-X(22))-8
783 P(233)=ABS(X(4)-X(23))-5
784 P(234)=ABS(X(4)-X(24))-7
785 P(235)=ABS(X(4)-X(25))-5
786 P(236)=ABS(X(4)-X(26))-7
787 P(237)=ABS(X(4)-X(27))-6
788 P(238)=ABS(X(4)-X(28))-8
789 P(239)=ABS(X(4)-X(29))-6.5
790 P(240)=ABS(X(4)-X(30))-6
791 P(241)=ABS(X(5)-X(21))-3
792 P(242)=ABS(X(5)-X(22))-6
793 P(243)=ABS(X(5)-X(23))-3
794 P(244)=ABS(X(5)-X(24))-5
795 P(245)=ABS(X(5)-X(25))-3
796 P(246)=ABS(X(5)-X(26))-5
797 P(247)=ABS(X(5)-X(27))-4
798 P(248)=ABS(X(5)-X(28))-6
799 P(249)=ABS(X(5)-X(29))-4.5
800 P(250)=ABS(X(5)-X(30))-4
801 P(251)=ABS(X(6)-X(21))-5
802 P(252)=ABS(X(6)-X(22))-8
803 P(253)=ABS(X(6)-X(23))-5
804 P(254)=ABS(X(6)-X(24))-7
805 P(255)=ABS(X(6)-X(25))-5
806 P(256)=ABS(X(6)-X(26))-7
807 P(257)=ABS(X(6)-X(27))-6
808 P(258)=ABS(X(6)-X(28))-8
809 P(259)=ABS(X(6)-X(29))-6.5
810 P(260)=ABS(X(6)-X(30))-6
811 P(261)=ABS(X(7)-X(21))-4
812 P(262)=ABS(X(7)-X(22))-7
813 P(263)=ABS(X(7)-X(23))-4
814 P(264)=ABS(X(7)-X(24))-6
815 P(265)=ABS(X(7)-X(25))-4
816 P(266)=ABS(X(7)-X(26))-6
817 P(267)=ABS(X(7)-X(27))-5
818 P(268)=ABS(X(7)-X(28))-7
819 P(269)=ABS(X(7)-X(29))-5.5
820 P(270)=ABS(X(7)-X(30))-5
821 P(271)=ABS(X(8)-X(21))-6
822 P(272)=ABS(X(8)-X(22))-9
823 P(273)=ABS(X(8)-X(23))-6
824 P(274)=ABS(X(8)-X(24))-8
825 P(275)=ABS(X(8)-X(25))-6
826 P(276)=ABS(X(8)-X(26))-8
827 P(277)=ABS(X(8)-X(27))-7
828 P(278)=ABS(X(8)-X(28))-9
829 P(279)=ABS(X(8)-X(29))-7.5
830 P(280)=ABS(X(8)-X(30))-7
831 P(281)=ABS(X(9)-X(21))-4.5
832 P(282)=ABS(X(9)-X(22))-7.5
833 P(283)=ABS(X(9)-X(23))-4.5
834 P(284)=ABS(X(9)-X(24))-6.5
835 P(285)=ABS(X(9)-X(25))-4.5
836 P(286)=ABS(X(9)-X(26))-6.5
837 P(287)=ABS(X(9)-X(27))-5.5
838 P(288)=ABS(X(9)-X(28))-7.5
839 P(289)=ABS(X(9)-X(29))-6
840 P(290)=ABS(X(9)-X(30))-5.5
841 P(291)=ABS(X(10)-X(21))-4
842 P(292)=ABS(X(10)-X(22))-7
843 P(293)=ABS(X(10)-X(23))-4
844 P(294)=ABS(X(10)-X(24))-6
845 P(295)=ABS(X(10)-X(25))-4
846 P(296)=ABS(X(10)-X(26))-6
847 P(297)=ABS(X(10)-X(27))-5
848 P(298)=ABS(X(10)-X(28))-7
849 P(299)=ABS(X(10)-X(29))-5.5
850 P(300)=ABS(X(10)-X(30))-5
851 P(301)=ABS(X(11)-X(21))-3
852 P(302)=ABS(X(11)-X(22))-6
853 P(303)=ABS(X(11)-X(23))-3
854 P(304)=ABS(X(11)-X(24))-5
855 P(305)=ABS(X(11)-X(25))-3
856 P(306)=ABS(X(11)-X(26))-5
857 P(307)=ABS(X(11)-X(27))-4
858 P(308)=ABS(X(11)-X(28))-6
859 P(309)=ABS(X(11)-X(29))-4.5
860 P(310)=ABS(X(11)-X(30))-4
861 P(311)=ABS(X(12)-X(21))-6
862 P(312)=ABS(X(12)-X(22))-9
863 P(313)=ABS(X(12)-X(23))-6
864 P(314)=ABS(X(12)-X(24))-8
865 P(315)=ABS(X(12)-X(25))-6
866 P(316)=ABS(X(12)-X(26))-8
867 P(317)=ABS(X(12)-X(27))-7
868 P(318)=ABS(X(12)-X(28))-9
869 P(319)=ABS(X(12)-X(29))-7.5
870 P(320)=ABS(X(12)-X(30))-7
871 P(321)=ABS(X(13)-X(21))-3
872 P(322)=ABS(X(13)-X(22))-6
873 P(323)=ABS(X(13)-X(23))-3
874 P(324)=ABS(X(13)-X(24))-5
875 P(325)=ABS(X(13)-X(25))-3
876 P(326)=ABS(X(13)-X(26))-5
877 P(327)=ABS(X(13)-X(27))-4
878 P(328)=ABS(X(13)-X(28))-6
879 P(329)=ABS(X(13)-X(29))-4.5
880 P(330)=ABS(X(13)-X(30))-4
881 P(331)=ABS(X(14)-X(21))-5
882 P(332)=ABS(X(14)-X(22))-8
883 P(333)=ABS(X(14)-X(23))-5
884 P(334)=ABS(X(14)-X(24))-7
885 P(335)=ABS(X(14)-X(25))-5
886 P(336)=ABS(X(14)-X(26))-7
887 P(337)=ABS(X(14)-X(27))-6
888 P(338)=ABS(X(14)-X(28))-8
889 P(339)=ABS(X(14)-X(29))-6.5
890 P(340)=ABS(X(14)-X(30))-6
891 P(341)=ABS(X(15)-X(21))-3
892 P(342)=ABS(X(15)-X(22))-6
893 P(343)=ABS(X(15)-X(23))-3
894 P(344)=ABS(X(15)-X(24))-5
895 P(345)=ABS(X(15)-X(25))-3
896 P(346)=ABS(X(15)-X(26))-5
897 P(347)=ABS(X(15)-X(27))-4
898 P(348)=ABS(X(15)-X(28))-6
899 P(349)=ABS(X(15)-X(29))-4.5
900 P(350)=ABS(X(15)-X(30))-4
901 P(351)=ABS(X(16)-X(21))-5
902 P(352)=ABS(X(16)-X(22))-8
903 P(353)=ABS(X(16)-X(23))-5
904 P(354)=ABS(X(16)-X(24))-7
905 P(355)=ABS(X(16)-X(25))-5
906 P(356)=ABS(X(16)-X(26))-7
907 P(357)=ABS(X(16)-X(27))-6
908 P(358)=ABS(X(16)-X(28))-8
909 P(359)=ABS(X(16)-X(29))-6.5
910 P(360)=ABS(X(16)-X(30))-6
911 P(361)=ABS(X(17)-X(21))-4
912 P(362)=ABS(X(17)-X(22))-7
913 P(363)=ABS(X(17)-X(23))-4
914 P(364)=ABS(X(17)-X(24))-6
915 P(365)=ABS(X(17)-X(25))-4
916 P(366)=ABS(X(17)-X(26))-6
917 P(367)=ABS(X(17)-X(27))-5
918 P(368)=ABS(X(17)-X(28))-7
919 P(369)=ABS(X(17)-X(29))-5.5
920 P(370)=ABS(X(17)-X(30))-5
921 P(371)=ABS(X(18)-X(21))-6
922 P(372)=ABS(X(18)-X(22))-9
923 P(373)=ABS(X(18)-X(23))-6
924 P(374)=ABS(X(18)-X(24))-8
925 P(375)=ABS(X(18)-X(25))-6
926 P(376)=ABS(X(18)-X(26))-8
927 P(377)=ABS(X(18)-X(27))-7
928 P(378)=ABS(X(18)-X(28))-9
929 P(379)=ABS(X(18)-X(29))-7.5
930 P(380)=ABS(X(18)-X(30))-7
931 P(381)=ABS(X(19)-X(21))-4.5
932 P(382)=ABS(X(19)-X(22))-7.5
933 P(383)=ABS(X(19)-X(23))-4.5
934 P(384)=ABS(X(19)-X(24))-6.5
935 P(385)=ABS(X(19)-X(25))-4.5
936 P(386)=ABS(X(19)-X(26))-6.5
937 P(387)=ABS(X(19)-X(27))-5.5
938 P(388)=ABS(X(19)-X(28))-7.5
939 P(389)=ABS(X(19)-X(29))-6
940 P(390)=ABS(X(19)-X(30))-5.5
941 P(391)=ABS(X(20)-X(21))-4
942 P(392)=ABS(X(20)-X(22))-7
943 P(393)=ABS(X(20)-X(23))-4
944 P(394)=ABS(X(20)-X(24))-6
945 P(395)=ABS(X(20)-X(25))-4
946 P(396)=ABS(X(20)-X(26))-6
947 P(397)=ABS(X(20)-X(27))-5
948 P(398)=ABS(X(20)-X(28))-7
949 P(399)=ABS(X(20)-X(29))-5.5
950 P(400)=ABS(X(20)-X(30))-5
951 P(401)=ABS(X(21)-X(22))-6
952 P(402)=ABS(X(21)-X(23))-3
953 P(403)=ABS(X(21)-X(24))-5
954 P(404)=ABS(X(21)-X(25))-3
955 P(405)=ABS(X(21)-X(26))-5
956 P(406)=ABS(X(21)-X(27))-4
957 P(407)=ABS(X(21)-X(28))-6
958 P(408)=ABS(X(21)-X(29))-4.5
959 P(409)=ABS(X(21)-X(30))-4
960 P(410)=ABS(X(22)-X(23))-6
961 P(411)=ABS(X(22)-X(24))-8
962 P(412)=ABS(X(22)-X(25))-6
963 P(413)=ABS(X(22)-X(26))-8
964 P(414)=ABS(X(22)-X(27))-7
965 P(415)=ABS(X(22)-X(28))-9
966 P(416)=ABS(X(22)-X(29))-7.5
967 P(417)=ABS(X(22)-X(30))-7
968 P(418)=ABS(X(23)-X(24))-5
969 P(419)=ABS(X(23)-X(25))-3
970 P(420)=ABS(X(23)-X(26))-5
971 P(421)=ABS(X(23)-X(27))-4
972 P(422)=ABS(X(23)-X(28))-6
973 P(423)=ABS(X(23)-X(29))-4.5
974 P(424)=ABS(X(23)-X(30))-4
975 P(425)=ABS(X(24)-X(25))-5
976 P(426)=ABS(X(24)-X(26))-7
977 P(427)=ABS(X(24)-X(27))-6
978 P(428)=ABS(X(24)-X(28))-8
979 P(429)=ABS(X(24)-X(29))-6.5
980 P(430)=ABS(X(24)-X(30))-6
981 P(431)=ABS(X(25)-X(26))-5
982 P(432)=ABS(X(25)-X(27))-4
983 P(433)=ABS(X(25)-X(28))-6
984 P(434)=ABS(X(25)-X(29))-4.5
985 P(435)=ABS(X(25)-X(30))-4
986 P(436)=ABS(X(26)-X(27))-6
987 P(437)=ABS(X(26)-X(28))-8
988 P(438)=ABS(X(26)-X(29))-6.5
989 P(439)=ABS(X(26)-X(30))-6
990 P(440)=ABS(X(27)-X(28))-7
991 P(441)=ABS(X(27)-X(29))-5.5
992 P(442)=ABS(X(27)-X(30))-5
993 P(443)=ABS(X(28)-X(29))-7.5
994 P(444)=ABS(X(28)-X(30))-7
995 P(445)=ABS(X(29)-X(30))-5.5
996 REM
997 REM
998 REM
999 REM
1088 FOR INSI=11 TO 445
1091 IF P(INSI)<0 THEN P(INSI)=P(INSI) ELSE P(INSI)=0
1095 NEXT INSI
1111 REM
1115 REM
1118 REM
1121 REM
1131 REM
1211 PSUM=0
1215 FOR IOCT=11 TO 445
1218 PSUM=PSUM+P(IOCT)
1221 NEXT IOCT
1231 PS1=99*555555!*PSUM
1238 P11B=3*ABS(X(1)-X(2))+2*ABS(X(1)-X(3))+0*ABS(X(1)-X(4))+0*ABS(X(1)-X(5))+2*ABS(X(1)-X(6))+10*ABS(X(1)-X(7))+5*ABS(X(1)-X(8))+0*ABS(X(1)-X(9))+5*ABS(X(1)-X(10))+2*ABS(X(1)-X(11))
1239 P12B=5*ABS(X(1)-X(12))+0*ABS(X(1)-X(13))+0*ABS(X(1)-X(14))+2*ABS(X(1)-X(15))
1240 P13B=4*ABS(X(2)-X(3))+0*ABS(X(2)-X(4))+10*ABS(X(2)-X(5))+4*ABS(X(2)-X(6))+0*ABS(X(2)-X(7))+0*ABS(X(2)-X(8))+2*ABS(X(2)-X(9))+2*ABS(X(2)-X(10))+1*ABS(X(2)-X(11))
1241 P14B=0*ABS(X(2)-X(12))+5*ABS(X(2)-X(13))+0*ABS(X(2)-X(14))+0*ABS(X(2)-X(15))
1243 P15B=3*ABS(X(3)-X(4))+4*ABS(X(3)-X(5))+0*ABS(X(3)-X(6))+5*ABS(X(3)-X(7))+5*ABS(X(3)-X(8))+5*ABS(X(3)-X(9))+1*ABS(X(3)-X(10))+4*ABS(X(3)-X(11))
1244 P16B=1*ABS(X(3)-X(12))+0*ABS(X(3)-X(13))+4*ABS(X(3)-X(14))+0*ABS(X(3)-X(15))
1245 P17B=0*ABS(X(4)-X(5))+0*ABS(X(4)-X(6))+0*ABS(X(4)-X(7))+2*ABS(X(4)-X(8))+2*ABS(X(4)-X(9))+0*ABS(X(4)-X(10))+6*ABS(X(4)-X(11))
1246 P18B=0*ABS(X(4)-X(12))+2*ABS(X(4)-X(13))+5*ABS(X(4)-X(14))+2*ABS(X(4)-X(15))
1248 P19B=5*ABS(X(5)-X(6))+2*ABS(X(5)-X(7))+0*ABS(X(5)-X(8))+0*ABS(X(5)-X(9))+0*ABS(X(5)-X(10))+0*ABS(X(5)-X(11))+2*ABS(X(5)-X(12))+0*ABS(X(5)-X(13))+0*ABS(X(5)-X(14))+0*ABS(X(5)-X(15))
1249 P20B=1*ABS(X(6)-X(7))+2*ABS(X(6)-X(8))+2*ABS(X(6)-X(9))+1*ABS(X(6)-X(10))+4*ABS(X(6)-X(11))+10*ABS(X(6)-X(12))+10*ABS(X(6)-X(13))+2*ABS(X(6)-X(14))+5*ABS(X(6)-X(15))
1250 P21B=10*ABS(X(7)-X(8))+10*ABS(X(7)-X(9))+5*ABS(X(7)-X(10))+10*ABS(X(7)-X(11))+10*ABS(X(7)-X(12))+6*ABS(X(7)-X(13))+0*ABS(X(7)-X(14))+0*ABS(X(7)-X(15))
1251 P22B=1*ABS(X(8)-X(9))+3*ABS(X(8)-X(10))+5*ABS(X(8)-X(11))+0*ABS(X(8)-X(12))+0*ABS(X(8)-X(13))+0*ABS(X(8)-X(14))+2*ABS(X(8)-X(15))
1255 P23B=10*ABS(X(9)-X(10))+2*ABS(X(9)-X(11))+1*ABS(X(9)-X(12))+5*ABS(X(9)-X(13))+2*ABS(X(9)-X(14))+0*ABS(X(9)-X(15))
1257 P24B=5*ABS(X(10)-X(11))+5*ABS(X(10)-X(12))+6*ABS(X(10)-X(13))+0*ABS(X(10)-X(14))+1*ABS(X(10)-X(15))
1258 P25B=0*ABS(X(11)-X(12))+0*ABS(X(11)-X(13))+1*ABS(X(11)-X(14))+2*ABS(X(11)-X(15))
1259 P26B=5*ABS(X(12)-X(13))+5*ABS(X(12)-X(14))+2*ABS(X(12)-X(15))+2*ABS(X(13)-X(14))+0*ABS(X(13)-X(15))+2*ABS(X(14)-X(15))
1301 P27B=0*ABS(X(1)-X(16))+5*ABS(X(1)-X(17))+6*ABS(X(1)-X(18))+3*ABS(X(1)-X(18))+0*ABS(X(1)-X(20))
1302 P28B=0*ABS(X(2)-X(16))+0*ABS(X(2)-X(17))+2*ABS(X(2)-X(18))+0*ABS(X(2)-X(19))+1*ABS(X(2)-X(20))
1303 P29B=4*ABS(X(3)-X(16))+0*ABS(X(3)-X(17))+6*ABS(X(3)-X(18))+3*ABS(X(3)-X(19))+2*ABS(X(3)-X(20))
1304 P30B=5*ABS(X(4)-X(16))+1*ABS(X(4)-X(17))+1*ABS(X(4)-X(18))+1*ABS(X(4)-X(19))+1*ABS(X(4)-X(20))
1305 P31B=0*ABS(X(5)-X(16))+2*ABS(X(5)-X(17))+1*ABS(X(5)-X(18))+0*ABS(X(5)-X(19))+0*ABS(X(5)-X(20))
1306 P32B=5*ABS(X(6)-X(16))+0*ABS(X(6)-X(17))+5*ABS(X(6)-X(18))+0*ABS(X(6)-X(19))+0*ABS(X(6)-X(20))
1307 P33B=10*ABS(X(7)-X(16))+2*ABS(X(7)-X(17))+1*ABS(X(7)-X(18))+10*ABS(X(7)-X(19))+1*ABS(X(7)-X(20))
1308 P34B=4*ABS(X(8)-X(16))+5*ABS(X(8)-X(17))+2*ABS(X(8)-X(18))+10*ABS(X(8)-X(19))+6*ABS(X(8)-X(20))
1309 P35B=3*ABS(X(9)-X(16))+0*ABS(X(9)-X(17))+2*ABS(X(9)-X(18))+0*ABS(X(9)-X(19))+0*ABS(X(9)-X(20))
1310 P36B=5*ABS(X(10)-X(16))+5*ABS(X(10)-X(17))+0*ABS(X(10)-X(18))+5*ABS(X(10)-X(19))+2*ABS(X(10)-X(20))
1311 P37B=1*ABS(X(11)-X(16))+0*ABS(X(11)-X(17))+2*ABS(X(11)-X(18))+0*ABS(X(11)-X(19))+0*ABS(X(11)-X(20))
1312 P38B=0*ABS(X(12)-X(16))+0*ABS(X(12)-X(17))+0*ABS(X(12)-X(18))+0*ABS(X(12)-X(19))+2*ABS(X(12)-X(20))
1313 P39B=4*ABS(X(13)-X(16))+2*ABS(X(13)-X(17))+2*ABS(X(13)-X(18))+1*ABS(X(13)-X(19))+0*ABS(X(13)-X(20))
1314 P40B=1*ABS(X(14)-X(16))+0*ABS(X(14)-X(17))+5*ABS(X(14)-X(18))+3*ABS(X(14)-X(19))+10*ABS(X(14)-X(20))
1315 P41B=4*ABS(X(15)-X(16))+5*ABS(X(15)-X(17))+1*ABS(X(15)-X(18))+0*ABS(X(15)-X(19))+1*ABS(X(15)-X(20))
1316 P42B=0*ABS(X(16)-X(17))+3*ABS(X(16)-X(18))+0*ABS(X(16)-X(19))+2*ABS(X(16)-X(20))
1317 P43B=2*ABS(X(17)-X(18))+2*ABS(X(17)-X(19))+0*ABS(X(17)-X(20))+5*ABS(X(18)-X(19))+1*ABS(X(18)-X(20))+0*ABS(X(19)-X(20))
1321 P44B=1*ABS(X(1)-X(21))+10*ABS(X(1)-X(22))+0*ABS(X(1)-X(23))+10*ABS(X(1)-X(24))+2*ABS(X(1)-X(25))+1*ABS(X(1)-X(26))+1*ABS(X(1)-X(27))+1*ABS(X(1)-X(28))+0*ABS(X(1)-X(29))+1*ABS(X(1)-X(30))
1322 P45B=6*ABS(X(2)-X(21))+1*ABS(X(2)-X(22))+0*ABS(X(2)-X(23))+1*ABS(X(2)-X(24))+2*ABS(X(2)-X(25))+2*ABS(X(2)-X(26))+5*ABS(X(2)-X(27))+1*ABS(X(2)-X(28))+10*ABS(X(2)-X(29))+5*ABS(X(2)-X(30))
1323 P46B=5*ABS(X(3)-X(21))+5*ABS(X(3)-X(22))+2*ABS(X(3)-X(23))+1*ABS(X(3)-X(24))+0*ABS(X(3)-X(25))+0*ABS(X(3)-X(26))+3*ABS(X(3)-X(27))+1*ABS(X(3)-X(28))+0*ABS(X(3)-X(29))+2*ABS(X(3)-X(30))
1324 P47B=2*ABS(X(4)-X(21))+2*ABS(X(4)-X(22))+4*ABS(X(4)-X(23))+0*ABS(X(4)-X(24))+2*ABS(X(4)-X(25))+0*ABS(X(4)-X(26))+2*ABS(X(4)-X(27))+2*ABS(X(4)-X(28))+5*ABS(X(4)-X(29))+5*ABS(X(4)-X(30))
1325 P48B=2*ABS(X(5)-X(21))+0*ABS(X(5)-X(22))+5*ABS(X(5)-X(23))+1*ABS(X(5)-X(24))+0*ABS(X(5)-X(25))+2*ABS(X(5)-X(26))+1*ABS(X(5)-X(27))+0*ABS(X(5)-X(28))+2*ABS(X(5)-X(29))+1*ABS(X(5)-X(30))
1326 P49B=0*ABS(X(6)-X(21))+10*ABS(X(6)-X(22))+0*ABS(X(6)-X(23))+0*ABS(X(6)-X(24))+0*ABS(X(6)-X(25))+4*ABS(X(6)-X(26))+0*ABS(X(6)-X(27))+10*ABS(X(6)-X(28))+1*ABS(X(6)-X(29))+1*ABS(X(6)-X(30))
1327 P50B=5*ABS(X(7)-X(21))+5*ABS(X(7)-X(22))+2*ABS(X(7)-X(23))+3*ABS(X(7)-X(24))+5*ABS(X(7)-X(25))+0*ABS(X(7)-X(26))+2*ABS(X(7)-X(27))+0*ABS(X(7)-X(28))+1*ABS(X(7)-X(29))+3*ABS(X(7)-X(30))
1328 P51B=0*ABS(X(8)-X(21))+5*ABS(X(8)-X(22))+5*ABS(X(8)-X(23))+2*ABS(X(8)-X(24))+5*ABS(X(8)-X(25))+0*ABS(X(8)-X(26))+5*ABS(X(8)-X(27))+5*ABS(X(8)-X(28))+0*ABS(X(8)-X(29))+2*ABS(X(8)-X(30))
1329 P52B=4*ABS(X(9)-X(21))+0*ABS(X(9)-X(22))+5*ABS(X(9)-X(23))+2*ABS(X(9)-X(24))+0*ABS(X(9)-X(25))+5*ABS(X(9)-X(26))+2*ABS(X(9)-X(27))+2*ABS(X(9)-X(28))+5*ABS(X(9)-X(29))+2*ABS(X(9)-X(30))
1330 P53B=3*ABS(X(10)-X(21))+5*ABS(X(10)-X(22))+0*ABS(X(10)-X(23))+5*ABS(X(10)-X(24))+2*ABS(X(10)-X(25))+10*ABS(X(10)-X(26))+10*ABS(X(10)-X(27))+1*ABS(X(10)-X(28))+5*ABS(X(10)-X(29))+2*ABS(X(10)-X(30))
1331 P54B=0*ABS(X(11)-X(21))+6*ABS(X(11)-X(22))+6*ABS(X(11)-X(23))+0*ABS(X(11)-X(24))+4*ABS(X(11)-X(25))+5*ABS(X(11)-X(26))+3*ABS(X(11)-X(27))+2*ABS(X(11)-X(28))+2*ABS(X(11)-X(29))+10*ABS(X(11)-X(30))
1332 P55B=0*ABS(X(12)-X(21))+4*ABS(X(12)-X(22))+5*ABS(X(12)-X(23))+10*ABS(X(12)-X(24))+1*ABS(X(12)-X(25))+0*ABS(X(12)-X(26))+0*ABS(X(12)-X(27))+0*ABS(X(12)-X(28))+0*ABS(X(12)-X(29))+1*ABS(X(12)-X(30))
1333 P56B=6*ABS(X(13)-X(21))+2*ABS(X(13)-X(22))+1*ABS(X(13)-X(23))+5*ABS(X(13)-X(24))+5*ABS(X(13)-X(25))+0*ABS(X(13)-X(26))+0*ABS(X(13)-X(27))+1*ABS(X(13)-X(28))+5*ABS(X(13)-X(29))+5*ABS(X(13)-X(30))
1334 P57B=0*ABS(X(14)-X(21))+0*ABS(X(14)-X(22))+4*ABS(X(14)-X(23))+2*ABS(X(14)-X(24))+0*ABS(X(14)-X(25))+0*ABS(X(14)-X(26))+4*ABS(X(14)-X(27))+2*ABS(X(14)-X(28))+5*ABS(X(14)-X(29))+5*ABS(X(14)-X(30))
1335 P58B=0*ABS(X(15)-X(21))+5*ABS(X(15)-X(22))+0*ABS(X(15)-X(23))+2*ABS(X(15)-X(24))+0*ABS(X(15)-X(25))+0*ABS(X(15)-X(26))+5*ABS(X(15)-X(27))+1*ABS(X(15)-X(28))+1*ABS(X(15)-X(29))+0*ABS(X(15)-X(30))
1336 P59B=2*ABS(X(16)-X(21))+0*ABS(X(16)-X(22))+2*ABS(X(16)-X(23))+0*ABS(X(16)-X(24))+5*ABS(X(16)-X(25))+0*ABS(X(16)-X(26))+5*ABS(X(16)-X(27))+2*ABS(X(16)-X(28))+5*ABS(X(16)-X(29))+10*ABS(X(16)-X(30))
1337 P60B=0*ABS(X(17)-X(21))+0*ABS(X(17)-X(22))+6*ABS(X(17)-X(23))+5*ABS(X(17)-X(24))+3*ABS(X(17)-X(25))+5*ABS(X(17)-X(26))+0*ABS(X(17)-X(27))+0*ABS(X(17)-X(28))+5*ABS(X(17)-X(29))+1*ABS(X(17)-X(30))
1338 P61B=2*ABS(X(18)-X(21))+10*ABS(X(18)-X(22))+10*ABS(X(18)-X(23))+4*ABS(X(18)-X(24))+0*ABS(X(18)-X(25))+0*ABS(X(18)-X(26))+5*ABS(X(18)-X(27))+0*ABS(X(18)-X(28))+0*ABS(X(18)-X(29))+0*ABS(X(18)-X(30))
1339 P62B=5*ABS(X(19)-X(21))+5*ABS(X(19)-X(22))+1*ABS(X(19)-X(23))+0*ABS(X(19)-X(24))+5*ABS(X(19)-X(25))+2*ABS(X(19)-X(26))+1*ABS(X(19)-X(27))+2*ABS(X(19)-X(28))+10*ABS(X(19)-X(29))+10*ABS(X(19)-X(30))
1340 P63B=5*ABS(X(20)-X(21))+2*ABS(X(20)-X(22))+1*ABS(X(20)-X(23))+3*ABS(X(20)-X(24))+1*ABS(X(20)-X(25))+5*ABS(X(20)-X(26))+6*ABS(X(20)-X(27))+5*ABS(X(20)-X(28))+5*ABS(X(20)-X(29))+3*ABS(X(20)-X(30))
1341 P64B=4*ABS(X(21)-X(22))+0*ABS(X(21)-X(23))+1*ABS(X(21)-X(24))+0*ABS(X(21)-X(25))+0*ABS(X(21)-X(26))+0*ABS(X(21)-X(27))+5*ABS(X(21)-X(28))+0*ABS(X(21)-X(29))+0*ABS(X(21)-X(30))
1342 P65B=5*ABS(X(22)-X(23))+0*ABS(X(22)-X(24))+4*ABS(X(22)-X(25))+4*ABS(X(22)-X(26))+5*ABS(X(22)-X(27))+0*ABS(X(22)-X(28))+2*ABS(X(22)-X(29))+5*ABS(X(22)-X(30))
1343 P66B=0*ABS(X(23)-X(24))+4*ABS(X(23)-X(25))+4*ABS(X(23)-X(26))+1*ABS(X(23)-X(27))+0*ABS(X(23)-X(28))+2*ABS(X(23)-X(29))+2*ABS(X(23)-X(30))
1344 P67B=5*ABS(X(24)-X(25))+5*ABS(X(24)-X(26))+0*ABS(X(24)-X(27))+1*ABS(X(24)-X(28))+0*ABS(X(24)-X(29))+0*ABS(X(24)-X(30))
1345 P68B=1*ABS(X(25)-X(26))+0*ABS(X(25)-X(27))+10*ABS(X(25)-X(28))+1*ABS(X(25)-X(29))+0*ABS(X(25)-X(30))
1346 P69B=0*ABS(X(26)-X(27))+0*ABS(X(26)-X(28))+0*ABS(X(26)-X(29))+0*ABS(X(26)-X(30))
1347 P70B=0*ABS(X(27)-X(28))+0*ABS(X(27)-X(29))+10*ABS(X(27)-X(30))+2*ABS(X(28)-X(29))+2*ABS(X(28)-X(30))+2*ABS(X(29)-X(30))
1388 P1=P11B+P12B+P13B+P14B+P15B+P16B+P17B+P18B+P19B+P20B+P21B+P22B+P23B+P24B+P25B+P26B
1395 P2=P27B+P28B+P29B+P30B+P31B+P32B+P33B+P34B+P35B+P36B+P37B+P38B+P39B+P40B+P41B+P42B+P43B
1397 P3=P44B+P45B+P46B+P47B+P48B+P49B+P50B+P51B+P52B+P53B+P54B+P55B+P56B+P57B+P58B+P59B+P60B
1398 P4=P61B+P62B+P63B+P64B+P65B+P66B+P67B+P68B+P69B+P70B
1448 P5=P1+P2+P3+P4
1449 REM
1450 P=-P5+PS1
1451 IF P<=M THEN 1670
1452 M=P
1453 PP1=P5
1454 FOR KLX=1 TO 30
1455 A(KLX)=X(KLX)
1456 NEXT KLX
1457 GOTO 128
1670 NEXT I
1702 NEXT JJ
1706 NEXT J
1777 NEXT INEW
1778 REM
1888 NEXT JJJ
1890 IF M>-47000! THEN 1916 ELSE 1999
1916 PRINT JJJJ,A(1),A(2),A(3),A(4),A(5),M,PP1
1918 PRINT A(6),A(7),A(8),A(9),A(10)
1920 PRINT A(11),A(12),A(13),A(14),A(15)
1922 PRINT A(16),A(17),A(18),A(19),A(20)
1924 PRINT A(21),A(22),A(23),A(24),A(25)
1926 PRINT A(26),A(27),A(28),A(29),A(30)
1999 NEXT JJJJ
This BASIC computer program was run with the IBM basica/D interpreter, and its candidate solutions from JJJJ=-32000 through JJJJ=-31963 (in compressed form and to be interpreted in accordance with line 1916 through line 1926; copied manually from the computer monitor) are presented below.
-31963 39963.88 39924.88 40009.88 40071.88
39937.88 -45961.73 45961.73
39971.88 40005.88 40045.88 40022.38 39986.88
40001.88 39943.88 39998.88 40059.88 39966.88
40034.88 39959.88 39952.88 40028.38 40065.88
40012.88 39979.88 39990.88 39932.88 40039.88
39916.88 39994.88 40079.88 40053.38 40016.88
And this same computer run also produced candidate solutions with objective function values of 46428.14, 46955.42, 46684.12, and 46540.5 at JJJJ=-31988, -31975, -31972, and -31965, respectively.
The total time for the output above was obtained in 55 hours on a personal computer with an Intel 2.66 GHz. chip and the IBM interpreter.
If one makes changes to the computer program above, such as a change of changing line 128 to, for example, 128 FOR I=1 TO 3, one may or may not obtain this objective function value of 45961.73. For this problem, the optimal objective function value is 44965.0, Amaral [1]. Generally speaking, in order to increase the probability of getting an optimal solution, one preferably should run several computers simultaneously, each with a different line 128, for example. This multicomputer approach can mitigate the possible danger of missed optimality and can produce a usable solution faster than just running one computer.
References
[1] A.R.S. Amaral, A new lower bound for the single row facility layout problem, Discrete Applied Mathematics, in press.
[2] W.C. Conley, Optimization: A Simplified Approach, Petrocelli, Princeton, New Jersey, 1981.
[3] S.S. Heragu, A. Kusiak, Efficient models for the facility layut problem, European Journal of Operational Research 53 (1991) 1-13.
[4] S.S. Heragu, Recent models and techniques for solving the layout problem, European Journal of Operational Research 57 (1992) 136-144.
[5] C.E. Nugent, T.E. Vollmann, J. Ruml, An experimental comparison of techniques for the assignment of facilities to locations, Operations Research 16 (1968) 150-173.
[6] D.M. Simmons, One-dimensional space allocation: An ordering algorithm, Operations Research 17 (1969) 812-826.
