コード例
Dim mah(4, 4) As Byte, x(16) As Byte, y(16) As Byte
Dim n As Byte, cn As Integer
Private Sub CommandButton1_Click()
n = Cells(3, 8)
cn = 0
zhy
ms (0)
Cells(4, 12) = cn
End Sub
Sub zhy()
Dim i As Byte
For i = 0 To n * n - 1
x(i) = i Mod n
y(i) = Int(i / n)
Next
End Sub
Sub ms(g As Byte)
Dim i As Byte, j As Byte, k As Byte, a As Byte, b As Byte, w As Byte
a = x(g)
b = y(g)
For i = 0 To n * n - 1
mah(b, a) = i + 1
If g > 0 Then
For j = 0 To g - 1
If mah(b, a) = mah(y(j), x(j)) Then GoTo tobi
Next
End If
If a = n - 1 Then
w = 0
For j = 0 To n - 1
w = w + mah(b, j)
Next
If w <> Int(n * (n * n + 1) / 2) Then GoTo tobi
End If
If b = n - 1 Then
w = 0
For j = 0 To n - 1
w = w + mah(j, a)
Next
If w <> Int(n * (n * n + 1) / 2) Then GoTo tobi
End If
If (a = n - 1) And (b = n - 1) Then
w = 0
For j = 0 To n - 1
w = w + mah(j, j)
Next
If w <> Int(n * (n * n + 1) / 2) Then GoTo tobi
End If
If (a = 0) And (b = n - 1) Then
w = 0
For j = 0 To n - 1
w = w + mah(j, n - 1 - j)
Next
If w <> Int(n * (n * n + 1) / 2) Then GoTo tobi
End If
If g + 1 < n * n Then
ms (g + 1)
Else
For j = 0 To n - 1
For k = 0 To n - 1
Cells(6 + j + Int(cn / 10) * (n + 1), 1 + k + (cn Mod 10) * (n + 1)) = mah(j, k)
Next
Next
cn = cn + 1
End If
tobi:
Next
End Sub
Private Sub CommandButton2_Click()
Rows("6:2000").Select
Selection.ClearContents
Cells(4, 12).Select
Selection.ClearContents
Cells(1, 1).Select
End Sub
実行画面
解説は次話でかなり詳しく行います。
今回のプログラムは力業で４次魔方陣が限界です。
もちろん理論的には、Dim mah(4, 4) As Byte, x(16) As Byte, y(16) As Byteの要素数を大きくとれば、
何次魔方陣でも可能ですが、
４次魔方陣でさえ、すべての計算が終わるまで５分程度要してしまいます。
５次魔方陣になると、作成を１００個に限定したとしても２０分は軽くかかってしまうでしょう。
私の作った最速の魔方陣作成ソフト魔方陣（細胞構成法）と比べると２６次魔方陣あたりなら作成速度は、
１０００兆分の１より遙かに低速です。
最速の魔方陣作成ソフト魔方陣（細胞構成法）なら、２６次魔方陣でも１秒で数百という単位で魔方陣を作り出してしまいますが、
今回のソフトでは、１００年コンピュータを走らせたとしても１個も見つけられないでしょう。
理由は簡単です。２６次魔方陣なら、
２６！＝１×２×３×・・・×２６＝4.03291461126606E+26（約40000000000000000000000000）通りのすべての場合を調べているからです。
これだとさすがのコンピュータも歯が立たないわけです。
従いまして、今回作ったプログラムは講義の進行とともに改良が加えられ、最終的には30次魔方陣を1秒で数百作り出すソフトまで変身していきます。
今回は本当にはじめの一歩です。
ダウンロード用参考ファイル
第２話へ 第４話へ
