階優化

8k 階優化 _Franklin 型幻 方的 快速構作方法

梁培基 · ^張航輔 · ^張俠輔

摘要: 優化 Franklin 型幻方是一種具有許多特殊性質的幻方。 富蘭克林幻方是迄今 為止奇妙性質最多的幻方, 《有趣的數論》 一書 (潘承彪譯, 北大出版社出版) 稱為

「最神奇的幻方」而享譽國際, 開 「曲線幻方」 研究之先河, 深受幻方愛好者所崇敬。

本文給出一種快速構作這類幻方的方法, 並將其公式化, 用之能計算出這類幻方 任意位置處的元素, 因此也給使用電腦來構作這類幻方提供一條便利的途徑。

關鍵字: 幻方, 全對稱幻方, Franklin 型幻方, 優化 Franklin 型幻方。

一、 引言

本傑明．富蘭克林 (Benjamin Franklin, 1706∼1790), 著名政治活動家, 發明家, 他對 幻方, 幻圓都有濃厚的興趣, 進行過潛心研究, 留傳下兩個後人稱之為 「Franklin 幻方」 的 8 階和 16 階的方陣 [1]。 這兩個幻方具有普通幻方所沒有的許多奇妙的特性。 但遺憾的是這兩個 幻方的兩條對角線上元素之和不等於幻方常數, 因此嚴格說來不能算作幻方, 這點白璧之瑕影 響了富蘭克林幻方的完美性, 但比起創造一個幻方 「新品種」 的豐功偉績, 還是微乎其微的。 雖 然, 他沒有把構作這類方陣的方法傳給後人。

1990 年, 筆者在 《數學傳播》 十四卷二期 「多功能組合幻方」 一文中, 發表了一個具有

「Franklin 幻方」 的全部性質, 並且兩條對角線及各左、 右折斷對角線上元素之和也都等於該幻 方常數的 8 階幻方。 我們稱這類幻方為 「優化 Franklin 幻方」, 簡稱為 「優化 F 型幻方」[2]。

二、 8k 階優化 F 型幻方的快速構作方法

定義: 如果由連續自然數 1, 2, . . . , n² 排列一個 n = 8k (k = 1, 2, . . .) 階方陣具有下列各條 性質, 則稱為 8k 階優化 F 型幻方:

1. 每行、 每列 n 個元素之和都等於定值 C;

2. 每條左、 右對角線及折斷對角線上 n 個元素之和都等於定值 C;

3. 將每個行、 列等分為兩半, 每一半上 n/2 個元素之和都等於 C/2;

4. 將兩條左、 右對角線分別等分為兩半, 每一半上 n/2 個元素之和都等於 C/2;

5. 在任意 m (2 ≤ m ≤ 2k) 個行上, 沿左右對稱的任一橫向波狀折線 「 」 或 「 」 線上 的 n 個元素之和都等於 C, 稱為 「行 W 性質」(因其形狀像字母 W, 故名);

以 8 階優化佛蘭克林幻方的 「行 W 線」 為例:

(1) 從第 1 行第 1 列開始向右下方斜行前進, 到第 2 行第 2 列。 之後, 平移至第 2 行的第 3 列位置上;

(2) 向右上方斜行前進到第 1 行的第 4 列。 之後, 平移至第 1 行的第 5 列位置上;

(3) 向右下方斜行前進到第 2 行的第 6 列, 之後, 平移到第 2 行第 7 列位置上;

(4) 向右上方斜行前進到第 1 行的第8 列, 完成了 「行 W 線」 的過程。

需要注意的是, 到達轉彎處, 兩個元素並列。 對於 「列 W 線」, 可參照圖 1。

行、 列四通 V 形線 行、 列四通 W 形線 行、 列與半行, 半列及 2 × 2 小方

S⁸ = 260 S⁸ = 260 S⁸ = 260, S⁴ = 130 圖 1

6. 在任意 m (2 ≤ m ≤ 2k) 個列上, 沿左右對稱的任一橫向波狀折線 或 線上的 n 個 元素之和都等於 C, 稱為 「列 W 性質」;

上述各性質中的 C 均為幻方常數 n(n²+ 1)/2。 以上 6 條性質將在後面給出證明。

另外還有三條有趣的性質: 證明從略。

