原文链接: 求解四阶幻方
上一篇: 汇编实验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 =