魔方矩阵

看到魔方矩阵，好奇，好玩儿，正好赶上周五，就来放松一下，总结一下几种魔方矩阵的规律，并写一下C++实现过程

定义：平面魔方的一般定义：将自然数 1 到 N^2， 排列 N 行 N 列的方阵，使每行、每列及两条主对角线上的 N 个数的和都等于 (N^2+1)/2，这样的方阵称为 N 阶幻方。

此定义摘自百度百科，话说百度百科的这个定义还靠谱，那后面的例子就不行了，有错误。

1、N为奇数

（1）将1放在第一行中间一列。
（2）从2开始直到n×n止各数依次按下列规则存放：
 按右上45°方向行走，也就是说下一个数的行数比前一个数的行数减1，列数加1
<1>.如果行范围超出矩阵范围，则行数变为n-1行，列不变。
<2>.如果列范围超出矩阵范围，则列数变为0列，行不变。
<3>.如果按上面规则确定的位置上已有数，或上一个数是第0行第n-1列时，则把下一个数放在上一个数的下面。

（3）按照上述规则(1)(2)将N^2 个数都放入矩阵中，得到的矩阵即为N阶魔方矩阵。

demo:

8	1	6
3	5	7
4	9	2

code:

奇数的实现
#include <stdio.h>
#include <math.h>
int main()
{
	int a[16][16],i,j,n,k;
	printf("Plase input n(1-15,it must be odd.): ");
	scanf("%d",&n);
	while(!(n>=1&&n<=15)||n%2==0)
	{
		printf("The number is invalid.Plase insert again:");
		scanf("%d",&n);
	}
	for(i=1;i<=n;i++)
		for(j=1;j<=n;j++)
		a[i][j]=0;
	j=n/2+1;
	a[1][j]=1;
	i=1;
	for(k=2;k<=n*n;k++)
	{
		i=i-1;j=j+1;
		if(i==0&&j==n+1)
		{
			i=i+2;j=j-1;
		}
		else
		{
			if(i==0)i=n;
			if(j==n+1)j=1;
		}
		if(a[i][j]==0)a[i][j]=k;
		else
		{
			i=i+2;
			j=j-1;
			a[i][j]=k;
		}
	}
	for(i=1;i<=n;i++)
	{
		for(j=1;j<=n;j++)
			printf("%5d",a[i][j]);
		printf("\n");
	}
	return 0;
}

2、N为4的倍数

采用对称元素交换法：


（1）把数1到n×n按从上至下，从左到右顺序填入矩阵。 
（2）将方阵划分为4×4子方阵，每个子方阵中画出子主对角线、子副对角线，然后按照下面步骤操作：
<1>.主对角线上的元素关于原点对称交换。
<2>.副对角线上的元素不动。
<3>.子对角线（子主对角线、子副对角线）关于副对角线对称交换。
<4>.各子方阵中其余元素关于主对角线对称交换。
（3）按照上述规则(1)(2)操作元素，得到的矩阵即为N阶魔方矩阵。

 最后得到的8阶魔方矩阵为：

 最后得到的8阶魔方矩阵为：

64	2	3	61	60	6	7	57
9	55	54	12	13	51	50	16
17	47	46	20	21	43	42	24
40	26	27	37	36	30	31	33
32	34	35	29	28	38	39	25
41	23	22	44	45	19	18	48
49	15	14	52	53	11	10	56
8	58	59	5	4	62	63	1

3、N为其它数（4N+2）

（1）把大方阵分解为4个奇数(2m+1阶)子方阵。
（2）按照第1节中奇数阶魔方矩阵填数的方法给分解的4个子方阵赋值，四个子方阵赋值规则为：
上左子阵最小(i)，下右子阵次小(i+v)，下左子阵最大(i+3v)，上右子阵次大(i+2v)，即4个子方阵对应元素相差v，其中v=n*n/4。
四个子矩阵由小到大排列方式为 ① ③ ④ ②。
（3）按照（2）中描述赋值之后，作相应的元素交换：
<1>a(i,j)与a(i+u,j)在同一列做交换(j<t-1或j>n-t+1，矩阵行列编号从0开始)。
<2>a(t-1,0)与a(t+u-1,0)；a(t-1,t-1)与a(t+u-1,t-1)两对元素交换。
其中u=n/2，t=(n+2)/4。
(4)到此完成。

参考：http://www.cnblogs.com/yaochc/p/3404179.html