(1) 「四通 V」 性質。 例如: 從幻方左上角第一行、 第一列 (i = 1, j = 1) 開始向 i + 1, j + 1 (右下方) 的方向斜行前進, 斜行到第 n/2 行、 第 n/2 列; 之後, 平移到第 n/2 + 1列; 開 始向 i − 1, j + 1 (右上方) 斜行前進, 斜行到第 1 行的第 n 列為止。 所經過的路線形狀 像字母 「V」, 其 「V」 形線上 n 個元素之和等於幻和, 且滿足四個方向: 「 , , , 」 故 稱為 「四通 V」 性質。 可以參照圖 1 畫出其餘的 「四通 V 形線」。 與 「W 形線」 不同的是,

「V 形線」 只有 1 個轉彎處, 而 「W 形線」 有 3 個轉彎處。

(2) 在任意偶數×偶數小方陣中, 偶數 2 個元素之和都等於定值的性質。

(3) 在 8 階優化 F 型幻方中, 某個 2 × 2 小方陣 4 個數的平方和與另外一個 2 × 2 小方陣 4 個數的平方和之和等於 8 階幻方平方幻和。 即 S8² = 11180。 (詳見文尾, 16 個 2 × 2 小方 陣平方和表)。

當 k = 1 時, 8 階優化 F 型幻方的上述性質用圖形表示如圖 1。

我們知道, 普通幻方僅具有上述性質 1 與性質 2 的前半部分, 全對稱幻方也只具有性質 1 與性質 2 的全部。 可見, 優化 F 型幻方是一種具有更多特殊性質的全對稱幻方。 然而, 這也使 得人們對它的構作顯得更加困難。 在本文中, 我們給出一種構作這類幻方的方法, 利用它可以迅 速作出任意 8k (k = 1, 2, . . .) 階的優化 F 幻方。 這種方法建立在下述定理的基礎上:

定理1: 設 A = (ai,j), B = (bi,j) 為自然數集 {1, 2, . . . , n} (n = 8k, k = 1, 2, . . .) 上的 兩個相互正交 (即 A、 B 的相同位置處的元素對構成的 n² 個有序偶兩兩互異) 的方陣, 如果 A、 B 分別具有優化 F 型幻方的各條性質 [對於 A、 B 各性質中的定值 C = n(n + 1)/2 = 4k(8k + 1)], 令

hi,j = n(ai,j − 1) + bi,j (i, j = 1, 2, . . . , n), 則 H = (h_i,j) 為 8k 階優化 F 型幻方。

證明: 根據定理條件, H 中對應於 8k 階優化 F 型幻方的性質 1, 2, 5, 6 中的 n = 8k 個元 素的和分別為

X

(n)

hi,j= n X

(n)

ai,j − n +X

(n)

bi,j

= nhn(n + 1) 2 − ni

+n(n + 1)

2 = n(n²+ 1) 2 . H = (hi,j) 中對應於性質 3, 4 中的 n/2 個元素的和為

X

(n/2)

hi,j= n X

(n/2)

ai,j − n/2 + X

(n/2)

bi,j

= nhn(n + 1)

4 − n/2i

+ n(n + 1)

4 = n(n²+ 1) 4 . 從而, H = (h_i,j) 滿足 8k 階優化 F 型幻方的各條性質。

另外, 筆者在文 [3] 的定理 1 中已證, 由 A、B 的正交性可知 H 中的元素恰由 n² 個互 異的連續自然數 1, 2, . . . , n² 組成, 因此, 根據定義, H 為 8k 階優化 F 型幻方。

由定理 1 可知, 我們只要找到滿足該定理條件的兩個 8k 階方陣 A 與 B, 即可構作出一 個 8k 階優化 F 型幻方 H.

下面, 給出一種構作這樣的 A, B 方陣的方法, 它可分為如下幾個步驟 [在下文中, 用 A_p,q (p, q = 1, 2, . . . , 2k) 表示 A 中的 4 階子方陣, 用 a^p,q_r,s (r, s = 1, 2, 3, 4) 表示 A_p,q 中的元 素];

1. 按照下列順序依次填寫出 A 的前 4 列上的 k 個 4 階子方陣 A¹,1, A²,1, . . . , Ak,1 中的前 2 行:

A

1 1 4k 1 4k 4k +1 1 4k 4k +1 8k

2 4k −1 2 4k +2 4k −1 2 8k −1 4k +2 4k −1 2

A

3 3 4k −2 3 4k −2 4k +3 3 4k −2 4k +3 8k −2

4 4k −3 4 4k +4 4k −3 4 8k −3 4k +4 4k −3 4

. . .

. .

. ⇒

. .

. ⇒

. .

. ⇒

. . .

A

2k −3 2k −3 2k +4 2k −3 2k +4 6k −3 2k −3 2k +4 6k −3 6k +4

2k −2 2k +3 2k −2 6k −2 2k +3 2k −2 6k +3 6k −2 2k +3 2k −2 A

