魔方阵

import java.util.Scanner; public class magicSquare2 { public static int MAX=100; public static int[][] matrix=new int[MAX][MAX]; /**===================================================================** * [##private void init(){}]: * 参数：无 * 返回值：无 * 修饰符：public * 功能：初始化 matrix矩阵,初始化为0 **===================================================================** */ private void init() { for(int i=0;i<MAX;i++) { for(int j=0;j<MAX;j++) { matrix[i][j]=0; } } } /**====================================================================** * [##private void odd(int digit,int startrow,int startcol,int rowcol){}]: * 参数： int digit 表示开始数字 (测试为非负) * int startrow 表示开始行的值 * int startcol 表示开始列的值 * int rowcol 表示每行每列有多少值 * startrow,startcol,rowcol都表示非负 * 返回值：无 * 修饰符：private * 功能：制造奇数魔方阵 * **=====================================================================** */ private void odd(int startdigit,int startrow,int startcol,int rowcol) { int row,col; int row0,col0; row0=row=startrow; col0=col=startcol+rowcol/2; row=startrow; matrix[row][col]=startdigit; for(int i=startdigit+1;i<=startdigit+rowcol*rowcol-1;i++) { if((row-1)<startrow) { row=startrow+rowcol-1; }else { row--; } if((col+1)>startcol+rowcol-1) { col=startcol; }else{ col++; } //System.out.println(row+"/t"+col); if(matrix[row][col]==0) { //System.out.println("*********************"); row0=row; col0=col; matrix[row][col]=i; //System.out.println(i); }else{ row=row0; if((row+1)>startrow+rowcol-1) { row=startrow; }else{ row++; } col=col0; matrix[row][col]=i; //System.out.println(i); } } } /**=======================================================================** * [##private void odd_even(int rowcol){}]: * 参数：int rowcol 表示生成魔方阵的行列 * 返回值 ：无 * 修饰符：private * 功能：生成4n+2阶魔方阵 **========================================================================** */ private void odd_even(int rowcol){ int k=rowcol/2; odd(1,k,k,k); odd(k*k+1,0,0,k); odd(2*k*k+1,0,k,k); odd(3*k*k+1,k,0,k); //System.out.println("交换前(一次生成数据):"); //print(rowcol); int m=k/2; /*保存临时数据，便于交换*/ int t; for(int i=0;i<m-1;i++) { for(int j=0;j<k;j++) { t=matrix[i][j+k]; matrix[i][j+k]=matrix[i][j]; matrix[i][j]=t; } } //System.out.println("一次交换(交换前m-1行):"); //print(rowcol); for(int i=k+m+1;i<rowcol;i++) { for(int j=0;j<k;j++) { t=matrix[i][j+k]; matrix[i][j+k]=matrix[i][j]; matrix[i][j]=t; } } //System.out.println("二次交换(交换后m行):"); //print(rowcol); /**交换矩阵的中心值*/ t=matrix[k+m][m+k]; matrix[k+m][m+k]=matrix[k+m][m]; matrix[k+m][m]=t; //System.out.println("三次交换(交换中心值):"); //print(rowcol); /**交换矩阵的最大值*/ t=matrix[rowcol-1][m+k]; matrix[rowcol-1][m+k]=matrix[rowcol-1][m]; matrix[rowcol-1][m]=t; //System.out.println("四次交换(交换最大值):"); //print(rowcol); } /**=======================================================================** * [##private void even_even(int rowcol){}]: * 参数：int rowcol 代表生成魔方阵的行列数 * 修饰符：private * 返回值：无 * 功能：生成4n阶魔方阵 **========================================================================** */ private void even_even(int rowcol){ int location=0; for(int i=0;i<rowcol;i++) { for(int j=0;j<rowcol;j++) { /**从1开始计数*/ location=i*rowcol+j; if((i%4)==(j%4)||(i%4)+(j%4)==3) { matrix[i][j]=rowcol*rowcol-location; }else{ matrix[i][j]=location+1; } } } //print(rowcol); } /**=======================================================================** * [##private void print(int rowcol){}]: * 参数：int rowcol 表示生成的魔方阵行列 * 返回值：无 * 修饰符：private * 功能：打印魔方阵 **=======================================================================** */ private void print(int rowcol) { //System.out.println(); for(int i=0;i<rowcol;i++) { for(int j=0;j<rowcol;j++) { System.out.print(matrix[i][j]+"/t"); } System.out.println(); } } /**对生成的魔方阵进行校验,检查生成的魔方阵是否正确*/ private void verify(int rowcol){ int sum=0; int sum1=0; for(int i=0;i<rowcol;i++) { sum1=0; for(int j=0;j<rowcol;j++) { sum1+=matrix[i][j]; } if(sum==0){ sum=sum1; }else if(sum!=sum1){ System.out.println("Error"); return; } } for(int i=0;i<rowcol;i++) { sum1=0; for(int j=0;j<rowcol;j++) { sum1+=matrix[j][i]; } if(sum==0){ sum=sum1; }else if(sum1!=sum){ System.out.println("Error"); return; } } System.out.println("Verity OK! The number is "+sum1); } public static void main(String args[]) { magicSquare2 mg=new magicSquare2(); mg.init(); int rowcol=0; Scanner input=new Scanner(System.in); System.out.print("请输入一个正整数:"); while(true) { rowcol=input.nextInt(); if(rowcol<=0||rowcol>100) continue; else break; } //System.out.println(); if(rowcol%2==1) mg.odd(1,0,0,rowcol); else if(rowcol%4==2) mg.odd_even(rowcol); else mg.even_even(rowcol); mg.print(rowcol); mg.verify(rowcol); } }     先说奇数阶的算法，这是最容易的算法：
n=2*m+1,m为自然数
1)将数字1填在(0,(n+1)/2)   ；要注意c中是从下标0开始
2)从左上往右下依次填。
3)由2),列的下标出界(超过n-1)时，行加1，以n为摸的余数为应填的列数；
4)由2),行的下标出界(超过n-1)时，列加1，以n为摸的余数为应填的行数；
5)由2),行列都未出界，但已添上其他数，应在当前位置左横移一个位置进行填数。

