奇数阶始元臻美和幻方的构造方法 The Method of Constructing Odd-Order Shi Yuan Perfect Sum Magic Square

doi:10.12677/AAM.2021.108279

Advances in Applied Mathematics
Vol. 10 No. 08 ( 2021 ), Article ID: 44391 , 9 pages
10.12677/AAM.2021.108279

奇数阶始元臻美和幻方的构造方法

惠晶^*，刘兴祥^#

●How to Cite this Article

延安大学数学与计算机科学学院，陕西延安

收稿日期：2021年7月4日；录用日期：2021年7月23日；发布日期：2021年8月6日

摘要

目的：研究奇数阶始元臻美和幻方的构造方法。方法：通过构造特殊函数及固定向量与矩阵表示相结合的思路进行研究。结果与结论：给出了奇数阶始元臻美和幻方的构造方法并提出关于偶数阶臻美和幻方不存在的猜想。

关键词

矩阵，奇数阶始元臻美和幻方，臻美和幻方的存在性，性质

The Method of Constructing Odd-Order Shi Yuan Perfect Sum Magic Square

Jing Hui^*, Xingxiang Liu^#

College of Mathematics and Computer Science, Yan’an University, Yan’an Shaanxi

Received: Jul. 4^th, 2021; accepted: Jul. 23^rd, 2021; published: Aug. 6^th, 2021

ABSTRACT

Object: To study the construction methods of odd-order Shi Yuan Perfect sum magic square. Methods: The research is carried out by constructing special function and combining fixed vector and matrix representation. Result and Conclusion: The construction methods of odd-order Shi Yuan Perfect sum magic square are given and conjecture about the existence of even-order perfect sum magic square is put forward.

Keywords:Matrix, Shi Yuan Perfect Sum Magic Square, The Existence of Perfect Sum Magic Square, Property

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

幻方 [1] [2] 又称“河图”、“洛书”或“纵横图”，是使得m行m列矩阵每行每列以及对角线上的数字和都相等的方法。幻方作为一类特殊的整数矩阵，不仅可应用于大量的美术设计和西方建筑中，还可应用于封面包装，其重要性不言而喻，自然而然得到众多学者的青睐。而臻美和幻方作为一类特殊的整数矩阵及和幻方，其构造问题也就有必要进行深入研究。

吴硕辛老师在中国文化出版社出版的《 $m i (q)$ 语言的理论与应用》 [3] 专著中提出了臻美和幻方，并说

$(\begin{matrix} 0 & 13 & 21 & 9 & 17 \\ 6 & 19 & 2 & 10 & 23 \\ 12 & 20 & 8 & 16 & 4 \\ 18 & 1 & 14 & 22 & 5 \\ 24 & 7 & 15 & 3 & 11 \end{matrix})$

既是臻美和幻方又是完美和幻方，而

$(\begin{matrix} 0 & 7 & 9 & 14 \\ 13 & 10 & 4 & 3 \\ 6 & 1 & 15 & 8 \\ 11 & 12 & 2 & 5 \end{matrix})$

是完美和幻方但不是臻美和幻方。

虽然他对臻美和幻方按其理论给出了定义，但由于其理论较为复杂且繁琐，对于大多数读者来说难以读懂，故需要对该理论进行解释补充。在通过对他提出的定义进行研究后，本文结合现今和幻方的规范定义形式，给出了臻美和幻方的相关定义、性质及奇数阶始元臻美和幻方的构造方法 [4] [5] [6]。

注：为了叙述方便，现对文中的符号进行说明。

1) 下文中规定 $m (\mod m) = m$ 。

2) 本文中 $E_{i j} (m, m)$ 表示第i行第j列的元素是1而其余元素都是零的m阶矩阵。

3) $(i + k) (\mod m)$ 表示 $i + k$ 对m做除法所得的余数， $(k - i + 1) (m o d m)$ 表示 $k - i + 1$ 对m做除法所得的余数。

