魔术方阵 的算法

以N(没有偶阶魔方)阶魔方为例,首先第一行的中间数为1,下一个数

如果不是N的倍数加1,列加1,行减1(如果列大于N,则列减N,如果行小于0则行加N)
如果是N的倍数加1,则行加1(如果行大于N,则行减N)

程序如下：
#include "stdio.h"
void magicsquare(int n)
{
int i,j,k,a[16][16];
for(i=0;i<16;i++) //初始化方阵
for(j=0;j<16;j++)
a[i][j]=0;

i=0;j=(n-1)/2; //第一行中间的数为1
a[i][j]=1;

for(k=2;k<=n*n;k++)
{
if((k-1)%n==0) //是N的倍数加1
{
i++; //行加1
if(i>n-1)
i-=n;
a[i][j]=k;
}
else //不是N的倍数加1
{
i--;j++; //行减1,列加1
if(i<0)
i+=n;
if(j>n-1)j-=n;
a[i][j]=k;
}
}
for(i=0;i{
for(j=0;jprintf("%5d",a[i][j]);
printf("\n");
}
}

int main()
{
int n;
printf("Enter n:");
while(scanf("%d",&n)&&n!=0)
if(n%2) //没有偶价魔方
magicsquare(n);
return 0;
}