各阶幻方的构造算法
各阶幻方的构造算法
奇数阶幻方
拉-卢贝尔算法
这个算法又称“阶梯法”。算法如下:
- 将1置于第一行中间。
- 将下一个数字置于当前数字的右上角。如果已经处于方阵的边界,则放在方阵的对边(如图1中的2和4)。
- 若出现下面两种情况,则将下一个数字放于当前数字的下方:
- 当前位置的右上角已经有一个数字(如图2中的6和11)。
- 当前位置已经是方阵的右上方(如图2中的16)。
- 结束,如图3.
拉-卢贝尔算法
菱形算法
另一种由康韦(J.H.Conway)建立的算法被称为“菱形算法”,步骤如下(以5x5为例):
- 从左边中间开始,将奇数在方阵内填成一个菱形。
- 将方阵分成5条对角线,每条对角线上有5个方格。如果图1所示。
- 从第一条对角线开始将偶数填入剩余的空格内,图2中填满了前两条对角线。
- 结束,如图3。
菱形算法
单偶数阶幻方
侓克斯算法
这个算法也是由康韦给出的。思想是将方阵分成多个2x2的小方阵,小方阵按照位置分成L、U、X三种类型。然后在大体上按照卢-拉贝尔算法来走,在每个小方阵中根据小方阵的类型来填数。具体算法如下:
- 将方阵分成(2i+1)个(2x2)的小方阵,小方阵的分类这样确定:前i+1行是L类型,后面一行是U类型,最后的i-1行是X类型,然后交换第i+1行和第i+2行中间小方阵的类型。对于10x10的方阵如图1。
- L、U、X的填法如图2。
- 最终结果如图3。
LUX算法
加法算法
将一个幻方加上另外一个幻方所得的和仍然具有幻方的特性,只是可能会有重复项,这是幻方的加法特性。下面的方法就是根据这个特性设计的,首先建立两个方阵A、B,具有幻方的特性(横、纵、斜和相同),然后让A加上B的i倍,就得到一个幻方。假如我们要作一个4i+2阶幻方(此处以14为例)。具体算法如下:
- 先作一个14(4i+2)阶的方阵A,这个方阵分成4个7(2i+1)阶小方阵,每个小方阵是一个奇数阶的幻方,奇数阶幻方构造方法已经有了。如图1。
- 再作一个14(4i+2)阶的方阵B,这个方阵只由0、1、2、3构成,具体作法如下:
- 第一列:3(i)个3,4(i+1)个0,5(i+2)个2,2(i-1)个1
- 前7(2i+1)列与都与第一列相同,只有第4(i+1)列例外,该列第一个3和第一个0交换位置。
- 后7(2i+1)列与前7(2i+1)列相同,不过3和0交换,1和2交换。
- 结果如图2。
- 构造幻方C=A+i2B。如图3,C即所求。
加法算法
替代算法
这个是最难的一个。对于2i+2的幻方(以6为例):
- 将幻方分成A、B、C、D、E、F、G、H和I几个区,如图1。
- 类似于双偶数阶幻方的分割算法,将处于上述分区中的格子填数,如图2。
- 然后在空格上(i+1和3i+1行列)填上数,如图3。
替代算法
好吧,太难了,中译本翻译得又很不清楚,不写了,等有机会看到更好的版本再加上。
双偶数阶幻方
分割算法
- 将方阵分成16个小方阵,如图1。
- 先在A、C、E、G、I方阵中填入数字,其他方阵跳过,如图2。
- 再逆序(从右下往左上)赶往余下的数字,如图3。
- 结束,如图3
分割算法
下面是一个8次的方阵:
8次方阵
对角线算法
- 将数字顺序填入方阵内,如图1。
- 将方阵分成四个相同大小的方阵。并找出每个小方阵的对角线,如图1阴影部分。
- 将阴影部分旋转180度,如图2。
对角线算法