2. 预备知识

为给出臻美和幻方的定义及其构造定理，我们需要以下定义。

定义1 [7] 若矩阵， $S_{1} \in R$ ，(其中 $a_{i j} \in F$ $(i, j = 1, 2, \dots, m)$ )满足：

①；

②；

③；

④；

则称矩阵A为数域F上的m阶和幻方，并称S₁为m阶和幻方A的幻和。

定义2 [8] 若矩阵， $S_{2} \in R$ ，(其中 $a_{i j} \in F$ $(i, j = 1, 2, \dots, m)$ )满足：

①；

②；

③；

④；

⑤ 对于 $\forall k (k = 0, 1, 2, \dots, m - 1)$ ，有 $\sum_{i = 1}^{m} E_{1 i} (1, m) A E_{(k - i + 1) (m o d m), 1} (m, 1) = S_{2}$ ；

⑥ 对于 $\forall k (k = 0, 1, 2, \dots, m - 1)$ ，有 $\sum_{i = 1}^{m} E_{1 i} (1, m) A E_{(i + k) (m o d m), 1} (m, 1) = S_{2}$ ；

则称矩阵A为m阶完美和幻方，并称S₂为m阶完美和幻方A的幻和。

定义3 [9] 设矩阵 $A = {(a_{i j})}_{m \times m} \in S^{m \times m}$ ( $S = {1, 2, \dots, m}$ )， $\forall i, j, k \in {1, 2, \dots, m}$ ，当 $i \neq j$ 有 $a_{i k} \neq a_{j k}, a_{k i} \neq a_{k j}$ ，则称矩阵A为m阶拉丁方。

定义4 [10] 设矩阵 $A = {(a_{i j})}_{m \times m} \in S^{m \times m}$ ( $S = {1, 2, \dots, m}$ )，满足：

① 对 $\forall i, j, k \in {1, 2, \dots, m}$ ，当 $j \neq k$ 有 $a_{i j} \neq a_{i k}$ 恒成立；

② 对 $\forall i, j, k \in {1, 2, \dots, m}$ ，当 $j \neq k$ 有 $a_{i k} \neq a_{j k}$ 恒成立；

③ 对 $\forall i, j, k \in {1, 2, \dots, m}$ ，当 $j \neq k$ 有 $a_{i i} \neq a_{j j}$ 恒成立；

④ 对 $\forall i, j, k \in {1, 2, \dots, m}$ ，当 $j \neq k$ 有 $a_{m + 1 - i} \neq a_{m + 1 - j}$ 恒成立；

则称矩阵A是m阶拉丁幻方。

定义5 [8] 设矩阵 $A = {(a_{i j})}_{m \times m} \in S^{m \times m}$ ( $S = {1, 2, \dots, m}$ )，满足：

① 矩阵A的每一行，每一列都是 $1, 2, \dots, m$ 的一个全排列；

② ${a_{i, (i + k) (\mod m)}} (k = 0, 1, 2, \dots, m - 1)$ 都是 $1, 2, \dots, m$ 的一个全排列；

③ ${a_{i, (k - i + 1) (\mod m)}} (k = 0, 1, 2, \dots, m - 1)$ 都是 $1, 2, \dots, m$ 的一个全排列；

则称矩阵A是m阶完美拉丁幻方。

3. 主要结果

3.1. 定义

首先给出奇数阶始元完美和幻方的定义。

定义6 若奇数阶矩阵， $S_{3} \in R$ (其中 $a_{i j} \in {1, 2, \dots, m^{2}}$ $(i, j = 1, 2, \dots, m)$ )，满足：

①；

②；

③；

④；

⑤ 对于 $\forall k (k = 0, 1, 2, \dots, m - 1)$ ，有 $\sum_{i = 1}^{m} E_{1 i} (1, m) A E_{(k - i + 1) (m o d m), 1} (m, 1) = S_{3}$ ；

