OFFSET
0,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: Product_{n>0} (1-x^n)^6/(1-x^(6*n)).
Euler transform of period 6 sequence [ -6, -6, -6, -6, -6, -5, ...].
Empirical: Sum_{n>=0} a(n) / exp(n*Pi) = (1/32) * Pi^(5/4) * 2^(5/6) * (sqrt(2) * (3^(1/2)-1))^(3/2) * (sqrt(2) * (1+3^(1/2)))^(5/3) / Gamma(3/4)^2 / Gamma(11/12)^(3/2) / Gamma(7/12)^(3/2) = A389010. - Simon Plouffe, Sep 22 2025
MAPLE
with(numtheory):
a:= proc(n) option remember; `if`(n=0, 1, add(add(d*
`if`(irem(d, 6)=0, -5, -6), d=divisors(j))*a(n-j), j=1..n)/n)
end:
seq(a(n), n=0..70); # Alois P. Heinz, Jan 07 2017
MATHEMATICA
QP = QPochhammer; QP[x]^6/QP[x^6] + O[x]^70 // CoefficientList[#, x]& (* Jean-François Alcover, Mar 25 2017 *)
PROG
(PARI) q='q+O('q^66); Vec( eta(q)^6/eta(q^6) ) \\ Joerg Arndt, Mar 25 2017
CROSSREFS
KEYWORD
sign
AUTHOR
Seiichi Manyama, Jan 07 2017
STATUS
approved
