Vorige Pagina Volgende Pagina About the Author

' Constructs 10 x 10 Borders for 6 x 6 Center Squares (Prime Numbers)
' Split Border Lines

' Tested with Office 2007 under Windows 7

Sub Priem10e1()

    Dim a1(2448), a(169), a10(100), b1(43300), b(43300), c(100)
 
    y = MsgBox("Locked", vbCritical, "Routine Priem10e1")
    End

    Sheets("Klad1").Select

    n5 = 0: n9 = 0: k1 = 1: k2 = 1
    ShtNm1 = "Pairs8"
     
    t1 = Timer

    For j100 = 362 To 7958

        Rcrd1a = j100

'       Read Prime Numbers From Sheet ShtNm1

        Pr3 = Sheets(ShtNm1).Cells(Rcrd1a, 1).Value
        Cntr3 = Pr3 / 2                                  'Center Elment (Hypothetical)
        s1 = 3 * Cntr3                                   'MC3
        s4 = 2 * Pr3                                     'MC4
        s10 = 5 * Pr3                                    'MC10
        nVar = Sheets(ShtNm1).Cells(Rcrd1a, 5).Value

        If nVar < 100 Then GoTo 1000

        Erase b1
        For j1 = 1 To nVar
            x = Sheets(ShtNm1).Cells(Rcrd1a, 6 + j1).Value
            b1(x) = x
        Next j1
       
'       Restore available pairs in a1()

        n10 = 0
        For j1 = 1 To pMax
            If b1(j1) <> 0 Then
                n10 = n10 + 1
                a1(n10) = b1(j1)
            End If
        Next j1
        m1 = 1: m2 = n10: n10 = 0
        If a1(1) = 1 Then m1 = 2: m2 = m2 - 1

'       Determine Corner Squares and Border Sections 2 x 4

        GoSub 2000                          'Determine 4 Corner Squares
        If n10 < 4 Then GoTo 910

        GoSub 3000                          'Determine 4 Border Sections
        If n10 < 8 Then GoTo 910
 
                             
        GoSub 800                           'Double Check Identical Integers a9()
        If fl1 = 1 Then
           n9 = n9 + 1: GoSub 1640          'Print results (lines)
'          n9 = n9 + 1: GoSub 1650          'Print results (squares)
        End If

910  Erase b1, b, c
1000 Next j100

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

End

'   Determine Corner Squares 3 x 3

2000 fl1 = 1

For j9 = m1 To m2                                                     'a(9)
If b1(a1(j9)) = 0 Then GoTo 90
If b(a1(j9)) = 0 Then b(a1(j9)) = a1(j9): c(9) = a1(j9) Else GoTo 90
a(9) = a1(j9)

