魔方矩阵 (维数为奇数的情况)

魔方矩阵 (维数为奇数的情况)

最近换了OS，很多以前的资料丢失了...这个还在...

实现N*N的奇数维数的魔方矩阵，偶数的魔方矩阵可以实现的， 但是觉得麻烦 ，以后再去折腾偶数的情况吧 

算法介绍：

N 为奇数时

(1) 将1放在第一行中间一列;

(2) 从2开始直到n×n止各数依次按下列规则存放：

按 45°方向行走，如向右上

每一个数存放的行比前一个数的行数减1，列数减1

(3) 如果行列范围超出矩阵范围，则回绕。

例如1在第1行，则2应放在最下一行，列数同样减1;

(4) 如果按上面规则确定的位置上已有数，或上一个数是第1行第n列时，

则把下一个数放在上一个数的下面。

N为4的倍数时

采用对称元素交换法。

首先把数1到n×n按从上至下，从左到右顺序填入矩阵

然后将方阵的所有4×4子方阵中的两对角线上的数关于大方阵中心作中心对称交换（注意是各各子矩阵对角线上面的数）， 即a(i,j)与a(n+1-i,n+1-j)交换，

所有其它位置上的数不变。(或者将对角线不变，其它位置对称交换也可)

N 为其它偶数时

当n为非4倍数的偶数(即4n+2形)时：首先把大方阵分解为4个奇数(2m+1阶)子方阵。

按上述奇数阶魔方给分解的4个子方阵对应赋值

上左子阵最小(i)，下右子阵次小(i+v)，下左子阵最大(i+3v)，上右子阵次大(i+2v)

即4个子方阵对应元素相差v，其中v=n*n/4

四个子矩阵由小到大排列方式为 ① ③ ④ ②

然后作相应的元素交换：a(i,j)与a(i+u,j)在同一列做对应交换(j<t-1或j>n-t+1),

注意其中j可以去零。

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 上述交换使每行每列与两对角线上元素之和相等。

*/

#include<stdio.h>
#include<stdlib.h>

int main()
{
	int n=0;//矩阵的维数
	int num=0;//赋值变量
	int start_coordinate;
	int row=0,rank=0;
	printf("Please enter the rank number of the cube number array \
	and \nmake sure the number is odd number or the anwers is worng:\n");

	while((!scanf("%d",&n))||(n%2==0))
	{
		printf("The data you entered is wrong ! Please enter again!\n");
	}

	//***********************申请一个N*N的内存空间*********************************************
	int **pp_num=(int**)malloc(n*sizeof(int*));

	for(row=0;row<n;row++)
	{
		*(pp_num+row)=(int*)malloc(n*sizeof(int));
	}
	//*****************************************************************************************

	for(row=0;row<n;row++)//初始化
	{
		for(rank=0;rank<n;rank++)
		{
			pp_num[row][rank]=0;
		}
	}

	start_coordinate=n/2;//计算开始填充时,1的坐标点，1始终放在第一行！算法要求
	row=0;
	rank=start_coordinate;

	for(num=1;num<=n*n;num++)
	{
		pp_num[row][rank]=num;
		rank++;
		row--;

		if(rank==n)
		{
			rank=0;
		}

		if(row<0)
		{
			row=n-1;
		}

		if(pp_num[row][rank]!=0&&!(row==n-1&&rank==0))
		{
		    row=row+2;
		    rank--;
		}
		else if(row==n-1&&rank==0)
		{
		     row=1;
		     rank=n-1;
		}

	}

	for(row=0;row<n;row++)
	{
	     for(rank=0;rank<n;rank++)
	     {
		printf("\t%d",pp_num[row][rank]);
	     }

	     if(rank==n)
	     {
	      	printf("\n\n");
	     }
	}

	for(row=0;row<n;row++)//释放掉申请的内存空间！！！
	{
		free(*(pp_num+row));
	}
	free(pp_num);

	pp_num=NULL;

	return 0;
}

应该是13年6月份写的东东~ Just a demo.