Vorige Pagina About the Author

' Generates Most Perfect Magic Squares of order 12, magic sum 870
' 1/3 Rows and 1/3 Columns sum to 290

' Tested with Office 365 under Windows 10

Sub MgcSqr12e()

Dim a(144), b(144), c(144)

y = MsgBox("Locked", vbCritical, "Routine MgcSqr12e")
End

n2 = 0: n9 = 0: k1 = 1: k2 = 1
m1 = 1: m2 = 144: s1 = 870: Pr12 = 145

'   Generate data
    
    Sheets("Klad1").Select
    
    t1 = Timer

For j144 = 1 To 1 ''m1 To m2                                               'a(144)
If b(j144) = 0 Then b(j144) = j144: c(144) = j144 Else GoTo 1440
a(144) = j144
a(66) = 145 - a(144): If b(a(66)) = 0 Then b(a(66)) = a(66): c(66) = a(66) Else GoTo 660

For j143 = 138 To 138 ''m1 To m2                                           'a(143)
If b(j143) = 0 Then b(j143) = j143: c(143) = j143 Else GoTo 1430
a(143) = j143
a(65) = 145 - a(143): If b(a(65)) = 0 Then b(a(65)) = a(65): c(65) = a(65) Else GoTo 650

For j142 = m1 To m2                                               'a(142)
If b(j142) = 0 Then b(j142) = j142: c(142) = j142 Else GoTo 1420
a(142) = j142
a(64) = 145 - a(142): If b(a(64)) = 0 Then b(a(64)) = a(64): c(64) = a(64) Else GoTo 640

Cells(2, 1).Value = j142

a(141) = 4 * s1 / 12 - a(142) - a(143) - a(144)
If a(141) < 1 Or a(141) > 144 Then GoTo 1410
If b(a(141)) = 0 Then b(a(141)) = a(141): c(141) = a(141) Else GoTo 1410
a(63) = 145 - a(141): If b(a(63)) = 0 Then b(a(63)) = a(63): c(63) = a(63) Else GoTo 630

For j140 = m1 To m2                                               'a(140)
If b(j140) = 0 Then b(j140) = j140: c(140) = j140 Else GoTo 1400
a(140) = j140
a(62) = 145 - a(140): If b(a(62)) = 0 Then b(a(62)) = a(62): c(62) = a(62) Else GoTo 620

Cells(3, 1).Value = j140

a(139) = -a(140) + a(143) + a(144)
If a(139) < 1 Or a(139) > 144 Then GoTo 1390
If b(a(139)) = 0 Then b(a(139)) = a(139): c(139) = a(139) Else GoTo 1390
a(61) = 145 - a(139): If b(a(61)) = 0 Then b(a(61)) = a(61): c(61) = a(61) Else GoTo 610

For j138 = m1 To m2                                               'a(138)
If b(j138) = 0 Then b(j138) = j138: c(138) = j138 Else GoTo 1380
a(138) = j138
a(72) = 145 - a(138): If b(a(72)) = 0 Then b(a(72)) = a(72): c(72) = a(72) Else GoTo 720

a(137) = 4 * s1 / 12 - a(138) - a(143) - a(144)
If a(137) < 1 Or a(137) > 144 Then GoTo 1370
If b(a(137)) = 0 Then b(a(137)) = a(137): c(137) = a(137) Else GoTo 1370
a(71) = 145 - a(137): If b(a(71)) = 0 Then b(a(71)) = a(71): c(71) = a(71) Else GoTo 710

a(136) = a(138) - a(142) + a(144)
If a(136) < 1 Or a(136) > 144 Then GoTo 1360
If b(a(136)) = 0 Then b(a(136)) = a(136): c(136) = a(136) Else GoTo 1360
a(70) = 145 - a(136): If b(a(70)) = 0 Then b(a(70)) = a(70): c(70) = a(70) Else GoTo 700

a(135) = -a(138) + a(142) + a(143)
If a(135) < 1 Or a(135) > 144 Then GoTo 1350
If b(a(135)) = 0 Then b(a(135)) = a(135): c(135) = a(135) Else GoTo 1350
a(69) = 145 - a(135): If b(a(69)) = 0 Then b(a(69)) = a(69): c(69) = a(69) Else GoTo 690

a(134) = a(138) - a(140) + a(144)
If a(134) < 1 Or a(134) > 144 Then GoTo 1340
If b(a(134)) = 0 Then b(a(134)) = a(134): c(134) = a(134) Else GoTo 1340
a(68) = 145 - a(134): If b(a(68)) = 0 Then b(a(68)) = a(68): c(68) = a(68) Else GoTo 680

