单偶数阶和幻方的一种构造方法
A Method of Constructing a Sum Magic Square of Single Even Order
DOI: 10.12677/AAM.2022.1112940, PDF, HTML, XML,  被引量 下载: 195  浏览: 295  科研立项经费支持
作者: 张宇婷*, 刘兴祥#, 朱雅妮#:延安大学数学与计算机科学学院,陕西 延安
关键词: 和幻方分块矩阵奇数阶和幻方单偶数阶和幻方Sum Magic Square Chunked Matrix Odd-Order Sum Magic Square Sum Magic Square of Single Even Order
摘要: 根据和幻方的定义及其性质,并通过观察研究,结合已知奇数阶和幻方及以分块矩阵为工具,给出单偶数阶和幻方的分块矩阵构造法及证明。
Abstract: According to the definitions and properties of sum magic square, and adopted observational study, combined with the odd-order sum magic square known to us and used chunked matrix as a tool, the method of chunked matrix constructed about sum magic square of single even order is given and proved.
文章引用:张宇婷, 刘兴祥, 朱雅妮. 单偶数阶和幻方的一种构造方法[J]. 应用数学进展, 2022, 11(12): 8917-8922. https://doi.org/10.12677/AAM.2022.1112940

1. 引言

从古至今,幻方作为一种特殊的矩阵深受广大研究者的喜爱,引发了许多人的研究。经过一代代数学家与数学爱好者的共同努力,关于幻方的研究成果已经相当丰富。文献 [1] [2] [3] 是关于和幻方的部分研究成果。但在深入研究的过程中,发现在和幻方构造方法中仍存在着一些问题:在已有和幻方的构造方法中,利用特殊矩阵构造和幻方的方法并不多,大部分都是具体阶数和幻方的构造方法。由此可以了解到给出一般阶数和幻方的构造方法较为困难。为了解决这类问题,可以结合分块矩阵分别找出奇数阶、单偶数阶、双偶数阶和幻方构造的方法。文中将给出单偶数阶和幻方的一种构造方法,即利用分块矩阵构造单偶数阶和幻方。

2. 预备知识

定义1 [4] [5] 设矩阵 A = ( a i j ) n × n Z n × n , n N * ,若矩阵A满足以下条件:

i { 1 , 2 , , n } E 1 j ( 1, n ) ( j = 1 n ( A E i j ( n , n ) ) ) i = 1 n E i 1 ( n , 1 ) = S r

j { 1 , 2 , , n } j = 1 n E 1 j ( 1 , n ) ( i = 1 n ( A E i j ( n , n ) ) ) E i 1 ( n , 1 ) = S c

S r = S c = j = 1 n E 1 j ( 1 , n ) ( i = 1 n ( A E i i ) ) i = 1 n E i 1 ( n , 1 ) = j = 1 n E 1 j ( 1 , n ) ( i = 1 n A E i , n + 1 i ) i = 1 n E i 1 ( n , 1 ) = S

则称矩阵 A = ( a i j ) n × n 为n阶和幻方,并称S为n阶和幻方A的幻和。

定义2 [4] [5] 设矩阵 A = ( a i j ) n × n Z n × n , n N * ,若矩阵A满足以下条件:

i { 1 , 2 , , n } E 1 j ( 1, n ) ( j = 1 n ( A E i j ( n , n ) ) ) i = 1 n E i 1 ( n , 1 ) = S r

j { 1 , 2 , , n } j = 1 n E 1 j ( 1 , n ) ( i = 1 n ( A E i j ( n , n ) ) ) E i 1 ( n , 1 ) = S c

S r = S c = j = 1 n E 1 j ( 1 , n ) ( i = 1 n ( A E i i ) ) i = 1 n E i 1 ( n , 1 ) = j = 1 n E 1 j ( 1 , n ) ( i = 1 n A E i , n + 1 i ) i = 1 n E i 1 ( n , 1 ) = S

④ 当 i k j l 时, i , j , k , l { 1 , 2 , , n } 均有 a i j a k l

则称矩阵 A = ( a i j ) n × n 为n阶异元和幻方,并称S为n阶异元和幻方A的幻和。

定义3 [6] [7] 设矩阵 A = ( a i j ) n × n a i j { 1 , 2 , , n 2 } ( i , j = 1 , 2 , , n ) , n N * ,若矩阵A满足以下条件:

i { 1 , 2 , , n } E 1 j ( 1, n ) ( j = 1 n ( A E i j ( n , n ) ) ) i = 1 n E i 1 ( n , 1 ) = S r

j { 1 , 2 , , n } j = 1 n E 1 j ( 1 , n ) ( i = 1 n ( A E i j ( n , n ) ) ) E i 1 ( n , 1 ) = S c

S r = S c = j = 1 n E 1 j ( 1 , n ) ( i = 1 n ( A E i i ) ) i = 1 n E i 1 ( n , 1 ) = j = 1 n E 1 j ( 1 , n ) ( i = 1 n A E i , n + 1 i ) i = 1 n E i 1 ( n , 1 ) = S

④ 当 i k j l 时, i , j , k , l { 1 , 2 , , n } 均有 a i j a k l

则称矩阵 A = ( a i j ) n × n 为n阶始元和幻方,并称S为n阶始元和幻方A的幻和。

定义4 [6] [7] 设矩阵 A = ( a i j ) n × n a i j { a + 1 , a + 2 , a + n 2 } ( i , j = 1 , 2 , , n 2 ) , n N * ,若矩阵A满足以下条件:

i { 1 , 2 , , n } E 1 j ( 1, n ) ( j = 1 n ( A E i j ( n , n ) ) ) i = 1 n E i 1 ( n , 1 ) = S r

j { 1 , 2 , , n } j = 1 n E 1 j ( 1 , n ) ( i = 1 n ( A E i j ( n , n ) ) ) E i 1 ( n , 1 ) = S c

S r = S c = j = 1 n E 1 j ( 1 , n ) ( i = 1 n ( A E i i ) ) i = 1 n E i 1 ( n , 1 ) = j = 1 n E 1 j ( 1 , n ) ( i = 1 n A E i , n + 1 i ) i = 1 n E i 1 ( n , 1 ) = S

④ 当 i k j l 时, i , j , k , l { 1 , 2 , , n } 均有 a i j a k l

则称矩阵 A = ( a i j ) n × n 为n阶连元和幻方,并称S为n阶连元和幻方A的幻和。

定义5 [8] 把一个 m × n 矩阵A,在行的方向分成 s p 块,在列的方向分成 t q 块,称为A的 s p × t q 分块矩

阵,记作 A = [ A s i × t j ] s p × t q i = 1 p s i = m j = 1 q t j = n 其中 A s i × t j ( i = 1 , 2 , , p ; j = 1 , 2 , , q ) ,称为A的子块,它

们是各种类型的小矩阵。

定义6 [9] 设AB都是 m × n 矩阵,并且对AB用同样的方法进行分块: A = [ A s 1 × t 1 A s 1 × t 2 A s 1 × t q A s 2 × t 1 A s 2 × t 2 A s 2 × t q A s p × t 1 A s p × t 2 A s p × t q ] B = [ B s 1 × t 1 B s 1 × t 2 B s 1 × t q B s 2 × t 1 B s 2 × t 2 B s 2 × t q B s p × t 1 B s p × t 2 B s p × t q ]

i = 1 p s i = m j = 1 q t j = n ,其中 A s i × t j B s i × t j ( i = 1 , 2 , , p ; j = 1 , 2 , , q ) 都是 s i × t j 矩阵,即 A s i × t j B s i × t j 是同型矩阵,那么

A + B = [ A s 1 × t 1 + B s 1 × t 1 A s 1 × t 2 + B s 1 × t 2 A s 1 × t q + B s 1 × t q A s 2 × t 1 + B s 2 × t 1 A s 2 × t 2 + B s 2 × t 2 A s 2 × t q + B s 2 × t q A s p × t 1 + B s p × t 1 A s p × t 2 + B s p × t 2 A s p × t q + B s p × t q ]

定义7 [10] 设A m × n 矩阵,把A进行分块:

A = [ A s 1 × t 1 A s 1 × t 2 A s 1 × t q A s 2 × t 1 A s 2 × t 2 A s 2 × t q A s p × t 1 A s p × t 2 A s p × t q ]

其中 A s i × t j ( i = 1 , 2 , , p ; t = 1 , 2 , , q ; i = 1 p s i = m ; j = 1 q t j = n ) a为任意数,则

