n 阶魔方阵超简单算法
#include<stdio.h>#define M 507
main()
{
unsigned int a[M][M]={0},H[M]={0},V[M]={0};
unsigned int i,j,k,n;
unsigned int q,p;
p=0;q=0;
i=0;j=M/2;a[i][j]=1;//第一行中间列为1
for(k=2;k<=M*M;k++)
{
i=(i+(M-1))%M;j=(j+1)%M;//从2开始后一个数在前一个数右上方，超过则回绕
if(a[i][j]==0)//初始状态（0状态）
a[i][j]=k;
else//非0状态，填到该数字正下方
{
i=(i+1)%M+1;
j=(j+(M-1))%M;
a[i][j]=k;
}
}
printf("%d阶魔方阵为：\n",M);
for(i=0;i<M;i++)
for(j=0;j<M;j++)
{
printf("%5d",a[i][j]);
if((j+1)%M==0) printf("\n\n");
}
printf("\n");
/////////////////验证
printf("%d阶魔方验证：\n",M);
for(n=0;n<M;n++)
for(j=0;j<M;j++)
{
V[n]=V[n]+a[n][j];//每行之和
H[n]=H[n]+a[j][n];//每列只和
}
for(i=0;i<M;i++)
{ p=p+a[i][i];//对角线只和
q=q+a[M-1-i][i];
printf("V[%d]=%d,H[%d]=%d\n",i,V[i],i,H[i]);
}
printf("p=%d\t,q=%d",p,q);
printf("\n");
}