⑥ 对于 $\forall k (k = 0, 1, 2, \dots, m - 1)$ ，有 $\sum_{i = 1}^{m} E_{1 i} (1, m) A E_{(i + k) (m o d m), 1} (m, 1) = S_{3}$ ；

当 $i \neq k$ 或 $j \neq l$ 时有 $a_{i j} \neq a_{k l}$ $(i, j, k, l = 1, 2, \dots, m)$ ，则称矩阵A为m阶始元完美和幻方，并称S₃为m阶始元完美和幻方A的幻和。

由于奇数阶始元完美和幻方是特殊的奇数阶连元完美和幻方，所以就有如下定义。

定义7 若奇数阶矩阵， $S_{4} \in R$ (其中 $a_{i j} \in {a + 1, a + 2, \dots, a + m^{2}}$ $(a \in Z)$

$(i, j = 1, 2, \dots, m)$ )，满足以下条件：

①；

②；

③；

④；

⑤ 对于 $\forall k (k = 0, 1, 2, \dots, m - 1)$ ，有 $\sum_{i = 1}^{m} E_{1 i} (1, m) A E_{(k - i + 1) (m o d m), 1} (m, 1) = S_{4}$ ；

⑥ 对于 $\forall k (k = 0, 1, 2, \dots, m - 1)$ ，有 $\sum_{i = 1}^{m} E_{1 i} (1, m) A E_{(i + k) (m o d m), 1} (m, 1) = S_{4}$ ；

当 $i \neq k$ 或 $j \neq l$ 时有 $a_{i j} \neq a_{k l}$ ， $(i, j, k, l = 1, 2, \dots, m)$ ，则称矩阵A为m阶连元完美和幻方，并称S₄为m阶连元完美和幻方A的幻和。

又由于奇数阶连元完美和幻方是特殊的奇数阶异元完美和幻方，所以也就有如下定义。

定义8 若奇数阶矩阵， $S_{5} \in R$ ，(其中 $a_{i j} \in {a_{1}, a_{2}, \dots, a_{n^{2}}}$ $(i, j = 1, 2, \dots, m)$ )且 $\forall i, j \in {1, 2, \dots, m^{2}}$ 当 $i \neq j$ 时 $a_{i} \neq a_{j}$ ，满足：

①；

②；

③；

④；

⑤ 对于 $\forall k (k = 0, 1, 2, \dots, m - 1)$ ，有 $\sum_{i = 1}^{m} E_{1 i} (1, m) A E_{(k - i + 1) (m o d m), 1} (m, 1) = S_{5}$ ；

⑥ 对于 $\forall k (k = 0, 1, 2, \dots, m - 1)$ ，有 $\sum_{i = 1}^{m} E_{1 i} (1, m) A E_{(i + k) (m o d m), 1} (m, 1) = S_{5}$ ；

当 $i \neq k$ 或 $j \neq l$ 时有 $a_{i j} \neq a_{k l}$ ， $(i, j, k, l = 1, 2, \dots, m)$ ，则称矩阵A为m阶异元完美和幻方，并称S₅为m阶异元完美和幻方A的幻和。

接下来，我们给出奇数阶始元臻美和幻方的定义。

定义9 若奇数阶矩阵， $S_{6} \in R$ ，(其中 $a_{i j} \in {1, 2, \dots, m^{2}}$ $(i, j = 1, 2, \dots, m)$ )，满足：

① $a_{11}, a_{12}, \dots, a_{1 m}, a_{21}, \dots, a_{2 m}, \dots, a_{m 1}, \dots, a_{m m}$ 是 $1, 2, \dots, m^{2}$ 的一个全排列；

②；

③ $\sum_{i = 1}^{m} a_{i, (i + k) (\mod m)} = \sum_{i = 1}^{m} a_{i, (k - i + 1) (\mod m)} = S_{6}$ ，其中 $(k = 0, 1, 2, \dots, m - 1)$ ；

