OFFSET
0,4
COMMENTS
Compare g.f. to the product:
_ exp(x) = Product_{n>=1} [Product_{k=0..n} (1 + k*x)^(-(-1)^k*binomial(n,k))]^(1/n),
more explicitly:
_ exp(x) = [(1+x)^1/(1+0*x)^1] * [(1+x)^2/((1+0*x)^1*(1+2*x)^1)]^(1/2) * [(1+x)^3*(1+3*x)^1/((1+0*x)^1*(1+2*x)^3)]^(1/3) * [(1+x)^4*(1+3*x)^4/((1+0*x)^1*(1+2*x)^6*(1+4*x)^1)]^(1/4) * ...
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..320
FORMULA
G.f.: A(x) = exp( Sum_{n>=1} 2*A000670(n)*x^(n+1)/(n+1) ) where Sum_{n>=0} A000670(n)*x^n = Sum_{n>=0} n!*x^n/Product_{k=0..n} (1-k*x).
a(n) ~ n! / (n^2 * (log(2))^n). - Vaclav Kotesovec, Nov 19 2014
EXAMPLE
G.f.: A(x) = 1 + x^2 + 2*x^3 + 7*x^4 + 32*x^5 + 189*x^6 + 1382*x^7 +...
where
A(x) = [(1-x)^2/((1-0*x)^1*(1-2*x)^1)] * [(1-x)^4*(1-3*x)^4/((1-0*x)^1*(1-2*x)^6*(1-4*x)^1)]^(1/2) * [(1-x)^6*(1-3*x)^20*(1-5*x)^6/((1-0*x)^1*(1-2*x)^15*(1-4*x)^15*(1-6*x)^1)]^(1/3) * ...
Also,
log(A(x)) = 2*x^2/2 + 6*x^3/3 + 26*x^4/4 + 150*x^5/5 + 1082*x^6/6 + 9366*x^7/7 +...+ 2*A000670(n)*x^(n+1)/(n+1) +...
thus, the logarithmic derivative equals the series:
A'(x)/A(x) = 2*x/(1-x) + 2*2!*x^2/((1-x)*(1-2*x)) + 2*3!*x^3/((1-x)*(1-2*x)*(1-3*x)) + 2*4!*x^4/((1-x)*(1-2*x)*(1-3*x)*(1-4*x)) +...
PROG
(PARI) {a(n)=polcoeff(prod(m=1, n, prod(k=1, 2*m, (1-k*x+x*O(x^n))^(-(-1)^k*binomial(2*m, k)/m))), n)}
(PARI) /* Using o.g.f. exp( Sum_{n>=1} 2*A000670(n)*x^(n+1)/(n+1) ): */
{a(n)=polcoeff(exp(intformal(2*sum(m=1, n+1, m!*x^m/prod(k=1, m, 1-k*x+x*O(x^n))))), n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Sep 25 2011
STATUS
approved