a(133) = 4 * s1 / 12 - a(138) + a(140) - a(143) - 2 * a(144)
If a(133) < 1 Or a(133) > 144 Then GoTo 1330
If b(a(133)) = 0 Then b(a(133)) = a(133): c(133) = a(133) Else GoTo 1330
a(67) = 145 - a(133): If b(a(67)) = 0 Then b(a(67)) = a(67): c(67) = a(67) Else GoTo 670

For j132 = m1 To m2                                               'a(132)
If b(j132) = 0 Then b(j132) = j132: c(132) = j132 Else GoTo 1320
a(132) = j132
a(54) = 145 - a(132): If b(a(54)) = 0 Then b(a(54)) = a(54): c(54) = a(54) Else GoTo 540

a(131) = 4 * s1 / 12 - a(132) - a(143) - a(144)
If a(131) < 1 Or a(131) > 144 Then GoTo 1310
If b(a(131)) = 0 Then b(a(131)) = a(131): c(131) = a(131) Else GoTo 1310
a(53) = 145 - a(131): If b(a(53)) = 0 Then b(a(53)) = a(53): c(53) = a(53) Else GoTo 530

a(130) = a(132) - a(142) + a(144)
If a(130) < 1 Or a(130) > 144 Then GoTo 1300
If b(a(130)) = 0 Then b(a(130)) = a(130): c(130) = a(130) Else GoTo 1300
a(52) = 145 - a(130): If b(a(52)) = 0 Then b(a(52)) = a(52): c(52) = a(52) Else GoTo 520

a(129) = -a(132) + a(142) + a(143)
If a(129) < 1 Or a(129) > 144 Then GoTo 1290
If b(a(129)) = 0 Then b(a(129)) = a(129): c(129) = a(129) Else GoTo 1290
a(51) = 145 - a(129): If b(a(51)) = 0 Then b(a(51)) = a(51): c(51) = a(51) Else GoTo 510

a(128) = a(132) - a(140) + a(144)
If a(128) < 1 Or a(128) > 144 Then GoTo 1280
If b(a(128)) = 0 Then b(a(128)) = a(128): c(128) = a(128) Else GoTo 1280
a(50) = 145 - a(128): If b(a(50)) = 0 Then b(a(50)) = a(50): c(50) = a(50) Else GoTo 500

a(127) = 4 * s1 / 12 - a(132) + a(140) - a(143) - 2 * a(144)
If a(127) < 1 Or a(127) > 144 Then GoTo 1270
If b(a(127)) = 0 Then b(a(127)) = a(127): c(127) = a(127) Else GoTo 1270
a(49) = 145 - a(127): If b(a(49)) = 0 Then b(a(49)) = a(49): c(49) = a(49) Else GoTo 490

a(126) = a(132) - a(138) + a(144)
If a(126) < 1 Or a(126) > 144 Then GoTo 1260
If b(a(126)) = 0 Then b(a(126)) = a(126): c(126) = a(126) Else GoTo 1260
a(60) = 145 - a(126): If b(a(60)) = 0 Then b(a(60)) = a(60): c(60) = a(60) Else GoTo 600

a(125) = -a(132) + a(138) + a(143)
If a(125) < 1 Or a(125) > 144 Then GoTo 1250
If b(a(125)) = 0 Then b(a(125)) = a(125): c(125) = a(125) Else GoTo 1250
a(59) = 145 - a(125): If b(a(59)) = 0 Then b(a(59)) = a(59): c(59) = a(59) Else GoTo 590

a(124) = a(132) - a(138) + a(142)
If a(124) < 1 Or a(124) > 144 Then GoTo 1240
If b(a(124)) = 0 Then b(a(124)) = a(124): c(124) = a(124) Else GoTo 1240
a(58) = 145 - a(124): If b(a(58)) = 0 Then b(a(58)) = a(58): c(58) = a(58) Else GoTo 580

a(123) = 4 * s1 / 12 - a(132) + a(138) - a(142) - a(143) - a(144):
If a(123) < 1 Or a(123) > 144 Then GoTo 1230
If b(a(123)) = 0 Then b(a(123)) = a(123): c(123) = a(123) Else GoTo 1230
a(57) = 145 - a(123): If b(a(57)) = 0 Then b(a(57)) = a(57): c(57) = a(57) Else GoTo 570

a(122) = a(132) - a(138) + a(140)
If a(122) < 1 Or a(122) > 144 Then GoTo 1220
If b(a(122)) = 0 Then b(a(122)) = a(122): c(122) = a(122) Else GoTo 1220
a(56) = 145 - a(122): If b(a(56)) = 0 Then b(a(56)) = a(56): c(56) = a(56) Else GoTo 560

a(121) = -a(132) + a(138) - a(140) + a(143) + a(144)
If a(121) < 1 Or a(121) > 144 Then GoTo 1210
If b(a(121)) = 0 Then b(a(121)) = a(121): c(121) = a(121) Else GoTo 1210
a(55) = 145 - a(121): If b(a(55)) = 0 Then b(a(55)) = a(55): c(55) = a(55) Else GoTo 550

