editing
approved
editing
approved
for(n=0, 25, print1(round(A[n+1]), ", "))
approved
editing
editing
approved
O.g.f.: Sum_{n>=0} 2^n / Product_{k=0..n} (3 - k*x). - Paul D. Hanna, Oct 27 2014
...
O.g.f.: G(x) = 1 + 2*x + 14*x^2 + 162*x^3 + 2622*x^4 + 54546*x^5 +...
where
G(x) = 1/3 + 2/(3*(3-x)) + 2^2/(3*(3-x)*(3-2*x)) + 2^3/(3*(3-x)*(3-2*x)*(3-3*x)) + 2^4/(3*(3-x)*(3-2*x)*(3-3*x)*(3-4*x)) + 2^5/(3*(3-x)*(3-2*x)*(3-3*x)*(3-4*x)*(3-5*x)) +...
(PARI) \p100 \\ set precision
{A=Vec(sum(n=0, 600, 1.*2^n/prod(k=0, n, 3 - k*x + O(x^31))))}
for(n=0, 25, print1(round(A[n+1]), ", ")) \\ Paul D. Hanna, Oct 27 2014
approved
editing
proposed
approved
editing
proposed
E.g.f.: (x-2)/3 - LambertW(-2/3*exp((x-2)/3)). - Vaclav Kotesovec, Dec 26 2013
a(n) ~ n^(n-1) / (sqrt(3) * exp(n) * (3*log(3)-3*log(2)-1)^(n-1/2)). - Vaclav Kotesovec, Dec 26 2013
Rest[CoefficientList[1 + InverseSeries[Series[2 + 3*x - 2*Exp[x], {x, 0, 20}], x], x]* Range[0, 20]!] (* Vaclav Kotesovec, Dec 26 2013 *)
approved
editing
proposed
approved
editing
proposed
G.f.: 1/Q(0), where Q(k)= 1 - k*x - 2*x*(k+1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 01 2013
approved
editing