OFFSET
0,6
COMMENTS
a(n) is the number of partitions of n into parts 3*k-1 of k kinds (k>=1).
In general, if s>0, t>0, GCD(s,t)=1 and g.f. = Product_{k>=1} 1/(1 - x^(s*k-t))^k then a(n) ~ s^(t^2/(3*s^2) - 7/18) * n^(t^2/(6*s^2) - 25/36) * exp(d(s,t) - Pi^4 * t^2 / (432*s^2 * Zeta(3)) + Pi^2 * t * 2^(2/3) * s^(2/3) * n^(1/3) / (12 * s^2 * Zeta(3)^(1/3)) + 3*Zeta(3)^(1/3) * n^(2/3) / (2^(2/3)*s^(2/3))) / (2^(t^2/(6*s^2) + 11/36) * sqrt(3*Pi) * Zeta(3)^(t^2/(6*s^2) - 7/36)), where d(s,t) = Integral_{x=0..infinity} 1/x * (exp(-(s-t)*x)/(1 - exp(-s*x))^2 - 1/(s^2*x^2) - t/(s^2*x) + exp(-x)*(1/12 - t^2/(2*s^2))) dx. - Vaclav Kotesovec, Oct 12 2015
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..10000
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
Vaclav Kotesovec, Graph - The asymptotic ratio (70000 terms)
FORMULA
a(n) ~ Zeta(3)^(19/108) * exp(d1 - Pi^4 / (3888*Zeta(3)) + Pi^2 * n^(1/3) / (2^(4/3)*3^(7/3) * Zeta(3)^(1/3)) + 3^(1/3) * Zeta(3)^(1/3) * n^(2/3) / 2^(2/3)) / (2^(35/108) * 3^(23/27) * sqrt(Pi) * n^(73/108)), where d1 = A263030 = Integral_{x=0..infinity} 1/x*(exp(-2*x)/(1 - exp(-3*x))^2 - 1/(9*x^2) - 1/(9*x) + exp(-x)/36) = -0.188708191979528532376410098649207973592114467268429221509... . - Vaclav Kotesovec, Oct 08 2015
MAPLE
with(numtheory):
a:= proc(n) option remember; `if`(n=0, 1, add(add(d*
`if`(irem(d+3, 3, 'r')=2, r, 0), d=divisors(j))*a(n-j), j=1..n)/n)
end:
seq(a(n), n=0..60); # Alois P. Heinz, Oct 05 2015
MATHEMATICA
nmax=100; CoefficientList[Series[Product[1/(1-x^(3k-1))^k, {k, 1, nmax}], {x, 0, nmax}], x]
nmax=100; CoefficientList[Series[E^Sum[1/j*x^(2*j)/(1-x^(3j))^2, {j, 1, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Oct 04 2015
STATUS
approved