魔方

魔方是一个从1到n2个整数的n×n矩阵，一个从1到n2个整数的n×n魔方中的每一行、每一列及对角线上的整数的和都是相同的。下面是一个n=5的魔方(和是75):

15	8	1	24	17
16	14	7	5	23
22	20	13	6	4
3	21	19	12	10
9	2	25	18	11

当n为奇数时，H. Coxeter给出了一个生成魔方的简单规则：

、在最上一行的中间位置填入1，然后不断寻找下一个位置，顺次填入2，3，…。

2、寻找下一个位置的方向是左上；即如当前位置是i,j,则下一个位置就是i-1,j-1。

2、把平面看成是环状的封闭体（即第一行的上方是最后一行，第一列的左方是最右一列）；

3、当前位置是i,j，如下一个预填入的位置不空，则以i+1,j作为下一个预填入的位置。

void magic(int n)

// create a magic square of size n, n is odd

{

   const int MaxSize = 51; // maximum square size

   int square[MaxSize][MaxSize], k, l;

    // check correctness of n

   if((n > (MaxSize)) || (n < 1)){

      cerr << “Error!…n out of range” << range” << endl;

      return;

   }

   else if(!(n%2)){

      cerr << “Error!…n is even” << range” << endl;

      return;

   }

   // n is odd. Coxeter’s rule can be used

   for(int i = 0; i < n; i++) // initialize square to 0

       for(int j = 0; j < n; j++)

           square[i][j] = 0;

    square[0][(n-1)/2] = 1; // middle of first row

    // i and j are current position

   int key = 2; i = 0; int j = (n-1)/2;

while (key <= n * n){

   // move up and left

      if (i – 1 < 0) k = n –1; else k = i –1;

      if (j – 1 < 0) l = n –1; else l = j –1;

      if (square[k][l]) i = (i+1)%n; // square occupied, move down

      else{ // square[k][l] is unoccupied

           i = k;

           j = l;

      }

      square[i][j] = key;

      key++;

   } // end of while

    //output the magic square

   cout << “magic square of size” << n << endl;

   for(i = 0; i < n; i++){

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

            cout << square[i][j] << ““;

       cout << endl;

   }

}

Program ： Magic square