数学模型第六章2.docVIP

幻方趣谈 一、幻方的概念 幻方在我国古代也称为纵横图，在西方也称为魔方(magic). 定义 若一个 n阶由1~n2的正整数组成，且每行、每列与两对角线上的n个元素之和都相等. 则称此矩阵为n阶幻方. 每行的n个元素之和称为幻和，并记为Sn. 例如，下面分别是3阶幻方和4阶幻方 8 1 6 3 5 7 4 9 2 1 15 14 4 12 6 7 9 8 10 11 5 13 3 2 16 显然，一个幻方经旋转或转置后，仍为同阶幻方. 由定义可知，幻和的计算公式 S3=3(1+9)/2=15, S4=4(1+16)/2=34, S5=5(1+25)/2=65 注：不存在2阶幻方 二、幻方的起源 传说,我国远在夏禹治水时(公元前23世纪), 陕西的洛河常常泛滥成灾，威胁着两岸人们的生活与生产. 于是，大禹日夜奔忙，三过家门而不入，带领人们开沟挖渠，疏通河道，驯服了河水，感动了上天. 事后，一只神龟从河中跃出, 背上有一个九种花纹的图,后人把这个图称为“洛书”. 它就是从1到9连续自然数排成3行的图. 4 9 2 3 5 7 8 1 6 此图我国古代也称为九宫图. 最早见于记载的4阶纵横图,是在印度卡俱拉霍地方发现的一个11世纪的碑文上. 它是一个极不平凡的4阶纵横图,有着十分玄妙的性质. 7 12 1 14 2 13 8 11 16 3 10 5 9 6 15 4 它除了一般四阶幻方的通有的性质外，还有如下特性： （1）任一“折断的对角线”上4个数之和也等于幻和34（2）任一2阶子阵的4个数之和也等于幻和34（3）任一3阶子阵的4角4个数之和也等于幻和34（4）任一3阶子阵的2对角数之和恰是幻和34的一半17 顺便提一下，1977年美国发射寻求星外文明的宇宙飞船旅行者1号、2号上除了携带向宇宙人问候的“地球之声”(古今音乐、近六十种语言的问候话,三十五种自然界的各种声响唱片)外,还带了一些图片,其中有这张四阶幻方图. 1980年，上海博物馆在整理明代古墓的出土文物时，发现了一块玉佩上有一个四阶幻方，它也有上述玄妙的性质： 8 11 14 1 13 2 7 12 3 16 9 6 10 5 4 15 三、自然顺序方阵及其性质 定义 把自然数1~n2从小到大排成n阶方阵: ， ，把A称为自然顺序n阶方阵. 把行（列）号之和等于?n+1的两行（列）称为对称行（列），当n为奇数时，设n=2k-1，称为中心数，它位于A的中央. 位于对称行（列）同列（行）的两个数与（）称为行（列）对称数，而关于中心对称的两个数与称为对称数. 位于不同行，不同列的数称为独立数. 对于矩阵A,有如下重要的性质： 性质1. 任两个对称数与之和都是. 证. . 性质2. 任意n个独立数之和为幻和Sn. 证. 设是A的n个独立数，则与都是从1到n各取值一次，故. 因此 推论.主(次)对角线上n个数之和为Sn.任一折断的对角线上n个数之和也为Sn. 性质3 任两个对称行（列）的2n个数之和都等于2 Sn. 证.第i行的行和为 ，（i=1,2,…,n）. 从而 第j列的列和为 , （j=1,2,…,n）. 从而 性质4 当n=2k-1时，第k(中间)行(列)的n个数之和为Sn. 证.第k行的行和 第k列的列和 四、奇数阶幻方的构造方法及原理 在这里，设n=2k?1，（k=1,2,…）. (一) 构造方法1----连续摆数法 只需按以下步骤填写，即可得到一个n阶幻方. (1) 先画一个n×n方格表； (2) 把1填写在第一行中间； (3) 当m填好后，若m的右上方空，则把m+1填在此格，否则，把m+1填在m的下方.(注意，这里我们把最左列视作在最右列的右方，把最底行视作在第一行的上方) 17 24 1 8 15 23 5 7 14 16 4 6 13 20 22 10 12 19 21 3 11 18 25 2 9 例如 填写一个3阶幻方和5阶幻方 8 1 6 3 5 7 4 9 2 可验证其幻和分别为15和65. 设是按以上方法构造的n阶方阵，是自然顺序方阵，仔细观察知，A的每一行对应B的一条折断对角线或对角线，从而推得： ， （1） 注意，上式，右边的下标是按模n取值的，即大于n时就减n,小于1时就加n. (二)．原理 1．B的列和 从（1）式可见，若给定j, 让i取1，2，…，n时, i每增加1，（1）式右边的行号与列号也分别增加1(mod n)，最终都取遍1~n, 不会有重复，即B的每列数都是 A的n个独立数. 故其和是Sn. 2．B的行和 从（1）式可见，若给定i, 让j取1，2，…，n时, j每增加1，（1）式右边的行号增加1(mod n), 取遍1~n,不重复; 列号分别增加2(mod n)，一旦大于n就减去n(奇数),