a A = [ a A s 1 × t 1 a A s 1 × t 2 a A s 1 × t q a A s 2 × t 1 a A s 2 × t 2 a A s 2 × t q a A s p × t 1 a A s p × t 2 a A s p × t q ]

3. 主要结果

定义6 设矩阵 A = ( a i j ) n × n ( N * ) n × n ,若 i = 1 n a i j = j = 1 n a i j = i = 1 n a i i = i = 1 n a i , n + 1 i = S ,则称矩阵A是具有和幻性的n阶和幻阵,Sn阶和幻阵A的幻和。

性质1 若矩阵 A = ( a i j ) n × n ( N * ) n × n ,则 a N * ,存在aA也是具有和幻性的n阶和幻阵。

性质2 若矩阵AB为具有和幻性的n阶和矩阵,则 a , b N * ,存在 a A + b B 也是具有和幻性的n阶和幻阵。

定理1 设矩阵T是一个2n+1阶始元和幻方,

矩阵 A 2 ( 2 n + 1 ) × 2 ( 2 n + 1 ) = ( B ( 2 n + 1 ) × 2 F ( 2 n + 1 ) × 2 C ( 2 n + 1 ) × ( n 1 ) G ( 2 n + 1 ) × ( n 1 ) D ( 2 n + 1 ) × n H ( 2 n + 1 ) × n J ( 2 n + 1 ) × ( n 1 ) M ( 2 n + 1 ) × ( n 1 ) K ( 2 n + 1 ) × ( n + 2 ) N ( 2 n + 1 ) × ( n + 2 ) )

其中

B = 3 × ( b i j ) ( 2 n + 1 ) × 2 × E 2 × 2 ( b i j = { 1 , ( i + j ) ( mod 2 ) = 0 0 , ( i + j ) ( mod 2 ) = 1 )

C = 3 × [ 1 ] ( 2 n + 1 ) × ( n 1 )

D = 0 × [ 1 ] ( 2 n + 1 ) × n

F = 3 × ( f i j ) ( 2 n + 1 ) × 2 × E 2 × 2 ( f i j = { 1 , ( i + j ) ( mod 2 ) = 1 0 , ( i + j ) ( mod 2 ) = 0 )

G = 0 × [ 1 ] ( 2 n + 1 ) × ( n 1 )

H = 3 × [ 1 ] ( 2 n + 1 ) × n

J = [ 1 ] ( 2 n + 1 ) × ( n 1 )

K = 2 × [ 1 ] ( 2 n + 1 ) × ( n + 2 )

M = 2 × [ 1 ] ( 2 n + 1 ) × ( n 1 )

N = [ 1 ] ( 2 n + 1 ) × ( n + 2 )

则矩阵 P = ( 2 n + 1 ) 2 × A + ( T T T T ) 为4n + 2阶单偶数阶和幻方。

证明:要证明P为和幻方,则需证明矩阵P的每行每列的和均相等,由于单偶数阶和幻方是利用奇数阶和幻方并通过分块矩阵构造而来。奇数阶和幻方每行、每列、主对角线及副对角线的幻和是相同的,所以只要证明矩阵A的每行、每列、主对角线及副对角线的和均为 3 ( 2 n + 1 ) 即可。

1) 由于矩阵A的构造形式,在计算它的行和时可以将矩阵A分块为两部分来进行计算,即前 2 n + 1 行与后 2 n + 1 行,并分别计算这两部分每一行的行和。

2 n + 1 行:矩阵B的每行的和为3,矩阵C的每行的和为 3 ( n 1 ) ,矩阵D的每行的和为0,矩阵E的每行的和为 ( n 1 ) ,矩阵F的每行的和为 2 ( n + 2 ) ,所以 S r ( i ) = 3 + 3 × ( n 1 ) + 0 + ( n 1 ) + 2 ( n + 2 ) = 6 n + 3 ( i = 1 , 2 , , 2 n + 1 )

2 n + 1 行:矩阵F的每行的和为3,矩阵G的每行的和为0,矩阵H的每行的和为3n,矩阵M的每行的和为 2 ( n 1 ) ,矩阵N的每行的和为(n+2),所以 S r ( i ) = 3 + 0 × ( n 1 ) + 3 n + 2 ( n 1 ) + ( n + 2 ) = 6 n + 3 ( i = 2 n + 2 , 2 n + 3 , , 4 n + 2 )

因此在矩阵A中, S r ( i ) = 6 n + 3 ( i = 1 , 2 , , 4 n + 2 )

2) 由于矩阵A的构造形式,在计算它的列和时,可以将矩阵A分为五部分来进行计算,在每一部分中两个分块矩阵是同型的,则可分别计算每一部分的每一列的列和,即:

① 第1列及第2列: B + F = 3 × ( c i j ) ( 4 n + 2 ) , 2 × E 2 × 2 , ( c i j = { 1 , ( i + j ) ( mod 2 ) = 0 0 , ( i + j ) ( mod 2 ) = 1 ) ,所以 S c ( j ) = 3 ( 2 n + 1 ) + 0 × ( 2 n + 1 ) = 6 n + 3 ( j = 1 , 2 )

② 第3列至第 ( n 1 ) 列: C + G = [ 3 ( 2 n + 1 ) × ( n 1 ) 0 ( 2 n + 1 ) × ( n 1 ) ] ,所以 S c ( j ) = 3 ( 2 n + 1 ) + 0 × ( 2 n + 1 ) = 6 n + 3 ( j = 3 , 4 , , n 1 )

③ 第n列至第2n列: D + H = [ 0 ( 2 n + 1 ) × n 3 ( 2 n + 1 ) × n ] ,所以 S c ( j ) = 0 × ( 2 n + 1 ) + 3 ( 2 n + 1 ) = 6 n + 3 ( j = n , n + 1 , , 2 n )

④ 第 ( 2 n + 1 ) 列至第3n列: J + M = [ 1 ( 2 n + 1 ) × ( n 1 ) 2 ( 2 n + 1 ) × ( n 1 ) ] ,所以 S c ( j ) = ( 2 n + 1 ) + 2 ( 2 n + 1 ) = 6 n + 3 ( j = 2 n + 1 , 2 n + 2 , , 3 n )

⑤ 第 ( 3 n + 1 ) 列到第 ( 4 n + 2 ) 列: K + N = [ 2 ( 2 n + 1 ) × ( n + 2 ) 1 ( 2 n + 1 ) × ( n + 2 ) ] ,所以 S c ( j ) = ( 2 n + 1 ) + 2 ( 2 n + 1 ) = 6 n + 3 ( j = 3 n + 1 , 3 n + 2 , , 4 n + 2 )

因此矩阵A中, S c ( j ) = 6 n + 3 ( j = 1 , 2 , , 4 n + 2 )

3) S m d = 3 ( n + 1 ) + 0 × n + 2 ( n 1 ) + ( n + 2 ) = 6 n + 3

4) S c d = 0 × ( n + 1 ) + 3 n + ( n 1 ) + 2 ( n + 2 ) = 6 n + 3

所以在矩阵A中,存在 S r = S c = S m d = S c d = 6 n + 3 ,那么在 ( 2 n + 1 ) 2 × A 中也有 S r = S c = S m d = S c d = ( 2 n + 1 ) 2 ( 6 n + 3 )

又由于T是奇数阶和幻方,必存在 S r = S c = S m d = S c d = S T

所以在 P = ( 2 n + 1 ) 2 × A + ( T T T T ) 中,有 S r = S c = S m d = S c d = S P ,同时可以计算得出 S P = 2 ( 4 n 3 + 6 n 2 + 4 n + 1 ) + ( 6 n + 3 ) ( 2 n + 1 ) 2

根据构造的定理可以得出以下性质:

推论1 如果一个4n + 2阶矩阵P是一个始元和幻方,则 P T 也是一个始元和幻方。

证明:矩阵P的第i行第j列元素就是矩阵 P T 的第j行第i列元素,即有 [ P ] i j = [ P T ] j i 。通过矩阵的转置并未改变幻方中原有的元素,所以在矩阵 P T 中的元素也是 1 [ 2 ( 2 n + 1 ) ] 2 中的互不相同的数。由于矩阵的转置只是对行与列的元素进行了整体交换,并未改变矩阵中每行每列的总和,即原矩阵的行和等于现矩阵的列和,则在 P T 中有 S r = S c = S m d = S c d = 2 ( 4 n 3 + 6 n 2 + 4 n + 1 ) + ( 6 n + 3 ) ( 2 n + 1 ) 2 。因此 P T 也是一个始元和幻方,且 S P = 2 ( 4 n 3 + 6 n 2 + 4 n + 1 ) + ( 6 n + 3 ) ( 2 n + 1 ) 2