For j120 = m1 To m2                                               'a(120)
If b(j120) = 0 Then b(j120) = j120: c(120) = j120 Else GoTo 1200
a(120) = j120
a(42) = 145 - a(120): If b(a(42)) = 0 Then b(a(42)) = a(42): c(42) = a(42) Else GoTo 420

a(119) = -a(120) + a(143) + a(144)
If a(119) < 1 Or a(119) > 144 Then GoTo 1190
If b(a(119)) = 0 Then b(a(119)) = a(119): c(119) = a(119) Else GoTo 1190
a(41) = 145 - a(119): If b(a(41)) = 0 Then b(a(41)) = a(41): c(41) = a(41) Else GoTo 410

a(118) = a(120) + a(142) - a(144)
If a(118) < 1 Or a(118) > 144 Then GoTo 1180
If b(a(118)) = 0 Then b(a(118)) = a(118): c(118) = a(118) Else GoTo 1180
a(40) = 145 - a(118): If b(a(40)) = 0 Then b(a(40)) = a(40): c(40) = a(40) Else GoTo 400

a(117) = 4 * s1 / 12 - a(120) - a(142) - a(143)
If a(117) < 1 Or a(117) > 144 Then GoTo 1170
If b(a(117)) = 0 Then b(a(117)) = a(117): c(117) = a(117) Else GoTo 1170
a(39) = 145 - a(117): If b(a(39)) = 0 Then b(a(39)) = a(39): c(39) = a(39) Else GoTo 390

a(116) = a(120) + a(140) - a(144)
If a(116) < 1 Or a(116) > 144 Then GoTo 1160
If b(a(116)) = 0 Then b(a(116)) = a(116): c(116) = a(116) Else GoTo 1160
a(38) = 145 - a(116): If b(a(38)) = 0 Then b(a(38)) = a(38): c(38) = a(38) Else GoTo 380

a(115) = -a(120) - a(140) + a(143) + 2 * a(144)
If a(115) < 1 Or a(115) > 144 Then GoTo 1150
If b(a(115)) = 0 Then b(a(115)) = a(115): c(115) = a(115) Else GoTo 1150
a(37) = 145 - a(115): If b(a(37)) = 0 Then b(a(37)) = a(37): c(37) = a(37) Else GoTo 370

a(114) = a(120) + a(138) - a(144)
If a(114) < 1 Or a(114) > 144 Then GoTo 1140
If b(a(114)) = 0 Then b(a(114)) = a(114): c(114) = a(114) Else GoTo 1140
a(48) = 145 - a(114): If b(a(48)) = 0 Then b(a(48)) = a(48): c(48) = a(48) Else GoTo 480

a(113) = 4 * s1 / 12 - a(120) - a(138) - a(143)
If a(113) < 1 Or a(113) > 144 Then GoTo 1130
If b(a(113)) = 0 Then b(a(113)) = a(113): c(113) = a(113) Else GoTo 1130
a(47) = 145 - a(113): If b(a(47)) = 0 Then b(a(47)) = a(47): c(47) = a(47) Else GoTo 470

a(112) = a(120) + a(138) - a(142)
If a(112) < 1 Or a(112) > 144 Then GoTo 1120
If b(a(112)) = 0 Then b(a(112)) = a(112): c(112) = a(112) Else GoTo 1120
a(46) = 145 - a(112): If b(a(46)) = 0 Then b(a(46)) = a(46): c(46) = a(46) Else GoTo 460

a(111) = -a(120) - a(138) + a(142) + a(143) + a(144)
If a(111) < 1 Or a(111) > 144 Then GoTo 1110
If b(a(111)) = 0 Then b(a(111)) = a(111): c(111) = a(111) Else GoTo 1110
a(45) = 145 - a(111): If b(a(45)) = 0 Then b(a(45)) = a(45): c(45) = a(45) Else GoTo 450

a(110) = a(120) + a(138) - a(140)
If a(110) < 1 Or a(110) > 144 Then GoTo 1100
If b(a(110)) = 0 Then b(a(110)) = a(110): c(110) = a(110) Else GoTo 1100
a(44) = 145 - a(110): If b(a(44)) = 0 Then b(a(44)) = a(44): c(44) = a(44) Else GoTo 440

a(109) = 4 * s1 / 12 - a(120) - a(138) + a(140) - a(143) - a(144):
If a(109) < 1 Or a(109) > 144 Then GoTo 1090
If b(a(109)) = 0 Then b(a(109)) = a(109): c(109) = a(109) Else GoTo 1090
a(43) = 145 - a(109): If b(a(43)) = 0 Then b(a(43)) = a(43): c(43) = a(43) Else GoTo 430

