### Python 实现 4n+2 阶偶数幻方算法 以下是基于引用内容[^2]中描述的 (4M+2) 阶魔方阵算法实现的一个完整的 Python 程序： #### 步骤解析 1. 构造一个 \( n \times n \) 的基础魔方阵，其元素为等差数列 {4k | k = 0, 1, ..., n²-1}。 2. 将该基础魔方阵扩展为 \( 2n \times 2n \)，通过复制每个单元格形成 \( 2 \times 2 \) 子矩阵。 3. 对扩展后的矩阵进行上下和左右对称扩展操作。 4. 最终将扩展后的两个矩阵叠加，生成目标 \( 4n+2 \) 阶幻方。 下面是具体的代码实现： ```python import numpy as np def generate_base_magic_square(n): """生成基础 n x n 魔方阵""" magic_square = np.zeros((n, n), dtype=int) num_list = [4 * k for k in range(n ** 2)] index = 0 for row in range(n): for col in range(n): magic_square[row][col] = num_list[index] index += 1 # 调整使得每行、列、对角线之和满足条件 total_sum = 2 * n ** 3 - 2 * n correction_matrix = np.full((n, n), int(total_sum / n)) magic_square += correction_matrix return magic_square def expand_to_double_size(magic_square): """将 n x n 扩展为 2n x 2n""" n = magic_square.shape[0] expanded_square = np.block([[magic_square[i//2,j//2] for j in range(2*n)] for i in range(2*n)]) return expanded_square def add_symmetric_extension(matrix): """添加上下和左右对称扩展""" top_extension = matrix[::-1, :] # 上下翻转作为顶部扩展 bottom_part = np.vstack([matrix, top_extension]) left_extension = bottom_part[:, ::-1] # 左右翻转作为左侧扩展 final_matrix = np.hstack([left_extension, bottom_part]) return final_matrix def combine_matrices(base_matrix, extension_matrix): """组合基础矩阵与扩展矩阵""" combined_matrix = base_matrix + extension_matrix return combined_matrix def create_even_order_magic_square(order): """创建指定阶次的偶数幻方""" if order % 2 != 0 or order < 6: raise ValueError("Order must be an even number greater than or equal to 6.") n = order // 2 base_square = generate_base_magic_square(n) double_expanded_square = expand_to_double_size(base_square) symmetric_extended_square = add_symmetric_extension(double_expanded_square) final_magic_square = combine_matrices(double_expanded_square, symmetric_extended_square) return final_magic_square[:order, :order] # 测试程序 if __name__ == "__main__": N = 10 # 创建 10x10 幻方（即 4*2+2） magic_square = create_even_order_magic_square(N) print(f"{N}-阶幻方:") print(magic_square) # 验证每一行、列以及对角线之和是否相同 expected_sum = sum(range(1, N*N + 1)) // N rows_sums = np.sum(magic_square, axis=1) cols_sums = np.sum(magic_square, axis=0) diag_sum1 = np.trace(magic_square) diag_sum2 = np.trace(np.fliplr(magic_square)) is_valid = all(s == expected_sum for s in [*rows_sums, *cols_sums, diag_sum1, diag_sum2]) print(f"\n验证结果: {'有效' if is_valid else '无效'}") ``` --- #### 关键点解释 1. `generate_base_magic_square` 函数用于构建初始的基础魔方阵，并调整使其满足特定和的要求。 2. `expand_to_double_size` 函数负责将 \( n \times n \) 的矩阵扩展成 \( 2n \times 2n \) 形式的子矩阵结构。 3. `add_symmetric_extension` 添加必要的上下和左右对称扩展部分。 4. `combine_matrices` 组合多个阶段的结果以完成最终的幻方构造。 ---