判断魔方阵c语言程序设计_魔方阵原理及十种解法（C语言）

一、魔方阵的简介1.何谓矩阵？矩阵就是由方程组的系数及常数所构成的方阵。把用在解线性方程组上既方便，又直观。2.何谓n阶方阵？若一个矩阵是由n个横列与n个纵行所构成，共有个小方格，则称这个方阵是一个n阶方阵。3.何谓魔方阵？4 9 2 3 5 7 8 1 6定义：由n*n个数字所组成的n阶方阵，具有各对角线，各横列与纵行的数字和都相等的性质，称为魔方阵。而这个相等的和称为魔术数字。若填入的数字是从1到n*n，称此种魔方阵为n阶正规魔方阵。4.最早的魔方阵相传古时为了帮助治水专家大禹统治天下，由水中浮出两只庞大动物背上各负有一图，只有大禹才可指挥其中之由龙马负出的为河图，出自黄河；另一由理龟负出的洛书出自洛河。洛书5.最早的四阶魔方阵最早的四阶方阵刻在印度一所庙宇石上，年代大约是十一世纪。古代印度人十分崇拜这种幻方，至今从古神殿的遗址，墓碑上常常还可以发现四阶幻方的遗迹。6.欧洲最早的魔方阵欧洲最早的方阵是公元1514年德国画家Albrecht Dure在他著名的铜板画Melencolia上的4×4幻方，有趣的是，他连创造年代(1514)也镶在这个方阵中，而且上下左右，四个小方阵的和皆为34，是欧洲最古老的幻方。二、魔方阵的变形1.相异魔方阵对应方格所填的数字不相同的两个魔方阵，如下图16 02 03 13 01 15 04 14 05 11 10 08 12 06 09 07 09 07 06 12 13 03 16 02 04 1415 01 08 10 05 11

2.相等魔方阵每一个对应方格所填的数字都相同的两个魔方阵，如下图4 9 2 4 9 2 3 5 7 3 5 7 8 1 6 8 1 6

3.全等魔方阵若一个魔方阵能够经过旋转或镜射的方式变成和另一个魔方阵相等，称这两个魔方阵全等。在计算魔方阵个数时，我们把全等的魔方阵视为同一种。2 7 6 8 3 4 9 5 11 5 9 4 3 8 6 7 2

4,刚性变形法(1)顺时针方向旋转90 (2)顺时针方向旋转180 (3)顺时针方向旋转270 (4)左右翻转(绕铅直对称轴镜射) (5)上下翻转(绕水平对称轴镜射) (6)左上右下翻转(绕右上至左下对角线镜射) (7)右上左下翻转(绕左上至右下对角线镜射)原始方阵旋转90。旋转180。旋转270。01 03 16 14 12 08 13 01 07 05 10 12 14 04 09 07 13 15 02 04 10 06 15 0309 11 06 08 16 02 11 05 08 06 11 09 05 11 02 16 04 02 15 13 03 15 06 10 12 1005 07 07 09 04 14 14 16 03 01 01 13 08 12左右翻转上下翻转左上右下翻转右上左下翻转14 16 03 01 12 10 05 07 07 09 04 14 01 13 0812 04 02 15 13 08 06 11 09 05 11 02 16 03 15 06 10 09 11 06 08 13 15 02 04 1006 15 03 16 02 11 05 07 05 10 12 01 03 16 14 12 08 13 01 14 04 09 07 5.加值变形法正规魔方阵：首项是1，公差是1加值变形后的魔方阵：首项是a，公差是r原始魔方首项为5公差为1首项为5公差为2 14 16 03 01 18 20 07 05 31 35 09 05 04 02 15 13 08 06 19 17 11 07 33 2909 11 06 08 13 15 10 12 21 25 15 19 07 05 10 12 11 09 14 16 17 13 23 27

6.互补变形法就是将魔方阵中的每一个数字都替换成互补数的变形方式。在n阶魔方阵中，数字k的互补数=(1+n*n)–k原始魔方阵互补变形魔方阵14 16 03 01 03 01 14 16 04 02 15 13 13 15 02 04 09 11 0608 08 06 11 09 07 05 10 12 10 12 07 05

