求解四阶幻方

最新推荐文章于 2023-02-13 11:26:41 发布

阿豪boy

最新推荐文章于 2023-02-13 11:26:41 发布

阅读量1k

点赞数

文章标签：剪枝算法 dfs python 数据结构

本文链接：https://blog.csdn.net/qq_35516360/article/details/122066544

版权

本文探讨了解决四阶幻方问题的算法，通过剪枝技巧减少搜索量，介绍了幻方的特殊性质，并对比了C、Python、C++、Java和JavaScript在实现该算法时的性能差异，强调了Java和C++的高效性。

摘要由CSDN通过智能技术生成

原文链接: 求解四阶幻方

上一篇: 汇编实验4.3

下一篇: 汇编实验 4.2

由于 16！ = 20922789888000 规模在10^13次方级别，需要进行大量剪枝操作

利用四阶幻方的特殊性质可以进行剪枝

幻方的一般性质为：幻方每一行之和、每一列之和、两条对角线之和都相等，都等于幻和(四阶幻和为34)。

四阶幻方还有一些性质：中心四个数之和、四个顶角数字之和、边上两对行中间四个数之和、边上两对列中间四个数之和也都等于幻和(34)。

 0  1  2  3
 4  5  6  7
 8  9 10 11
12 13 14 15

12 13  3  6
 1  4 14 15
16  9  7  2
 5  8 10 11


相比三阶幻方，四阶幻方多了下面几个约束条件
1+2+13+14
4+7+8+11
5+6+9+10
0+3+12+15

利用上面性质对排列局面剪枝，可以极大减少搜索量。

另外，每搜索到一个幻方，就可以得到与其旋转反射相同的8个幻方，搜索量可以缩减为1/8

C语言跑了十几秒

#include <stdio.h>
int a[17], b[17], m;
void s(int i){/*四阶幻方全解搜索程序，C代码,运行时间7秒*/
    int n = 0, j = 0;
    while (++j < 17)
        if (!a[j]){
            a[b[i] = j] = 1;
            switch (i){
            case 1: case 2: case 3: case 5: case 6: case 7: case 9: case 10: s(i + 1); break;
            case 11: if (b[6] + b[7] + b[10] + b[11] == 34) s(12); break;
            case 4: case 8: case 12: if (b[i - 3] + b[i - 2] + b[i - 1] + b[i] == 34) s(i + 1); break;
            case 13: if (b[1] + b[5] + b[9] + b[13] == 34 && b[4] + b[7] + b[10] + b[13] == 34) s(14); break;
            case 14: case 15: if (b[i - 12] + b[i - 8] + b[i - 4] + b[i] == 34) s(i + 1); break;
            case 16: for (printf("\n"), ++m; ++n < 17; n % 4 ? 0 : printf("\n") ) printf("%2d ", b[n]);
            }
            a[j] = 0;
        }
}
int main(void){
    s(1);
    printf("四阶幻方总数:%d个(含旋转反射相同)", m);
    return 0;
}

四级幻方共7040种

Python跑了30分钟。。。

'''
三阶幻方8个
四阶幻方7040个
'''

# 幻方阶数
n = 4
# 幻和，即每行，每列和对角线的和
s = n * n * (n * n + 1) // 2 // n
# 数组
a = []
# 使用过的数字
used = set()
# 答案数组
ans = []


# 检查是否满足幻方条件
def check():
    flag = True
    for i in range(n):
        #  行的幻和
        r = sum([a[i * n + j] for j in range(n)])

        # 列的幻和
        c = sum([a[i + j * n] for j in range(n)])
        flag = flag and r == c == s
    # 主对角线和次对角线的幻和
    s1 = sum([a[i * n + i] for i in range(n)])
    s2 =