④ $A = B + m C$ (B、C都是m阶完美拉丁幻方)；

则称A为m阶始元臻美和幻方，并称S₆为m阶始元臻美和幻方A的幻和。

定义10 若奇数阶矩阵， $S_{7} \in R$ (其中 $a_{i j} \in {a + 1, a + 2, \dots, a + m^{2}}$ $(a \in Z)$ $(i, j = 1, 2, \dots, m)$ )，满足：

① $a_{11}, a_{12}, \dots, a_{1 m}, a_{21}, \dots, a_{2 m}, \dots, a_{m 1}, \dots, a_{m m}$ 是 $a + 1, a + 2, \dots, a + m^{2}$ 的一个全排列；

②；

③ $\sum_{i = 1}^{m} a_{i, (i + k) (\mod m)} = \sum_{i = 1}^{m} a_{i, (k - i + 1) (\mod m)} = S_{7}$ ，其中 $(k = 0, 1, 2, \dots, m - 1)$ ；

④ $A = B + m C$ (B、C都是m阶完美拉丁幻方)；

则称A为m阶连元臻美和幻方，并称S₇为m阶连元臻美和幻方A的幻和。

定义11 若奇数阶矩阵， $S_{8} \in R$ (其中 $a_{i j} \in {a_{1}, a_{2}, \dots, a_{n^{2}}}$ $(i, j = 1, 2, \dots, m)$ )且 $\forall i, j \in {1, 2, \dots, m^{2}}$ 当 $i \neq j$ 时 $a_{i} \neq a_{j}$ ，满足：

① $a_{11}, a_{12}, \dots, a_{1 m}, a_{21}, \dots, a_{2 m}, \dots, a_{m 1}, \dots, a_{m m}$ 是 $a_{1}, a_{2}, \dots, a_{n^{2}}$ 的一个全排列；

②；

③ $\sum_{i = 1}^{m} a_{i, (i + k) (\mod m)} = \sum_{i = 1}^{m} a_{i, (k - i + 1) (\mod m)} = S_{8}$ ，其中 $(k = 0, 1, 2, \dots, m - 1)$ ；

④ $A = B + m C$ (B，C都是m阶完美拉丁幻方)；

则称A为m阶异元臻美和幻方，并称S₈为m阶异元臻美和幻方A的幻和。

为了证明奇数阶臻美和幻方的构造定理，下面我们给出具有和幻性矩阵的定义。

定义12 设奇数阶矩阵 $A = {(a_{i j})}_{m \times m}$ ，如果，则称A是具有和幻性的m阶矩阵，S₉为m阶矩阵A的幻和。

3.2. 性质

根据臻美和幻方的定义可知，有以下结论成立。

性质1. 设奇数阶矩阵 $A = {(a_{i j})}_{n \times n}$ 是异元臻美和幻方，并且规定 $n (\mod n) = n$ ，则 $\forall k = 0, 1, 2, \dots, n - 1$ 有 ${(a_{(i + k) (\mod n), j})}_{n \times n}$ 及 ${(a_{i, (j + k) (\mod n)})}_{n \times n}$ 也是异元臻美和幻方。

由于连元臻美和幻方是特殊的异元臻美和幻方，始元臻美和幻方是特殊的连元臻美和幻方，所以有以下推论成立。

推论1. 设奇数阶矩阵 $A = {(a_{i j})}_{n \times n}$ 是连元臻美和幻方，并且规定 $n (\mod n) = n$ ，则 $\forall k = 0, 1, 2, \dots, n - 1$ 有 ${(a_{(i + k) (\mod n), j})}_{n \times n}$ 及 ${(a_{i, (j + k) (\mod n)})}_{n \times n}$ 也是连元臻美和幻方。

推论2. 设奇数阶矩阵 $A = {(a_{i j})}_{n \times n}$ 是始元臻美和幻方，并且规定 $n (\mod n) = n$ ，则 $\forall k = 0, 1, 2, \dots, n - 1$ 有 ${(a_{(i + k) (\mod n), j})}_{n \times n}$ 及 ${(a_{i, (j + k) (\mod n)})}_{n \times n}$ 也是始元臻美和幻方。

