OFFSET
1,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000
FORMULA
L.g.f.: -log( (1-x + sqrt((1-x)^2 - 4*x^3*(1-x)))/2 ).
a(n) ~ 2^n/sqrt(Pi*n). - Vaclav Kotesovec, Mar 06 2014
EXAMPLE
L.g.f.: L(x) = x + x^2/2 + 4*x^3/3 + 5*x^4/4 + 6*x^5/5 + 16*x^6/6 + 29*x^7/7 +...
such that the l.g.f. equals the series:
L(x) = ((1-x) + x)*x +
((1-x)^2 + 2^2*x*(1-x) + x^2)*x^2/2 +
((1-x)^3 + 3^2*x*(1-x)^2 + 3^2*x^2*(1-x) + x^3)*x^3/3 +
((1-x)^4 + 4^2*x*(1-x)^3 + 6^2*x^2*(1-x)^2 + 4^2*x^3*(1-x) + x^4)*x^4/4 +
((1-x)^5 + 5^2*x*(1-x)^4 + 10^2*x^2*(1-x)^3 + 10^2*x^3*(1-x)^2 + 5^2*x^4*(1-x) + x^5)*x^5/5 +...
where exponentiation yields the g.f. of A216604:
exp(L(x)) = 1 + x + x^2 + 2*x^3 + 3*x^4 + 4*x^5 + 7*x^6 + 12*x^7 + 19*x^8 + 33*x^9 + 59*x^10 + 102*x^11 + 181*x^12 +...+ A216604(n)*x^n +...
MATHEMATICA
Rest[CoefficientList[Series[-Log[(1-x+Sqrt[(1-x)^2-4*x^3*(1-x)])/2], {x, 0, 20}], x]*Range[0, 20]] (* Vaclav Kotesovec, Mar 06 2014 *)
PROG
(PARI) {a(n)=n*polcoeff(sum(m=1, n+1, x^m/m*sum(k=0, m, binomial(m, k)^2*x^k*(1-x)^(m-k) + x*O(x^n))), n)}
(PARI) {a(n)=n*polcoeff(log(2/(1-x+sqrt((1-x)^2-4*x^3*(1-x) +x*O(x^n)))), n)}
for(n=1, 40, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Oct 03 2012
STATUS
approved