Vorige Pagina About the Author

' Generates Quaternary Pantriagonal Magic Cubes (Complete, Hor. Magic Planes) of order 4 for integers 0 thru 3
' Sudoku Comparable (Optional)

' Tested with Office 2007 under Windows 7

Sub SudCube47()

Dim a(64), b(4)

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

n2 = 0: n9 = 0: k1 = 1: k2 = 1
m1 = 0: m2 = 3: s1 = 6

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

For j64 = m1 To m2                                                   'a(64)
    a(64) = j64
For j63 = m1 To m2                                                   'a(63)
    a(63) = j63
For j62 = m1 To m2                                                   'a(62)
    a(62) = j62

    a(61) = s1 - a(62) - a(63) - a(64): If a(61) < m1 Or a(61) > m2 Then GoTo 620
    
For j60 = m1 To m2                                                   'a(60)
    a(60) = j60
For j59 = m1 To m2                                                   'a(59)
    a(59) = j59
For j58 = m1 To m2                                                   'a(58)
    a(58) = j58

    a(57) = s1 - a(58) - a(59) - a(60): If a(57) < m1 Or a(57) > m2 Then GoTo 580

For j56 = m1 To m2                                                   'a(56)
    a(56) = j56

    a(55) =      a(56) - a(58) + a(60) - a(61) + a(64): 
    If a(55) < m1 Or a(55) > m2 Then GoTo 560
    a(54) = s1 - a(55) - a(58) - a(59): If a(54) < m1 Or a(54) > m2 Then GoTo 560
    a(53) =    - a(56) + a(58) + a(59): If a(53) < m1 Or a(53) > m2 Then GoTo 560
    a(52) = s1 - a(56) - a(60) - a(64): If a(52) < m1 Or a(52) > m2 Then GoTo 560
    a(51) =      a(54) + a(58) - a(63): If a(51) < m1 Or a(51) > m2 Then GoTo 560
    a(50) = s1 - a(51) - a(62) - a(63): If a(50) < m1 Or a(50) > m2 Then GoTo 560
    a(49) =      a(56) + a(60) - a(61): If a(49) < m1 Or a(49) > m2 Then GoTo 560

For j48 = m1 To m2                                                   'a(48)
    a(48) = j48
    
    a(43) = s1 - a(48) - a(59) - a(64): If a(43) < m1 Or a(43) > m2 Then GoTo 480
    a(38) =      a(48) - a(54) + a(64): If a(38) < m1 Or a(38) > m2 Then GoTo 480
    
For j47 = m1 To m2                                                   'a(47)
    a(47) = j47
    
    a(37) = a(47) + a(56) - a(58) - a(59) + a(63): If a(37) < m1 Or a(37) > m2 Then GoTo 470
    
For j46 = m1 To m2                                                   'a(46)
    a(46) = j46
    
    a(45) = s1 - a(46) - a(47) - a(48): If a(45) < m1 Or a(45) > m2 Then GoTo 460
    a(44) =      a(45) - a(46) - a(48) - a(60) + 2 * a(61) + a(63): 
    If a(44) < m1 Or a(44) > m2 Then GoTo 460
    a(42) = s1 - a(45) - a(58) - a(61): If a(42) < m1 Or a(42) > m2 Then GoTo 460
    a(41) =      a(42) - a(43) - a(47) + 2 * a(58) + a(60) - a(63): 
    If a(41) < m1 Or a(41) > m2 Then GoTo 460
    a(40) =      a(46) - a(56) + a(62): If a(40) < m1 Or a(40) > m2 Then GoTo 460
    a(39) =      a(45) - a(56) + a(58) - a(60) + 2 * a(61) - a(64): 
    If a(39) < m1 Or a(39) > m2 Then GoTo 460
    a(36) = s1 - a(45) - a(54) - a(59): If a(36) < m1 Or a(36) > m2 Then GoTo 460
    a(35) =    - a(39) - a(47) + a(48) + a(59) + a(64):             
    If a(35) < m1 Or a(35) > m2 Then GoTo 460
    a(34) = s1 - a(35) - a(46) - a(47): If a(34) < m1 Or a(34) > m2 Then GoTo 460
    a(33) =    - a(36) + a(46) + a(47): If a(33) < m1 Or a(33) > m2 Then GoTo 460
        
    a(1)  = s1 / 2 - a(43): a(2)  = s1 / 2 - a(44): a(3)  = s1 / 2 - a(41): a(4)  = s1 / 2 - a(42):
    a(5)  = s1 / 2 - a(47): a(6)  = s1 / 2 - a(48): a(7)  = s1 / 2 - a(45): a(8)  = s1 / 2 - a(46):
    a(9)  = s1 / 2 - a(35): a(10) = s1 / 2 - a(36): a(11) = s1 / 2 - a(33): a(12) = s1 / 2 - a(34):
    a(13) = s1 / 2 - a(39): a(14) = s1 / 2 - a(40): a(15) = s1 / 2 - a(37): a(16) = s1 / 2 - a(38):
    a(17) = s1 / 2 - a(59): a(18) = s1 / 2 - a(60): a(19) = s1 / 2 - a(57): a(20) = s1 / 2 - a(58):
    a(21) = s1 / 2 - a(63): a(22) = s1 / 2 - a(64): a(23) = s1 / 2 - a(61): a(24) = s1 / 2 - a(62):
    a(25) = s1 / 2 - a(51): a(26) = s1 / 2 - a(52): a(27) = s1 / 2 - a(49): a(28) = s1 / 2 - a(50):
    a(29) = s1 / 2 - a(55): a(30) = s1 / 2 - a(56): a(31) = s1 / 2 - a(53): a(32) = s1 / 2 - a(54):
    