7.田字变形法以中心点为准将魔方阵分成四个相等的小方阵。当n为奇数时，中央的行列要独立出来。n是偶数时n是奇数时ABA A1 BDC? D1 E B1 D C1 C将魔方阵如下重组: CD ? C C1 D? BA B1E D1 B A1 A可参考下面的范例:原始4阶魔方阵4阶田字变形魔方阵原始5阶魔方阵5阶田字变形魔方阵01 03 16 14 11 09 08 06 18 22 01 10 14 13 17 09 21 05 13 15 02 04 05 07 1210 04 08 12 16 25 24 03 20 07 11 08 06 11 09 16 14 01 03 15 19 23 02 06 02 0623 15 19 12 10 05 07 02 04 13 15 21 05 09 13 17 10 14 01 18 22 07 11 20 24 0316 25 12 04 08 8.井字对换变形法(1)任选一数k，1<=k<=n。(2)将方阵的第k行和其互补行(第n+1-k行)对换。(3)将方阵的第k列和其互补列(第n+1-k列)对换。(4)为方便称呼，此时姑且命名为k值井字对换变形。原始魔方阵1值井字对换魔方阵2值井字对换魔方阵1,2值井字对换魔方阵01 03 16 14 07 10 05 12 01 16 03 14 07 05 10 12 13 15 02 04 04 15 02 13 0811 06 09 09 11 06 08 08 06 11 09 09 06 11 08 13 02 15 04 04 02 15 13 12 10 0507 14 03 16 01 12 05 10 07 14 16 03 01 9.拓朴变形法( 1)任选不相等的两数1<=k1,k2<=n，但当(n+1)/2为奇数时，k1及k2不能等于(n+1)/2。(2)将座标含k1值的全改成k2。含k2值的全改成k1。含n+1-k1值的全改成n+1-k2。含n+1-k2值的全改成n+1-k1。(3)为方便称呼，此时姑且命名为k1,k2拓朴变形。可参考下面的3阶1，3拓朴变形示意图(1,1) (1,2) (1,3) (3,3) (3,2) (3,1) (2,1) (2 ,2) (2,3) (2,3) (2,2) (2,1)(3,1) (3,2) (3,3) (1,3) (1,2) (1 ,1)可参考下面的范例原始魔方阵1,2拓朴变形魔方阵18 22 01 10 14 08 0412 25 16 04 08 12 16 25 22 18 01 14 10 15 19 23 02 06 19 15 23 06 02 21 05 0913 17 11 07 20 03 24 07 11 20 24 03 05 21 09 17 13注：其实对换变形、田字变形可包含于拓朴变形中，只因想法不同故列出参考三、奇数阶魔方阵的建构法1.杨辉法发明者：杨辉适用：三阶魔方阵方法：九子斜排，上下对易，左右相更，四维挺进***1*** **4*2** 4 9 2 *7*5*3*  3 5 7 **8*6**  8 1 6 ***9***杨辉法推广－菱形法发明者：Bachet de Meziriac适用：奇数阶魔方阵方法：数字斜排，上下对易，左右相更，四维挺进********01******** ******06**02****** ? ? ? ****11**07**03*** * 11 24 07 20 03**16**12**08**04** 04 12 25 08 16 21**17**13**09**05 17 05 13 21 09 **22**18**14**10** 10 18 01 14 22 ****23**19**15**** ? 23 06 19 02 15 ******24**20******* *******25********

