用两种方法求解九宫算问题

问题：

在一个3*3的方格内填入9个数字，数字的范围从1~9，每个数字仅出现一次，使得每行、每列以及两个对角线的数字之和为15。

 一、VBA程序

算法思路：

1、将9个数字分成3组，每组3个数字，使得每组数字之和为15。

2、将第二个数字为5个分组单独选出，作为方格的第二行数组。

3、循环抽取其它的分组，与第二分组进行重复性判断。

4、如果没有重复，记录这3种分组。

5、将这3种分组传人判断函数，核实是否满足条件。

6、如果满足条件，将结果输出。

 VBA源程序如下：

运行结果

6  7  2

Sub jiugongsuan_All()

Dim cs(1 To 9) As Integer
Dim jgs1(1 To 3) As Integer
Dim Ot() As Integer
'第一、三行数组
Dim Sec() As Integer
'第二行数组
Dim temp1(1 To 100, 1 To 3) As Integer
Dim temp2(1 To 100, 1 To 3) As Integer
'temp1、temp2分别为临时数组
Dim Totali As Integer
'第二个数为5，三个数之和为15的数组个数
Dim Totaln As Integer
'不包括第二个数为5，三个数之和为15的数组个数
Dim at(1 To 3) As Integer
Dim bt(1 To 3) As Integer
Dim ct(1 To 3) As Integer
'at、bt、ct分别为中间变量数组

Dim scjg(1 To 3, 1 To 3) As Integer
'scjg为结果数组


Dim i As Integer

Totali = 0
Totaln = 0

For i = 1 To 9
'初始数组赋值
cs(i) = i

Next

For m = 1 To 9
 If cs(m) <> 0 Then
'   判断该数字是否已被占用
   jgs1(1) = cs(m)
   cs(m) = 0
'   占用该数字
    
   For n = 1 To 9
    If cs(n) <> 0 Then
       jgs1(2) = cs(n)
        cs(n) = 0
        
      For q = 1 To 9
        If cs(q) <> 0 Then
         jgs1(3) = cs(q)
          cs(q) = 0
     
       
            If jgs1(1) + jgs1(2) + jgs1(3) = 15 Then
      
             ' Debug.Print jgs1(1); jgs1(2); jgs1(3)
                 
               If jgs1(2) = 5 Then
               '第二行数据
                 Totali = Totali + 1
                temp1(Totali, 1) = jgs1(1)
                temp1(Totali, 2) = jgs1(2)
                temp1(Totali, 3) = jgs1(3)
               Else
'                第一行，第三行数据
                 Totaln = Totaln + 1
                temp2(Totaln, 1) = jgs1(1)
                temp2(Totaln, 2) = jgs1(2)
                temp2(Totaln, 3) = jgs1(3)
              
                     
              
              End If
              
              
                 
             End If
      
            cs(q) = q
            End If
       
       Next
       
         cs(n) = n
       End If
       
       
       Next
     
     cs(m) = m
 End If
 

 
Next

ReDim Sec(1 To Totali, 1 To 3)
ReDim Ot(1 To Totaln, 1 To 3)


For m1 = 1 To Totali
Sec(m1, 1) = temp1(m1, 1)
Sec(m1, 2) = temp1(m1, 2)
Sec(m1, 3) = temp1(m1, 3)

'第二行数组赋值

Next

For n1 = 1 To Totaln
Ot(n1, 1) = temp2(n1, 1)
Ot(n1, 2) = temp2(n1, 2)
Ot(n1, 3) = temp2(n1, 3)

Next

'第一、三行数组赋值

For m1 = 1 To Totali

    bt(1) = Sec(m1, 1)
    bt(2) = Sec(m1, 2)
    bt(3) = Sec(m1, 3)
   
     For n1 = 1 To Totaln
       at(1) = Ot(n1, 1)
       at(2) = Ot(n1, 2)
       at(3) = Ot(n1, 3)
       
     If cfpd(at, bt) Then
 
           
      Else
    
      For k1 = 1 To Totaln
       ct(1) = Ot(k1, 1)
       ct(2) = Ot(k1, 2)
       ct(3) = Ot(k1, 3)
       
      
      If Not cfpd(at, ct) And Not cfpd(bt, ct) Then
'          第一、三行，第二、三行重复情况判断
     
      scjg(1, 1) = at(1)
      scjg(1, 2) = at(2)
      scjg(1, 3) = at(3)
      scjg(2, 1) = bt(1)
      scjg(2, 2) = bt(2)
      scjg(2, 3) = bt(3)
      scjg(3, 1) = ct(1)
      scjg(3, 2) = ct(2)
      scjg(3, 3) = ct(3)

'       调用函数
       finalpd (scjg)
    
    Else
     
   
    End If
 
        
     
     Next
     
     End If
       
    Next

Next


End Sub

Function cfpd(a() As Integer, b() As Integer) As Boolean
'2个数组重复情况判断


 If b(1) <> a(1) And b(1) <> a(2) And b(1) <> a(3) Then
  cfpd = False
 Else
  cfpd = True
  Exit Function
  
 End If
 
 If b(2) <> a(1) And b(2) <> a(2) And b(2) <> a(3) Then
 cfpd = False
 Else
 cfpd = True
 Exit Function
 
 End If
 
 If b(3) <> a(1) And b(3) <> a(2) And b(3) <> a(3) Then
 cfpd = False
Else
 cfpd = True
 Exit Function
 
 
End If

End Function


Sub finalpd(sz As Variant)
'列求和、对角线求和判断

Dim i%, j%
Dim myjgs(1 To 3, 1 To 3)
Dim Hpd As Boolean
Dim Lpd As Boolean
Dim DJXpd As Boolean


For i = 1 To 3
 For j = 1 To 3
 myjgs(i, j) = sz(i, j)

 Next
Next


'列判断
For i = 1 To 3

If myjgs(1, i) + myjgs(2, i) + myjgs(3, i) = 15 Then
Lpd = True
Else
Lpd = False
Exit For
End If

Next

'对角线判断

If myjgs(1, 1) + myjgs(2, 2) + myjgs(3, 3) = 15 And myjgs(1, 3) + myjgs(2, 2) + myjgs(3, 1) = 15 Then
 DJXpd = True

Else
 DJXpd = False
End If


If Lpd And DJXpd Then

'结果输出

Debug.Print myjgs(1, 1); myjgs(1, 2); myjgs(1, 3)

Debug.Print myjgs(2, 1); myjgs(2, 2); myjgs(2, 3)

Debug.Print myjgs(3, 1); myjgs(3, 2); myjgs(3, 3)

Debug.Print "-------------------------"

End If


End Sub

运行结果：

6  7  2

1  5  9

8  3  4

-------------------------

 8 3  4

 1 5  9

 6 7  2

-------------------------

 4 9  2

 3 5  7

 8 1  6

-------------------------

 8 1  6

 3 5  7

 4 9  2

-------------------------

 2 9  4

 7 5  3

 6 1  8

-------------------------

 6  1  8

 7 5  3

 2 9  4

-------------------------

 2 7  6

 9 5  1

 4 3  8

-------------------------

 4 3  8

 9 5  1

 2 7  6 

  

二、EXCEL规划求解

1、选择数据菜单-à“规划求解”

2、设置可变单元格及约束条件

约束条件

Ø  行之和为15

Ø  列之和为15

Ø  所有数字均为整数

Ø  所有数字>=1

Ø  所有数字<=9

Ø  所有数字连乘的积为9的阶乘，这种方式可以限制每个数字出现1次，且仅出现1次。

 

图一

图二

这种方法只能求出一个解。