For j8 = m1 To m2                                                     'a(8)
If b1(a1(j8)) = 0 Then GoTo 80
If b(a1(j8)) = 0 Then b(a1(j8)) = a1(j8): c(8) = a1(j8) Else GoTo 80
a(8) = a1(j8)

    a(7) = s1 - a(8) - a(9):
    If a(7) < a1(m1) Or a(7) > a1(m2) Then GoTo 70:
    If b1(a(7)) = 0 Then GoTo 70
    If b(a(7)) = 0 Then b(a(7)) = a(7): c(7) = a(7) Else GoTo 70
    
    a(6) = 4 * s1 / 3 - a(8) - 2 * a(9):
    If a(6) < a1(m1) Or a(6) > a1(m2) Then GoTo 60:
    If b1(a(6)) = 0 Then GoTo 60
    If b(a(6)) = 0 Then b(a(6)) = a(6): c(6) = a(6) Else GoTo 60

    a(5) = s1 / 3:
    If a(5) < a1(m1) Or a(5) > a1(m2) Then GoTo 50:
    ''If b1(a(5)) = 0 Then GoTo 50                                    'Might be Non Prime
    If b(a(5)) = 0 Then b(a(5)) = a(5): c(5) = a(5) Else GoTo 50
    
    a(4) = -2 * s1 / 3 + a(8) + 2 * a(9):
    If a(4) < a1(m1) Or a(4) > a1(m2) Then GoTo 40:
    If b1(a(4)) = 0 Then GoTo 40
    If b(a(4)) = 0 Then b(a(4)) = a(4): c(4) = a(4) Else GoTo 40
    
    a(3) = -s1 / 3 + a(8) + a(9):
    If a(3) < a1(m1) Or a(3) > a1(m2) Then GoTo 30:
    If b1(a(3)) = 0 Then GoTo 30
    If b(a(3)) = 0 Then b(a(3)) = a(3): c(3) = a(3) Else GoTo 30
    
    a(2) = 2 * s1 / 3 - a(8):
    If a(2) < a1(m1) Or a(2) > a1(m2) Then GoTo 20:
    If b1(a(2)) = 0 Then GoTo 20
    If b(a(2)) = 0 Then b(a(2)) = a(2): c(2) = a(2) Else GoTo 20
    
    a(1) = 2 * s1 / 3 - a(9):
    If a(1) < a1(m1) Or a(1) > a1(m2) Then GoTo 10:
    If b1(a(1)) = 0 Then GoTo 10
    If b(a(1)) = 0 Then b(a(1)) = a(1): c(1) = a(1) Else GoTo 10

        n10 = n10 + 1
        Select Case n10
        
            Case 1:
            
                 a10(1) = a(9):  a10(2) = a(7):  a10(3) = a(8):
                 a10(11) = a(3): a10(12) = a(1): a10(13) = a(2):
                 a10(21) = a(6): a10(22) = a(4):
                 
                 n32 = 9: GoSub 950                     'Remove used pairs from b1()
                 Erase b, c: GoTo 90
                 
            Case 2:
            
                 a10(8) = a(8):  a10(9) = a(9):  a10(10) = a(7):
                 a10(18) = a(2): a10(19) = a(3): a10(20) = a(1):
                                 a10(29) = a(6): a10(30) = a(4):
            
                 n32 = 9: GoSub 950                     'Remove used pairs from b1()
                 Erase b, c: GoTo 90
                 
            Case 3:
            
                                 a10(79) = a(6): a10(80) = a(4):
                 a10(88) = a(8): a10(89) = a(9): a10(90) = a(7):
                 a10(98) = a(2): a10(99) = a(3): a10(100) = a(1):
                 
                 n32 = 9: GoSub 950                     'Remove used pairs from b1()
                 Erase b, c: GoTo 90
                 
            Case 4:
            
                 a10(71) = a(6): a10(72) = a(4):
                 a10(81) = a(9): a10(82) = a(7): a10(83) = a(8):
                 a10(91) = a(3): a10(92) = a(1): a10(93) = a(2):
            
                 n32 = 9: GoSub 950                     'Remove used pairs from b1()
                 Erase b, c: Return
                 
        End Select
                          
5  b(c(1)) = 0: c(1) = 0
10 b(c(2)) = 0: c(2) = 0
20 b(c(3)) = 0: c(2) = 0
30 b(c(4)) = 0: c(4) = 0
40 b(c(5)) = 0: c(5) = 0
50 b(c(6)) = 0: c(6) = 0
60 b(c(7)) = 0: c(7) = 0
70 b(c(8)) = 0: c(8) = 0
80 Next j8
    
    b(c(9)) = 0: c(9) = 0
90 Next j9

   fl1 = 0
   Return
   
'   Determine Border Sections 2 x 4

