魔方矩阵问题

输出一个N*N的矩阵，矩阵的数字从1~N2，其每一行、每一列和对角线的各数之和均相等。这样的矩阵成为魔方矩阵。

分析：

要求每一行、每一列和对角线的各数之和均相等，很显然N行都相等，则每一行的和应该等于(1+2+3+…+N2)/N。

魔方阵计算规律(行，列以1开始)：

1.将1放在第一行，中间一列；

2.从2开始至N*N各数按如下规律： 每一个数存放的行比上一个数的行减1；每一个数存放的列比上一个数的列加1；

3.当一个数行为1，下一个数行为N；

4.当一个数列数为N,下一个数列数为1，行数减1；

5.若按上述规则确定的位置有数字，或上一个数位第1行第N列，下一个数字位置为上一个数的正下方（即行数减1，列数不变）。

例如N等于5时的数字排列如下图：

 

代码如下：

#include <stdio.h>

#define N 5

int main()

{

         int a[N][N] = {0};

         int i,j,count = 1;

         int row,col;

         row = 0;

         col = N / 2;

         while (count <= N*N)

         {

                   a[row][col] = count;

                   i = row;

                   j = col;

                   if (i == 0)

                            i = N - 1;

                   else

                            i = i - 1;

                   j = (j + 1) % N;

                   if (a[i][j]!=0||(row==0&&col==N-1))

                   {

                            i = row + 1;

                            j = col;

                   }

                   row = i;

                   col = j;

                   count++;

         }

         for (i = 0; i < N; i++)

         {

                   for (j = 0; j < N; j++)

                            printf("%3d",a[i][j]);

                   printf("\n");

         }

}

运行结果（N=5、7、9）：

       

 

 

 

进一步分析：

 

观察n和行列的关系，可以发现，当每N一组时，每组的第一个数字都是放在前面一个数字的下方，如：6、11、16、21，修改程序如下：

//文件:魔方矩阵.c

//作者:Ding Yatao

//日期:2020年12月

#include <stdio.h>

#define N 5

#define HDirection       -1    //1 下一个数在右侧 ，-1 下一个数在左侧

int main()

{

         int a[N][N] = {0};

         int i,row=0,col=N/2;

         for(i=1;i<=N*N;i++)

         {

                   a[row][col] = i;

                   if((i+1)%N==1)

                            row=(row+1)%N;

                   else

                   {

                            row=(row-1+N)%N;     //加N后对N求余数，可以自动修正超出范围的情况

                            col=(col+HDirection+N)%N; // HDirection用于列数的变化方向

                   }

         }

         for (row = 0; row < N; row++)

         {

                   for (col= 0; col < N; col++)

                            printf("%3d",a[row][col]);

                   printf("\n");

         }

}

当HDirection等于1，输出结果同上，当HDirection等于-1，相当于上图的水平翻转，如下图：

   

          同样，读者可以考虑垂直翻转图形的输出。