2k −1 2k −1 2k +2 2k −1 2k +2 6k −1 2k −1 2k +2 6k −1 6k +2

2k 2k +1 2k 6k 2k +1 2k 6k +1 6k 2k +1 2k

2. 在上述各子陣內, 將第 1, 2 行的 1, 3 列分別與 2, 4 列交換後作為第 3, 4 行, 即

Ap,1 =













a^p,11,1 a^p,11,2 a^p,11,3 a^p,11,4

a^p,12,1 a^p,12,2 a^p,12,3 a^p,12,4

a^p,11,2 a^p,11,1 a^p,11,4 a^p,11,3

a^p,12,2 a^p,12,1 a^p,12,4 a^p,12,3













(p = 1, 2, . . . , k).

3. 將 Ap,1 (p = 1, 2, . . . , k) 的第 1, 2, 3, 4 行分別為 Ap+k, 1 的第 3, 4, 1, 2 行, 得 出 A 的前 4 列上其餘 k 個 4 階子方陣 A_k+1,1, Ak+2,1, . . . , A²k,1 即對於 p = 1, 2, . . . , k, s = 1, 2, 3, 4, 有

a^p+k,1_r,s = a^p,1_r+2,s (r = 1, 2), a^p+k,1_r,s = a^p,1_r−2,s (r = 3, 4).

4. 將 A 的第 1 ∼ 4 列完全不變的複製成第 5 ∼ 8, 9 ∼ 12, . . ., 8k − 3 ∼ 8k 列, 即 A_p,q = A_p,1 (p = 1, 2, . . . , 2k; q = 2, 3, . . . , 2k),

這樣就造出了 8k 階方陣 A。

5. 將 A 轉置, 就得到 B 陣。 即 B = A^T。

利用此方法作出的 A、B 方陣, 就能按照定理 1 中的公式 h_i,j = 8k(ai,j − 1) + bi,j

(i, j = 1, 2, . . . , 8k) 計算出 H = (hi,j) 的每個元素。 也可用一種更直觀更簡便的方法, 即利

用自然數方陣 Z = z(i, j) (i, j = 1, 2, . . . , 8k):

Z =













1 2 · · · 8k 8k + 1 8k + 2 · · · 16k

... ... . .. ...

· · · (8k)²













及 h_i,j = Zai,jbi,j (i, j = 1, 2, . . . , 8k), 由 A、B 迅速填寫出 8k 階優化 F 型幻方 H = (hi,j)。 這種方法可由定理 1 推出 [3], 我們稱之為 「方陣定位法」。

下面, 給出按照上述方法構作的 8 階和 16 階優化 F 型幻方的實例及它們的 A 陣。

8 階:

16 階 A 陣

在 16 階優化 F 型幻方 H 陣中, 我們用實線標示行 V 形線與列 V 形線, 用虛線標示行 W 形線與列 W 形線。

根據上面曲線的畫法, 可以方便地畫出其餘的行、 列四通 V 形線及四通 W 形線, 從略。

事實上, 上述構作 A 方陣的方法可以用下列計算公式表示出來:

A_p,q=



























































2p−1 4k−2p+2 4k+2p−1 8k−2p+2 8k−2p+1 4k+2p 4k−2p+1 2p

4k−2p+2 2p−1 8k−2p+2 4k+2p−1 4k+2p 8k−2p+1 2p 4k−2p+1













p = 1, 2, . . . , k q = 1, 2, . . . , 2k

!













4k−2(p−k)+2 2(p−k)−1 8k−2(p−k)+2 4k+2(p−k)−1 4k+2(p−k) 8k−2(p−k)+1 2(p−k) 4k−2(p−k)+1 2(p−k)−1 4k−2(p−k)+2 4k+2(p−k)−1 8k−2(p−k)+2 8k−2(p−k)+1 4k+2(p−k) 4k−2(p−k)+1 2(p−k)











 p = k+1, k+2, . . . , 2k

q = 1, 2, . . . , 2k

!

(1)

或者等價地用下列計算公式表示 A_i,j 出來:

j = 1, 2, . . . , 8k

j ≡ 1(mod 4) j ≡ 2(mod 4) j ≡ 3(mod 4) j ≡ 1(mod 4) i+1

2 4k− i−1

2 4k+ i+1

2 8k− i−1

2 i ≡ 1(mod 4)

i = 1, 2, . . . , 4k 8k− i

2 4k+ i

2 +1 4k− i 2

i+2

2 i ≡ 2(mod 4) 4k− i−3

2

i−1

2 8k− i−3

