魔方阵实现的方法.docxVIP

（1）当n为奇数时，即n=2*m+1时，算法为：首先，把1填在第一行的正中间。其次，若数k填在第i行第j列的格子中，那么k+1应填在它的左上方,即i-1,j-1；如果左上方没有格子，若i-1=0,那么k+1填在第n行第j-1列的格子中;若j-1=0,那么k+1填在第i-1行第n列的格子中 ;若i-1=0,j-1=0,那么k+1填在第n行第n列的格子中。再次，若按上述方法找到的格子都已经填过了数，那么，数k+1填在第k个数的正下方。具体实现时：用i0,j0表明格子不存在。找到的格子都已经填过了数的这类数有个特点，那就是可以整除n，即k%n=0。（2）当n为双偶数时，即n=2*2*m时，算法为：将所要求的魔方阵分为上下左右4个n/2*n/2的小方阵。首先，在左上角的方阵中做记号，每行每列，各取一半打上记号。然后，将其向其余3个方阵中映像。接着，从左上角格子开始,按从左到右,从上到下的次序将1到n*n的值往方阵中填写,但遇到作了标记的格子时跳过。再从右下角格子开始,按从右到左,从下到上的次序将1到n*n的值往方阵中填写,但遇到已填过值的格子时跳过，既可以构成一个魔方阵。具体实现时：左上角和右下角的n/2*n/2的小矩阵,偶行偶列,奇行奇列赋值为0,其余赋值为1；右上角和左上角的n/2*n/2的小矩阵,偶行奇列,奇行偶列赋值为0,其余赋值为1。然后，从左上角格子开始,按从左到右,从上到下的次序将1到n*n的值往方阵中赋,但遇到值为0时跳过；再从右下角格子开始,按从右到左,从下到上的次序将1到n*n的值往方阵中赋,但遇到值为1时跳过。（3）当n为单偶数时，即n=2*（2*m+1）时，算法为：把n阶模仿阵先均分为A,B,C,D四个同样的小方阵，先按先前奇数阶魔方的形成方法在A,B,C,D中构成四个奇阶魔方，其中，A用1到n/2*n/2填写；B用n/2*n/2+1到2*n/2*n/2填写；C用2*n/2*n/2+1到3*n/2*n/2填写；D用3*n/2*n/2+1到n*n填写。先在A的中间一行从左侧的第二列起取m个格子，把这些格子中的数字同D相应格子中的数字对调；然后，在C中从最后一列起在各行中取m-1个格子，把这些格子中的数字同B相应格子中的数字对调。具体实现时：设A,B,C,D四个二维数组，对于A调用形成奇阶魔方的函数array_pointerjishu(intn,int matrix1[][100])，将A用1到n/2*n/2排成魔方阵。B,C,D魔方阵有个特点就是分别在前一个魔方阵中各个数字上加n/2*n/2所构成的新魔方阵，则设了一个实现此功能的函数array_pointer first(intn,int A2[][100],int B2[][100])。然后，采用引用的方法设了一个实现相应位置上函数交换的函数void swap(int w,int v)。最后，将四个小方阵联结起来构成所要求的魔方阵方法一；/*此程序使我按照上面介绍的魔方阵的规律编写的,不过只能求奇数魔方阵,经过测试可以算到508阶*/#define N 7#includestdio.hvoid main(){int a[N][N],i,j,k; for(i=0;iN;i++) /*先处理第一行*/for(j=0;jN;j++) { a[i][j]=0; /*先令所有元素都为0*/ } j=(N-1)/2; /*判断j的位置*/ a[0][j]=1; /*将1放在第一行中间一列*/ for(k=2;k=N*N;k++) /*再从2开始处理*/ {i=i-1; /*存放的行比前一个数的行数减1*/ j=j+1; /*存放的列比前一个数的列数加1*/ if((i0)(j==N)) /*前一个数是第一行第N列时,把下一个数放在上一个数的下面*/ {i=i+2; j=j-1; }else { if(i0) /*当行数减到第一行，返回到最后一行*/i=N-1; if(jN-1) /*当列数加到最后一行，返回到第一行*/ j=0; } if(a[i][j]==0) /*如果该元素为0,继续执行程序*/a[i][j]=k; else /*如果该元素不为0,就说明要填的数的位置已经被占,则该数放在上一个数的下面*/ {i=i+2; j=j-1;a[i][j]=k; } } for(i=0;iN;i++) /*输出数组*/ {for(j=0;jN;j++)pr