A208237 G.f.: Sum_{n>=0} n! * x^n * Product_{k=1..n} (1 + k*x) / (1 + k*x + k^2*x^2).

1, 1, 2, 5, 15, 54, 223, 1045, 5474, 31685, 200895, 1384470, 10304431, 82376101, 703949762, 6403761365, 61784985615, 630180031734, 6775001385343, 76572619018165, 907658144193314, 11259399965148005, 145879271404693215, 1970471655222795990, 27702625497930064591

Offset

0,3

Comments

Compare to the identity:

Sum_{n>=0} n! * x^n * Product_{k=1..n} (1 + x) / (1 + k*x + k*x^2) = 1/(1-x-x^2).

Compare also to the g.f. of A136127:

x*Sum_{n>=0} n! * x^n * Product_{k=1..n} (2 + k*x) / (1 + 2*k*x + k^2*x^2).

Links

Vaclav Kotesovec, Table of n, a(n) for n = 0..270

Formula

a(n) ~ 2 * 3^(n/2 + 5/4) * n^(n+2) / (exp(n) * Pi^(n+3/2)). - Vaclav Kotesovec, Nov 02 2014

Example

G.f.: A(x) = 1 + x + 2*x^2 + 5*x^3 + 15*x^4 + 54*x^5 + 223*x^6 + 1045*x^7 +...
where
A(x) = 1 + x*(1+x)/(1+x+x^2) + 2!*x^2*(1+x)*(1+2*x)/((1+x+x^2)*(1+2*x+4*x^2)) + 3!*x^3*(1+x)*(1+2*x)*(1+3*x)/((1+x+x^2)*(1+2*x+4*x^2)*(1+3*x+9*x^2)) + 4!*x^4*(1+x)*(1+2*x)*(1+3*x)*(1+4*x)/((1+x+x^2)*(1+2*x+4*x^2)*(1+3*x+9*x^2)*(1+4*x+16*x^2)) +...

Prog

(PARI) {a(n)=polcoeff( sum(m=0, n, m!*x^m*prod(k=1, m, (1+k*x)/(1+k*x+k^2*x^2 +x*O(x^n))) ), n)}
for(n=0, 30, print1(a(n), ", "))

Crossrefs

Cf. A221370, A204064, A204066, A208236, A136127.

Sequence in context: A056841 A185040 A369599 * A321958 A107112 A193318

Adjacent sequences: A208234 A208235 A208236 * A208238 A208239 A208240

Keyword

nonn

Author

Paul D. Hanna, Jan 11 2013

Status

approved