OFFSET
0,3
COMMENTS
Euler transform of the generalized octagonal numbers (A001082).
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..5000
N. J. A. Sloane, Transforms
Eric Weisstein's World of Mathematics, Octagonal Number
FORMULA
G.f.: Product_{k>=1} 1/(1 - x^k)^A001082(k+1).
a(n) ~ exp(Pi * 2^(3/2) * n^(3/4) / (3*5^(1/4)) + 3*Zeta(3) * sqrt(5*n) / (2*Pi^2) - (45*Zeta(3)^2 / Pi^5 + Pi/6) * 5^(1/4) * (n^(1/4) / 2^(5/2)) + 225 * Zeta(3)^3 / (4*Pi^8) - Zeta(3) / (32*Pi^2) + 1/8) * Pi^(1/8) / (A^(3/2) * 2^(77/48) * 5^(5/32) * n^(21/32)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Nov 07 2017
MATHEMATICA
nmax = 32; CoefficientList[Series[Product[1/((1 - x^(2 k - 1))^(k (3 k - 2)) (1 - x^(2 k))^(k (3 k + 2))), {k, 1, nmax}], {x, 0, nmax}], x]
a[n_] := a[n] = If[n == 0, 1, Sum[Sum[d (d^2 + d - Ceiling[d/2]^2), {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[a[n], {n, 0, 32}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 06 2017
STATUS
approved