求N奇数阶幻方

1. 如果矩阵满足条件，那么对任意,也满足条件。
证明显然。

设为奇数，我们现在构造一个n阶幻方包含0到所有数

这里x,y满足同余式

待确定。

由于该方程组的系数矩阵的行列式为1，所以对任意i,j有唯一解。我们接下来确定a,b：

首先验证每行每列的和均相等，即

由于对任意i,，当x取遍模n的剩余类时，也会取遍所有的剩余类

故每行的和为：
同理验证每列。

对于对角线：令b = n - 1, a = (n-1)/2，可以验证y = (n - 1)/2, x = i
所以其和为：n*(0 + ... n - 1) + n(n - 1)/2 = n(n-1)(n+1)/2 = n(n^2 - 1)/2

偶数不成立是因为， (n-1)/2不能整除。
==========
这种用同余式构造组合结构的方法很常见，比如拉丁方的构造、n皇后问题的构造等。

作者：Changlong Wu
链接：https://www.zhihu.com/question/30498489/answer/48393702
来源：知乎
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。

 

import java.util.Scanner;
//奇数阶幻方的实现
public class Main {
    //n为用户输入的n阶奇数幻方   n为奇数
    public static int[][] magicOdd(int n) {
        //构造一个（n+2）*（n+2）矩阵
        int[][] square = new int[n + 1][n + 1];
        int i = 0;
        int j = (n + 1) / 2;
        //从第一行的中间那个数字（是1）开始填幻方
        //n阶幻方一共有n*n个数字（从1~n*n）
        //奇数阶幻方的实现算法
        for (int key = 1; key <= n * n; key++) {
            if ((key % n) == 1)
                i++;
            else {      // //填充当前数的右上角那个数
                i--;
                j++;
            }
            if (i == 0) {      ////判断条件：若是在(n+2)*(n+2)阶方阵的第一行
                i = n;
            }
            if (j > n) {
                j = 1;
            }
            square[i][j] = key;
        }

        //对(n+2)*(n+2)阶的方阵进行筛选出中间的n*n阶幻方

        int[][] matrix = new int[n][n];

        for (int k = 0; k < matrix.length; k++) {

            for (int l = 0; l < matrix[0].length; l++) {
                matrix[k][l] = square[k + 1][l + 1];
            }
        }
        return matrix;
    }
    public static void main(String[] args) {
        Scanner a = new Scanner(System.in);
        int b = a.nextInt();
        //b为用户输入的奇数
        System.out.println();
        int[][] magic = Main.magicOdd(b);
        for (int k = 0; k < magic.length; k++) {
            for (int l = 0; l < magic[0].length; l++) {
                System.out.print(magic[k][l] + "    ");
            }
            System.out.println();
        }
    }
}