2 4k− i−1

2 i ≡ 3(mod 4) 4k+ i

2 8k− i

2 +1 i

2 4k− i

2 +1 i ≡ 4(mod 4) 6k− i−1

2

i+1

2 −2k 10k− i−1

2 2k+ i+1

2 i ≡ 1(mod 4) 2k+ i+2

2 10k− i 2

i+2

2 −2k 6k− i

2 i ≡ 2(mod 4) i = 4k+1, 4k+2, . . . , 8k i−1

2 −2k 6k− i−3

2 2k+ i−1

2 10k− i−3

2 i ≡ 3(mod 4) 10k− i

2 +1 2k+ i

2 6k− i

2 +1 i

2 −2k i ≡ 4(mod 4)

(2) 公式 (2) 提供了一種計算 8k 階優化 F 型幻方在任意指定行列位置處的元素的方法, 利 用它可以很方便地編制程式在電腦上構作出任意 8k 階的優化 F 型幻方。

三、 8k 階優化 F 型幻方構作方法的證明

按照上述方法構作的 A、B 方陣顯然在自然數集 {1, 2, . . . , n} 上, 根據定理 1, 只須證明 A、B 相互正交, 並具有 8k 階優化 F 型幻方的各條性質即可。

證明 A、B 相互正交。

設 (a^m,n_i,j , b^m,n_i,j ) = (a^p,q_k,l, b^p,q_k,l), 其中 a^m,n_i,j 與 a^p,q_k,l 分别位於 A 之 4 階子方陣 A_m,n 與 Ap,q 中, b^m,n_i,j 與 b^p,q_k,l 分別位於 B 之 4 階子方陣 B_m,n 與 B_p,q 中, i, j, k, l ∈ {1, 2, 3, 4}, p, q, m, n ∈ {1, 2, . . . , 2k}。

由 A 的構作步驟 4 知 A_α,β = Aα,1, 其中 α = 1, 2, . . . , 2k; β = 2, 3, . . . , 2k。 因此由 假設知

a^m,1_i,j = a^m,n_i,j = a^p,q_k,l = a^p,1_k,l, (3) 由 A 的構作步驟 1、2、3 知, 在列標為 1 的 2k 個 4 階子方陣中恰有 k 對子方陣 A_p1,1

與 A_p2,1 (p¹ = 1, 2, . . . , k; p² = p¹+ k) 具有相同的元素, 而其餘任意兩個子方陣內均無相 同元素。 在每對具有相同元素的子陣 A_p1,1 與 A_p2,1 內恰存在下列 8 組相同元素

a^p1¹,1^,1 = a^p3,2^1,1= a^p1²,2^,1= a^p3,1^2,1, a^p1¹,2^,1 = a^p3¹,1^,1 = a^p1²,1^,1 = a^p3²,2^,1, a^p1¹,3^,1 = a^p3,4^1,1= a^p1²,4^,1= a^p3,3^2,1, a^p1¹,4^,1 = a^p3¹,3^,1 = a^p1²,3^,1 = a^p3²,4^,1, a^p2¹,1^,1 = a^p4,2^1,1= a^p2²,2^,1= a^p4,1^2,1, a^p2¹,2^,1 = a^p4¹,1^,1 = a^p2²,1^,1 = a^p4²,1^,1, a^p2¹,3^,1 = a^p4,4^1,1= a^p2²,4^,1= a^p4,3^2,1, a^p2¹,4^,1 = a^p4¹,3^,1 = a^p2²,3^,1 = a^p4²,4^,1, 而其餘任意兩個元素必定互不相同。 上面各組元素之下標所成集合如下:

{(1, 1), (3, 2), (1, 2), (3, 1)}, {(1, 3), (3, 4), (1, 4), (3, 3)}, {(2, 1), (4, 2), (2, 2), (4, 1)}, {(2, 3), (4, 4), (2, 4), (4, 3)}。

(4)

由 (3) 式知, m ≡ p(mod k), 且 (i, j) 與 (k, l) 必同屬於 (4) 其中一組集合。 因為 B = A^T, 故 b^m,n_i,j = a^n,m_j,i , b^p,q_k,l = b^q,p_l,k, 因此由假設知

a^n,1_j,i = a^n,m_j,i = b^m,n_i,j = b^p,q_k,l = a^q,p_l,k = a^p,1_l,k, 故 n ≡ q(mod k), 且 (j, i) 與 (l, k) 也必同屬於 (4) 其中一組集合。

然而 (4) 之 4 組集合中, 滿足 (i, j) 與 (k, l) 同屬於其中一組集合, 且 (j, i) 與 (l, k) 也 同屬於其中一組集合, 唯有 i = k, j = l。

