扩展四阶幻方

hujunhua

用与前述相同的方法进行搜索，结果如下

1. e = 1;
   
2. octagon[abcd_List] := Complement[Range[23], e + (Plus @@@ Subsets@abcd)]
   
3. powersetQ[oct_List] := MemberQ[oct, oct[[1]] + oct[[2]]]
   
4. bcf[oct_List] := {oct[[1]], oct[[2]], 23 - oct[[1]] - oct[[2]]}
   
5. abcd = Select[IntegerPartitions[23 - 2 e, {4}], DuplicateFreeQ[Plus @@@ Subsets[#]] && powersetQ@octagon@# &&Length[Intersection[#, bcf@octagon@#]] > 1 &];
   
6. abcdef = {#, bcf@octagon@#} & /@ abcd

复制代码

输出结果：{{a,b,c,d},{b,c,f}}, 无序，b,c为前后组的交。

1. {{12, 6, 2, 1}, {5, 6, 12}},
   
2. {{12, 5, 3, 1}, {3, 8, 12}},
   
3. {{12, 4, 3, 2}, {2, 9, 12}},
   
4. {{11, 6, 3, 1}, {3, 6, 14}},
   
5. {{8, 7, 4, 2}, {2, 4, 17}}

复制代码

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

hujunhua

将上楼程序中的首行改成e=2就行了。输出格式同上。

1. {{12, 4, 2, 1}, {1, 10, 12}},
   
2. {{10, 6, 2, 1}, {1, 6, 16}},
   
3. {{8, 6, 4, 1}, {1, 4, 18}}

复制代码

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

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

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

hujunhua

统计：

        0在外围实为144解。
加上 0在中间4×4部分的288解，共432解。

hujunhua

以上各楼中，都有提到参数分组，有的是{a,b,c}{d}，有的是{a,d}{b,c},为什么会这样呢？
这是由内外参数的交决定的。

mathe

中心16格采用8#方案，而对于边上8个数字，只要满足上面俩数和=$\frac{\Sigma}2-c$, 左边俩数和=$\frac{\Sigma}2-b$，右边俩数和=$\frac{\Sigma}2+b$,下面俩数和=$\frac{\Sigma}2+c$即可。而每个满足条件的方案中，上下左右俩数各自可以任意交换位置，也就是一个结果就对应 16 组不同的解。
另外如果将中心16格划分为左上左下右上右下4个2×2方块，并且对角交换，

如图可以看出同样需要分别加入和标准值差b或c的值，也就是说同样的搜索结果还可以对应另外8组。
于是最后我们只要搜索满足条件的a,b,c,d,e以及四个边上较小的数字就可以了，其中每组对应16个结果。
附件中搜索结果为100组，每组给出了一个例子结果，所以总共结果应该1600组（已经去除对称情况）
e4j.zip (3.1 KB, 下载次数: 1)

2020-10-8 11:27 上传

点击文件名下载附件

hujunhua

mathe 发表于 2020-10-8 11:28
无额外约束的原扩展幻方的解数
中心16格采用8#方案，而对于边上8个数字，只要满足上面俩数和比标志值($\frac{\Sigma}2=a+b+c+d+2e$)小c,左 ...

对扩展部分增加的额外约束是要求如下图中底色相同的两格之和等于∑/2.
前文中有时称之为强约扩展幻方。前文中所有参数图和数字图都是强约扩展幻方。

增加扩展对径格约束后，扩展幻方多了以下特点：
1、3×3的倾斜45°正方形四角之和等于∑。（mathe在群中指出的）
2、扩展对角线，如上图中的22+7+5+12=∑。

这个原始扩展幻方的搜索结果好像与chyanog在群中透露的数字不同。
chyanog说不去重有117760解，那么去重后应该是920解，见6#。
不知道出入在哪里。

mathe

我统计了一下，扩展部分如同hujunhua这样邻边对偶的有52种，对边对偶的有48种。
如果继续考虑我找出的方案中一种只对应hujunhua方案中两种（可以同时将所有扩展部分数字两两交换），那么应该只对应52*2=104种，如果继续添加2e格子在中心的情况，那么有208种。好像数目偏少？

mathe

前面有bug,应该是230组不等价的解，满足扩展约束的108组。
这下好像和chyanog的匹配了。
每组中心16个格子内外置换可以乘上系数$2$，扩展部分交换位置又可以乘上$2^4$, 旋转和翻转又可以乘上$2^3$，得出结果怎么同chyanog正好差了一倍？还有哪种对称情况没有计数？

hujunhua

mathe 发表于 2020-10-8 18:12
前面有bug,应该是230组不等价的解，满足扩展约束的108组。
这下好像和chyanog的匹配了。
每组中心16个 ...

中间4×4的格子内外置换是5#所说的取景平移的一种（+(2,2)），取景平移还有上两行与下两行对易、左两列与右两列对易两种，所以还要乘以2——这就对上了。

hujunhua

==