a(108) = 4 * s1 / 12 - a(120) - a(132) - a(144)
If a(108) < 1 Or a(108) > 144 Then GoTo 1080
If b(a(108)) = 0 Then b(a(108)) = a(108): c(108) = a(108) Else GoTo 1080
a(30) = 145 - a(108): If b(a(30)) = 0 Then b(a(30)) = a(30): c(30) = a(30) Else GoTo 300

a(107) = a(120) + a(132) - a(143)
If a(107) < 1 Or a(107) > 144 Then GoTo 1070
If b(a(107)) = 0 Then b(a(107)) = a(107): c(107) = a(107) Else GoTo 1070
a(29) = 145 - a(107): If b(a(29)) = 0 Then b(a(29)) = a(29): c(29) = a(29) Else GoTo 290

a(106) = 4 * s1 / 12 - a(120) - a(132) - a(142)
If a(106) < 1 Or a(106) > 144 Then GoTo 1060
If b(a(106)) = 0 Then b(a(106)) = a(106): c(106) = a(106) Else GoTo 1060
a(28) = 145 - a(106): If b(a(28)) = 0 Then b(a(28)) = a(28): c(28) = a(28) Else GoTo 280

a(105) = -4 * s1 / 12 + a(120) + a(132) + a(142) + a(143) + a(144)
If a(105) < 1 Or a(105) > 144 Then GoTo 1050
If b(a(105)) = 0 Then b(a(105)) = a(105): c(105) = a(105) Else GoTo 1050
a(27) = 145 - a(105): If b(a(27)) = 0 Then b(a(27)) = a(27): c(27) = a(27) Else GoTo 270

a(104) = 4 * s1 / 12 - a(120) - a(132) - a(140)
If a(104) < 1 Or a(104) > 144 Then GoTo 1040
If b(a(104)) = 0 Then b(a(104)) = a(104): c(104) = a(104) Else GoTo 1040
a(26) = 145 - a(104): If b(a(26)) = 0 Then b(a(26)) = a(26): c(26) = a(26) Else GoTo 260

a(103) = a(120) + a(132) + a(140) - a(143) - a(144)
If a(103) < 1 Or a(103) > 144 Then GoTo 1030
If b(a(103)) = 0 Then b(a(103)) = a(103): c(103) = a(103) Else GoTo 1030
a(25) = 145 - a(103): If b(a(25)) = 0 Then b(a(25)) = a(25): c(25) = a(25) Else GoTo 250

a(102) = 4 * s1 / 12 - a(120) - a(132) - a(138)
If a(102) < 1 Or a(102) > 144 Then GoTo 1020
If b(a(102)) = 0 Then b(a(102)) = a(102): c(102) = a(102) Else GoTo 1020
a(36) = 145 - a(102): If b(a(36)) = 0 Then b(a(36)) = a(36): c(36) = a(36) Else GoTo 360

a(101) = -4 * s1 / 12 + a(120) + a(132) + a(138) + a(143) + a(144)
If a(101) < 1 Or a(101) > 144 Then GoTo 1010
If b(a(101)) = 0 Then b(a(101)) = a(101): c(101) = a(101) Else GoTo 1010
a(35) = 145 - a(101): If b(a(35)) = 0 Then b(a(35)) = a(35): c(35) = a(35) Else GoTo 350

a(100) = 4 * s1 / 12 - a(120) - a(132) - a(138) + a(142) - a(144)
If a(100) < 1 Or a(100) > 144 Then GoTo 1000
If b(a(100)) = 0 Then b(a(100)) = a(100): c(100) = a(100) Else GoTo 1000
a(34) = 145 - a(100): If b(a(34)) = 0 Then b(a(34)) = a(34): c(34) = a(34) Else GoTo 340

a(99) = a(120) + a(132) + a(138) - a(142) - a(143)
If a(99) < 1 Or a(99) > 144 Then GoTo 990
If b(a(99)) = 0 Then b(a(99)) = a(99): c(99) = a(99) Else GoTo 990
a(33) = 145 - a(99): If b(a(33)) = 0 Then b(a(33)) = a(33): c(33) = a(33) Else GoTo 330

a(98) = 4 * s1 / 12 - a(120) - a(132) - a(138) + a(140) - a(144)
If a(98) < 1 Or a(98) > 144 Then GoTo 980
If b(a(98)) = 0 Then b(a(98)) = a(98): c(98) = a(98) Else GoTo 980
a(32) = 145 - a(98): If b(a(32)) = 0 Then b(a(32)) = a(32): c(32) = a(32) Else GoTo 320