性质2. 如果奇数阶矩阵 $A = {(a_{i j})}_{n \times n}$ 是异元臻美和幻方，且幻和为 $Δ$ ( $Δ$ is u-digit number，即 $Δ$ 是一个u位数)，规定 $n (\mod n) = n$ ，将 $a_{i j}$ 变为 $c_{i j} a_{i j} = d_{i j}$ (其中 $a_{i j}$ is $k_{i j}$ -digit number (即 $a_{i j}$ 是一个 $k_{i j}$ 位数)， $c_{i j}$ is $u - k_{i j}$ -digit number (即 $c_{i j}$ 是一个 $u - k_{i j}$ 位数) $(i, j = 0, 1, 2, \dots, n)$ )，则有

$\begin{matrix} \sum_{i = 1}^{n} d_{i 1} d_{i 2} \dots d_{i n} = \sum_{j = 1}^{n} d_{1 j} d_{2 j} \dots d_{n j} \\ = \sum_{j = 1}^{n} d_{1, (1 + j) (m o d n)} d_{2, (2 + j) (m o d n)} \dots d_{n, (n + j) (m o d n)} \\ = \sum_{j = 1}^{n} d_{n, (1 + j) (m o d n)} d_{n - 1, (2 + j) (m o d n)} \dots d_{1, (n + j) (m o d n)} \end{matrix}$

3.3. 奇数阶臻美和幻方的构造方法

定理1 设 $α = {(0, 1, 2, \dots, 2 n)}^{T}$ ， $β = α + {(1, 1, \dots, 1)}^{T}$ ，对 $\forall x = {(x_{1}, x_{2}, \dots, x_{2 n + 1})}^{T}$ 规定：

$f (x) = (x_{n + 1}, x_{n + 2}, \dots, x_{2 n + 1}, x_{1}, \dots, x_{n})$

$g (x) = (x_{n + 2}, x_{n +3}, \dots, x_{2 n + 1}, x_{1}, \dots, x_{n + 1})$ ;

利用 $f (x)$ 及 $g (x)$ 构造 $2 n + 1$ 阶矩阵B及C如下：

$B = (β, f (β), f^{2} (β), f^{3} (β), \dots, f^{2 n} ( β ))$

$C = (α, g (α), g^{2} (α), g^{3} (α), \dots, g^{2 n} (α))$ ;

如果 $2 n + 1$ 阶矩阵B，C为 $2 n + 1$ 阶完美拉丁幻方，则称矩阵 $A = B + (2 n + 1) C$ 是 $2 n + 1$ 阶始元臻美和幻方。

证明：要证明矩阵A是 $2 n + 1$ 阶始元臻美和幻方，就需要完成以下四步：

① A是具有和幻性的矩阵；

② A中的元素两两互不相等；

③ A中的元素取自集合 ${1, 2, \dots, {(2 n + 1)}^{2}}$ ；

④ B，C是完美拉丁幻方。

第一步：

由于B，C是完美拉丁幻方，则B，C是具有和幻性的矩阵。又由于 $A = B + (2 n + 1) C$ 且具有和幻性矩阵的线性组合仍然是具有和幻性的矩阵。则矩阵A就是具有和幻性的矩阵。

第二步：

由于B是关于 $1, 2, \dots, 2 n + 1$ 的拉丁幻方，而C是关于 $0, 1, \dots, 2 n$ 的拉丁幻方，且 $A = B + (2 n + 1) C$ ，则B与C正交。则 $a_{i j} = b_{i j} + (2 n + 1) c_{i j}$ ，对 $\forall i, j \in {1, 2, \dots, 2 n + 1}$ 两两互不相等，则 $a_{i j}$ $(i, j = 1, 2, \dots, 2 n + 1)$ 两两互不相等，也即A是 $2 n + 1$ 阶异元臻美和幻方。

