' Generates Semi-Latin Squares of order 16
' Most Perfect Franklin Pan Magic Squares
' Tested with Office 365 under Windows 10
Sub MstPrfct16a()
Dim a(256), b(256), C(256), s(128)
Dim a2(256), b2(256)
''y = MsgBox("Locked", vbCritical, "Routine MstPrfct16a")
''End
n2 = 0: n9 = 0: k1 = 1: k2 = 1
m1 = 0: m2 = 15: s1 = 120: p1 = s1 / 8
' Generate Semi Latin Squares
Sheets("Klad1").Select
t1 = Timer
For j256 = m1 To m2 'a(256)
a(256) = j256
For j255 = m1 To m2 'a(255)
a(255) = j255
If a(255) = a(256) Then GoTo 2550
For j254 = m1 To m2 'a(254)
a(254) = j254
If a(254) = a(255) Or a(254) = a(256) Then GoTo 2540
a(253) = s1 / 4 - a(254) - a(255) - a(256): If a(253) < m1 Or a(253) > m2 Then GoTo 2540:
If a(253) = a(254) Or a(253) = a(255) Or a(253) = a(256) Then GoTo 2540
For j252 = m1 To m2 'a(252)
a(252) = j252
If a(252) = a(253) Or a(252) = a(254) Or a(252) = a(255) Or a(252) = a(256) Then GoTo 2520
For j251 = m1 To m2 'a(251)
a(251) = j251
If a(251) = a(252) Or a(251) = a(253) Or a(251) = a(254) Or a(251) = a(255) Or a(251) = a(256) Then GoTo 2510
a(250) = -a(252) + a(254) + a(256): If a(250) < m1 Or a(250) > m2 Then GoTo 2510:
If a(250) = a(252) Or a(250) = a(253) Or a(250) = a(254) Or a(250) = a(255) Or a(250) = a(256) Then GoTo 2510
If a(250) = a(251) Then GoTo 2510
a(249) = s1 / 4 - a(251) - a(254) - a(256): If a(249) < m1 Or a(249) > m2 Then GoTo 2510:
If a(249) = a(252) Or a(249) = a(253) Or a(249) = a(254) Or a(249) = a(255) Or a(249) = a(256) Then GoTo 2510
If a(249) = a(250) Or a(249) = a(251) Then GoTo 2510
For j248 = m1 To m2 'a(248)
a(248) = j248
If a(248) = a(252) Or a(248) = a(253) Or a(248) = a(254) Or a(248) = a(255) Or a(248) = a(256) Then GoTo 2480
If a(248) = a(249) Or a(248) = a(250) Or a(248) = a(251) Then GoTo 2480
For j247 = m1 To m2 'a(247)
a(247) = j247
If a(247) = a(252) Or a(247) = a(253) Or a(247) = a(254) Or a(247) = a(255) Or a(247) = a(256) Then GoTo 2470
If a(247) = a(248) Or a(247) = a(249) Or a(247) = a(250) Or a(247) = a(251) Then GoTo 2470
a(246) = -a(248) + a(254) + a(256): If a(246) < m1 Or a(246) > m2 Then GoTo 2470:
If a(246) = a(252) Or a(246) = a(253) Or a(246) = a(254) Or a(246) = a(255) Or a(246) = a(256) Then GoTo 2470
If a(246) = a(247) Or a(246) = a(248) Or a(246) = a(249) Or a(246) = a(250) Or a(246) = a(251) Then GoTo 2470
a(245) = s1 / 4 - a(247) - a(254) - a(256): If a(245) < m1 Or a(245) > m2 Then GoTo 2470:
If a(245) = a(252) Or a(245) = a(253) Or a(245) = a(254) Or a(245) = a(255) Or a(245) = a(256) Then GoTo 2470
If a(245) = a(247) Or a(245) = a(248) Or a(245) = a(249) Or a(245) = a(250) Or a(245) = a(251) Then GoTo 2470
If a(245) = a(246) Then GoTo 2470
For j244 = m1 To m2 'a(244)
a(244) = j244
If a(244) = a(252) Or a(244) = a(253) Or a(244) = a(254) Or a(244) = a(255) Or a(244) = a(256) Then GoTo 2440
If a(244) = a(247) Or a(244) = a(248) Or a(244) = a(249) Or a(244) = a(250) Or a(244) = a(251) Then GoTo 2440
If a(244) = a(245) Or a(244) = a(246) Then GoTo 2440
For j243 = m1 To m2 'a(243)
a(243) = j243
a(242) = -a(244) + a(254) + a(256): If a(242) < m1 Or a(242) > m2 Then GoTo 2430:
a(241) = s1 / 4 - a(243) - a(254) - a(256): If a(241) < m1 Or a(241) > m2 Then GoTo 2430:
m11 = 241: m21 = 256: GoSub 2700: If fl1 = 0 Then GoTo 2430 'Check row 1
For j240 = m1 To m2 'a(240)
a(240) = j240
a(239) = s1 / 4 - a(240) - a(255) - a(256): If a(239) < m1 Or a(239) > m2 Then GoTo 2400:
a(238) = a(240) - a(254) + a(256): If a(238) < m1 Or a(238) > m2 Then GoTo 2400:
a(237) = -a(240) + a(254) + a(255): If a(237) < m1 Or a(237) > m2 Then GoTo 2400:
a(236) = a(240) - a(252) + a(256): If a(236) < m1 Or a(236) > m2 Then GoTo 2400:
a(235) = s1 / 4 - a(240) - a(251) - a(256): If a(235) < m1 Or a(235) > m2 Then GoTo 2400:
a(234) = a(240) + a(252) - a(254): If a(234) < m1 Or a(234) > m2 Then GoTo 2400:
a(233) = -a(240) + a(251) + a(254): If a(233) < m1 Or a(233) > m2 Then GoTo 2400:
a(232) = a(240) - a(248) + a(256): If a(232) < m1 Or a(232) > m2 Then GoTo 2400:
a(231) = s1 / 4 - a(240) - a(247) - a(256): If a(231) < m1 Or a(231) > m2 Then GoTo 2400:
a(230) = a(240) + a(248) - a(254): If a(230) < m1 Or a(230) > m2 Then GoTo 2400:
a(229) = -a(240) + a(247) + a(254): If a(229) < m1 Or a(229) > m2 Then GoTo 2400:
a(228) = a(240) - a(244) + a(256): If a(228) < m1 Or a(228) > m2 Then GoTo 2400:
a(227) = s1 / 4 - a(240) - a(243) - a(256): If a(227) < m1 Or a(227) > m2 Then GoTo 2400:
a(226) = a(240) + a(244) - a(254): If a(226) < m1 Or a(226) > m2 Then GoTo 2400:
a(225) = -a(240) + a(243) + a(254): If a(225) < m1 Or a(225) > m2 Then GoTo 2400:
m11 = 225: m21 = 240: GoSub 2700: If fl1 = 0 Then GoTo 2400 'Check row 2
a(224) = s1 / 8 - a(254): If a(224) < m1 Or a(224) > m2 Then GoTo 2400:
a(223) = -s1 / 8 + a(254) + a(255) + a(256): If a(223) < m1 Or a(223) > m2 Then GoTo 2400:
a(222) = s1 / 8 - a(256): If a(222) < m1 Or a(222) > m2 Then GoTo 2400:
a(221) = s1 / 8 - a(255): If a(221) < m1 Or a(221) > m2 Then GoTo 2400:
a(220) = s1 / 8 + a(252) - a(254) - a(256): If a(220) < m1 Or a(220) > m2 Then GoTo 2400:
a(219) = -s1 / 8 + a(251) + a(254) + a(256): If a(219) < m1 Or a(219) > m2 Then GoTo 2400:
a(218) = s1 / 8 - a(252): If a(218) < m1 Or a(218) > m2 Then GoTo 2400:
a(217) = s1 / 8 - a(251): If a(217) < m1 Or a(217) > m2 Then GoTo 2400:
a(216) = s1 / 8 + a(248) - a(254) - a(256): If a(216) < m1 Or a(216) > m2 Then GoTo 2400:
a(215) = -s1 / 8 + a(247) + a(254) + a(256): If a(215) < m1 Or a(215) > m2 Then GoTo 2400:
a(214) = s1 / 8 - a(248): If a(214) < m1 Or a(214) > m2 Then GoTo 2400:
a(213) = s1 / 8 - a(247): If a(213) < m1 Or a(213) > m2 Then GoTo 2400:
a(212) = s1 / 8 + a(244) - a(254) - a(256): If a(212) < m1 Or a(212) > m2 Then GoTo 2400:
a(211) = -s1 / 8 + a(243) + a(254) + a(256): If a(211) < m1 Or a(211) > m2 Then GoTo 2400:
a(210) = s1 / 8 - a(244): If a(210) < m1 Or a(210) > m2 Then GoTo 2400:
a(209) = s1 / 8 - a(243): If a(209) < m1 Or a(209) > m2 Then GoTo 2400:
m11 = 209: m21 = 224: GoSub 2700: If fl1 = 0 Then GoTo 2400 'Check row 3
a(208) = s1 / 8 - a(240) + a(254) - a(256): If a(208) < m1 Or a(208) > m2 Then GoTo 2400:
a(207) = s1 / 8 + a(240) - a(254) - a(255): If a(207) < m1 Or a(207) > m2 Then GoTo 2400:
a(206) = s1 / 8 - a(240): If a(206) < m1 Or a(206) > m2 Then GoTo 2400:
a(205) = -s1 / 8 + a(240) + a(255) + a(256): If a(205) < m1 Or a(205) > m2 Then GoTo 2400:
a(204) = s1 / 8 - a(240) - a(252) + a(254): If a(204) < m1 Or a(204) > m2 Then GoTo 2400:
a(203) = s1 / 8 + a(240) - a(251) - a(254): If a(203) < m1 Or a(203) > m2 Then GoTo 2400:
a(202) = s1 / 8 - a(240) + a(252) - a(256): If a(202) < m1 Or a(202) > m2 Then GoTo 2400:
a(201) = -s1 / 8 + a(240) + a(251) + a(256): If a(201) < m1 Or a(201) > m2 Then GoTo 2400:
a(200) = s1 / 8 - a(240) - a(248) + a(254): If a(200) < m1 Or a(200) > m2 Then GoTo 2400:
a(199) = s1 / 8 + a(240) - a(247) - a(254): If a(199) < m1 Or a(199) > m2 Then GoTo 2400:
a(198) = s1 / 8 - a(240) + a(248) - a(256): If a(198) < m1 Or a(198) > m2 Then GoTo 2400:
a(197) = -s1 / 8 + a(240) + a(247) + a(256): If a(197) < m1 Or a(197) > m2 Then GoTo 2400:
a(196) = s1 / 8 - a(240) - a(244) + a(254): If a(196) < m1 Or a(196) > m2 Then GoTo 2400:
a(195) = s1 / 8 + a(240) - a(243) - a(254): If a(195) < m1 Or a(195) > m2 Then GoTo 2400:
a(194) = s1 / 8 - a(240) + a(244) - a(256): If a(194) < m1 Or a(194) > m2 Then GoTo 2400:
a(193) = -s1 / 8 + a(240) + a(243) + a(256): If a(193) < m1 Or a(193) > m2 Then GoTo 2400:
m11 = 193: m21 = 208: GoSub 2700: If fl1 = 0 Then GoTo 2480 'Check row 4
For j192 = m1 To m2 'a(192)
a(192) = j192
a(191) = -a(192) + a(255) + a(256): If a(191) < m1 Or a(191) > m2 Then GoTo 1920:
a(190) = a(192) + a(254) - a(256): If a(190) < m1 Or a(190) > m2 Then GoTo 1920:
a(189) = s1 / 4 - a(192) - a(254) - a(255): If a(189) < m1 Or a(189) > m2 Then GoTo 1920:
a(188) = a(192) + a(252) - a(256): If a(188) < m1 Or a(188) > m2 Then GoTo 1920:
a(187) = -a(192) + a(251) + a(256): If a(187) < m1 Or a(187) > m2 Then GoTo 1920:
a(186) = a(192) - a(252) + a(254): If a(186) < m1 Or a(186) > m2 Then GoTo 1920:
a(185) = s1 / 4 - a(192) - a(251) - a(254): If a(185) < m1 Or a(185) > m2 Then GoTo 1920:
a(184) = a(192) + a(248) - a(256): If a(184) < m1 Or a(184) > m2 Then GoTo 1920:
a(183) = -a(192) + a(247) + a(256): If a(183) < m1 Or a(183) > m2 Then GoTo 1920:
a(182) = a(192) - a(248) + a(254): If a(182) < m1 Or a(182) > m2 Then GoTo 1920:
a(181) = s1 / 4 - a(192) - a(247) - a(254): If a(181) < m1 Or a(181) > m2 Then GoTo 1920:
a(180) = a(192) + a(244) - a(256): If a(180) < m1 Or a(180) > m2 Then GoTo 1920:
a(179) = -a(192) + a(243) + a(256): If a(179) < m1 Or a(179) > m2 Then GoTo 1920:
a(178) = a(192) - a(244) + a(254): If a(178) < m1 Or a(178) > m2 Then GoTo 1920:
a(177) = s1 / 4 - a(192) - a(243) - a(254): If a(177) < m1 Or a(177) > m2 Then GoTo 1920:
m11 = 177: m21 = 192: GoSub 2700: If fl1 = 0 Then GoTo 1920 'Check row 5
For j176 = m1 To m2 'a(176)
a(176) = j176
a(175) = s1 / 4 - a(176) - a(255) - a(256): If a(175) < m1 Or a(175) > m2 Then GoTo 1760:
a(174) = a(176) - a(254) + a(256): If a(174) < m1 Or a(174) > m2 Then GoTo 1760:
a(173) = -a(176) + a(254) + a(255): If a(173) < m1 Or a(173) > m2 Then GoTo 1760:
a(172) = a(176) - a(252) + a(256): If a(172) < m1 Or a(172) > m2 Then GoTo 1760:
a(171) = s1 / 4 - a(176) - a(251) - a(256): If a(171) < m1 Or a(171) > m2 Then GoTo 1760:
a(170) = a(176) + a(252) - a(254): If a(170) < m1 Or a(170) > m2 Then GoTo 1760:
a(169) = -a(176) + a(251) + a(254): If a(169) < m1 Or a(169) > m2 Then GoTo 1760:
a(168) = a(176) - a(248) + a(256): If a(168) < m1 Or a(168) > m2 Then GoTo 1760:
a(167) = s1 / 4 - a(176) - a(247) - a(256): If a(167) < m1 Or a(167) > m2 Then GoTo 1760:
a(166) = a(176) + a(248) - a(254): If a(166) < m1 Or a(166) > m2 Then GoTo 1760:
a(165) = -a(176) + a(247) + a(254): If a(165) < m1 Or a(165) > m2 Then GoTo 1760:
a(164) = a(176) - a(244) + a(256): If a(164) < m1 Or a(164) > m2 Then GoTo 1760:
a(163) = s1 / 4 - a(176) - a(243) - a(256): If a(163) < m1 Or a(163) > m2 Then GoTo 1760:
a(162) = a(176) + a(244) - a(254): If a(162) < m1 Or a(162) > m2 Then GoTo 1760:
a(161) = -a(176) + a(243) + a(254): If a(161) < m1 Or a(161) > m2 Then GoTo 1760:
m11 = 161: m21 = 176: GoSub 2700: If fl1 = 0 Then GoTo 1760 'Check row 6
a(160) = s1 / 8 - a(192) - a(254) + a(256): If a(160) < m1 Or a(160) > m2 Then GoTo 1760:
a(159) = -s1 / 8 + a(192) + a(254) + a(255): If a(159) < m1 Or a(159) > m2 Then GoTo 1760:
a(158) = s1 / 8 - a(192): If a(158) < m1 Or a(158) > m2 Then GoTo 1760:
a(157) = s1 / 8 + a(192) - a(255) - a(256): If a(157) < m1 Or a(157) > m2 Then GoTo 1760:
a(156) = s1 / 8 - a(192) + a(252) - a(254): If a(156) < m1 Or a(156) > m2 Then GoTo 1760:
a(155) = -s1 / 8 + a(192) + a(251) + a(254): If a(155) < m1 Or a(155) > m2 Then GoTo 1760:
a(154) = s1 / 8 - a(192) - a(252) + a(256): If a(154) < m1 Or a(154) > m2 Then GoTo 1760:
a(153) = s1 / 8 + a(192) - a(251) - a(256): If a(153) < m1 Or a(153) > m2 Then GoTo 1760:
a(152) = s1 / 8 - a(192) + a(248) - a(254): If a(152) < m1 Or a(152) > m2 Then GoTo 1760:
a(151) = -s1 / 8 + a(192) + a(247) + a(254): If a(151) < m1 Or a(151) > m2 Then GoTo 1760:
a(150) = s1 / 8 - a(192) - a(248) + a(256): If a(150) < m1 Or a(150) > m2 Then GoTo 1760:
a(149) = s1 / 8 + a(192) - a(247) - a(256): If a(149) < m1 Or a(149) > m2 Then GoTo 1760:
a(148) = s1 / 8 - a(192) + a(244) - a(254): If a(148) < m1 Or a(148) > m2 Then GoTo 1760:
a(147) = -s1 / 8 + a(192) + a(243) + a(254): If a(147) < m1 Or a(147) > m2 Then GoTo 1760:
a(146) = s1 / 8 - a(192) - a(244) + a(256): If a(146) < m1 Or a(146) > m2 Then GoTo 1760:
a(145) = s1 / 8 + a(192) - a(243) - a(256): If a(145) < m1 Or a(145) > m2 Then GoTo 1760:
m11 = 145: m21 = 160: GoSub 2700: If fl1 = 0 Then GoTo 1760 'Check row 7
a(144) = s1 / 8 - a(176) + a(254) - a(256): If a(144) < m1 Or a(144) > m2 Then GoTo 1760:
a(143) = s1 / 8 + a(176) - a(254) - a(255): If a(143) < m1 Or a(143) > m2 Then GoTo 1760:
a(142) = s1 / 8 - a(176): If a(142) < m1 Or a(142) > m2 Then GoTo 1760:
a(141) = -s1 / 8 + a(176) + a(255) + a(256): If a(141) < m1 Or a(141) > m2 Then GoTo 1760:
a(140) = s1 / 8 - a(176) - a(252) + a(254): If a(140) < m1 Or a(140) > m2 Then GoTo 1760:
a(139) = s1 / 8 + a(176) - a(251) - a(254): If a(139) < m1 Or a(139) > m2 Then GoTo 1760:
a(138) = s1 / 8 - a(176) + a(252) - a(256): If a(138) < m1 Or a(138) > m2 Then GoTo 1760:
a(137) = -s1 / 8 + a(176) + a(251) + a(256): If a(137) < m1 Or a(137) > m2 Then GoTo 1760:
a(136) = s1 / 8 - a(176) - a(248) + a(254): If a(136) < m1 Or a(136) > m2 Then GoTo 1760:
a(135) = s1 / 8 + a(176) - a(247) - a(254): If a(135) < m1 Or a(135) > m2 Then GoTo 1760:
a(134) = s1 / 8 - a(176) + a(248) - a(256): If a(134) < m1 Or a(134) > m2 Then GoTo 1760:
a(133) = -s1 / 8 + a(176) + a(247) + a(256): If a(133) < m1 Or a(133) > m2 Then GoTo 1760:
a(132) = s1 / 8 - a(176) - a(244) + a(254): If a(132) < m1 Or a(132) > m2 Then GoTo 1760:
a(131) = s1 / 8 + a(176) - a(243) - a(254): If a(131) < m1 Or a(131) > m2 Then GoTo 1760:
a(130) = s1 / 8 - a(176) + a(244) - a(256): If a(130) < m1 Or a(130) > m2 Then GoTo 1760:
a(129) = -s1 / 8 + a(176) + a(243) + a(256): If a(129) < m1 Or a(129) > m2 Then GoTo 1760:
m11 = 129: m21 = 144: GoSub 2700: If fl1 = 0 Then GoTo 1760 'Check row 8
For j128 = m1 To m2 'a(128)
a(128) = j128
a(127) = -a(128) + a(255) + a(256): If a(127) < m1 Or a(127) > m2 Then GoTo 1280:
a(126) = a(128) + a(254) - a(256): If a(126) < m1 Or a(126) > m2 Then GoTo 1280:
a(125) = s1 / 4 - a(128) - a(254) - a(255): If a(125) < m1 Or a(125) > m2 Then GoTo 1280:
a(124) = a(128) + a(252) - a(256): If a(124) < m1 Or a(124) > m2 Then GoTo 1280:
a(123) = -a(128) + a(251) + a(256): If a(123) < m1 Or a(123) > m2 Then GoTo 1280:
a(122) = a(128) - a(252) + a(254): If a(122) < m1 Or a(122) > m2 Then GoTo 1280:
a(121) = s1 / 4 - a(128) - a(251) - a(254): If a(121) < m1 Or a(121) > m2 Then GoTo 1280:
a(120) = a(128) + a(248) - a(256): If a(120) < m1 Or a(120) > m2 Then GoTo 1280:
a(119) = -a(128) + a(247) + a(256): If a(119) < m1 Or a(119) > m2 Then GoTo 1280:
a(118) = a(128) - a(248) + a(254): If a(118) < m1 Or a(118) > m2 Then GoTo 1280:
a(117) = s1 / 4 - a(128) - a(247) - a(254): If a(117) < m1 Or a(117) > m2 Then GoTo 1280:
a(116) = a(128) + a(244) - a(256): If a(116) < m1 Or a(116) > m2 Then GoTo 1280:
a(115) = -a(128) + a(243) + a(256): If a(115) < m1 Or a(115) > m2 Then GoTo 1280:
a(114) = a(128) - a(244) + a(254): If a(114) < m1 Or a(114) > m2 Then GoTo 1280:
a(113) = s1 / 4 - a(128) - a(243) - a(254): If a(113) < m1 Or a(113) > m2 Then GoTo 1280:
m11 = 113: m21 = 128: GoSub 2700: If fl1 = 0 Then GoTo 1280 'Check row 9
For j112 = m1 To m2 'a(112)
a(112) = j112
a(111) = s1 / 4 - a(112) - a(255) - a(256): If a(111) < m1 Or a(111) > m2 Then GoTo 1120:
a(110) = a(112) - a(254) + a(256): If a(110) < m1 Or a(110) > m2 Then GoTo 1120:
a(109) = -a(112) + a(254) + a(255): If a(109) < m1 Or a(109) > m2 Then GoTo 1120:
a(108) = a(112) - a(252) + a(256): If a(108) < m1 Or a(108) > m2 Then GoTo 1120:
a(107) = s1 / 4 - a(112) - a(251) - a(256): If a(107) < m1 Or a(107) > m2 Then GoTo 1120:
a(106) = a(112) + a(252) - a(254): If a(106) < m1 Or a(106) > m2 Then GoTo 1120:
a(105) = -a(112) + a(251) + a(254): If a(105) < m1 Or a(105) > m2 Then GoTo 1120:
a(104) = a(112) - a(248) + a(256): If a(104) < m1 Or a(104) > m2 Then GoTo 1120:
a(103) = s1 / 4 - a(112) - a(247) - a(256): If a(103) < m1 Or a(103) > m2 Then GoTo 1120:
a(102) = a(112) + a(248) - a(254): If a(102) < m1 Or a(102) > m2 Then GoTo 1120:
a(101) = -a(112) + a(247) + a(254): If a(101) < m1 Or a(101) > m2 Then GoTo 1120:
a(100) = a(112) - a(244) + a(256): If a(100) < m1 Or a(100) > m2 Then GoTo 1120:
a(99) = s1 / 4 - a(112) - a(243) - a(256): If a(99) < m1 Or a(99) > m2 Then GoTo 1120:
a(98) = a(112) + a(244) - a(254): If a(98) < m1 Or a(98) > m2 Then GoTo 1120:
a(97) = -a(112) + a(243) + a(254): If a(97) < m1 Or a(97) > m2 Then GoTo 1120:
m11 = 97: m21 = 112: GoSub 2700: If fl1 = 0 Then GoTo 1120 'Check row 10
a(96) = s1 / 8 - a(128) - a(254) + a(256): If a(96) < m1 Or a(96) > m2 Then GoTo 1120:
a(95) = -s1 / 8 + a(128) + a(254) + a(255): If a(95) < m1 Or a(95) > m2 Then GoTo 1120:
a(94) = s1 / 8 - a(128): If a(94) < m1 Or a(94) > m2 Then GoTo 1120:
a(93) = s1 / 8 + a(128) - a(255) - a(256): If a(93) < m1 Or a(93) > m2 Then GoTo 1120:
a(92) = s1 / 8 - a(128) + a(252) - a(254): If a(92) < m1 Or a(92) > m2 Then GoTo 1120:
a(91) = -s1 / 8 + a(128) + a(251) + a(254): If a(91) < m1 Or a(91) > m2 Then GoTo 1120:
a(90) = s1 / 8 - a(128) - a(252) + a(256): If a(90) < m1 Or a(90) > m2 Then GoTo 1120:
a(89) = s1 / 8 + a(128) - a(251) - a(256): If a(89) < m1 Or a(89) > m2 Then GoTo 1120:
a(88) = s1 / 8 - a(128) + a(248) - a(254): If a(88) < m1 Or a(88) > m2 Then GoTo 1120:
a(87) = -s1 / 8 + a(128) + a(247) + a(254): If a(87) < m1 Or a(87) > m2 Then GoTo 1120:
a(86) = s1 / 8 - a(128) - a(248) + a(256): If a(86) < m1 Or a(86) > m2 Then GoTo 1120:
a(85) = s1 / 8 + a(128) - a(247) - a(256): If a(85) < m1 Or a(85) > m2 Then GoTo 1120:
a(84) = s1 / 8 - a(128) + a(244) - a(254): If a(84) < m1 Or a(84) > m2 Then GoTo 1120:
a(83) = -s1 / 8 + a(128) + a(243) + a(254): If a(83) < m1 Or a(83) > m2 Then GoTo 1120:
a(82) = s1 / 8 - a(128) - a(244) + a(256): If a(82) < m1 Or a(82) > m2 Then GoTo 1120:
a(81) = s1 / 8 + a(128) - a(243) - a(256): If a(81) < m1 Or a(81) > m2 Then GoTo 1120:
m11 = 81: m21 = 96: GoSub 2700: If fl1 = 0 Then GoTo 1120 'Check row 11
a(80) = s1 / 8 - a(112) + a(254) - a(256): If a(80) < m1 Or a(80) > m2 Then GoTo 1120:
a(79) = s1 / 8 + a(112) - a(254) - a(255): If a(79) < m1 Or a(79) > m2 Then GoTo 1120:
a(78) = s1 / 8 - a(112): If a(78) < m1 Or a(78) > m2 Then GoTo 1120:
a(77) = -s1 / 8 + a(112) + a(255) + a(256): If a(77) < m1 Or a(77) > m2 Then GoTo 1120:
a(76) = s1 / 8 - a(112) - a(252) + a(254): If a(76) < m1 Or a(76) > m2 Then GoTo 1120:
a(75) = s1 / 8 + a(112) - a(251) - a(254): If a(75) < m1 Or a(75) > m2 Then GoTo 1120:
a(74) = s1 / 8 - a(112) + a(252) - a(256): If a(74) < m1 Or a(74) > m2 Then GoTo 1120:
a(73) = -s1 / 8 + a(112) + a(251) + a(256): If a(73) < m1 Or a(73) > m2 Then GoTo 1120:
a(72) = s1 / 8 - a(112) - a(248) + a(254): If a(72) < m1 Or a(72) > m2 Then GoTo 1120:
a(71) = s1 / 8 + a(112) - a(247) - a(254): If a(71) < m1 Or a(71) > m2 Then GoTo 1120:
a(70) = s1 / 8 - a(112) + a(248) - a(256): If a(70) < m1 Or a(70) > m2 Then GoTo 1120:
a(69) = -s1 / 8 + a(112) + a(247) + a(256): If a(69) < m1 Or a(69) > m2 Then GoTo 1120:
a(68) = s1 / 8 - a(112) - a(244) + a(254): If a(68) < m1 Or a(68) > m2 Then GoTo 1120:
a(67) = s1 / 8 + a(112) - a(243) - a(254): If a(67) < m1 Or a(67) > m2 Then GoTo 1120:
a(66) = s1 / 8 - a(112) + a(244) - a(256): If a(66) < m1 Or a(66) > m2 Then GoTo 1120:
a(65) = -s1 / 8 + a(112) + a(243) + a(256): If a(65) < m1 Or a(65) > m2 Then GoTo 1120:
m11 = 65: m21 = 80: GoSub 2700: If fl1 = 0 Then GoTo 1120 'Check row 12
For j64 = m1 To m2 'a(64)
a(64) = j64
a(63) = -a(64) + a(255) + a(256): If a(63) < m1 Or a(63) > m2 Then GoTo 640:
a(62) = a(64) + a(254) - a(256): If a(62) < m1 Or a(62) > m2 Then GoTo 640:
a(61) = s1 / 4 - a(64) - a(254) - a(255): If a(61) < m1 Or a(61) > m2 Then GoTo 640:
a(60) = a(64) + a(252) - a(256): If a(60) < m1 Or a(60) > m2 Then GoTo 640:
a(59) = -a(64) + a(251) + a(256): If a(59) < m1 Or a(59) > m2 Then GoTo 640:
a(58) = a(64) - a(252) + a(254): If a(58) < m1 Or a(58) > m2 Then GoTo 640:
a(57) = s1 / 4 - a(64) - a(251) - a(254): If a(57) < m1 Or a(57) > m2 Then GoTo 640:
a(56) = a(64) + a(248) - a(256): If a(56) < m1 Or a(56) > m2 Then GoTo 640:
a(55) = -a(64) + a(247) + a(256): If a(55) < m1 Or a(55) > m2 Then GoTo 640:
a(54) = a(64) - a(248) + a(254): If a(54) < m1 Or a(54) > m2 Then GoTo 640:
a(53) = s1 / 4 - a(64) - a(247) - a(254): If a(53) < m1 Or a(53) > m2 Then GoTo 640:
a(52) = a(64) + a(244) - a(256): If a(52) < m1 Or a(52) > m2 Then GoTo 640:
a(51) = -a(64) + a(243) + a(256): If a(51) < m1 Or a(51) > m2 Then GoTo 640:
a(50) = a(64) - a(244) + a(254): If a(50) < m1 Or a(50) > m2 Then GoTo 640:
a(49) = s1 / 4 - a(64) - a(243) - a(254): If a(49) < m1 Or a(49) > m2 Then GoTo 640:
m11 = 49: m21 = 64: GoSub 2700: If fl1 = 0 Then GoTo 640 'Check row 13
For j48 = m1 To m2 'a(48)
a(48) = j48
a(47) = s1 / 4 - a(48) - a(255) - a(256): If a(47) < m1 Or a(47) > m2 Then GoTo 480:
a(46) = a(48) - a(254) + a(256): If a(46) < m1 Or a(46) > m2 Then GoTo 480:
a(45) = -a(48) + a(254) + a(255): If a(45) < m1 Or a(45) > m2 Then GoTo 480:
a(44) = a(48) - a(252) + a(256): If a(44) < m1 Or a(44) > m2 Then GoTo 480:
a(43) = s1 / 4 - a(48) - a(251) - a(256): If a(43) < m1 Or a(43) > m2 Then GoTo 480:
a(42) = a(48) + a(252) - a(254): If a(42) < m1 Or a(42) > m2 Then GoTo 480:
a(41) = -a(48) + a(251) + a(254): If a(41) < m1 Or a(41) > m2 Then GoTo 480:
a(40) = a(48) - a(248) + a(256): If a(40) < m1 Or a(40) > m2 Then GoTo 480:
a(39) = s1 / 4 - a(48) - a(247) - a(256): If a(39) < m1 Or a(39) > m2 Then GoTo 480:
a(38) = a(48) + a(248) - a(254): If a(38) < m1 Or a(38) > m2 Then GoTo 480:
a(37) = -a(48) + a(247) + a(254): If a(37) < m1 Or a(37) > m2 Then GoTo 480:
a(36) = a(48) - a(244) + a(256): If a(36) < m1 Or a(36) > m2 Then GoTo 480:
a(35) = s1 / 4 - a(48) - a(243) - a(256): If a(35) < m1 Or a(35) > m2 Then GoTo 480:
a(34) = a(48) + a(244) - a(254): If a(34) < m1 Or a(34) > m2 Then GoTo 480:
a(33) = -a(48) + a(243) + a(254): If a(33) < m1 Or a(33) > m2 Then GoTo 480:
m11 = 33: m21 = 48: GoSub 2700: If fl1 = 0 Then GoTo 480 'Check row 14
a(32) = s1 / 8 - a(64) - a(254) + a(256): If a(32) < m1 Or a(32) > m2 Then GoTo 480:
a(31) = -s1 / 8 + a(64) + a(254) + a(255): If a(31) < m1 Or a(31) > m2 Then GoTo 480:
a(30) = s1 / 8 - a(64): If a(30) < m1 Or a(30) > m2 Then GoTo 480:
a(29) = s1 / 8 + a(64) - a(255) - a(256): If a(29) < m1 Or a(29) > m2 Then GoTo 480:
a(28) = s1 / 8 - a(64) + a(252) - a(254): If a(28) < m1 Or a(28) > m2 Then GoTo 480:
a(27) = -s1 / 8 + a(64) + a(251) + a(254): If a(27) < m1 Or a(27) > m2 Then GoTo 480:
a(26) = s1 / 8 - a(64) - a(252) + a(256): If a(26) < m1 Or a(26) > m2 Then GoTo 480:
a(25) = s1 / 8 + a(64) - a(251) - a(256): If a(25) < m1 Or a(25) > m2 Then GoTo 480:
a(24) = s1 / 8 - a(64) + a(248) - a(254): If a(24) < m1 Or a(24) > m2 Then GoTo 480:
a(23) = -s1 / 8 + a(64) + a(247) + a(254): If a(23) < m1 Or a(23) > m2 Then GoTo 480:
a(22) = s1 / 8 - a(64) - a(248) + a(256): If a(22) < m1 Or a(22) > m2 Then GoTo 480:
a(21) = s1 / 8 + a(64) - a(247) - a(256): If a(21) < m1 Or a(21) > m2 Then GoTo 480:
a(20) = s1 / 8 - a(64) + a(244) - a(254): If a(20) < m1 Or a(20) > m2 Then GoTo 480:
a(19) = -s1 / 8 + a(64) + a(243) + a(254): If a(19) < m1 Or a(19) > m2 Then GoTo 480:
a(18) = s1 / 8 - a(64) - a(244) + a(256): If a(18) < m1 Or a(18) > m2 Then GoTo 480:
a(17) = s1 / 8 + a(64) - a(243) - a(256): If a(17) < m1 Or a(17) > m2 Then GoTo 480:
m11 = 17: m21 = 32: GoSub 2700: If fl1 = 0 Then GoTo 480 'Check row 15
a(16) = s1 / 8 - a(48) + a(254) - a(256): If a(16) < m1 Or a(16) > m2 Then GoTo 480:
a(15) = s1 / 8 + a(48) - a(254) - a(255): If a(15) < m1 Or a(15) > m2 Then GoTo 480:
a(14) = s1 / 8 - a(48): If a(14) < m1 Or a(14) > m2 Then GoTo 480:
a(13) = -s1 / 8 + a(48) + a(255) + a(256): If a(13) < m1 Or a(13) > m2 Then GoTo 480:
a(12) = s1 / 8 - a(48) - a(252) + a(254): If a(12) < m1 Or a(12) > m2 Then GoTo 480:
a(11) = s1 / 8 + a(48) - a(251) - a(254): If a(11) < m1 Or a(11) > m2 Then GoTo 480:
a(10) = s1 / 8 - a(48) + a(252) - a(256): If a(10) < m1 Or a(10) > m2 Then GoTo 480:
a(9) = -s1 / 8 + a(48) + a(251) + a(256): If a(9) < m1 Or a(9) > m2 Then GoTo 480:
a(8) = s1 / 8 - a(48) - a(248) + a(254): If a(8) < m1 Or a(8) > m2 Then GoTo 480:
a(7) = s1 / 8 + a(48) - a(247) - a(254): If a(7) < m1 Or a(7) > m2 Then GoTo 480:
a(6) = s1 / 8 - a(48) + a(248) - a(256): If a(6) < m1 Or a(6) > m2 Then GoTo 480:
a(5) = -s1 / 8 + a(48) + a(247) + a(256): If a(5) < m1 Or a(5) > m2 Then GoTo 480:
a(4) = s1 / 8 - a(48) - a(244) + a(254): If a(4) < m1 Or a(4) > m2 Then GoTo 480:
a(3) = s1 / 8 + a(48) - a(243) - a(254): If a(3) < m1 Or a(3) > m2 Then GoTo 480:
a(2) = s1 / 8 - a(48) + a(244) - a(256): If a(2) < m1 Or a(2) > m2 Then GoTo 480:
a(1) = -s1 / 8 + a(48) + a(243) + a(256): If a(1) < m1 Or a(1) > m2 Then GoTo 480:
m11 = 1: m21 = 16: GoSub 2700: If fl1 = 0 Then GoTo 480 'Check row 16
GoSub 2750: If fl1 = 0 Then GoTo 480 'Check diagonal 1, 2
' Calculate c() = 16 * a2() + b2() + 1
Erase a2
For i1 = 1 To 256: a2(i1) = a(i1): Next i1
GoSub 1500: If fl1 = 0 Then GoTo 480
n9 = n9 + 1 'Total
' GoSub 2600 'Print results (Lines)
GoSub 2650 'Print results (Squares)
' Cells(2, 1).Value = n9 'Counting
480 Next j48
640 Next j64
1120 Next j112
1280 Next j128
1760 Next j176
1920 Next j192
2400 Next j240
2430 Next j243
2440 Next j244
2470 Next j247
2480 Next j248
2510 Next j251
2520 Next j252
2540 Next j254
2550 Next j255
2560 Next j256
t2 = Timer
t10 = Str(t2 - t1) + " sec., " + Str(n9) + " Solutions for sum" + Str(s1)
y = MsgBox(t10, 0, "Routine MstPrfct16a")
End
' Construct Most Perfect Franklin Pan Magic (16 x 16)
1500 fl1 = 1
' Transpose a2()
b2(1) = a2(1): b2(2) = a2(17): b2(3) = a2(33): b2(4) = a2(49): b2(5) = a2(65): b2(6) = a2(81):
b2(7) = a2(97): b2(8) = a2(113): b2(9) = a2(129): b2(10) = a2(145): b2(11) = a2(161): b2(12) = a2(177):
b2(13) = a2(193): b2(14) = a2(209): b2(15) = a2(225): b2(16) = a2(241):
b2(17) = a2(2): b2(18) = a2(18): b2(19) = a2(34): b2(20) = a2(50): b2(21) = a2(66): b2(22) = a2(82):
b2(23) = a2(98): b2(24) = a2(114): b2(25) = a2(130): b2(26) = a2(146): b2(27) = a2(162): b2(28) = a2(178):
b2(29) = a2(194): b2(30) = a2(210): b2(31) = a2(226): b2(32) = a2(242):
b2(33) = a2(3): b2(34) = a2(19): b2(35) = a2(35): b2(36) = a2(51): b2(37) = a2(67): b2(38) = a2(83):
b2(39) = a2(99): b2(40) = a2(115): b2(41) = a2(131): b2(42) = a2(147): b2(43) = a2(163): b2(44) = a2(179):
b2(45) = a2(195): b2(46) = a2(211): b2(47) = a2(227): b2(48) = a2(243):
b2(49) = a2(4): b2(50) = a2(20): b2(51) = a2(36): b2(52) = a2(52): b2(53) = a2(68): b2(54) = a2(84):
b2(55) = a2(100): b2(56) = a2(116): b2(57) = a2(132): b2(58) = a2(148): b2(59) = a2(164): b2(60) = a2(180):
b2(61) = a2(196): b2(62) = a2(212): b2(63) = a2(228): b2(64) = a2(244):
b2(65) = a2(5): b2(66) = a2(21): b2(67) = a2(37): b2(68) = a2(53): b2(69) = a2(69): b2(70) = a2(85):
b2(71) = a2(101): b2(72) = a2(117): b2(73) = a2(133): b2(74) = a2(149): b2(75) = a2(165): b2(76) = a2(181):
b2(77) = a2(197): b2(78) = a2(213): b2(79) = a2(229): b2(80) = a2(245):
b2(81) = a2(6): b2(82) = a2(22): b2(83) = a2(38): b2(84) = a2(54): b2(85) = a2(70): b2(86) = a2(86):
b2(87) = a2(102): b2(88) = a2(118): b2(89) = a2(134): b2(90) = a2(150): b2(91) = a2(166): b2(92) = a2(182):
b2(93) = a2(198): b2(94) = a2(214): b2(95) = a2(230): b2(96) = a2(246):
b2(97) = a2(7): b2(98) = a2(23): b2(99) = a2(39): b2(100) = a2(55): b2(101) = a2(71): b2(102) = a2(87):
b2(103) = a2(103): b2(104) = a2(119): b2(105) = a2(135): b2(106) = a2(151): b2(107) = a2(167): b2(108) = a2(183):
b2(109) = a2(199): b2(110) = a2(215): b2(111) = a2(231): b2(112) = a2(247):
b2(113) = a2(8): b2(114) = a2(24): b2(115) = a2(40): b2(116) = a2(56): b2(117) = a2(72): b2(118) = a2(88):
b2(119) = a2(104): b2(120) = a2(120): b2(121) = a2(136): b2(122) = a2(152): b2(123) = a2(168): b2(124) = a2(184):
b2(125) = a2(200): b2(126) = a2(216): b2(127) = a2(232): b2(128) = a2(248):
b2(129) = a2(9): b2(130) = a2(25): b2(131) = a2(41): b2(132) = a2(57): b2(133) = a2(73): b2(134) = a2(89):
b2(135) = a2(105): b2(136) = a2(121): b2(137) = a2(137): b2(138) = a2(153): b2(139) = a2(169): b2(140) = a2(185):
b2(141) = a2(201): b2(142) = a2(217): b2(143) = a2(233): b2(144) = a2(249):
b2(145) = a2(10): b2(146) = a2(26): b2(147) = a2(42): b2(148) = a2(58): b2(149) = a2(74): b2(150) = a2(90):
b2(151) = a2(106): b2(152) = a2(122): b2(153) = a2(138): b2(154) = a2(154): b2(155) = a2(170): b2(156) = a2(186):
b2(157) = a2(202): b2(158) = a2(218): b2(159) = a2(234): b2(160) = a2(250):
b2(161) = a2(11): b2(162) = a2(27): b2(163) = a2(43): b2(164) = a2(59): b2(165) = a2(75): b2(166) = a2(91):
b2(167) = a2(107): b2(168) = a2(123): b2(169) = a2(139): b2(170) = a2(155): b2(171) = a2(171): b2(172) = a2(187):
b2(173) = a2(203): b2(174) = a2(219): b2(175) = a2(235): b2(176) = a2(251):
b2(177) = a2(12): b2(178) = a2(28): b2(179) = a2(44): b2(180) = a2(60): b2(181) = a2(76): b2(182) = a2(92):
b2(183) = a2(108): b2(184) = a2(124): b2(185) = a2(140): b2(186) = a2(156): b2(187) = a2(172): b2(188) = a2(188):
b2(189) = a2(204): b2(190) = a2(220): b2(191) = a2(236): b2(192) = a2(252):
b2(193) = a2(13): b2(194) = a2(29): b2(195) = a2(45): b2(196) = a2(61): b2(197) = a2(77): b2(198) = a2(93):
b2(199) = a2(109): b2(200) = a2(125): b2(201) = a2(141): b2(202) = a2(157): b2(203) = a2(173): b2(204) = a2(189):
b2(205) = a2(205): b2(206) = a2(221): b2(207) = a2(237): b2(208) = a2(253):
b2(209) = a2(14): b2(210) = a2(30): b2(211) = a2(46): b2(212) = a2(62): b2(213) = a2(78): b2(214) = a2(94):
b2(215) = a2(110): b2(216) = a2(126): b2(217) = a2(142): b2(218) = a2(158): b2(219) = a2(174): b2(220) = a2(190):
b2(221) = a2(206): b2(222) = a2(222): b2(223) = a2(238): b2(224) = a2(254):
b2(225) = a2(15): b2(226) = a2(31): b2(227) = a2(47): b2(228) = a2(63): b2(229) = a2(79): b2(230) = a2(95):
b2(231) = a2(111): b2(232) = a2(127): b2(233) = a2(143): b2(234) = a2(159): b2(235) = a2(175): b2(236) = a2(191):
b2(237) = a2(207): b2(238) = a2(223): b2(239) = a2(239): b2(240) = a2(255):
b2(241) = a2(16): b2(242) = a2(32): b2(243) = a2(48): b2(244) = a2(64): b2(245) = a2(80): b2(246) = a2(96):
b2(247) = a2(112): b2(248) = a2(128): b2(249) = a2(144): b2(250) = a2(160): b2(251) = a2(176): b2(252) = a2(192):
b2(253) = a2(208): b2(254) = a2(224): b2(255) = a2(240): b2(256) = a2(256):
Erase C
For i1 = 1 To 256
C(i1) = 16 * a2(i1) + b2(i1) + 1
Next i1
fl1 = 1: n20 = 0
For j1 = 1 To 256
a20 = C(j1):
For j2 = (1 + j1) To 256
If a20 = C(j2) Then fl1 = 0: Return
Next j2
Next j1
Return
' Print Results (lines)
2600
For i1 = 1 To 256
Cells(n9, i1).Value = a(i1)
Next i1
Cells(n9, 257).Value = n9
Cells(1, 259).Value = n9
Return
' Print results (squares)
2650 n2 = n2 + 1
If n2 = 3 Then
n2 = 1: k1 = k1 + 17: k2 = 1
Else
If n9 > 1 Then k2 = k2 + 17
End If
Cells(k1, k2 + 1).Font.Color = -4165632
Cells(k1, k2 + 1).Value = n9
i3 = 0
For i1 = 1 To 16
For i2 = 1 To 16
i3 = i3 + 1
Cells(k1 + i1, k2 + i2).Value = C(i3) ''a(i3)
Next i2
Next i1
Return
' Check Row i (i = 1 ... 16)
2700 fl1 = 1
For i1 = m11 To m21: b(i1 - m11 + 1) = a(i1): Next i1
For j1 = 1 To 16
a20 = b(j1):
For j2 = (1 + j1) To 16
If a20 = b(j2) Then fl1 = 0: Return
Next j2
Next j1
Return
' Check Diagonals 1, 2
2750 fl1 = 1
i2 = 1
For i1 = 1 To 16:
b(i1) = a(i2): i2 = i2 + 17
Next i1
For j1 = 1 To 16
a20 = b(j1):
For j2 = (1 + j1) To 16
If a20 = b(j2) Then fl1 = 0: Return
Next j2
Next j1
i2 = 16
For i1 = 1 To 16:
b(i1) = a(i2): i2 = i2 + 15
Next i1
For j1 = 1 To 16
a20 = b(j1):
For j2 = (1 + j1) To 16
If a20 = b(j2) Then fl1 = 0: Return
Next j2
Next j1
Return
End Sub
