OFFSET
0,3
COMMENTS
In general, if m > 1 and g.f. = Sum_{k>=0} p(k)^m * x^k / Sum_{k>=0} p(k)*x^k, then a(n, m) ~ exp(Pi*sqrt(2*(m^2 - 1)*n/3)) * ((m^2 - 1)^(m - 3/4) / (2^(2*m - 3/4) * 3^(m/2 - 1/4) * m^(2*m - 1) * n^(m - 1/4))).
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1000
FORMULA
a(n) ~ 2^(3/4) * 3^(3/2) * 5^(13/4) * exp(Pi*sqrt(10*n)) / (2^22 * n^(15/4)).
MATHEMATICA
nmax = 25; CoefficientList[Series[Sum[PartitionsP[k]^4*x^k, {k, 0, nmax}] / Sum[PartitionsP[k]*x^k, {k, 0, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, May 20 2018
STATUS
approved