Vorige Pagina About the Author

' 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

Vorige Pagina About the Author