然后是偶数阶：
分两种情况，一种是n%4==2,一种是n%4==0
前一种：n=2*(2*m+1),m为自然数
1)将n阶方阵分为四个小魔方阵ABCD如下排列：

B   C
D   A

因为n*n=4*(2*m+1)*(2*m+1),
记u=n/2=2*m+1,分为1~u*u,u*u+1~2*u*u,2*u*u+1~3*u*u,3*u*u+1~4*u*u
即在调用子函数的时候分别如下面传递参数:
A(0),B(u*u),C(2*u*u),D(3*u*u)
分别在ABCD中按照前面的填法把奇数阶填好(注意加上所传参数作为基数，每一个元素都要加上这个值)，最后做如下交换：
(1)B中第0~(m-1)-1行中元素与C中相对应元素交换
(2)D中第(n-1)-m+1~(n-1)共m行的每行中的元素与A中相对应元素交换
(3)交换D:(u+m,m)与A中对应元素(矩阵中心值)
(4)交换D:(n-1,m)与A中对应元素(实际为矩阵最大值n*n)

所谓对应位置，指相对于小魔方阵的左顶角的相对的行列位置
上面的这些你可以用数学进行证明，利用魔方阵常数(注意n阶的和u阶的关系)
 
后一种：n=4*m,m为自然数
因为行列都是4的倍数，因而可以将整个矩阵分为每4*4的小矩阵。
先判断一个数是否在划为4*4小矩阵的对角线上，
如果在，则填该位置的数为n*n-i+1(i为该元素的相对位置，从1开始，比如n阶的第s行第t个元素则其i=s*n+t)
如果不在，则填上i。

 

 

