奇阶幻方的算术构造

题记

这是到目前为止，我所发现的最简奇阶幻方构造公式。为什么凭一条简洁的初等数学通项公式，就能填充任意奇阶幻方中的所有数字，且满足幻方的所有要求？这样看似无甚理由的上下、左右连加或连减，就不会出现重复数字或掉了一个数字？

——求质疑！

本法构造的幻方，先根据一条通项公式，求得坐标轴上的1至n²数列中的部分数后，再依据本法的构建规律，便可填充该幻方其余所有的数字。

先建立两两垂直的两对坐标B、D轴与A、C轴（图一，坐标上的点，用格子代替，以便填数），

通过如下坐标轴通项公式

Nₘ=n₀±[Jm+2m(m-1) ]
(m= 0，1，2，3…)

确定各坐标轴格子上的数字（坐标轴正方向，式中取+号，负方向取-号），再参照图例所示的排列规律，依次填写，便可求得任意奇阶幻方。

式中，Nₘ表示各个坐标轴要填的数。n₀是1至n²自然数数列的正中间数（n₀=（1+ n²）/2） 。J称初进数，对应A、B、C、D轴，J分别为1、2、3、4。m是各坐标格子序号，原点序号为0，填n₀。

可以证明，由此方法得到的幻方，每一行、每一列及两条对角线上的数字之和（称幻和）

∑= n（1+ n²）/2

均相等，满足幻方的要求。

一、无重复数的证明

证明备用公式 (由通项公式导出)

Aₘ =n₀+2m²-m， Aₘ＇＝n₀-2m²+m
Bₘ =n₀+2m² ， Bₘ＇＝n₀-2m²
Cₘ =n₀+2m²+m， Cₘ＇=n₀-2m²-m
Dₘ=n₀+2m²+2m， Dₘ＇=n₀-2m²-2m

我们先把构造分布图以四条坐标轴为界，分为8个区域（图一），根据通项公式和排列规律，导出8个区域的一般分布表式。要确定每个区域的一般分布表式，关键是确定每行（列）的空格数K与格子序号m的关系。由分布示意图容易看出，当m=1时，K=0；m=2时，K=1等等，即空格数与格子序号的关系，是K=m-1。

下面，我们先来求解各区域的一般分

区域Ⅰ₁.

由示意图可知，该区域是由横轴上的数字逐格减1得到的，也就是说，对于该区域来说，在格子序号为m的一般列上，分布着如下一组数字
Bₘ ，Bₘ -1，Bₘ -2，…，Bₘ -（m-1）。

这一组数字，构成了一个连续递减数列。为了比较的方便，我们把它转换成连续递增数列，并以如下的式子表达

Bₘ -（m-1）→Bₘ

（在本文的表式中，符号→，表示由符号两端的两个数确定的一个递增（递减）数列。另外，为了表述的方便，我们把相差为1的两个数，称为“连续”。）
将B轴通项公式Bₘ = n₀ +2m² 代入，求得该数列的一般表式为

n₀ +2m²-m+1→n₀ +2m² （1）

区域Ⅰ₂

由分布图可知，该区域的每行数字，是由C轴正方向上同序号的对应数字（Cₘ ）逐格减1得到的。所以，在格子序号为m的一般行上，分布的是如下一组数字

Cₘ -1，Cₘ -2，…， Cₘ -（m-1）。

同上，也把它转换成递增数列，则是

Cₘ -（m-1）→ Cₘ -1

将Cₘ 归在此数列中，

则是 Cₘ -（m-1）→ Cₘ

由此，按C轴通项公式Cₘ= n₀ +2m²+m，求得该数列的一般表式是

n₀+2m²+1→n₀+2m²+m（2）

依次类推，容易得到其余六个区域的表式（图二）。

这样，格子序号