因此　(a^m,n_i,j , b^m,n_i,j ) = (a^p,q_i,j, b^p,q_i,j), m ≡ p(mod k) 且 n ≡ q(mod k), 由 A 的構作步驟 3 知 m = p 且 n = q。 故 A 與 B 相互正交。

2. 證明 A、B 有與 8k 階優化 F 型幻方相同的各條性質 (對於 A、B 常數 C = 4k(8k + 1))。

1) 由 A 的構作步驟 4 知 Ap,1 = Ap,2 = · · · = Ap,2k (p = 1, 2, . . . , 2k), 從而, 當 i = 1, 2, . . . , 4k 時, 由 (1) 式可得

8k

X

j=1

ai,j =















2k(2p − 1 + 4k − 2p + 2 + 4k + 2p − 1 + 8k − 2p + 2) = 4k(8k + 1) (i ≡ 1(mod 4)), 2k(8k − 2p + 1 + 4k + 2p + 4k − 2p + 1 + 2p) = 4k(8k + 1)

(i ≡ 2(mod 4)).

由 A 的構作步驟 2 知, 當 i ≡ 3, 0(mod 4) 時, A 的各行中的元素集合分別與 i ≡ 1, 2(mod 4) 中的相同; 由 A 的構作步驟 3, 4 知, A 的第 4k + 1, 4k + 2, . . . , 8k 行中的每個

行內的元素集合分別與第 1, 2, . . . , 4k 行中一個行內的元素集合相同, 因而 A 的各行元素之和 均為 4k(8k + 1)。

另一方面, 對於 j = 1, 2, . . . , 8k, 由 (1) 式可得

4k

X

i=1

ai,j=























k

P

p=1

(2p−1+8k−2p+1+4k−2p+2+4k+2p) =

k

P

p=1

(16k + 2) = 2k(8k + 1) (j ≡ 1(mod 4)),

k

P

p=1

(4k+2p−1+4k−2p+1+8k−2p+2+2p) =

k

P

p=1

(16k+2) = 2k(8k+1) (j ≡ 3(mod 4)).

由 A 的構作步驟 2 知, 在 A 的第 1 ∼ 4k 行上, 當 j ≡ 2, 0 (mod 4) 時, A 的各列元素集 合分別與 j ≡ 1, 3 (mod 4) 的相同, 因而對於任意 j = 1, 2, . . . , 8k, 均有

4k

X

i=1

ai,j = 2k(8k + 1).

由 A 的構作步驟 3 知, A 的各列中的下半部分的元素集合與上半部分相同, 因此, A 的各個 半列元素之和均為 2k(8k + 1), 各列元素之和均為 4k(8k + 1)。

從而, 已經證明了 A 滿足性質 1, 3.

2) A 的左對角線及各左折斷對角線上 8k 個元素之和為

Lm =

8k−m

X

i=1

ai,i+m+

8k

X

i=8k+1−m

ai,i+m−8k (m = 0, 1, . . . , 8k − 1).

由 Ap,q = Ap1 (p, q = 1, 2, . . . , 2k), 從 (1) 式可知, 當 m ≡ 0 (mod 4) 時,

Lm=

2k

X

p=1 4

X

r=1

a^p,1_r,r

=

k

X

p=1

(2p−1+4k+2p+8k−2p+2+4k−2p+1)

+

2k

X

p=k+1

[4k−2(p−k)+2+8k−2(p−k)+1+4k+2(p−k)−1+2(p−k)]

=

k

X

p=1

(16k + 2) +

2k

X

p=k+1

(16k + 2) = 4k(8k + 1).

當 m ≡ 1 (mod 4) 時

Lm=

2k

X

p=1

(a^p,11,2+ a^p,12,3+ a^p,13,4+ a^p,14,1)

=

k

X

p=1

(4k−2p+2+4k−2p+1+4k+2p−1+4k+2p)

+

2k

X

p=k+1

[2(p−k)−1+2(p−k)+8k−2(p−k)+2+8k−2(p−k)+1]

=

k

X

p=1

(16k + 2) +

2k

X

p=k+1

(16k + 2) = 4k(8k + 1).

當 m ≡ 2 (mod 4) 時

Lm=

2k

X

p=1

(a^p,11,3+ a^p,12,4+ a^p,13,1+ a^p,14,2)

=

k

X

p=1

(4k+2p−1+2p+4k−2p+2+8k−2p+1)

+

2k

X

p=k+1

[8k−2(p−k)+2+4k−2(p−k)+1+2(p−k)−1+4k+2(p−k)]