#include <stdio.h> #define N 100 unsigned matrix[N][N]; void init(void); void odd(unsigned startdigit,unsigned char startrow,unsigned char startcol,unsigned char rowcol); void odd_even(unsigned char rowcol); void even_even(unsigned char rowcol); void print(unsigned char rowcol); void verify(unsigned char rowcol); void init(void){ unsigned char i,j; for(i=0;i<N;i++) { for(j=0;j<N;j++) { matrix[i][j]=0; } } } void odd(unsigned startdigit,unsigned char startrow,unsigned char startcol,unsigned char rowcol) { unsigned char row,col,row0,col0; unsigned i; //这一句可以省去吧 row0=row=startrow; col=startcol+(rowcol-1)/2; row=startrow; matrix[row][col]=startdigit; for(i=startdigit+1;i<=startdigit+rowcol*rowcol-1;i++) { if(row-1<startrow) { row=startrow+rowcol-1; //printf("test1/n"); }else { row--; //printf("test2/n"); } if(col+1>startcol+rowcol-1) { col=startcol; //printf("test3/n"); } else{ col++; } if(!(matrix[row][col])) { row0=row; col0=col; matrix[row][col]=i; //printf("%d/n",i); //printf("test4/n"); }else { if(row0+1>startrow+rowcol-1) { row=startrow; //printf("test5/n"); }else { row=row0+1; //printf("test6/n"); } col=col0; matrix[row][col]=i; //printf("%d/n",i); } } } void odd_even(unsigned char rowcol) { unsigned char k=rowcol/2; unsigned char row,col; unsigned char t,l; odd(1,0,0,k); odd(k*k+1,k,k,k); odd(2*k*k+1,0,k,k); odd(3*k*k+1,k,0,k); print(rowcol); for(col=0;col<=k/2-1;col++) { for(row=0;row<k;row++) { if(!(col==0&&row==k/2)) { t=matrix[row+k][col]; matrix[row+k][col]=matrix[row][col]; matrix[row][col]=t; } } } t=matrix[(k-1)/2+k][(k-1)/2]; matrix[(k-1)/2+k][(k-1)/2]=matrix[(k-1)/2][(k-1)/2]; matrix[(k-1)/2][(k-1)/2]=t; l=(k-1)/2; for(col=0;col<l-1;col++) { for(row=0;row<k;row++) { t=matrix[row+k][rowcol-col-1]; matrix[row+k][rowcol-col-1]=matrix[row][rowcol-col-1]; matrix[row][rowcol-col-1]=t; } } } void even_even(unsigned char rowcol) { unsigned order1=rowcol*rowcol; unsigned order2=1; unsigned char row,col,k; for(row=0;row<rowcol;row++) { for(col=0;col<rowcol;col++) { if((row%4==col%4)||(row%4+col%4==3)) { matrix[row][col]=order1--; order2++; } else{ matrix[row][col]=order2++; order1--; } } } } /* void print(unsigned char rowcol) { char *str="/xc4/xc4/xc4/xc4"; unsigned char i,j; for(i=0;i<rowcol;i++) { for(j=0;j<rowcol;j++) { if(i==0) { if(j==0) { printf("/xda"); } else { printf("/xc2"); } } else { if(j==0) { printf("/xc3"); } else { printf("/xc5"); } } printf("%s",str); } if(i==0) { printf("/xbf/n"); } else { printf("/xb4/n"); } for(j=0;j<rowcol;j++) { printf("/xb3%4d",matrix[i][j]); } printf("/xb3/n"); } printf("/xc0%s",str); for(j=0;j<rowcol;j++) { printf("/xc1%s",str); } printf("/xd9/n"); } */ void print(unsigned char rowcol) { unsigned i; unsigned j; printf("/n"); for(i=0;i<rowcol;i++) { for(j=0;j<rowcol;j++) { printf("%5d",matrix[i][j]); } printf("/n"); } } void verify(unsigned char rowcol) { unsigned char i,j; unsigned sum=0,sum1=0; for(i=0;i<rowcol;i++) { sum1=0; for(j=0;j<rowcol;j++) { sum1+=matrix[i][j]; } if(sum==0) { sum=sum1; } else if(sum1!=sum) { printf("Error!/n"); return; } } for(i=0;i<rowcol;i++) { sum1=0; for(j=0;j<rowcol;j++) { sum1+=matrix[j][i]; } if(sum==0) { sum=sum1; } else if(sum1!=sum) { printf("Error!/n"); return; } } printf("Verity OK! the sum is %d/n",sum); } /* this is a test,test the number if is right; void print() { unsigned i; unsigned j; for(i=0;i<5;i++) { for(j=0;j<5;j++) { printf("%3d",matrix[i][j]); } printf("/n"); } } */ int main() { unsigned char number; printf("Please input a number: "); scanf ("%d",&number); init(); if(number%2) { odd(1,0,0,number); } else if(number%4) { odd_even(number); } else { even_even(number); } print(number); //print(); verify(number); return 0; }