a(97) = -4 * s1 / 12 + a(120) + a(132) + a(138) - a(140) + a(143) + 2 * a(144)
If a(97) < 1 Or a(97) > 144 Then GoTo 970
If b(a(97)) = 0 Then b(a(97)) = a(97): c(97) = a(97) Else GoTo 970
a(31) = 145 - a(97): If b(a(31)) = 0 Then b(a(31)) = a(31): c(31) = a(31) Else GoTo 310

For j96 = m1 To m2                                                'a(96)
If b(j96) = 0 Then b(j96) = j96: c(96) = j96 Else GoTo 960
a(96) = j96
a(18) = 145 - a(96): If b(a(18)) = 0 Then b(a(18)) = a(18): c(18) = a(18) Else GoTo 180

a(95) = -a(96) + a(143) + a(144)
If a(95) < 1 Or a(95) > 144 Then GoTo 950
If b(a(95)) = 0 Then b(a(95)) = a(95): c(95) = a(95) Else GoTo 950
a(17) = 145 - a(95): If b(a(17)) = 0 Then b(a(17)) = a(17): c(17) = a(17) Else GoTo 170

a(94) = a(96) + a(142) - a(144)
If a(94) < 1 Or a(94) > 144 Then GoTo 940
If b(a(94)) = 0 Then b(a(94)) = a(94): c(94) = a(94) Else GoTo 940
a(16) = 145 - a(94): If b(a(16)) = 0 Then b(a(16)) = a(16): c(16) = a(16) Else GoTo 160

a(93) = 4 * s1 / 12 - a(96) - a(142) - a(143)
If a(93) < 1 Or a(93) > 144 Then GoTo 930
If b(a(93)) = 0 Then b(a(93)) = a(93): c(93) = a(93) Else GoTo 930
a(15) = 145 - a(93): If b(a(15)) = 0 Then b(a(15)) = a(15): c(15) = a(15) Else GoTo 150

a(92) = a(96) + a(140) - a(144)
If a(92) < 1 Or a(92) > 144 Then GoTo 920
If b(a(92)) = 0 Then b(a(92)) = a(92): c(92) = a(92) Else GoTo 920
a(14) = 145 - a(92): If b(a(14)) = 0 Then b(a(14)) = a(14): c(14) = a(14) Else GoTo 140

a(91) = -a(96) - a(140) + a(143) + 2 * a(144)
If a(91) < 1 Or a(91) > 144 Then GoTo 910
If b(a(91)) = 0 Then b(a(91)) = a(91): c(91) = a(91) Else GoTo 910
a(13) = 145 - a(91): If b(a(13)) = 0 Then b(a(13)) = a(13): c(13) = a(13) Else GoTo 130

a(90) = a(96) + a(138) - a(144)
If a(90) < 1 Or a(90) > 144 Then GoTo 900
If b(a(90)) = 0 Then b(a(90)) = a(90): c(90) = a(90) Else GoTo 900
a(24) = 145 - a(90): If b(a(24)) = 0 Then b(a(24)) = a(24): c(24) = a(24) Else GoTo 240

a(89) = 4 * s1 / 12 - a(96) - a(138) - a(143)
If a(89) < 1 Or a(89) > 144 Then GoTo 890
If b(a(89)) = 0 Then b(a(89)) = a(89): c(89) = a(89) Else GoTo 890
a(23) = 145 - a(89): If b(a(23)) = 0 Then b(a(23)) = a(23): c(23) = a(23) Else GoTo 230

a(88) = a(96) + a(138) - a(142)
If a(88) < 1 Or a(88) > 144 Then GoTo 880
If b(a(88)) = 0 Then b(a(88)) = a(88): c(88) = a(88) Else GoTo 880
a(22) = 145 - a(88): If b(a(22)) = 0 Then b(a(22)) = a(22): c(22) = a(22) Else GoTo 220

a(87) = -a(96) - a(138) + a(142) + a(143) + a(144)
If a(87) < 1 Or a(87) > 144 Then GoTo 870
If b(a(87)) = 0 Then b(a(87)) = a(87): c(87) = a(87) Else GoTo 870
a(21) = 145 - a(87): If b(a(21)) = 0 Then b(a(21)) = a(21): c(21) = a(21) Else GoTo 210

a(86) = a(96) + a(138) - a(140)
If a(86) < 1 Or a(86) > 144 Then GoTo 860
If b(a(86)) = 0 Then b(a(86)) = a(86): c(86) = a(86) Else GoTo 860
a(20) = 145 - a(86): If b(a(20)) = 0 Then b(a(20)) = a(20): c(20) = a(20) Else GoTo 200

a(85) = 4 * s1 / 12 - a(96) - a(138) + a(140) - a(143) - a(144)
If a(85) < 1 Or a(85) > 144 Then GoTo 850
If b(a(85)) = 0 Then b(a(85)) = a(85): c(85) = a(85) Else GoTo 850
a(19) = 145 - a(85): If b(a(19)) = 0 Then b(a(19)) = a(19): c(19) = a(19) Else GoTo 190

