|
PROG
|
(Sage)
# see Python in A159842 for the definition of dc, fin, per, u, N, N2
def fin_d(d):
return fin(*(d.get(n+1, 0) for n in range(max(d))))
def a(n): # see Hanany & Seong 2011, Table 1 row D=5 and Table 9
return (dc(u, N, N2, lambda n: n**3)(n) +
10 * dc(u, u, N, N2, fin(1, -1, 0, 8))(n) +
15 * dc(u, u, N, N, fin_d({1: 1, 2: -3, 4: 14, 8: -12, 16: 16}))(n) +
20 * dc(u, u, N, per(0, 1, -1), fin(1, 0, -1, 0, 0, 0, 0, 0, 9))(n) +
20 * dc(u, u, u, per(0, 1, -1), fin(1, -1, 0, 2), fin(1, 0, -1, 0, 0, 0, 0, 0, 3))(n) +
30 * dc(u, u, u, per(0, 1, 0, -1), fin_d({1: 1, 2: -2, 4: 3, 16: 6, 32: -8, 64: 8}))(n) +
24 * dc(u, per(0, 1, -1, -1, 1), per(0, 1, I, -I, -1), per(0, 1, -I, I, -1))(n)) / 120
print([a(n) for n in range(1, 100)])
|