1701005566
1701005567
因此可知，采用该速乘算法是可行的。
1701005568
1701005569
1701005570
1701005571
1701005572
我和数学有约：趣味数学及算法解析 [:1701004189]
1701005573
我和数学有约：趣味数学及算法解析 3.5 三阶幻方
1701005574
1701005575
【问题】什么是三阶幻方？
1701005576
1701005577
【分析】
1701005578
1701005579
一个“三阶幻方”是指把数字1～9填入3×3的方格，使得每一行、每一列和两条对角线的三个数之和正好都相同。
1701005580
1701005581
表3-1就是一个三阶幻方，每条直线上的三个数之和都等于15。
1701005582
1701005583
表3-1 幻方
1701005584
1701005585
8 1 6 3 5 7 4 9 2 大家或许都听说过幻方这玩意儿，但不知道幻方中的一些美妙的性质。
1701005586
1701005587
例如，任意一个三阶幻方都满足各行所组成的三位数的平方和，等于各行逆序所组成的三位数的平方和。
1701005588
1701005589
对于表3-1中的三阶幻方，有：
1701005590
1701005591
1701005592
1701005593
1701005594
进行前四阶幻方分析，采用程序设计如下：
1701005595
1701005596
clc,clear,close all %清屏和清除变量 warning off %消除警告 magic(1) %1阶魔方 magic(2) %2阶魔方 magic(3) %3阶魔方 magic(4) %4阶魔方
1701005597
1701005598
运行程序输出结果如下：
1701005599
1701005600
ans = 1 ans = 1 3 4 2 ans = 8 1 6 3 5 7 4 9 2 ans = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1
1701005601
1701005602
进行正逆数字平方和验证如下：
1701005603
1701005604
＞＞ 13^2+42^2-31^2-24^2 ans = 396 ＞＞ 816^2+357^2+492^2-618^2-753^2-294^2 ans = 0 ＞＞ 162313^2+511108^2+97612^2+414151^2-313261^2-801115^2-21679^2-151414^2 ans = -2.9469e+11
1701005605
1701005606
整理结果如下：
1701005607
1701005608
1701005609
1701005610
1701005611
从结果可看出，对于一阶和三阶幻方是满足正逆顺序数字平方和相等的，对于二阶和四阶幻方是不满足的，那么对于五阶、六阶和七阶等等是否满足三阶幻方的关系呢？留待读者自己推导。
1701005612
1701005613
1701005614
1701005615