2.简捷连续填制法发明者：De La Loubere适用：奇数阶魔方阵方法：1立首列中，右一上一，受阻下一* * 1 * * * * 1 * * * * 1 * * * * 1 * * * * 1 8 * 17 24 1 8 15 * * * * * * * ** * * 5 * * * * 5 * * * * 5 7 * * 23 5 7 14 16 * * * * * * * * * * 4 * * * * 46 * * * 4 6 * * * 4 6 13 20 22 * * * * * * * * * * * * * * 3 * * * * 3 10 * * *3 10 12 19 21 3 * * * * * * * * 2 * * * * 2 * * * * 2 * 11 * * 2 9 11 18 25 2 9简捷连续填制法推广适用：奇数阶魔方阵方法：(1)1立首列中，右1上1，受阻下1；(2)1立中央上，右1下2，受阻上2；(3 )1立首行中，右2下1，受阻右1；(4)1立首列中，右1下2，受阻下4；(5)1立中央上，右1上1，受阻上2；(6)1立首行中，左1上1，受阻右1。3.辅助方阵法适用：五阶以上奇数阶魔方阵方法：(1)制作辅助方阵一：在左上端填入1，按照走马的方式往右二下一的方向填入1至底，接着在每列的1的右边依序填入2、3、4…n 1 * * * * 1 23 4 5 1 2 3 4 5 * * 1 * * * * 1 * * 4 5 1 2 3 * * * * 1 * * * * 1 2 3 4 5 1 * 1* * * * 1 * * * 5 1 2 3 4 * * * 1 * * * * 1 * 3 4 5 1 2 (2)制作辅助方阵二：在左上端填入0，按照走马的方式往右一下二的方向填入0至底，接着在每列的0的下面依序填入n、2n、3n…n(n-1) 0 * * ** 0 15 5 20 10 1 17 8 24 15 * * * 0 * 5 20 10 0 15 9 25 11 2 18 * 0 * * * 10 015 5 20 12 3 19 10 21 * * * * 0 15 5 20 10 0 20 6 22 13 4 * * 0 * * 20 10 0 155 23 14 5 16 7 (3)将辅助方阵一与二的对应方格内之数相加填到一个新的方阵对应方格内，则新的方阵即为一个魔方阵。4.扩阶法适用：n阶魔方阵，n为大于5的正整数方法：(1)先建构出一个n-2阶的魔方阵(2)将(1)所建立的魔方阵每个数字再加上2n-2 (3)将步骤(2)的魔方阵外面再加上一圈，这一圈有4n-4个空格，填入的数为1~(2n-2)和( -2n+3)~n*n，先算出魔术数字会比较好填。8 1 6 16 09 14 23 01 02 20 19 3 5 7 11 13 15 22 16 09 14 04 4 9 2 12 17 10 0511 13 15 21 08 12 17 10 18 07 25 24 06 03 5.方阵合成法适用：当魔方阵的阶数n可分解成两个大于2之整数p、q的乘积方法：假设A=[(a)ij]与B=[(b)ij]分别代表m阶与n阶魔方阵，对每个t=1，2…m*m，令At=[(a)ij+m*m(t–1)]。接着我们将方阵B中的数t以m阶魔方阵(A)t代替，就可以得到一个m．n阶魔方阵。可参考下面的范例魔方阵A魔方阵B 35 28 33 26 19 24 71 64 69 8 1 6 4 3 830 32 34 21 23 25 66 68 70 3 5 7 9 5 1 31 36 29 22 27 20 67 72 65 4 9 2 2 7 680 73 78 44 37 42 08 01 06 75 77 79 39 41 43 03 05 07 A4 A3 A8 76 81 74 40 4538 04 09 02 A9 A5 A1 17 10 15 62 55 60 53 46 51 A2 A7 A6 12 14 16 57 59 61 4850 52 13 18 11 58 63 56 49 54 47四、偶数阶魔方阵的建构法1、杨辉法发明者：杨辉适用：4阶魔方阵方法：以十六子，依次递作四行排列，先以外四子对换，一换十六、四换十三，以四内角对换，六换十一、七换十，横直上下斜角，皆三十四数，对换止可施之于小。01 02 03 04 16 02 03 13 16 02 03 13 05 06 07 08 05 06 07 08 05 11 10 08 09 1011 12 09 10 11 12 09 07 06 12 13 14 15 16 04 14 15 01 04 14 15 01杨辉法推广－消去对角线法适用：四之倍数阶魔方阵方法：(1)先将整个方阵划分成k*k个4阶方阵，然后在每个4阶方阵的对角线上做记号( 2)由左而右、由上而下，遇到没有记号的位置才填数字，但不管是否填入数字，每移动一格数字都要加1 (3)自右下角开始，由右而左、由下而上，遇到没有数字的位置就填入数字，但每移动一格数字都要加1

2.井字法适用：四之倍数阶魔方阵方法：(1)把1~n*n从左上角依序填入方阵内(2)用两条铅直线和两条水平线将方阵分隔成四个角落各有一个n/4阶的子方阵，和中心位置有一个n/2阶的子方阵(3)以方阵中心为对称点，将五个子方阵的数字作对称交换，其它的数字不要动。这样的方阵会是一个魔方阵01 02 03 04 16 02 03 13 05 06 07 08 05 06 07 08 09 10 11 12 09 10 11 12 13 1415 16 04 14 15 01

3.辅助方阵法适用：偶数阶魔方阵(因非四的倍数作法相当复杂，在此只介绍四的倍数的作法)方法：(1)制作辅助方阵一：第一列由左向右排列1至n，第二列由右向左排列1至n，第三列同第二列，第四列同第一列，若超过四列，重复一至四列的作法，完成辅助方阵一(2)制作辅助方阵二：将第一个辅助方阵各方格内的数字x换成n(x-1)，再做行列互换，完成辅助方阵二(3)将辅助方阵一与二的对应方格内之数相加填到一个新的方阵对应方格内，则新的方阵即为一个魔方阵01 02 03 04 00 1212 00 01 14 15 04 04 03 02 01 04 08 08 04 08 11 10 05 04 03 02 01 08 04 04 0812 07 06 09 01 02 03 04 12 00 00 12 13 02 03 16

4.扩阶法同奇数阶之作法5.方阵合成法同奇数阶之作法