Dim cn, mah(5, 5), n As Integer
Sub Main() '私は社長だ。
Console.WriteLine ("これはすべてのn次方陣を求めるソフトです。")
Console.WriteLine ("何次の方陣を発生させるのかをnに入力し")
Console.WriteLine ("エンターするとn次方陣がすべて生成されます。")
Console.Write ("n=")
n = Console.ReadLine()
cn = 0
Dim hj, ow As Object
hj = Timer
f (0)
ow = Timer
Console.WriteLine("魔方陣生成にかかった時間は{0:f6}秒です。", ow - hj)
Console.WriteLine("生成された{0:d}次方陣は{1:d}個です。", n, cn)
End Sub
Sub f(g As Integer)
Dim i, j, k, h, w, x, y As Integer
y = Int(g / n)
x = g Mod n
For i = 0 To n * n - 1
If g = 0 Then mah(y, x) = i + 1
h = 1
If g > 0 Then
For j = 0 To g - 1
If mah(Int(j / n), j Mod n) = i + 1 Then
h = 0
Exit For
End If
Next
If h = 1 Then mah(y, x) = i + 1
End If
If h = 1 Then
If x = n - 1 Then
w = 0
For j = 0 To n - 1
w += mah(y, j)
Next
If w <> Int(n * (n * n + 1) / 2) Then h = 0
End If
End If
If h = 1 Then
If y = n - 1 Then
w = 0
For j = 0 To n - 1
w += mah(j, x)
Next
If w <> Int(n * (n * n + 1) / 2) Then h = 0
End If
End If
If h = 1 Then
If y = n - 1 And x = 0 Then
w = 0
For j = 0 To n - 1
w += mah(j, n - 1 - j)
Next
If w <> Int(n * (n * n + 1) / 2) Then h = 0
End If
End If
If h = 1 Then
If y = n - 1 And x = n - 1 Then
w = 0
For j = 0 To n - 1
w += mah(j, j)
Next
If w <> Int(n * (n * n + 1) / 2) Then h = 0
End If
End If
If h = 1 Then
If g + 1 < n * n Then
f (g + 1)
Else
cn += 1
For j = 0 To n - 1
For k = 0 To n - 1
’二桁への対応
If mah(j, k) < 10 Then Console.Write(" {0:d} ", mah(j, k))
If mah(j, k) >= 10 Then Console.Write("{0:d} ", mah(j, k))
Next
Console.WriteLine()
Next
Console.WriteLine()
End If
End If
Next
End Sub
