' Generates Semi-Latin Squares of order 16
' Most Perfect Franklin Pan Magic Squares, Barink Restrictions
' Tested with Office 365 under Windows 10
Sub MstPrfct16b()
Dim a(256), b(256), C(256), s(128)
Dim a2(256), b2(256)
y = MsgBox("Locked", vbCritical, "Routine MstPrfct16b")
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
a(251) = -a(252) + a(255) + a(256): If a(251) < m1 Or a(251) > m2 Then GoTo 2520:
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 2520
a(250) = -a(252) + a(254) + a(256): If a(250) < m1 Or a(250) > m2 Then GoTo 2520:
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 2520
If a(250) = a(251) Then GoTo 2520
a(249) = s1 / 4 + a(252) - a(254) - a(255) - 2 * a(256): If a(249) < m1 Or a(249) > m2 Then GoTo 2520:
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 2520
If a(249) = a(250) Or a(249) = a(251) Then GoTo 2520
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
a(247) = -a(248) + a(255) + a(256): If a(247) < m1 Or a(247) > m2 Then GoTo 2480:
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 2480
If a(247) = a(248) Or a(247) = a(249) Or a(247) = a(250) Or a(247) = a(251) Then GoTo 2480
a(246) = -a(248) + a(254) + a(256): If a(246) < m1 Or a(246) > m2 Then GoTo 2480:
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 2480
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 2480
a(245) = s1 / 4 + a(248) - a(254) - a(255) - 2 * a(256): If a(245) < m1 Or a(245) > m2 Then GoTo 2480:
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 2480
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 2480
If a(245) = a(246) Then GoTo 2480
For j244 = m1 To m2 'a(244)
a(244) = j244
a(243) = -a(244) + a(255) + a(256): If a(243) < m1 Or a(243) > m2 Then GoTo 2440:
a(242) = -a(244) + a(254) + a(256): If a(242) < m1 Or a(242) > m2 Then GoTo 2440:
a(241) = s1 / 4 + a(244) - a(254) - a(255) - 2 * a(256): If a(241) < m1 Or a(241) > m2 Then GoTo 2440:
m11 = 241: m21 = 256: GoSub 2700: If fl1 = 0 Then GoTo 2440 '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(252) - a(255) - 2 * 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(252) + a(254) + a(255) + a(256): 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(248) - a(255) - 2 * 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(248) + a(254) + a(255) + a(256): 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(244) - a(255) - 2 * 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(244) + a(254) + a(255) + a(256): 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(252) + a(254) + a(255) + 2 * 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(252) - a(255) - a(256): 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(248) + a(254) + a(255) + 2 * 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(248) - a(255) - a(256): 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(244) + a(254) + a(255) + 2 * 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(244) - a(255) - a(256): 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(252) - a(254) - a(255) - a(256): 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(252) + a(255) + 2 * 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(248) - a(254) - a(255) - a(256): 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(248) + a(255) + 2 * 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(244) - a(254) - a(255) - a(256): 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(244) + a(255) + 2 * a(256): If a(193) < m1 Or a(193) > m2 Then GoTo 2400:
m11 = 193: m21 = 208: GoSub 2700: If fl1 = 0 Then GoTo 2400 '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(252) + a(255) + 2 * 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(252) - a(254) - a(255) - a(256): 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(248) + a(255) + 2 * 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(248) - a(254) - a(255) - a(256): 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(244) + a(255) + 2 * 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(244) - a(254) - a(255) - a(256): 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
a(176) = -a(192) + a(240) + a(256): If a(176) < m1 Or a(176) > m2 Then GoTo 1920:
a(175) = s1 / 4 + a(192) - a(240) - a(255) - 2 * a(256): If a(175) < m1 Or a(175) > m2 Then GoTo 1920:
a(174) = -a(192) + a(240) - a(254) + 2 * a(256): If a(174) < m1 Or a(174) > m2 Then GoTo 1920:
a(173) = a(192) - a(240) + a(254) + a(255) - a(256): If a(173) < m1 Or a(173) > m2 Then GoTo 1920:
a(172) = -a(192) + a(240) - a(252) + 2 * a(256): If a(172) < m1 Or a(172) > m2 Then GoTo 1920:
a(171) = s1 / 4 + a(192) - a(240) + a(252) - a(255) - 3 * a(256):
If a(171) < m1 Or a(171) > m2 Then GoTo 1920:
a(170) = -a(192) + a(240) + a(252) - a(254) + a(256): If a(170) < m1 Or a(170) > m2 Then GoTo 1920:
a(169) = a(192) - a(240) - a(252) + a(254) + a(255): If a(169) < m1 Or a(169) > m2 Then GoTo 1920:
a(168) = -a(192) + a(240) - a(248) + 2 * a(256): If a(168) < m1 Or a(168) > m2 Then GoTo 1920:
a(167) = s1 / 4 + a(192) - a(240) + a(248) - a(255) - 3 * a(256):
If a(167) < m1 Or a(167) > m2 Then GoTo 1920:
a(166) = -a(192) + a(240) + a(248) - a(254) + a(256): If a(166) < m1 Or a(166) > m2 Then GoTo 1920:
a(165) = a(192) - a(240) - a(248) + a(254) + a(255): If a(165) < m1 Or a(165) > m2 Then GoTo 1920:
a(164) = -a(192) + a(240) - a(244) + 2 * a(256): If a(164) < m1 Or a(164) > m2 Then GoTo 1920:
a(163) = s1 / 4 + a(192) - a(240) + a(244) - a(255) - 3 * a(256):
If a(163) < m1 Or a(163) > m2 Then GoTo 1920:
a(162) = -a(192) + a(240) + a(244) - a(254) + a(256): If a(162) < m1 Or a(162) > m2 Then GoTo 1920:
a(161) = a(192) - a(240) - a(244) + a(254) + a(255): If a(161) < m1 Or a(161) > m2 Then GoTo 1920:
m11 = 161: m21 = 176: GoSub 2700: If fl1 = 0 Then GoTo 1920 'Check row 6
a(160) = s1 / 8 - a(192) - a(254) + a(256): If a(160) < m1 Or a(160) > m2 Then GoTo 1920:
a(159) = -s1 / 8 + a(192) + a(254) + a(255): If a(159) < m1 Or a(159) > m2 Then GoTo 1920:
a(158) = s1 / 8 - a(192): If a(158) < m1 Or a(158) > m2 Then GoTo 1920:
a(157) = s1 / 8 + a(192) - a(255) - a(256): If a(157) < m1 Or a(157) > m2 Then GoTo 1920:
a(156) = s1 / 8 - a(192) + a(252) - a(254): If a(156) < m1 Or a(156) > m2 Then GoTo 1920:
a(155) = -s1 / 8 + a(192) - a(252) + a(254) + a(255) + a(256): If a(155) < m1 Or a(155) > m2 Then GoTo 1920:
a(154) = s1 / 8 - a(192) - a(252) + a(256): If a(154) < m1 Or a(154) > m2 Then GoTo 1920:
a(153) = s1 / 8 + a(192) + a(252) - a(255) - 2 * a(256): If a(153) < m1 Or a(153) > m2 Then GoTo 1920:
a(152) = s1 / 8 - a(192) + a(248) - a(254): If a(152) < m1 Or a(152) > m2 Then GoTo 1920:
a(151) = -s1 / 8 + a(192) - a(248) + a(254) + a(255) + a(256): If a(151) < m1 Or a(151) > m2 Then GoTo 1920:
a(150) = s1 / 8 - a(192) - a(248) + a(256): If a(150) < m1 Or a(150) > m2 Then GoTo 1920:
a(149) = s1 / 8 + a(192) + a(248) - a(255) - 2 * a(256): If a(149) < m1 Or a(149) > m2 Then GoTo 1920:
a(148) = s1 / 8 - a(192) + a(244) - a(254): If a(148) < m1 Or a(148) > m2 Then GoTo 1920:
a(147) = -s1 / 8 + a(192) - a(244) + a(254) + a(255) + a(256): If a(147) < m1 Or a(147) > m2 Then GoTo 1920:
a(146) = s1 / 8 - a(192) - a(244) + a(256): If a(146) < m1 Or a(146) > m2 Then GoTo 1920:
a(145) = s1 / 8 + a(192) + a(244) - a(255) - 2 * a(256): If a(145) < m1 Or a(145) > m2 Then GoTo 1920:
m11 = 145: m21 = 160: GoSub 2700: If fl1 = 0 Then GoTo 1920 'Check row 7
a(144) = s1 / 8 + a(192) - a(240) + a(254) - 2 * a(256): If a(144) < m1 Or a(144) > m2 Then GoTo 1920:
a(143) = s1 / 8 - a(192) + a(240) - a(254) - a(255) + a(256): If a(143) < m1 Or a(143) > m2 Then GoTo 1920:
a(142) = s1 / 8 + a(192) - a(240) - a(256): If a(142) < m1 Or a(142) > m2 Then GoTo 1920:
a(141) = -s1 / 8 - a(192) + a(240) + a(255) + 2 * a(256): If a(141) < m1 Or a(141) > m2 Then GoTo 1920:
a(140) = s1 / 8 + a(192) - a(240) - a(252) + a(254) - a(256): If a(140) < m1 Or a(140) > m2 Then GoTo 1920:
a(139) = s1 / 8 - a(192) + a(240) + a(252) - a(254) - a(255): If a(139) < m1 Or a(139) > m2 Then GoTo 1920:
a(138) = s1 / 8 + a(192) - a(240) + a(252) - 2 * a(256): If a(138) < m1 Or a(138) > m2 Then GoTo 1920:
a(137) = -s1 / 8 - a(192) + a(240) - a(252) + a(255) + 3 * a(256):
If a(137) < m1 Or a(137) > m2 Then GoTo 1920:
a(136) = s1 / 8 + a(192) - a(240) - a(248) + a(254) - a(256): If a(136) < m1 Or a(136) > m2 Then GoTo 1920:
a(135) = s1 / 8 - a(192) + a(240) + a(248) - a(254) - a(255): If a(135) < m1 Or a(135) > m2 Then GoTo 1920:
a(134) = s1 / 8 + a(192) - a(240) + a(248) - 2 * a(256): If a(134) < m1 Or a(134) > m2 Then GoTo 1920:
a(133) = -s1 / 8 - a(192) + a(240) - a(248) + a(255) + 3 * a(256):
If a(133) < m1 Or a(133) > m2 Then GoTo 1920:
a(132) = s1 / 8 + a(192) - a(240) - a(244) + a(254) - a(256): If a(132) < m1 Or a(132) > m2 Then GoTo 1920:
a(131) = s1 / 8 - a(192) + a(240) + a(244) - a(254) - a(255): If a(131) < m1 Or a(131) > m2 Then GoTo 1920:
a(130) = s1 / 8 + a(192) - a(240) + a(244) - 2 * a(256): If a(130) < m1 Or a(130) > m2 Then GoTo 1920:
a(129) = -s1 / 8 - a(192) + a(240) - a(244) + a(255) + 3 * a(256):
If a(129) < m1 Or a(129) > m2 Then GoTo 1920:
m11 = 129: m21 = 144: GoSub 2700: If fl1 = 0 Then GoTo 1920 '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(252) + a(255) + 2 * 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(252) - a(254) - a(255) - a(256): 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(248) + a(255) + 2 * 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(248) - a(254) - a(255) - a(256): 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(244) + a(255) + 2 * 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(244) - a(254) - a(255) - a(256): 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
a(112) = -a(128) + a(240) + a(256): If a(112) < m1 Or a(112) > m2 Then GoTo 1280:
a(111) = s1 / 4 + a(128) - a(240) - a(255) - 2 * a(256): If a(111) < m1 Or a(111) > m2 Then GoTo 1280:
a(110) = -a(128) + a(240) - a(254) + 2 * a(256): If a(110) < m1 Or a(110) > m2 Then GoTo 1280:
a(109) = a(128) - a(240) + a(254) + a(255) - a(256): If a(109) < m1 Or a(109) > m2 Then GoTo 1280:
a(108) = -a(128) + a(240) - a(252) + 2 * a(256): If a(108) < m1 Or a(108) > m2 Then GoTo 1280:
a(107) = s1 / 4 + a(128) - a(240) + a(252) - a(255) - 3 * a(256):
If a(107) < m1 Or a(107) > m2 Then GoTo 1280:
a(106) = -a(128) + a(240) + a(252) - a(254) + a(256): If a(106) < m1 Or a(106) > m2 Then GoTo 1280:
a(105) = a(128) - a(240) - a(252) + a(254) + a(255): If a(105) < m1 Or a(105) > m2 Then GoTo 1280:
a(104) = -a(128) + a(240) - a(248) + 2 * a(256): If a(104) < m1 Or a(104) > m2 Then GoTo 1280:
a(103) = s1 / 4 + a(128) - a(240) + a(248) - a(255) - 3 * a(256):
If a(103) < m1 Or a(103) > m2 Then GoTo 1280:
a(102) = -a(128) + a(240) + a(248) - a(254) + a(256): If a(102) < m1 Or a(102) > m2 Then GoTo 1280:
a(101) = a(128) - a(240) - a(248) + a(254) + a(255): If a(101) < m1 Or a(101) > m2 Then GoTo 1280:
a(100) = -a(128) + a(240) - a(244) + 2 * a(256): If a(100) < m1 Or a(100) > m2 Then GoTo 1280:
a(99) = s1 / 4 + a(128) - a(240) + a(244) - a(255) - 3 * a(256):
If a(99) < m1 Or a(99) > m2 Then GoTo 1280:
a(98) = -a(128) + a(240) + a(244) - a(254) + a(256): If a(98) < m1 Or a(98) > m2 Then GoTo 1280:
a(97) = a(128) - a(240) - a(244) + a(254) + a(255): If a(97) < m1 Or a(97) > m2 Then GoTo 1280:
m11 = 97: m21 = 112: GoSub 2700: If fl1 = 0 Then GoTo 1280 'Check row 10
a(96) = s1 / 8 - a(128) - a(254) + a(256): If a(96) < m1 Or a(96) > m2 Then GoTo 1280:
a(95) = -s1 / 8 + a(128) + a(254) + a(255): If a(95) < m1 Or a(95) > m2 Then GoTo 1280:
a(94) = s1 / 8 - a(128): If a(94) < m1 Or a(94) > m2 Then GoTo 1280:
a(93) = s1 / 8 + a(128) - a(255) - a(256): If a(93) < m1 Or a(93) > m2 Then GoTo 1280:
a(92) = s1 / 8 - a(128) + a(252) - a(254): If a(92) < m1 Or a(92) > m2 Then GoTo 1280:
a(91) = -s1 / 8 + a(128) - a(252) + a(254) + a(255) + a(256): If a(91) < m1 Or a(91) > m2 Then GoTo 1280:
a(90) = s1 / 8 - a(128) - a(252) + a(256): If a(90) < m1 Or a(90) > m2 Then GoTo 1280:
a(89) = s1 / 8 + a(128) + a(252) - a(255) - 2 * a(256): If a(89) < m1 Or a(89) > m2 Then GoTo 1280:
a(88) = s1 / 8 - a(128) + a(248) - a(254): If a(88) < m1 Or a(88) > m2 Then GoTo 1280:
a(87) = -s1 / 8 + a(128) - a(248) + a(254) + a(255) + a(256): If a(87) < m1 Or a(87) > m2 Then GoTo 1280:
a(86) = s1 / 8 - a(128) - a(248) + a(256): If a(86) < m1 Or a(86) > m2 Then GoTo 1280:
a(85) = s1 / 8 + a(128) + a(248) - a(255) - 2 * a(256): If a(85) < m1 Or a(85) > m2 Then GoTo 1280:
a(84) = s1 / 8 - a(128) + a(244) - a(254): If a(84) < m1 Or a(84) > m2 Then GoTo 1280:
a(83) = -s1 / 8 + a(128) - a(244) + a(254) + a(255) + a(256): If a(83) < m1 Or a(83) > m2 Then GoTo 1280:
a(82) = s1 / 8 - a(128) - a(244) + a(256): If a(82) < m1 Or a(82) > m2 Then GoTo 1280:
a(81) = s1 / 8 + a(128) + a(244) - a(255) - 2 * a(256): If a(81) < m1 Or a(81) > m2 Then GoTo 1280:
m11 = 81: m21 = 96: GoSub 2700: If fl1 = 0 Then GoTo 1280 'Check row 11
a(80) = s1 / 8 + a(128) - a(240) + a(254) - 2 * a(256): If a(80) < m1 Or a(80) > m2 Then GoTo 1280:
a(79) = s1 / 8 - a(128) + a(240) - a(254) - a(255) + a(256): If a(79) < m1 Or a(79) > m2 Then GoTo 1280:
a(78) = s1 / 8 + a(128) - a(240) - a(256): If a(78) < m1 Or a(78) > m2 Then GoTo 1280:
a(77) = -s1 / 8 - a(128) + a(240) + a(255) + 2 * a(256): If a(77) < m1 Or a(77) > m2 Then GoTo 1280:
a(76) = s1 / 8 + a(128) - a(240) - a(252) + a(254) - a(256): If a(76) < m1 Or a(76) > m2 Then GoTo 1280:
a(75) = s1 / 8 - a(128) + a(240) + a(252) - a(254) - a(255): If a(75) < m1 Or a(75) > m2 Then GoTo 1280:
a(74) = s1 / 8 + a(128) - a(240) + a(252) - 2 * a(256): If a(74) < m1 Or a(74) > m2 Then GoTo 1280:
a(73) = -s1 / 8 - a(128) + a(240) - a(252) + a(255) + 3 * a(256):
If a(73) < m1 Or a(73) > m2 Then GoTo 1280:
a(72) = s1 / 8 + a(128) - a(240) - a(248) + a(254) - a(256): If a(72) < m1 Or a(72) > m2 Then GoTo 1280:
a(71) = s1 / 8 - a(128) + a(240) + a(248) - a(254) - a(255): If a(71) < m1 Or a(71) > m2 Then GoTo 1280:
a(70) = s1 / 8 + a(128) - a(240) + a(248) - 2 * a(256): If a(70) < m1 Or a(70) > m2 Then GoTo 1280:
a(69) = -s1 / 8 - a(128) + a(240) - a(248) + a(255) + 3 * a(256):
If a(69) < m1 Or a(69) > m2 Then GoTo 1280:
a(68) = s1 / 8 + a(128) - a(240) - a(244) + a(254) - a(256): If a(68) < m1 Or a(68) > m2 Then GoTo 1280:
a(67) = s1 / 8 - a(128) + a(240) + a(244) - a(254) - a(255): If a(67) < m1 Or a(67) > m2 Then GoTo 1280:
a(66) = s1 / 8 + a(128) - a(240) + a(244) - 2 * a(256): If a(66) < m1 Or a(66) > m2 Then GoTo 1280:
a(65) = -s1 / 8 - a(128) + a(240) - a(244) + a(255) + 3 * a(256):
If a(65) < m1 Or a(65) > m2 Then GoTo 1280:
m11 = 65: m21 = 80: GoSub 2700: If fl1 = 0 Then GoTo 1280 '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(252) + a(255) + 2 * 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(252) - a(254) - a(255) - a(256): 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(248) + a(255) + 2 * 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(248) - a(254) - a(255) - a(256): 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(244) + a(255) + 2 * 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(244) - a(254) - a(255) - a(256): 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
a(48) = -a(64) + a(240) + a(256): If a(48) < m1 Or a(48) > m2 Then GoTo 640:
a(47) = s1 / 4 + a(64) - a(240) - a(255) - 2 * a(256): If a(47) < m1 Or a(47) > m2 Then GoTo 640:
a(46) = -a(64) + a(240) - a(254) + 2 * a(256): If a(46) < m1 Or a(46) > m2 Then GoTo 640:
a(45) = a(64) - a(240) + a(254) + a(255) - a(256): If a(45) < m1 Or a(45) > m2 Then GoTo 640:
a(44) = -a(64) + a(240) - a(252) + 2 * a(256): If a(44) < m1 Or a(44) > m2 Then GoTo 640:
a(43) = s1 / 4 + a(64) - a(240) + a(252) - a(255) - 3 * a(256):
If a(43) < m1 Or a(43) > m2 Then GoTo 640:
a(42) = -a(64) + a(240) + a(252) - a(254) + a(256): If a(42) < m1 Or a(42) > m2 Then GoTo 640:
a(41) = a(64) - a(240) - a(252) + a(254) + a(255): If a(41) < m1 Or a(41) > m2 Then GoTo 640:
a(40) = -a(64) + a(240) - a(248) + 2 * a(256): If a(40) < m1 Or a(40) > m2 Then GoTo 640:
a(39) = s1 / 4 + a(64) - a(240) + a(248) - a(255) - 3 * a(256):
If a(39) < m1 Or a(39) > m2 Then GoTo 640:
a(38) = -a(64) + a(240) + a(248) - a(254) + a(256): If a(38) < m1 Or a(38) > m2 Then GoTo 640:
a(37) = a(64) - a(240) - a(248) + a(254) + a(255): If a(37) < m1 Or a(37) > m2 Then GoTo 640:
a(36) = -a(64) + a(240) - a(244) + 2 * a(256): If a(36) < m1 Or a(36) > m2 Then GoTo 640:
a(35) = s1 / 4 + a(64) - a(240) + a(244) - a(255) - 3 * a(256):
If a(35) < m1 Or a(35) > m2 Then GoTo 640:
a(34) = -a(64) + a(240) + a(244) - a(254) + a(256): If a(34) < m1 Or a(34) > m2 Then GoTo 640:
a(33) = a(64) - a(240) - a(244) + a(254) + a(255): If a(33) < m1 Or a(33) > m2 Then GoTo 640:
m11 = 33: m21 = 48: GoSub 2700: If fl1 = 0 Then GoTo 640 'Check row 14
a(32) = s1 / 8 - a(64) - a(254) + a(256): If a(32) < m1 Or a(32) > m2 Then GoTo 640:
a(31) = -s1 / 8 + a(64) + a(254) + a(255): If a(31) < m1 Or a(31) > m2 Then GoTo 640:
a(30) = s1 / 8 - a(64): If a(30) < m1 Or a(30) > m2 Then GoTo 640:
a(29) = s1 / 8 + a(64) - a(255) - a(256): If a(29) < m1 Or a(29) > m2 Then GoTo 640:
a(28) = s1 / 8 - a(64) + a(252) - a(254): If a(28) < m1 Or a(28) > m2 Then GoTo 640:
a(27) = -s1 / 8 + a(64) - a(252) + a(254) + a(255) + a(256):
If a(27) < m1 Or a(27) > m2 Then GoTo 640:
a(26) = s1 / 8 - a(64) - a(252) + a(256): If a(26) < m1 Or a(26) > m2 Then GoTo 640:
a(25) = s1 / 8 + a(64) + a(252) - a(255) - 2 * a(256): If a(25) < m1 Or a(25) > m2 Then GoTo 640:
a(24) = s1 / 8 - a(64) + a(248) - a(254): If a(24) < m1 Or a(24) > m2 Then GoTo 640:
a(23) = -s1 / 8 + a(64) - a(248) + a(254) + a(255) + a(256): If a(23) < m1 Or a(23) > m2 Then GoTo 640:
a(22) = s1 / 8 - a(64) - a(248) + a(256): If a(22) < m1 Or a(22) > m2 Then GoTo 640:
a(21) = s1 / 8 + a(64) + a(248) - a(255) - 2 * a(256): If a(21) < m1 Or a(21) > m2 Then GoTo 640:
a(20) = s1 / 8 - a(64) + a(244) - a(254): If a(20) < m1 Or a(20) > m2 Then GoTo 640:
a(19) = -s1 / 8 + a(64) - a(244) + a(254) + a(255) + a(256): If a(19) < m1 Or a(19) > m2 Then GoTo 640:
a(18) = s1 / 8 - a(64) - a(244) + a(256): If a(18) < m1 Or a(18) > m2 Then GoTo 640:
a(17) = s1 / 8 + a(64) + a(244) - a(255) - 2 * a(256): If a(17) < m1 Or a(17) > m2 Then GoTo 640:
m11 = 17: m21 = 32: GoSub 2700: If fl1 = 0 Then GoTo 640 'Check row 15
a(16) = s1 / 8 + a(64) - a(240) + a(254) - 2 * a(256): If a(16) < m1 Or a(16) > m2 Then GoTo 640:
a(15) = s1 / 8 - a(64) + a(240) - a(254) - a(255) + a(256): If a(15) < m1 Or a(15) > m2 Then GoTo 640:
a(14) = s1 / 8 + a(64) - a(240) - a(256): If a(14) < m1 Or a(14) > m2 Then GoTo 640:
a(13) = -s1 / 8 - a(64) + a(240) + a(255) + 2 * a(256): If a(13) < m1 Or a(13) > m2 Then GoTo 640:
a(12) = s1 / 8 + a(64) - a(240) - a(252) + a(254) - a(256): If a(12) < m1 Or a(12) > m2 Then GoTo 640:
a(11) = s1 / 8 - a(64) + a(240) + a(252) - a(254) - a(255): If a(11) < m1 Or a(11) > m2 Then GoTo 640:
a(10) = s1 / 8 + a(64) - a(240) + a(252) - 2 * a(256): If a(10) < m1 Or a(10) > m2 Then GoTo 640:
a(9) = -s1 / 8 - a(64) + a(240) - a(252) + a(255) + 3 * a(256):
If a(9) < m1 Or a(9) > m2 Then GoTo 640:
a(8) = s1 / 8 + a(64) - a(240) - a(248) + a(254) - a(256): If a(8) < m1 Or a(8) > m2 Then GoTo 640:
a(7) = s1 / 8 - a(64) + a(240) + a(248) - a(254) - a(255): If a(7) < m1 Or a(7) > m2 Then GoTo 640:
a(6) = s1 / 8 + a(64) - a(240) + a(248) - 2 * a(256): If a(6) < m1 Or a(6) > m2 Then GoTo 640:
a(5) = -s1 / 8 - a(64) + a(240) - a(248) + a(255) + 3 * a(256):
If a(5) < m1 Or a(5) > m2 Then GoTo 640:
a(4) = s1 / 8 + a(64) - a(240) - a(244) + a(254) - a(256): If a(4) < m1 Or a(4) > m2 Then GoTo 640:
a(3) = s1 / 8 - a(64) + a(240) + a(244) - a(254) - a(255): If a(3) < m1 Or a(3) > m2 Then GoTo 640:
a(2) = s1 / 8 + a(64) - a(240) + a(244) - 2 * a(256): If a(2) < m1 Or a(2) > m2 Then GoTo 640:
a(1) = -s1 / 8 - a(64) + a(240) - a(244) + a(255) + 3 * a(256):
If a(1) < m1 Or a(1) > m2 Then GoTo 640:
m11 = 1: m21 = 16: GoSub 2700: If fl1 = 0 Then GoTo 640 'Check row 16
GoSub 2750: If fl1 = 0 Then GoTo 640 '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 640
n9 = n9 + 1 'Total
' GoSub 2600 'Print results (Lines)
GoSub 2650 'Print results (Squares)
' Cells(2, 1).Value = n9 'Counting
640 Next j64
1280 Next j128
1920 Next j192
2400 Next j240
2440 Next j244
2480 Next j248
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 MstPrfct16b")
End
' Construct Most Perfect Franklin Pan Magic (16 x 16)
' (Barink Restrictions)
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