a(84) = -a(96) + a(132) + a(144)
If a(84) < 1 Or a(84) > 144 Then GoTo 840
If b(a(84)) = 0 Then b(a(84)) = a(84): c(84) = a(84) Else GoTo 840
a(6) = 145 - a(84): If b(a(6)) = 0 Then b(a(6)) = a(6): c(6) = a(6) Else GoTo 60

a(83) = 4 * s1 / 12 + a(96) - a(132) - a(143) - 2 * a(144)
If a(83) < 1 Or a(83) > 144 Then GoTo 830
If b(a(83)) = 0 Then b(a(83)) = a(83): c(83) = a(83) Else GoTo 830
a(5) = 145 - a(83): If b(a(5)) = 0 Then b(a(5)) = a(5): c(5) = a(5) Else GoTo 50

a(82) = -a(96) + a(132) - a(142) + 2 * a(144)
If a(82) < 1 Or a(82) > 144 Then GoTo 820
If b(a(82)) = 0 Then b(a(82)) = a(82): c(82) = a(82) Else GoTo 820
a(4) = 145 - a(82): If b(a(4)) = 0 Then b(a(4)) = a(4): c(4) = a(4) Else GoTo 40

a(81) = a(96) - a(132) + a(142) + a(143) - a(144)
If a(81) < 1 Or a(81) > 144 Then GoTo 810
If b(a(81)) = 0 Then b(a(81)) = a(81): c(81) = a(81) Else GoTo 810
a(3) = 145 - a(81): If b(a(3)) = 0 Then b(a(3)) = a(3): c(3) = a(3) Else GoTo 30

a(80) = -a(96) + a(132) - a(140) + 2 * a(144)
If a(80) < 1 Or a(80) > 144 Then GoTo 800
If b(a(80)) = 0 Then b(a(80)) = a(80): c(80) = a(80) Else GoTo 800
a(2) = 145 - a(80): If b(a(2)) = 0 Then b(a(2)) = a(2): c(2) = a(2) Else GoTo 20

a(79) = 4 * s1 / 12 + a(96) - a(132) + a(140) - a(143) - 3 * a(144)
If a(79) < 1 Or a(79) > 144 Then GoTo 790
If b(a(79)) = 0 Then b(a(79)) = a(79): c(79) = a(79) Else GoTo 790
a(1) = 145 - a(79): If b(a(1)) = 0 Then b(a(1)) = a(1): c(1) = a(1) Else GoTo 10

a(78) = -a(96) + a(132) - a(138) + 2 * a(144)
If a(78) < 1 Or a(78) > 144 Then GoTo 780
If b(a(78)) = 0 Then b(a(78)) = a(78): c(78) = a(78) Else GoTo 780
a(12) = 145 - a(78): If b(a(12)) = 0 Then b(a(12)) = a(12): c(12) = a(12) Else GoTo 120

a(77) = a(96) - a(132) + a(138) + a(143) - a(144)
If a(77) < 1 Or a(77) > 144 Then GoTo 770
If b(a(77)) = 0 Then b(a(77)) = a(77): c(77) = a(77) Else GoTo 770
a(11) = 145 - a(77): If b(a(11)) = 0 Then b(a(11)) = a(11): c(11) = a(11) Else GoTo 110

a(76) = -a(96) + a(132) - a(138) + a(142) + a(144)
If a(76) < 1 Or a(76) > 144 Then GoTo 760
If b(a(76)) = 0 Then b(a(76)) = a(76): c(76) = a(76) Else GoTo 760
a(10) = 145 - a(76): If b(a(10)) = 0 Then b(a(10)) = a(10): c(10) = a(10) Else GoTo 100

a(75) = 4 * s1 / 12 + a(96) - a(132) + a(138) - a(142) - a(143) - 2 * a(144)
If a(75) < 1 Or a(75) > 144 Then GoTo 750
If b(a(75)) = 0 Then b(a(75)) = a(75): c(75) = a(75) Else GoTo 750
a(9) = 145 - a(75): If b(a(9)) = 0 Then b(a(9)) = a(9): c(9) = a(9) Else GoTo 90

a(74) = -a(96) + a(132) - a(138) + a(140) + a(144)
If a(74) < 1 Or a(74) > 144 Then GoTo 740
If b(a(74)) = 0 Then b(a(74)) = a(74): c(74) = a(74) Else GoTo 740
a(8) = 145 - a(74): If b(a(8)) = 0 Then b(a(8)) = a(8): c(8) = a(8) Else GoTo 80

