（Python）
导入itertools
导入数学
定义a（n）：
res=0
对于itertools.permutation中的p（[i for i in range（1，2*n+1）]，n）：
i=0
s=0
对于范围（1,2*n+1）中的j：
如果j不在p中：
s+=p[i]**j
i+=1
如果math.isqrt**2==s：
res=最大值（res，s）
返回res
（Python）
从itertools导入组合、排列
从sympy.theory.primetest导入为平方
a、 m，f=集合（范围（1，2*n+1）），0，[[b**c代表范围（2*n+1]）中的c]代表范围（2%n+1）中的b]
对于组合（a，n）中的b：
clist=排序（a集（b））
对于置换中的d（范围（n））：
k=总和（f[b[i]][clist[d[i]]）
如果k>m且is_square（k）：
m=k