'       Exclude solutions with identical numbers in rows, colums or pillars
        
        GoSub 800: If fl1 = 0 Then GoTo 460     'rows, colums or pillars
'       GoSub 900: If fl1 = 0 Then GoTo 460     'pantriagonals
                           
        n9 = n9 + 1: GoSub 740  'Print results (selected numbers)
'       n9 = n9 + 1: GoSub 750  'Print results (planes 1, 2, 3, 4)
    
460 Next j46
470 Next j47
480 Next j48
    
560 Next j56
    
580 Next j58
590 Next j59
600 Next j60

620 Next j62
630 Next j63
640 Next j64

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

End

'   Print results (selected numbers)

740 For i1 = 1 To 64
        Cells(n9, i1).Value = a(i1)
    Next i1
    Cells(n9, 64).Select
    Return

'   Print results (planes 11, 12, 13 and 14)

750 n2 = n2 + 1
    If n2 = 9 Then
        n2 = 1: k1 = k1 + 20: k2 = 1
    Else
        If n9 > 1 Then k2 = k2 + 5
    End If
       
    Cells(k1, k2 + 1).Select
    Cells(k1, k2 + 1).Font.Color = -4165632
    Cells(k1, k2 + 1).Value = CStr(n9)
       
    For i0 = 1 To 4
        i3 = (4 - i0) * 16
        For i1 = 1 To 4
            For i2 = 1 To 4
                i3 = i3 + 1
                Cells(k1 + i1 + (i0 - 1) * 5, k2 + i2).Value = a(i3)
            Next i2
        Next i1
    Next i0
    
    Return

'   Exclude solutions with identical numbers in rows, colums or pillars

800 fl1 = 1

'   Rows
    
    i1 = -3
    For i0 = 1 To 16
        i1 = i1 + 4
        b(1) = a(i1): b(2) = a(i1 + 1): b(3) = a(i1 + 2): b(4) = a(i1 + 3)
        GoSub 860
        If fl1 = 0 Then Return
    Next i0
   
'   Columns
    
    i1 = 0: i2 = 0
    For i0 = 1 To 16
        i1 = i1 + 1
        b(1) = a(i1): b(2) = a(i1 + 4): b(3) = a(i1 + 8): b(4) = a(i1 + 12)
        i2 = i2 + 1: If i2 = 4 Then i2 = 0: i1 = i1 + 12
        GoSub 860
        If fl1 = 0 Then Return
    Next i0
    
'   Pillars
    
    i1 = 0: i2 = 0
    For i0 = 1 To 16
        i1 = i0
        b(1) = a(i1): b(2) = a(i1 + 16): b(3) = a(i1 + 32): b(4) = a(i1 + 48)
        GoSub 860
        If fl1 = 0 Then Return
    Next i0
    
    Return

860 fl1 = 1
    For j1 = 1 To 4
       b2 = b(j1)
       For j2 = (1 + j1) To 4
           If b2 = b(j2) Then fl1 = 0: Return
       Next j2
    Next j1
    Return
   
'    Check Pan Triagonals