第三步：

由于 $1 \leq b_{i j} \leq 2 n + 1$ ， $1 \leq c_{i j} \leq 2 n$ ，则 $1 \leq b_{i j} + (2 n + 1) c_{i j} \leq {(2 n + 1)}^{2}$ ，则 $a_{i j} \in {1, 2, \dots, {(2 n + 1)}^{2}}$ $(i, j = 1, 2, \dots, 2 n + 1)$ ，也即A是 $2 n + 1$ 阶始元臻美和幻方。

第四步：

由已知可知B，C是完美拉丁幻方。

综合以上可知，矩阵A即为 $2 n + 1$ 阶始元臻美和幻方。

注：

奇数阶始元臻美和幻方的构造方法也是奇数阶始元完美和幻方的构造方法，同样也是奇数阶始元和幻方的构造方法。但奇数阶始元和幻方的构造方法不一定能构造出奇数阶始元完美和幻方，同时奇数阶始元完美和幻方的构造方法不一定能构造出奇数阶始元臻美和幻方。

(由于奇数阶始元臻美和幻方是特殊的奇数阶始元完美和幻方，也是特殊的奇数阶始元和幻方，故这两类和幻方是存在的且对应的奇数阶定义是有意义的。)

推论3 设 $α = {(0 + a, 1 + a, 2 + a, \dots, 2 n + a)}^{T}$ $(a \in Z)$ ， $β = α + {(1, 1, \dots, 1)}^{T}$ ，对 $\forall x = {(x_{1}, x_{2}, \dots, x_{2 n + 1})}^{T}$ 规定

$f (x) = (x_{n + 1}, x_{n + 2}, \dots, x_{2 n + 1}, x_{1}, \dots, x_{n})$

$g (x) = (x_{n + 2}, x_{n +3}, \dots, x_{2 n + 1}, x_{1}, \dots, x_{n + 1})$ ;

利用 $f (x)$ 及 $g (x)$ 构造 $2 n + 1$ 阶矩阵B及C如下：

$B = (β, f (β), f^{2} (β), f^{3} (β), \dots, f^{2 n} ( β ))$

$C = (α, g (α), g^{2} (α), g^{3} (α), \dots, g^{2 n} (α))$ ;

如果 $2 n + 1$ 阶矩阵B、C为 $2 n + 1$ 阶完美拉丁幻方，则称矩阵 $A = B + (2 n + 1) C$ 是 $2 n + 1$ 阶连元臻美和幻方。

由于定理及推论中涉及到的这 $2 n + 1$ 个元素都是两两互不相等的，所以始元臻美和幻方的构造方法也是异元臻美和幻方的构造方法。

接下来通过例子验证定义9及定理1的正确性。

例1. 给出一个11阶始元臻美和幻方。

解：由奇数阶臻美和幻方的构造定理1知，n阶完美拉丁幻方B，C分别为：

$B = {(\begin{matrix} 1 & 6 & 11 & 5 & 10 & 4 & 9 & 3 & 8 & 2 & 7 \\ 2 & 7 & 1 & 6 & 11 & 5 & 10 & 4 & 9 & 3 & 8 \\ 3 & 8 & 2 & 7 & 1 & 6 & 11 & 5 & 10 & 4 & 9 \\ 4 & 9 & 3 & 8 & 2 & 7 & 1 & 6 & 11 & 5 & 10 \\ 5 & 10 & 4 & 9 & 3 & 8 & 2 & 7 & 1 & 6 & 11 \\ 6 & 11 & 5 & 10 & 4 & 9 & 3 & 8 & 2 & 7 & 1 \\ 7 & 1 & 6 & 11 & 5 & 10 & 4 & 9 & 3 & 8 & 2 \\ 8 & 2 & 7 & 1 & 6 & 11 & 5 & 10 & 4 & 9 & 3 \\ 9 & 3 & 8 & 2 & 7 & 1 & 6 & 11 & 5 & 10 & 4 \\ 10 & 4 & 9 & 3 & 8 & 2 & 7 & 1 & 6 & 11 & 5 \\ 11 & 5 & 10 & 4 & 9 & 3 & 8 & 2 & 7 & 1 & 6 \end{matrix})}_{11 \times 11}$