=

k

X

p=1

(16k + 2) +

2k

X

p=k+1

(16k + 2) = 4k(8k + 1).

當 m ≡ 3 (mod 4) 時

Lm=

2k

X

p=1

(a^p,11,4+ a^p,12,1+ a^p,13,2+ a^p,14,3)

=

k

X

p=1

(8k−2p+2+8k−2p+1+2p−1+2p)

+

2k

X

p=k+1

[4k+2(p−k)−1+4k+2(p−k)+4k−2(p−k)+2+4k−2(p−k)+1]

=

k

X

p=1

(16k + 2) +

2k

X

p=k+1

(16k + 2) = 4k(8k + 1).

A 的右對角線及各右折斷對角線上 8k 個元素之和為 Rm =

8k−m

X

i=1

ai,8k+1−i−m+

8k

X

i=8k+1−m

ai,16k+1−i−m (m = 0, 1, . . . , 8k − 1).

同樣的, 由 A_p,q = Ap (p, q = 1, 2, . . . , 2k), 及 (1) 式可得, 當 m ≡ 0 (mod 4) 時,

Rm=

2k

X

p=1 4

X

r=1

a^p_r,5−r

=

k

X

p=1

(8k−2p+2+4k−2p+1+2p−1+4k+2p)

+

2k

X

p=k+1

[4k+2(p−k)−1+2(p−k)+4k−2(p−k)+2+8k−2(p−k)+1]

=

k

X

p=1

(16k + 2) +

2k

X

p=k+1

(16k + 2) = 4k(8k + 1).

當 m ≡ 1 (mod 4) 時

Rm=

2k

X

p=1

(a^p,11,3+ a^p,12,2+ a^p,13,1+ a^p,14,4)

=

k

X

p=1

(4k+2p−1+4k+2p+4k−2p+2+4k−2p+1)

+

2k

X

p=k+1

[8k−2(p−k)+2+8k−2(p−k)+1+2(p−k)−1+2(p−k)]

=

k

X

p=1

(16k + 2) +

2k

X

p=k+1

(16k + 2) = 4k(8k + 1).

當 m ≡ 2 (mod 4) 時

Rm=

2k

X

p=1

(a^p,11,2+ a^p,12,1+ a^p,13,4+ a^p,14,3)

=

k

X

p=1

(4k−2p+2+8k−2p+1+4k+2p−1+2p)

+

2k

X

p=k+1

[2(p−k)−1+4k+2(p−k)+8k−2(p−k)+2+4k−2(p−k)+1]

=

k

X

p=1

(16k + 2) +

2k

X

p=k+1

(16k + 2) = 4k(8k + 1).

當 m ≡ 3 (mod 4) 時 R_m=

2k

X

p=1

(a^p,11,1+ a^p,12,4+ a^p,13,3+ a^p,14,2)

=

k

X

p=1

(2p−1+2p+8k−2p+2+8k−2p+1)

+

2k

X

p=k+1

[4k−2(p−k)+2+4k−2(p−k)+1+4k+2(p−k)−1+4k+2(p−k)]

=

k

X

p=1

(16k + 2) +

2k

X

p=k+1

(16k + 2) = 4k(8k + 1).

從而, 我們已證, A 的左右對角線及各左右折斷對角線上 8k 個元素之和均為常數 4k(8k+

1), 故而滿足性質 2。

3) 由 (1) 式可知

4k

X

i=1

ai,i=

k

X

p=1

(2p − 1 + 4k + 2p + 8k − 2p + 2 + 4k − 2p + 1)

=

k

X

p=1

(16k + 2) = 2k(8k + 1),

4k

X

i=1

ai,8k+1−i=

k

X

p=1

(8k − 2p + 2 + 4k − 2p + 1 + 2p − 1 + 4k + 2p)

=

k

X

p=1

(16k + 2) = 2k(8k + 1).

而 2) 中已證,

8k

P

i=1

ai,i = 4k(8k + 1),

8k

P

i=1

ai,8k+1−i = 4k(8k + 1).

故 A 的各左右對角線的一半上的 4k 個元素之和均為 2k(8k + 1) 即 A 滿足性質 4。

4) 由 (1) 式顯然可知, 在任意 4 階子方陣 Ap,q (p, q = 1, 2, . . . , 2k), 中, 均有 a^p,q_r,1+ a^p,q_r,4 = 8k + 1, a^p,q_r,2+ a^p,q_r,3 = 8k + 1 (r = 1, 2, 3, 4).