a(73) = a(96) - a(132) + a(138) - a(140) + a(143)
If a(73) < 1 Or a(73) > 144 Then GoTo 730
If b(a(73)) = 0 Then b(a(73)) = a(73): c(73) = a(73) Else GoTo 730
a(7) = 145 - a(73): If b(a(7)) = 0 Then b(a(7)) = a(7): c(7) = a(7) Else GoTo 70

                           GoSub 1800: If fl1 = 0 Then GoTo 5

                           n9 = n9 + 1
'                          GoSub 2650               'Print results (squares)
'                          GoSub 2645               'Print results (selected numbers)
                           Cells(1, 1).Value = n9   'Counting
5

    b(c(7)) = 0: c(7) = 0
70  b(c(73)) = 0: c(73) = 0
730 b(c(8)) = 0: c(8) = 0
80  b(c(74)) = 0: c(74) = 0
740 b(c(9)) = 0: c(9) = 0
90  b(c(75)) = 0: c(75) = 0
750 b(c(10)) = 0: c(10) = 0
100 b(c(76)) = 0: c(76) = 0
760 b(c(11)) = 0: c(11) = 0
110 b(c(77)) = 0: c(77) = 0
770 b(c(12)) = 0: c(12) = 0
120 b(c(78)) = 0: c(78) = 0
780 b(c(1)) = 0: c(1) = 0
10  b(c(79)) = 0: c(79) = 0
790 b(c(2)) = 0: c(2) = 0
20  b(c(80)) = 0: c(80) = 0
800 b(c(3)) = 0: c(3) = 0
30  b(c(81)) = 0: c(81) = 0
810 b(c(4)) = 0: c(4) = 0
40  b(c(82)) = 0: c(82) = 0
820 b(c(5)) = 0: c(5) = 0
50  b(c(83)) = 0: c(83) = 0
830 b(c(6)) = 0: c(6) = 0
60  b(c(84)) = 0: c(84) = 0
840 b(c(19)) = 0: c(19) = 0
190 b(c(85)) = 0: c(85) = 0
850 b(c(20)) = 0: c(20) = 0
200 b(c(86)) = 0: c(86) = 0
860 b(c(21)) = 0: c(21) = 0
210 b(c(87)) = 0: c(87) = 0
870 b(c(22)) = 0: c(22) = 0
220 b(c(88)) = 0: c(88) = 0
880 b(c(23)) = 0: c(23) = 0
230 b(c(89)) = 0: c(89) = 0
890 b(c(24)) = 0: c(24) = 0
240 b(c(90)) = 0: c(90) = 0
900 b(c(13)) = 0: c(13) = 0
130 b(c(91)) = 0: c(91) = 0
910 b(c(14)) = 0: c(14) = 0
140 b(c(92)) = 0: c(92) = 0
920 b(c(15)) = 0: c(15) = 0
150 b(c(93)) = 0: c(93) = 0
930 b(c(16)) = 0: c(16) = 0
160 b(c(94)) = 0: c(94) = 0
940 b(c(17)) = 0: c(17) = 0
170 b(c(95)) = 0: c(95) = 0
950 b(c(18)) = 0: c(18) = 0
180 b(c(96)) = 0: c(96) = 0
960 Next j96

     b(c(31)) = 0: c(31) = 0
310  b(c(97)) = 0: c(97) = 0
970  b(c(32)) = 0: c(32) = 0
320  b(c(98)) = 0: c(98) = 0
980  b(c(33)) = 0: c(33) = 0
330  b(c(99)) = 0: c(99) = 0
990  b(c(34)) = 0: c(34) = 0
340  b(c(100)) = 0: c(100) = 0
1000 b(c(35)) = 0: c(35) = 0
350  b(c(101)) = 0: c(101) = 0
1010 b(c(36)) = 0: c(36) = 0
360  b(c(102)) = 0: c(102) = 0
1020 b(c(25)) = 0: c(25) = 0
250  b(c(103)) = 0: c(103) = 0
1030 b(c(26)) = 0: c(26) = 0
260  b(c(104)) = 0: c(104) = 0
1040 b(c(27)) = 0: c(27) = 0
270  b(c(105)) = 0: c(105) = 0
1050 b(c(28)) = 0: c(28) = 0
280  b(c(106)) = 0: c(106) = 0
1060 b(c(29)) = 0: c(29) = 0
290  b(c(107)) = 0: c(107) = 0
1070 b(c(30)) = 0: c(30) = 0
300  b(c(108)) = 0: c(108) = 0
1080 b(c(43)) = 0: c(43) = 0
430  b(c(109)) = 0: c(109) = 0
1090 b(c(44)) = 0: c(44) = 0
440  b(c(110)) = 0: c(110) = 0
1100 b(c(45)) = 0: c(45) = 0
450  b(c(111)) = 0: c(111) = 0
1110 b(c(46)) = 0: c(46) = 0
460  b(c(112)) = 0: c(112) = 0
1120 b(c(47)) = 0: c(47) = 0
470  b(c(113)) = 0: c(113) = 0
1130 b(c(48)) = 0: c(48) = 0
480  b(c(114)) = 0: c(114) = 0
1140 b(c(37)) = 0: c(37) = 0
370  b(c(115)) = 0: c(115) = 0
1150 b(c(38)) = 0: c(38) = 0
380  b(c(116)) = 0: c(116) = 0
1160 b(c(39)) = 0: c(39) = 0
390  b(c(117)) = 0: c(117) = 0
1170 b(c(40)) = 0: c(40) = 0
400  b(c(118)) = 0: c(118) = 0
1180 b(c(41)) = 0: c(41) = 0
410  b(c(119)) = 0: c(119) = 0
1190 b(c(42)) = 0: c(42) = 0
420  b(c(120)) = 0: c(120) = 0
1200 Next j120

     b(c(55)) = 0: c(55) = 0