推论2 如果一个 4 n + 2 阶矩阵P是一个始元和幻方,令 A = ( a ) ( 4 n + 2 ) ( 4 n + 2 ) , a N * ,则 A + P 是一个连元和幻方。

证明:由于P是一个始元和幻方,则P中元素互不相同,所以 A + P 中的元素必然互不相同。在 A + P 中,只是给原矩阵P中每个元素都加一个相同的数a,所以 A + P 中的每行每列的和相比矩阵P中的每行每列的和都增加了 a ( 4 n + 2 ) ,则在 A + P 中有 S r = S c = S m d = S c d = 2 ( 4 n 3 + 6 n 2 + 4 n + 1 ) + ( 6 n + 3 ) ( 2 n + 1 ) 2 + a ( 4 n + 2 ) 。因此 A + P 是一个连元和幻方,且 S A + P = 2 ( 4 n 3 + 6 n 2 + 4 n + 1 ) + ( 6 n + 3 ) ( 2 n + 1 ) 2 + a ( 4 n + 2 )

推论3 如果一个 4 n + 2 阶矩阵P是一个始元和幻方,则 a N * ,aP是一个异元和幻方。

证明:由于P是一个始元和幻方,则P中元素互不相同,所以在aP中的元素必然互不相同。在aP中,只是给原矩阵P中每个元素乘以一个相同的数a,所以在aP中的每行每列的和都是矩阵P中每行每列的和的a倍,则在aP中有 S r = S c = S m d = S c d = a [ 2 ( 4 n 3 + 6 n 2 + 4 n + 1 ) + ( 6 n + 3 ) ( 2 n + 1 ) 2 ] 。因此aP是一个异元和幻方,且 S a P = a [ 2 ( 4 n 3 + 6 n 2 + 4 n + 1 ) + ( 6 n + 3 ) ( 2 n + 1 ) 2 ]

4. 小结

单偶数阶和幻方作为一类基本幻方,其构造方法并不唯一,本文仅仅为读者提供一种构造单偶数阶和幻方的新思路,即利用分块矩阵同时结合已知奇数阶和幻方构造单偶数阶和幻方。关于此方法是否可以构造奇数阶、双偶数阶和幻方,有待于日后再做进一步的研究。

基金项目

地区科学基金项目(12161086)。

NOTES

*第一作者。

#通讯作者。

参考文献

[1] 何敏梅, 刘兴祥. 2k+1阶连元幻方的矩阵构造法[J]. 延安大学学报(自然科学版), 2017, 36(2): 11-13+17.
[2] 郭萍, 刘兴祥, 何敏梅. 偶数阶行列和始元幻阵的构造方法[J]. 河南科学, 2018, 36(4): 482-485.
[3] Johnson, C.R. (2018) A Matrix Theoretic Construction of Magic Squares. The American Mathematical Monthly, 79, 1004-1006.
[4] 郭萍, 刘兴祥. 和幻阵的定义及代数性质[J]. 延安大学学报(自然科学版), 2017, 36(1): 21-27.
[5] 张婧, 刘兴祥, 施钊. 完美和幻方的定义及其构造方法[J]. 湖北民族大学学报(自然科学版), 2020, 38(4): 420-423.
[6] 董朦朦. 幻性整数矩阵的保持性与JAD猜想[D]: [硕士学位论文]. 延安: 延安大学, 2020.
[7] 董朦朦, 刘兴祥, 田雨禾. 幻方可以分解为两个正交拉丁方的线性组合[J]. 重庆理工大学学报(自然科学), 2018, 32(8): 181-188.
[8] Hon, R.A. and Johnson, C.R. Matrix Analysis (卷1) [M]. 北京: 人民邮电出版社, 2005.
[9] Hon, R.A. and Johnson, C.R. Matrix Analysis (卷2) [M]. 北京: 人民邮电出版社, 2005.
[10] Hon, R.A. and Johnson, C.R. Matrix Analysis (卷3) [M]. 北京: 人民邮电出版社, 2009.