' Generates (Pan) Magic Squares of order 8 based on Sudoku Comparable Squares
' Tested with Office 2007 under Windows 7
Sub CnstrSqrs8a()
Dim b(3, 64), a(64), c(64), s9(18)
Sheets("Klad1").Select
y = MsgBox("Locked", vbCritical, "Routine CnstrSqrs8a")
End
n2 = 0: n9 = 0: k1 = 1: k2 = 1
n20 = 0: n10 = 0: k10 = 1: k20 = 1
s1 = 260
t1 = Timer
For j1 = 1 To 1008
For j2 = 1 To 864
j10 = j1: j20 = 1: sht1 = "B1": GoSub 100 'Read Sudoku Comparable Square 1
j10 = j2: j20 = 2: sht1 = "B2": GoSub 100 'Read Sudoku Comparable Square 2
For j4 = 1 To 64
a(j4) = b(1, j4) + 8 * b(2, j4) + 1
Next j4
GoSub 300: If fl1 = 0 Then GoTo 20 'Check identical numbers
n9 = n9 + 1: GoSub 740 'Print results (selected numbers)
' n9 = n9 + 1: GoSub 750 'Print results (squares)
GoSub 400 'Construct Square c() (Squared Elements)
GoSub 500: If fl1 = 0 Then GoTo 20 'Check Magic Sum c()
Cells(n9, 67).Value = "Bimagic"
n10 = n10 + 1: GoSub 745 'Print results (squared numbers)
' n10 = n10 + 1: GoSub 755 'Print results (squares)
20 Next j2
Next j1
t2 = Timer
t10 = Str(t2 - t1) + " sec., " + Str(n9) + " Solutions for sum" + Str(s1)
y = MsgBox(t10, 0, "Routine CnstrSqrs8a")
End
' Read Sudoku Comparable Squares (Line Format)
100 For i1 = 1 To 64
b(j20, i1) = Sheets(sht1).Cells(j10, i1).Value
Next i1
Return
' Check identical numbers
300 fl1 = 1
For i1 = 1 To 64
a2 = a(i1)
For i2 = (1 + i1) To 64
If a2 = a(i2) Then fl1 = 0: Return
Next i2
Next i1
Return
' Construct Square c() (Squared Elements)
400 For i1 = 1 To 64
c(i1) = a(i1) ^ 2
Next i1
Return
' Check Magic Sum c()
500 fl1 = 1
s9(1) = c(1) + c(2) + c(3) + c(4) + c(5) + c(6) + c(7) + c(8)
s9(2) = c(9) + c(10) + c(11) + c(12) + c(13) + c(14) + c(15) + c(16)
s9(3) = c(17) + c(18) + c(19) + c(20) + c(21) + c(22) + c(23) + c(24)
s9(4) = c(25) + c(26) + c(27) + c(28) + c(29) + c(30) + c(31) + c(32)
s9(5) = c(33) + c(34) + c(35) + c(36) + c(37) + c(38) + c(39) + c(40)
s9(6) = c(41) + c(42) + c(43) + c(44) + c(45) + c(46) + c(47) + c(48)
s9(7) = c(49) + c(50) + c(51) + c(52) + c(53) + c(54) + c(55) + c(56)
s9(8) = c(57) + c(58) + c(59) + c(60) + c(61) + c(62) + c(63) + c(64)
s9(9) = c(1) + c(9) + c(17) + c(25) + c(33) + c(41) + c(49) + c(57)
s9(10) = c(2) + c(10) + c(18) + c(26) + c(34) + c(42) + c(50) + c(58)
s9(11) = c(3) + c(11) + c(19) + c(27) + c(35) + c(43) + c(51) + c(59)
s9(12) = c(4) + c(12) + c(20) + c(28) + c(36) + c(44) + c(52) + c(60)
s9(13) = c(5) + c(13) + c(21) + c(29) + c(37) + c(45) + c(53) + c(61)
s9(14) = c(6) + c(14) + c(22) + c(30) + c(38) + c(46) + c(54) + c(62)
s9(15) = c(7) + c(15) + c(23) + c(31) + c(39) + c(47) + c(55) + c(63)
s9(16) = c(8) + c(16) + c(24) + c(32) + c(40) + c(48) + c(56) + c(64)
s9(17) = c(1) + c(10) + c(19) + c(28) + c(37) + c(46) + c(55) + c(64)
s9(18) = c(8) + c(15) + c(22) + c(29) + c(36) + c(43) + c(50) + c(57)
For i1 = 1 To 18
If s9(i1) <> 11180 Then fl1 = 0: Return
Next i1
Return
' Print results (selected numbers)
740 Cells(n9, 64).Select
For i1 = 1 To 64
Cells(n9, i1).Value = a(i1)
Next i1
Cells(n9, 65).Value = j1
Cells(n9, 66).Value = j2
Return
' Print results (selected squared numbers)
745
For i1 = 1 To 64
Sheets("Klad2").Cells(n10, i1).Value = c(i1) 'Squared Elements
Next i1
Sheets("Klad2").Cells(n10, 65).Value = j1
Sheets("Klad2").Cells(n10, 66).Value = j2
Sheets("Klad2").Cells(n10, 67).Value = n9
Return
' Print results (squares)
750 n2 = n2 + 1
If n2 = 5 Then
n2 = 1: k1 = k1 + 9: k2 = 1
Else
If n9 > 1 Then k2 = k2 + 9
End If
Cells(k1, k2 + 1).Select
Cells(k1, k2 + 1).Font.Color = -4165632
Cells(k1, k2 + 1).Value = n9
i3 = 0
For i1 = 1 To 8
For i2 = 1 To 8
i3 = i3 + 1
Cells(k1 + i1, k2 + i2).Value = a(i3)
Next i2
Next i1
Return
' Print results (squared numbers, squares)
755 n20 = n20 + 1
If n20 = 5 Then
n20 = 1: k10 = k10 + 9: k20 = 1
Else
If n10 > 1 Then k20 = k20 + 9
End If
Sheets("Klad2").Cells(k10, k20 + 1).Font.Color = -4165632
Sheets("Klad2").Cells(k10, k20 + 1).Value = n9
i3 = 0
For i1 = 1 To 8
For i2 = 1 To 8
i3 = i3 + 1
Sheets("Klad2").Cells(k10 + i1, k20 + i2).Value = c(i3)
Next i2
Next i1
Return
End Sub