$C = {(\begin{matrix} 0 & 6 & 1 & 7 & 2 & 8 & 3 & 9 & 4 & 10 & 5 \\ 1 & 7 & 2 & 8 & 3 & 9 & 4 & 10 & 5 & 0 & 6 \\ 2 & 8 & 3 & 9 & 4 & 10 & 5 & 0 & 6 & 1 & 7 \\ 3 & 9 & 4 & 10 & 5 & 0 & 6 & 1 & 7 & 2 & 8 \\ 4 & 10 & 5 & 0 & 6 & 1 & 7 & 2 & 8 & 3 & 9 \\ 5 & 0 & 6 & 1 & 7 & 2 & 8 & 3 & 9 & 4 & 10 \\ 6 & 1 & 7 & 2 & 8 & 3 & 9 & 4 & 10 & 5 & 0 \\ 7 & 2 & 8 & 3 & 9 & 4 & 10 & 5 & 0 & 6 & 1 \\ 8 & 3 & 9 & 4 & 10 & 5 & 0 & 6 & 1 & 7 & 2 \\ 9 & 4 & 10 & 5 & 0 & 6 & 1 & 7 & 2 & 8 & 3 \\ 10 & 5 & 0 & 6 & 1 & 7 & 2 & 8 & 3 & 9 & 4 \end{matrix})}_{11 \times 11}$

则矩阵

$A = {(\begin{matrix} 1 & 72 & 22 & 82 & 32 & 92 & 42 & 102 & 52 & 112 & 62 \\ 13 & 84 & 23 & 94 & 44 & 104 & 54 & 114 & 64 & 3 & 74 \\ 25 & 96 & 35 & 106 & 45 & 116 & 66 & 5 & 76 & 15 & 86 \\ 37 & 108 & 47 & 118 & 57 & 7 & 67 & 17 & 88 & 27 & 98 \\ 49 & 120 & 59 & 9 & 69 & 19 & 79 & 29 & 89 & 39 & 110 \\ 61 & 11 & 71 & 21 & 81 & 31 & 91 & 41 & 101 & 51 & 111 \\ 73 & 12 & 83 & 33 & 93 & 43 & 103 & 53 & 113 & 63 & 2 \\ 85 & 24 & 95 & 34 & 105 & 55 & 115 & 65 & 4 & 75 & 14 \\ 97 & 36 & 107 & 46 & 117 & 56 & 6 & 77 & 16 & 87 & 26 \\ 109 & 48 & 119 & 58 & 8 & 68 & 18 & 78 & 28 & 99 & 38 \\ 121 & 60 & 10 & 70 & 20 & 80 & 30 & 90 & 40 & 100 & 50 \end{matrix})}_{11 \times 11}$

对于矩阵A中的任意元素 $a_{i j}$ ，都有 $a_{i j} \in {1, 2, \dots, 11^{2}}$ $(i, j = 1, 2, \dots, 11)$ )，满足：

① $a_{11}, a_{12}, \dots, a_{1 ， 11}, a_{21}, \dots, a_{2 ， 11}, \dots, a_{11 ， 1}, \dots, a_{11 ， 11}$ 是 $1, 2, \dots, 11^{2}$ 的一个全排列；

②；

③ $\sum_{i = 1}^{11} a_{i, (i + k) (\mod 11)} = \sum_{i = 1}^{11} a_{i, (k - i + 1) (\mod 11)} = S$ ，其中 $(k = 0, 1, 2, \dots, 11 - 1)$ ；