因 A 中任意 m 行 (2 ≤ m ≤ 2k) 上左右對稱的每條橫向波形線 「 」 或 「 」 上的 8k 個元素, 恰由 2k 組 a^p,q_r,1, a^p,q_r,4 與 2k 組 a^p,q_r,2, a^p,q_r,3 構成。

因而, 這樣的 8k 個元素之和均為 4k(8k + 1), 即 A 滿足性質 5。

5) 由 (1) 式易見, 在任意 4 階子方陣, Ap,q 與 Ap+k,q (p = 1, 2, . . . , k; q = 1, 2, . . . , 2k) 中, 均有

a^p,q1,s + a^p,q4,s + a^p+k,q1,s + a^p+k,q4,s = 2(8k + 1), a^p,q2,s + a^p,q3,s + a^p+k,q2,s + a^p+k,q3,s = 2(8k + 1),

(s = 1, 2, 3, 4).

因 A 中任意 m 列 (2 ≤ m ≤ 2k) 中, 上下對稱的每條縱向波形線 或 上前 8k 個 元素恰由 k 组 a^p,q1,s, a^p,q4,s, a^p+k,q1,s , a^p+k,q4,s 與 k 組 a^p,q2,s, a^p,q3,s, a^p+k,q2,s , a^p+k,q3,s 構成。 因而這樣的 8k 個元素之和為 4k(8k + 1), 即 A 滿足性質 6。

從而, 已經證明 A 滿足 8k 階優化 F 型幻方的 6 條性質 (對於 A, 其中常數 C = 4k(8k + 1)), 由性質 1, 3, 5, 6 具有行列互換不變性, 因而 B = A^T, 亦具有性質 1, 3, 5, 6; 由於 A 經過轉置後各條左右對角線及折斷對角線僅僅排列順序發生了一些變化, 而元素集 合沒有改變, 因而 B 亦具有性質 2, 4。

至此, 我們證明了按上述方法構造的 8k 階方陣 A、 B 滿足定理 1 的所有條件。 因之, 由 它們作出的方陣 H = (h_i,j), h_i,j = 8k(a_i,j− 1) + b_i,j (i, j = 1, 2, . . . , 8k) 必為 8k 階優化 F 型幻方。

四、 結語

事實上, 按本文方法作出的 8k 階優化 F 型幻方還具有另外一些奇妙的性質, 例如:

將 8k 階優化 F 型幻方 H 分為 2k × 2k 個 4 階子方陣 Hp,q (p, q = 1, 2, . . . , 2k), 由這些 4 階子方陣組成的所有 4m 階 (m = 1, 2, . . . , 2k) 子方陣均為廣義全對稱幻方。 這樣 8 階優化 F 型幻方含有 4 個 4 階和 1 個 8 階全對稱幻方; 16 階優化 F 型幻方含有 16 個 4 階, 9 個 8 階, 4 個 12 階和 1 個 16 階全對稱幻方。 一般的, 8k 階優化 F 型幻方內含有 (2k − m + 1)² 個 4m (m = 1, 2, . . . , 2k) 階全對稱幻方。 並且, 其中那些被 8k 階方陣的水 準中軸線穿過, 同時上下行數正好相等的 8e (e = 1, 2, . . . , k) 階子方陣還是廣義 8e 階優化 F 型幻方。

8k 階優化 F 型幻方的每個 4 階子方陣都是幻和為 C/2k 的全對稱幻方, 因此, 任意地從 這 2k × 2k 個 4 階子方陣中取出 2k 個行、 列、 對角線或折斷對角線, 得到的 8k 個元素之和 均等於 8k 階幻和 C。

8k 階優化 F 型幻方內的任意 m 階 (m = 2, 4, . . . , 8k − 2) 子方陣中的 m² 個元素之 和都相等。

另外, 把 8 階優化 F 型幻方分成 16 個 2 × 2 小方陣, 依次為 H11, H¹², . . . , H⁴⁴。 把每 個 2 × 2 小方陣的平方和計算出來, 如下 16 個 2 × 2 小方陣平方和表:

H11 H12 H13 H14

5782 5654 5526 5654

H21 H22 H23 H24

5654 5782 5654 5526

H31 H32 H33 H34

5526 5654 5526 5398

H41 H42 H43 H44

5654 5526 5398 5526

把 H¹¹, H¹², H¹³, H¹⁴

H²¹, H²², H²³, H²⁴ 旋轉 180^◦, 分別與 H³¹, H³², H³³, H³⁴

H⁴¹, H⁴², H⁴³, H⁴⁴ 重合, 其兩個小方陣 8 元素 的平方和都等於 11180。

在上表中:

S4² = 5782 的兩組, 可以與 S4² = 5398 的兩組匹配使得 S8² = 11180。

S4² = 5526 的六組, 可以與 S4² = 5654 的六組匹配使得 S8² = 11180。

(S4² 表示 2 × 2 子方陣中 4 個數的平方和, S8² = 11180 是 8 階平方幻方的平方和。) 在這個優化F 型幻方裏, 出現了平方和相等的性質, 並且以2×2子方陣的形式出現, 頗為少見。

致謝: 本文從初稿到定稿經歷 4 次修改, 經過審稿老師的多次認真審理, 匡正謬誤, 並且不辭勞 苦親自給出簡捷、 清晰的證明 (A, B 陣相互正交), 要說感謝未免太過輕薄, 待日後有機會到貴 處拜謁學習!

附錄 :

紀念偉大科學家富蘭克林誕辰 310周年 (1706∼2016)

為紀念富蘭克林 (1706.1.17∼1790.4.17) 在幻方方面的巨大貢獻, 在富蘭克林誕辰 310 周年之際, 著名數學家李學數教授提議設計一個 「紀念富蘭克林誕辰 310 周年幻方」, 來紀念這 位身兼多職的著名科學家。

在下面這個幻方中, 具有 「優化富蘭克林幻方的全部性質」, (僅僅是使用的數字不是連續 自然數)。 幻方的幻和 S⁸ = 2106, 我們把 「310」 放在幻方的第 1 行, 第 1 列的位置上。 其它 性質如下:

生年:1706 年 = 第 1 行的第 2 列至第 8 列上 7 個元素之和;

= 第 1 列的第 2 行至第 8 行上 7 個元素之和;

= 對角線上 212 + 194 + 292 + 283 + 203 + 221 + 301 之和

= 折斷對角線上: 297 + 197 + 207 + 283 + 288 + 218 + 216 之和

月: 1 月, 隱含於 212 之中, 2 + 1 − 2 = 1。

日: 17 日, 隱含於 296 之中, 2 + 9 + 6 = 17。 或者, 隱含於 287 之中, 2 + 8 + 7 = 17。

逝世時間: 1790 年 = 1008 + 280 + 196 + 206,

(其中的 1008 等於任意半行、 或者半列、 或者 2 × 2 小方陣上的 4 個數之和等等。) 4 月, 隱含於 206 之中, 6 − 0 − 2 = 4. 17 日, 同上。

紀念富蘭克林誕辰 310 周年 8 階優化 F 型幻方 310 280 222 196 307 281 225 195

206 212 294 296 209 211 291 297 282 308 194 224 279 309 197 223 210 208 298 292 213 207 295 293 286 304 198 220 283 305 201 219 214 204 302 288 217 203 299 289 306 284 218 200 303 285 221 199 202 216 290 300 205 215 287 301

一友贊富蘭克林幻方

富蘭克林幻方奇, 富蘭克林幻方妙, 神奇奧妙驚天地。 特性最多無人超。

開創曲線和相等, 兩個幻方示後人, 和幻方中數第一。 豐功偉績滿天耀!

富蘭克林還有一個美好的願望, 他希望: 構造幻方不要過多繁雜的計算, 也不要調換元素 位置, 就像書寫自然數一樣快捷, 直接填寫出來。 這個美好的願望在 300 年後被實現了, 筆者 在 《數學傳播》1996 年 4 期發表了 「偶數階幻方的快速構作」, 其速度之快就像在空白紙上書 寫自然數一樣, 所以稱為 「直接書寫法」。 著名美國數學家李學數教授對 「直接書寫法」 的評語 是: 「方法之簡單, 速度之快捷, 其他任何方法歎之莫及!」。 從而, 完成了大發明家富蘭克林的夙 願, 他留下的兩個幻方難題都已經解決, 倘若富蘭克林先生有知, 定然 「含笑天堂」 矣!

參考文獻

1. Oystein Ore, Invitation to Number Theory, Yale University, 1967.

2. 梁培基, 顧同新。 多功能組合幻方。 數學傳播季刊, 14(2), 82-89, 1990。

3. 梁培基, 張航輔, 張俠輔。 幻方的一種構作方法。 雲南大學學報, 11(4), 310-319, 1989。

4. 梁培基, 顧同新。 平方幻方與雙重幻方的構造。 數學傳播季刊, 13(3), 65-69, 1989。

5. 潘承彪 譯。 《有趣的數論》。 北京大學出版社, 1985 年。

6. 李學數。《數學與數學家的故事》。 第 7 集, 上海科技出版社, 2015 年出版。

—

—