最简单的奇数阶幻方（魔方）算法详解

幻方，就是一个N*N的矩阵中填入1至N^2,使得每行、每列，以及两条对角线上的数之和都相等。当N为奇数时，可用以下的方法来填数。
（1）将1填入第一行的中间；
（2）将数字下一个数填在右上方。
    （i)出现越界。若越过第一行，则折返回最后一行。若越过最后一列，则折返回第一列。
    （ii）没越界，但右上方已有数了，此时将其填在上一个数的正下方。
（3）对于矩阵的最右上方的格子（第一行最后一列）要特殊处理，因为下一个数同时越过第一行和最后一列，折返时的格子也已经有数了。因此特殊处理，将其填在第二行的最后一列即可。
具体参考程序如下：

#include<stdio.h>
int i,j,n,x,y;
int a[101][101];
int main()
{
 scanf("%d",&n);
 for(i=0;i<=n+1;i++)
 {
  for(j=0;j<=n+1;j++) 
  if(1<=i&&i<=n&&1<=j&&j<=n) a[i][j]=0;   
   else a[i][j]=1;
 }
 x=1;y=n/2+1;
 a[x][y]=1;   //将数字1填在第一行的中间处
 for(i=2;i<=n*n;i++)     //待填数字i  
  {
  if(x==1&&y==n) x=2;   //右上方格子要特殊处理
  else
  {
        x--;y++;        //将数字i填在右上方
        if(x<1) x=n;    //越界判断
        if(y>n) y=1;    //越界判断
     if(a[x][y]>0) {x+=2;y--;}  //填在上一个数的正下方
  }
     a[x][y]=i;   
 }
 for(i=1;i<=n;i++)   //输出幻方阵
 {
  for(j=1;j<n;j++) printf("%d ",a[i][j]);
  printf("%d\n",a[i][n]);
 }
 return 0;
}