OFFSET
0,3
COMMENTS
The sum of the first n cubes equals the square of the n-th triangular number: (1^3 + 2^3 + 3^3 + ... + n^3) = (1 + 2 + 3 + ... + n)^2 = (n*(n+1)/2)^2.
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..100
G. H. Hardy and S. Ramanujan, Asymptotic formulae in combinatory analysis, Proceedings of the London Mathematical Society, 2, XVI, 1917, p. 373.
FORMULA
a(n) equals the coefficient of x^A000537(n) in Product_{k>=1} 1/(1 - x^(k^3)) for n >= 0, where A000537(n) = (n*(n+1)/2)^2.
Limit a(n)^(1/n) = exp( (Gamma(1/3) * Zeta(4/3))^(3/4) * (2/3)^(3/2) ) = 19.6784468967446... - see formula by Vaclav Kotesovec in A003108.
a(n) ~ (Gamma(1/3)*zeta(4/3))^(3/4) * exp(sqrt(2) * 3^(-3/2) * (Gamma(1/3)*zeta(4/3))^(3/4) * (2*n+1)) / (3 * sqrt(2) * Pi^2 * n^5). - Vaclav Kotesovec, Jan 25 2024
MAPLE
b:= proc(n, i) option remember; `if`(n=0 or i=1, 1,
b(n, i-1) +`if`(i^3>n, 0, b(n-i^3, i)))
end:
a:= n-> (t-> b(t, iroot(t, 3)))((n*(n+1)/2)^2):
seq(a(n), n=0..23); # Alois P. Heinz, Jan 24 2024
MATHEMATICA
Table[SeriesCoefficient[Product[1/(1 - x^(k^3)), {k, 1, (n*(n + 1)/2)^(2/3) + 1}], {x, 0, (n*(n + 1)/2)^2}], {n, 0, 20}] (* Vaclav Kotesovec, Jan 24 2024 *)
PROG
(PARI) {a(n) = my(M=(n*(n+1)/2)^2); polcoeff(1/prod(k=1, ceil(M^(1/3)), 1-x^(k^3)+x*O(x^M)), M)}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 23 2024
STATUS
approved