说起幻方，大家应该在小学时候就已经接触过了，最简单的就是九宫格，射雕英雄传中的那段至今还记得：戴九履一，左三右七，二四为肩，六八为足。下面我们就来看看这个有趣的问题。
幻方可以分为：奇数阶幻方、双偶阶幻方、单偶阶幻方。
奇数阶幻方
上面所说的九宫格就是典型的奇数阶幻方，奇数阶幻方值得是阶数为奇数的幻方。其最经典的填法是罗伯法。首先 把1（或最小的数）放在第一行正中；按以下规律排列剩下的(n×n－1)个数，具体步骤为：
（1）每一个数放在前一个数的右上一格；
（2）如果这个数所要放的格已经超出了顶行那么就把它放在底行，仍然要放在右一列；
（3）如果这个数所要放的格已经超出了最右列那么就把它放在最左列，仍然要放在上一行；
（4）如果这个数所要放的格已经超出了顶行且超出了最右列，那么就把它放在底行且最左列；
（5）如果这个数所要放的格已经有数填入，那么就把它放在前一个数的下一行同一列的格内。
上述步骤可以总结为七言绝句：
奇幻七绝
先填上行正中央，
依次斜填切莫忘。
上格没有顶格填，
顶格没有底格放。
下面有人通过作图可以很好的解释这几句话，现借鉴如下：
从上面的图可以看出，该图与我们前面的九宫格口诀不相符，上下颠倒了。但是这都是对的，本质上没有区别。
双偶数阶幻方
所谓双偶阶幻方就是当n可以被4整除时的偶阶幻方，即4K阶幻方。其最经典的填法为海尔法，下面以8阶幻方为例，具体的填法为：
（1）先把数字按顺序填。然后，按4×4把它分割成4块（如图）
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
|
9 |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
|
17 |
18 |
19 |
20 |
21 |
22 |
23 |
24 |
|
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
|
33 |
34 |
35 |
36 |
37 |
38 |
39 |
40 |
|
41 |
42 |
43 |
44 |
45 |
46 |
47 |
48 |
|
49 |
50 |
51 |
52 |
53 |
54 |
55 |
56 |
|
57 |
58 |
59 |
60 |
61 |
62 |
63 |
64 |
（2）每个小方阵对角线上的数字（如左上角小方阵部分），换成和它互补的数。
|
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 |
单偶数阶幻方
所谓单偶阶幻方就是当n不可以被4整除时的偶阶幻方，即4K+2阶幻方。如(n=6，10，……)的幻方。其经典的填法为斯特拉兹法，以10阶幻方为例，具体的步骤如下：
1）把魔方阵分为A，B，C，D四个象限，这样每一个象限肯定是奇数阶。用罗伯法，依次在A象限，D象限，B象限，C象限按奇数阶幻方的填法填数。
（2）在A象限的中间行、中间格开始，按自左向右的方向，标出k格。A象限的其它行则标出最左边的k格。将这些格，和C象限相对位置上的数互换位置。
（3）在B象限所有行的中间格，自右向左，标出k－1格。(注：6阶幻方由于k－1=0，所以不用再作B、D象限的数据交换)，将这些格，和D象限相对位置上的数互换位置。
具体的程序实现如下：
#include<stdio.h>
//注意由于matrix大小(可以更改)给定，能最大生成10阶幻方
bool check(int matrix[10][10],int n)//判断是否为幻方
{
int sum=0;
int temp=0;
int i=0,j=0,k=0;
for( i=0;i<n;i++)
sum=sum+matrix[0][i];//得到一行或列的总和
for(j=1;j<n;j++)//检查行
{
temp=0;
for(k=0;k<n;k++)
temp=temp+matrix[j][k];
if(temp!=sum)
return false;
}
for( j=0;j<n;j++)//检查列
{
temp=0;
for(k=0;k<n;k++)
temp=temp+matrix[k][j];
if(temp!=sum)
return false;
}
temp=0;
for(i=0;i<n;i++)
temp=temp+matrix[i][i];//检查主对角线
if(temp!=sum)
return false;
temp=0;
for(i=0;i<n;i++)
temp=temp+matrix[i][n-1-i];//检测副对角线
if(temp!=sum)
return false;
printf("该方阵为幻方！\n");
return true;
}
void Odd(int n,int matrix[10][10])//奇数阶幻方
{
int i=0,j=n/2;
int number=1;
for(int k=0;k<n*n;k++)
{
matrix[i][j]=number;
i--;
j++;
number++;
if(i<0&&j<n)//出上界
{
i=n-1;
}
else if(i>=0&&j>=n)//出右界
{
j=0;
}
else if(i<0&&j>=n)//右、上出界
{
if(matrix[n-1][0]!=0)//底格放
{
i=i+2;
j=j-1;
}
else
{
i=n-1;
j=0;
}
}
if(matrix[i][j]!=0)//底格放
{
i=i+2;
j=j-1;
}
}
}
void DoubleEven(int n,int matrix[10][10])//双偶数阶幻方
{
int number=1;
int temp=0;
int i=0,j=0,k=0;
for(i=0;i<n;i++)
for(j=0;j<n;j++)
matrix[i][j]=number++;//初始化
for(i=0;i<n;i=i+4)
for(j=0;j<n;j=j+4)
for(k=0;k<4;k++)
{
matrix[i+k][j+k]=n*n+1-matrix[i+k][j+k];//每个对角线的数换成互补的数
matrix[i+k][j+3-k]=n*n+1-matrix[i+k][j+3-k];
}
}
void SingleEven(int n,int matrix[10][10])//单偶数阶幻方
{
int degree=n/2;
int flag=n/4;
int i=0,j=0,k=0;
int temp=0;
int matrix1[10][10]={0};//将大矩阵化为A B C D四个小矩阵
int matrix2[10][10]={0};
int matrix3[10][10]={0};
int matrix4[10][10]={0};
Odd(degree,matrix1);//对每一个矩阵进行奇数幻方算法
for(i=0;i<degree;i++)
for(j=0;j<degree;j++)
{
matrix2[i][j]=matrix1[i][j]+degree*degree;
matrix3[i][j]=matrix1[i][j]+degree*degree*2;
matrix4[i][j]=matrix1[i][j]+degree*degree*3;
}
for(i=0;i<degree;i++)//对A C矩阵按照规则进行数据交换
for(j=0;j<flag;j++)
if(i!=(degree/2))
{
temp=matrix1[i][j];
matrix1[i][j]=matrix4[i][j];
matrix4[i][j]=temp;
}
else
{
temp=matrix1[i][j+degree/2];
matrix1[i][j+degree/2]=matrix4[i][j+degree/2];
matrix4[i][j+degree/2]=temp;
}
for(i=0;i<degree;i++)//对B D矩阵按照规则进行数据交换
for(j=0;j<flag-1;j++)
{
temp=matrix2[i][j+degree/2];
matrix2[i][j+degree/2]=matrix3[i][j+degree/2];
matrix3[i][j+degree/2]=temp;
}
//将新的四个矩阵赋给幻方矩阵matrix
for(i=0;i<degree;i++)
{
for(j=0;j<degree;j++)
matrix[i][j]=matrix1[i][j];
for(k=0;k<degree;k++)
matrix[i][j+k]=matrix3[i][k];
}
for(i=0;i<degree;i++)
{
for(j=0;j<degree;j++)
matrix[i+degree][j]=matrix4[i][j];
for(k=0;k<degree;k++)
matrix[i+degree][j+k]=matrix2[i][k];
}
}
void main()
{
int matrix[10][10]={0};
int n; printf("%d",6%2);
printf("请输入幻方的阶数：");
scanf("%d",&n);
if(n%2!=0)
Odd(n,matrix);
else if(n%4!=0)
SingleEven(n,matrix);
else
DoubleEven(n,matrix);
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
printf("%d ",matrix[i][j]);
printf("\n");
}
check(matrix,n);//检测是否为幻方
}
注：如果程序出错，可能是使用的开发平台版本不同，请点击如下链接： 解释说明
原文：http://blog.csdn.net/tengweitw/article/details/17687377
作者：nineheadedbird