login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A023872 Expansion of Product_{k>=1} (1 - x^k)^(-k^3). 20
1, 1, 9, 36, 136, 477, 1703, 5746, 19099, 61622, 195366, 607069, 1856516, 5586870, 16579850, 48549116, 140438966, 401592524, 1136121837, 3181700219, 8825733603, 24261363403, 66124058839, 178757752892, 479513547399, 1276792213203, 3375707760306, 8864712158225 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..5422 (first 1001 terms from Alois P. Heinz)
G. Almkvist, Asymptotic formulas and generalized Dedekind sums, Exper. Math., 7 (No. 4, 1998), pp. 343-359.
Vaclav Kotesovec, A method of finding the asymptotics of q-series based on the convolution of generating functions, arXiv:1509.08708 [math.CO], Sep 30 2015, p. 21.
FORMULA
a(n) ~ (3*Zeta(5))^(59/600) * exp(5 * n^(4/5) * (3*Zeta(5))^(1/5) / 2^(7/5) + Zeta'(-3)) / (2^(41/200) * n^(359/600) * sqrt(5*Pi)), where Zeta(5) = A013663 = 1.036927755143369926..., Zeta'(-3) = ((gamma + log(2*Pi) - 11/6)/30 - 3*Zeta'(4)/Pi^4)/4 = 0.00537857635777430114441697421... . - Vaclav Kotesovec, Feb 27 2015
G.f.: exp( Sum_{n>=1} sigma_4(n)*x^n/n ). - Seiichi Manyama, Mar 04 2017
a(n) = (1/n)*Sum_{k=1..n} sigma_4(k)*a(n-k). - Seiichi Manyama, Mar 04 2017
MAPLE
with(numtheory):
a:= proc(n) option remember; `if`(n=0, 1,
add(add(d*d^3, d=divisors(j)) *a(n-j), j=1..n)/n)
end:
seq(a(n), n=0..30); # Alois P. Heinz, Nov 02 2012
MATHEMATICA
max = 27; Series[ Product[ 1/(1-x^k)^k^3, {k, 1, max}], {x, 0, max}] // CoefficientList[#, x]& (* Jean-François Alcover, Mar 05 2013 *)
PROG
(PARI) m=30; x='x+O('x^m); Vec(prod(k=1, m, 1/(1-x^k)^k^3)) \\ G. C. Greubel, Oct 30 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); Coefficients(R! ( (&*[1/(1-x^k)^k^3: k in [1..m]]) )); // G. C. Greubel, Oct 30 2018
(SageMath) # uses[EulerTransform from A166861]
b = EulerTransform(lambda n: n^3)
print([b(n) for n in range(30)]) # Peter Luschny, Nov 11 2020
CROSSREFS
Column k=3 of A144048.
Sequence in context: A213283 A188162 A353389 * A034557 A285241 A231431
KEYWORD
nonn
AUTHOR
EXTENSIONS
Definition corrected by Franklin T. Adams-Watters and R. J. Mathar, Dec 04 2006
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 18 22:18 EDT 2024. Contains 371782 sequences. (Running on oeis4.)