3000 fl1 = 1

     For j1 = m1 To m2                                          'a(1)
     If b1(a1(j1)) = 0 Then GoTo 3010
     If b(a1(j1)) = 0 Then b(a1(j1)) = a1(j1): c(1) = a1(j1) Else GoTo 3010
     a(1) = a1(j1)

     a(5) = Pr3 - a(1): If b(a(5)) = 0 Then b(a(5)) = a(5): c(5) = a(5) Else GoTo 3050

     For j2 = m1 To m2                                          'a(2)
     If b1(a1(j2)) = 0 Then GoTo 3020
     If b(a1(j2)) = 0 Then b(a1(j2)) = a1(j2): c(2) = a1(j2) Else GoTo 3020
     a(2) = a1(j2)

     a(6) = Pr3 - a(2): If b(a(6)) = 0 Then b(a(6)) = a(6): c(6) = a(6) Else GoTo 3060

     For j3 = m1 To m2                                          'a(3)
     If b1(a1(j3)) = 0 Then GoTo 3030
     If b(a1(j3)) = 0 Then b(a1(j3)) = a1(j3): c(3) = a1(j3) Else GoTo 3030
     a(3) = a1(j3)

     a(7) = Pr3 - a(3): If b(a(7)) = 0 Then b(a(7)) = a(7): c(7) = a(7) Else GoTo 3070

     a(4) = s4 - a(3) - a(2) - a(1)
     If a(4) < a1(m1) Or a(4) > a1(m2) Then GoTo 3040
     If b1(a(4)) = 0 Then GoTo 3040
     If b(a(4)) = 0 Then b(a(4)) = a(4): c(4) = a(4) Else GoTo 3040

     a(8) = Pr3 - a(4): If b(a(8)) = 0 Then b(a(8)) = a(8): c(8) = a(8) Else GoTo 3080

        n10 = n10 + 1
        Select Case n10
        
            Case 5:
            
                a10(4) = a(1):  a10(5) = a(2):  a10(6) = a(3):  a10(7) = a(4):
                a10(14) = a(5): a10(15) = a(6): a10(16) = a(7): a10(17) = a(8):
            
                n32 = 8: GoSub 950                     'Remove used pairs from b1()
                Erase b, c: GoTo 3010
            
            Case 6:
            
                a10(39) = a(1): a10(40) = a(5):
                a10(49) = a(2): a10(50) = a(6):
                a10(59) = a(3): a10(60) = a(7):
                a10(69) = a(4): a10(70) = a(8):
            
                n32 = 8: GoSub 950                     'Remove used pairs from b1()
                Erase b, c: GoTo 3010
            
            Case 7:
            
                a10(84) = a(1): a10(85) = a(2): a10(86) = a(3): a10(87) = a(4):
                a10(94) = a(5): a10(95) = a(6): a10(96) = a(7): a10(97) = a(8):
            
                n32 = 8: GoSub 950                     'Remove used pairs from b1()
                Erase b, c: GoTo 3010
            
            Case 8:
            
                a10(31) = a(1): a10(32) = a(5):
                a10(41) = a(2): a10(42) = a(6):
                a10(51) = a(3): a10(52) = a(7):
                a10(61) = a(4): a10(62) = a(8):
            
                n32 = 8: GoSub 950                     'Remove used pairs from b1()
                Erase b, c: Return
        
        End Select

     b(c(8)) = 0: c(8) = 0
3080 b(c(4)) = 0: c(4) = 0
3040 b(c(7)) = 0: c(7) = 0
3070 b(c(3)) = 0: c(3) = 0
3030 Next j3

     b(c(6)) = 0: c(6) = 0
3060 b(c(2)) = 0: c(2) = 0
3020 Next j2

     b(c(5)) = 0: c(5) = 0
3050 b(c(1)) = 0: c(1) = 0
3010 Next j1

    fl1 = 0
    Return

'   Double Check Identical Numbers a10()

800 fl1 = 1
    For i1 = 1 To 100
       a20 = a10(i1): If a20 = 0 Then GoTo 810
       For i2 = (1 + i1) To 100
           If a20 = a10(i2) Then fl1 = 0: Return
       Next i2
810 Next i1
    Return

'   Remove used pairs from b1()

950 For i1 = 1 To n32
        b1(a(i1)) = 0
    Next i1
    Return

'    Print results (lines)

1640 Cells(n9, 1).Select
     For i1 = 1 To 100
         Cells(n9, i1).Value = a10(i1)
     Next i1
     Cells(n9, 101).Value = s10     'MC10
     Cells(n9, 102).Value = Rcrd1a  'Record
     Return

'    Print results (squares)

1650 n2 = n2 + 1
     If n2 = 3 Then
         n2 = 1: k1 = k1 + 11: k2 = 1
     Else
         If n9 > 1 Then k2 = k2 + 11
     End If

     Cells(k1, k2 + 1).Select
     Cells(k1, k2 + 1).Font.Color = -4165632
     Cells(k1, k2 + 1).Value = "MC = " + CStr(s10)
    
     i3 = 0
     For i1 = 1 To 10
         For i2 = 1 To 10
             i3 = i3 + 1
             Cells(k1 + i1, k2 + i2).Value = a10(i3)
         Next i2
     Next i1

     Return

End Sub

Vorige Pagina Volgende Pagina About the Author