④ $A = B + 11 C$ (B、C都是m阶完美拉丁幻方)；

因此，矩阵 $A = {(a_{i j})}_{11 \times 11}$ 满足定义9的条件。

所以，矩阵A为11阶始元臻美和幻方。

4. 偶数阶臻美和幻方的存在性

由奇数阶臻美和幻方的构造定理知，奇数阶臻美和幻方是存在的，但我们通过计算机程序求解得到 $n (n \geq 50)$ 阶的所有和幻方，对其进行 $A = B + n C$ 形式的分解，分解后发现B和C只有一个是拉丁幻方，另一个不是拉丁幻方。又由任意的一个始元和幻方都可以分解为一个行拉丁方和一个列拉丁方的形式 [3]，所以我们猜想：偶数阶臻美和幻方不存在。

关于此猜想的证明，目前来说，我们没有找到一个合适的证明方法，只有通过列举法把一些情况进行否定，具体证明方式有待于日后再做进一步的研究。

5. 小结

本文利用了构造特殊函数及固定向量与矩阵表示相结合的思路，研究了奇数阶始元臻美和幻方的构造方法，并通过具体例子验证了定义9及定理1的正确性，最后提出关于偶数阶臻美和幻方不存在的猜想。

文章引用

惠晶,刘兴祥. 奇数阶始元臻美和幻方的构造方法
The Method of Constructing Odd-Order Shi Yuan Perfect Sum Magic Square[J]. 应用数学进展, 2021, 10(08): 2692-2700. https://doi.org/10.12677/AAM.2021.108279

参考文献

1. Denes, J. and Keedwell, A.D. (1988) A Conjecture Concerning Magic Squares. The American Mathematical Monthly, 95, 845-849. https://doi.org/10.1080/00029890.1988.11972098

2. Peyman, F., Toussi, C.A., Trump, W., Haddadnia, J. and Matta, C.F. (2021) Striking Patterns in Natural Magic Squares’ Associated Electrostatic Potentials: Matrices of the 4th and 5th Order. Discrete Mathematics, 344, 112229. https://doi.org/10.1016/j.disc.2020.112229

3. 吴硕辛. mi(q)语言的理论与应用[M]. 北京: 中国文化出版社, 2016.

4. Johnson, C.R. (2018) A Matrix Theoretic Construction of Magic Squares. The American Mathematical Monthly, 79, 1004-1006. https://doi.org/10.1080/00029890.1972.11993170

5. Xu, D., Mao, Z., Chen, B. and Huang, P. (2011) New Recursive Construction of Magic Squares Using Kronecker Compositional Operations and Its Application in Engineering Computation. Systems Engineering Procedia, 2, 331-337. https://doi.org/10.1016/j.sepro.2011.10.046

6. Baksalary, O.M., Trenkler, D. and Trenkler, G. (2020) On Most Perfect Magic Squares of Order Four. Linear and Multilinear Algebra, 68, 1411-1423. ttps://doi.org/10.1080/03081087.2018.1545005

7. 董朦朦, 刘兴祥, 田雨禾. 幻方可以分解为两个正交拉丁方的线性组合[J]. 重庆理工大学学报(自然科学), 2018, 32(8): 181-188.

8. 张婧, 刘兴祥, 刘娟娟. 完美积幻方的定义及其构造[J]. 延安大学学报(自然科学版), 2020, 39(4): 90-93.

9. 刘兴祥, 曹燕飞. 4k阶拉丁幻方的构造方法[J]. 河南科学, 2016, 34(11): 1777-1780.

10. 曹燕飞, 刘兴祥. 3k阶拉丁幻方的构造方法[J]. 延安大学学报(自然科学版), 2016, 35(4): 24-26.

NOTES

^*惠晶(1998-)，女，陕西榆林人，延安大学本科生。

^#通讯作者。刘兴祥(1964-)，男，陕西合阳人，副教授，硕士生导师，研究方向：矩阵理论。

期刊菜单