550  b(c(121)) = 0: c(121) = 0
1210 b(c(56)) = 0: c(56) = 0
560  b(c(122)) = 0: c(122) = 0
1220 b(c(57)) = 0: c(57) = 0
570  b(c(123)) = 0: c(123) = 0
1230 b(c(58)) = 0: c(58) = 0
580  b(c(124)) = 0: c(124) = 0
1240 b(c(59)) = 0: c(59) = 0
590  b(c(125)) = 0: c(125) = 0
1250 b(c(60)) = 0: c(60) = 0
600  b(c(126)) = 0: c(126) = 0
1260 b(c(49)) = 0: c(49) = 0
490  b(c(127)) = 0: c(127) = 0
1270 b(c(50)) = 0: c(50) = 0
500  b(c(128)) = 0: c(128) = 0
1280 b(c(51)) = 0: c(51) = 0
510  b(c(129)) = 0: c(129) = 0
1290 b(c(52)) = 0: c(52) = 0
520  b(c(130)) = 0: c(130) = 0
1300 b(c(53)) = 0: c(53) = 0
530  b(c(131)) = 0: c(131) = 0
1310 b(c(54)) = 0: c(54) = 0
540  b(c(132)) = 0: c(132) = 0
1320 Next j132

     b(c(67)) = 0: c(67) = 0
670  b(c(133)) = 0: c(133) = 0
1330 b(c(68)) = 0: c(68) = 0
680  b(c(134)) = 0: c(134) = 0
1340 b(c(69)) = 0: c(69) = 0
690  b(c(135)) = 0: c(135) = 0
1350 b(c(70)) = 0: c(70) = 0
700  b(c(136)) = 0: c(136) = 0
1360 b(c(71)) = 0: c(71) = 0
710  b(c(137)) = 0: c(137) = 0
1370 b(c(72)) = 0: c(72) = 0
720  b(c(138)) = 0: c(138) = 0
1380 Next j138

     b(c(61)) = 0: c(61) = 0
610  b(c(139)) = 0: c(139) = 0
1390 b(c(62)) = 0: c(62) = 0
620  b(c(140)) = 0: c(140) = 0
1400 Next j140

     b(c(63)) = 0: c(63) = 0
630  b(c(141)) = 0: c(141) = 0
1410 b(c(64)) = 0: c(64) = 0
640  b(c(142)) = 0: c(142) = 0
1420 Next j142

     b(c(65)) = 0: c(65) = 0
650  b(c(143)) = 0: c(143) = 0
1430 Next j143

     b(c(66)) = 0: c(66) = 0
660  b(c(144)) = 0: c(144) = 0
1440 Next j144

    t2 = Timer
    
    t10 = Str(t2 - t1) + " sec., " + Str(n9) + " Solutions for sum" + Str(s1)
    y = MsgBox(t10, 0, "Routine MgcSqr12e")

End

'   Exclude solutions with identical numbers Latin Lines Order 12

1800 fl1 = 1
     For j1 = 1 To 144
        a2 = a(j1):
        For j2 = (1 + j1) To 144
            If a2 = a(j2) Then fl1 = 0: Return
        Next j2
     Next j1
     Return

'    Print results (selected numbers)

2645 Cells(n9, 138).Select
     For i1 = 1 To 144
         Cells(n9, i1).Value = a(i1)
     Next i1
    
     Return

'    Print results (squares)

2650 n2 = n2 + 1
     If n2 = 3 Then
         n2 = 1: k1 = k1 + 13: k2 = 1
     Else
         If n9 > 1 Then k2 = k2 + 13
     End If
     Cells(k1 + 1, k2 + 1).Select
    
     i3 = 0
     For i1 = 1 To 12
         For i2 = 1 To 12
             i3 = i3 + 1
             Cells(k1 + i1, k2 + i2).Value = a(i3)
         Next i2
     Next i1
    
     Return

End Sub

Vorige Pagina About the Author