900  fl1 = 1
                
     b(1) = a(49): b(2) = a(38): b(3) = a(27): b(4) = a(16): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(53): b(2) = a(42): b(3) = a(31): b(4) = a(4): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(57): b(2) = a(46): b(3) = a(19): b(4) = a(8): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(61): b(2) = a(34): b(3) = a(23): b(4) = a(12): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(50): b(2) = a(39): b(3) = a(28): b(4) = a(13): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(54): b(2) = a(43): b(3) = a(32): b(4) = a(1): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(58): b(2) = a(47): b(3) = a(20): b(4) = a(5): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(62): b(2) = a(35): b(3) = a(24): b(4) = a(9): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(51): b(2) = a(40): b(3) = a(25): b(4) = a(14): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(55): b(2) = a(44): b(3) = a(29): b(4) = a(2): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(59): b(2) = a(48): b(3) = a(17): b(4) = a(6): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(63): b(2) = a(36): b(3) = a(21): b(4) = a(10): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(52): b(2) = a(37): b(3) = a(26): b(4) = a(15): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(56): b(2) = a(41): b(3) = a(30): b(4) = a(3): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(60): b(2) = a(45): b(3) = a(18): b(4) = a(7): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(64): b(2) = a(33): b(3) = a(22): b(4) = a(11): GoSub 860: If fl1 = 0 Then Return
                    
     b(1) = a(52): b(2) = a(39): b(3) = a(26): b(4) = a(13): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(56): b(2) = a(43): b(3) = a(30): b(4) = a(1): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(60): b(2) = a(47): b(3) = a(18): b(4) = a(5): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(64): b(2) = a(35): b(3) = a(22): b(4) = a(9): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(49): b(2) = a(40): b(3) = a(27): b(4) = a(14): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(53): b(2) = a(44): b(3) = a(31): b(4) = a(2): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(57): b(2) = a(48): b(3) = a(19): b(4) = a(6): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(61): b(2) = a(36): b(3) = a(23): b(4) = a(10): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(50): b(2) = a(37): b(3) = a(28): b(4) = a(15): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(54): b(2) = a(41): b(3) = a(32): b(4) = a(3): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(58): b(2) = a(45): b(3) = a(20): b(4) = a(7): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(62): b(2) = a(33): b(3) = a(24): b(4) = a(11): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(51): b(2) = a(38): b(3) = a(25): b(4) = a(16): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(55): b(2) = a(42): b(3) = a(29): b(4) = a(4): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(59): b(2) = a(46): b(3) = a(17): b(4) = a(8): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(63): b(2) = a(34): b(3) = a(21): b(4) = a(12): GoSub 860: If fl1 = 0 Then Return
                    
     b(1) = a(61): b(2) = a(42): b(3) = a(23): b(4) = a(4): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(49): b(2) = a(46): b(3) = a(27): b(4) = a(8): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(53): b(2) = a(34): b(3) = a(31): b(4) = a(12): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(57): b(2) = a(38): b(3) = a(19): b(4) = a(16): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(62): b(2) = a(43): b(3) = a(24): b(4) = a(1): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(50): b(2) = a(47): b(3) = a(28): b(4) = a(5): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(54): b(2) = a(35): b(3) = a(32): b(4) = a(9): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(58): b(2) = a(39): b(3) = a(20): b(4) = a(13): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(63): b(2) = a(44): b(3) = a(21): b(4) = a(2): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(51): b(2) = a(48): b(3) = a(25): b(4) = a(6): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(55): b(2) = a(36): b(3) = a(29): b(4) = a(10): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(59): b(2) = a(40): b(3) = a(17): b(4) = a(14): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(64): b(2) = a(41): b(3) = a(22): b(4) = a(3): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(52): b(2) = a(45): b(3) = a(26): b(4) = a(7): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(56): b(2) = a(33): b(3) = a(30): b(4) = a(11): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(60): b(2) = a(37): b(3) = a(18): b(4) = a(15): GoSub 860: If fl1 = 0 Then Return
                    
     b(1) = a(64): b(2) = a(43): b(3) = a(22): b(4) = a(1): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(52): b(2) = a(47): b(3) = a(26): b(4) = a(5): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(56): b(2) = a(35): b(3) = a(30): b(4) = a(9): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(60): b(2) = a(39): b(3) = a(18): b(4) = a(13): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(61): b(2) = a(44): b(3) = a(23): b(4) = a(2): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(49): b(2) = a(48): b(3) = a(27): b(4) = a(6): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(53): b(2) = a(36): b(3) = a(31): b(4) = a(10): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(57): b(2) = a(40): b(3) = a(19): b(4) = a(14): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(62): b(2) = a(41): b(3) = a(24): b(4) = a(3): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(50): b(2) = a(45): b(3) = a(28): b(4) = a(7): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(54): b(2) = a(33): b(3) = a(32): b(4) = a(11): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(58): b(2) = a(37): b(3) = a(20): b(4) = a(15): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(63): b(2) = a(42): b(3) = a(21): b(4) = a(4): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(51): b(2) = a(46): b(3) = a(25): b(4) = a(8): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(55): b(2) = a(34): b(3) = a(29): b(4) = a(12): GoSub 860: If fl1 = 0 Then Return
     b(1) = a(59): b(2) = a(38): b(3) = a(17): b(4) = a(16): GoSub 860: If fl1 = 0 Then Return
          
Return

End Sub


Vorige Pagina About the Author