A155926 G.f. satisfies: A(x) = B(x*A(x)) where A(x) = Sum_{n>=0} a(n)*x^n/[n!*(n+1)!/2^n] and B(x) = Sum_{n>=0} x^n/[n!*(n+1)!/2^n].

1, 1, 4, 37, 621, 16526, 640207, 34039027, 2379382609, 211619306134, 23337543447296, 3125553148981176, 499716551101393705, 94016487294245251308, 20561796731966531616954, 5172827581575899147920471

Offset

0,3

Links

Table of n, a(n) for n=0..15.

Formula

a(n) = A105558(n)/(n+1) = A105556(2n,n)/(n+1) = [N^(n+1)](n+1,1)/(n+1) for n>=0, where N^(n+1) is the (n+1)-th matrix power of the Narayana triangle N=A001263.

G.f.: A(x) = Series_Reversion[x/B(x)]/x where B(x) = A(x/B(x)) = Sum_{n>=0} x^n/[n!*(n+1)!/2^n].

G.f. satisfies: A(x) = C(x*A(x)^2) and C(x) = A(x/C(x)^2) where C(x) is the g.f. of A155927.

Example

G.f.: A(x) = 1 + x + 4*x^2/3 + 37*x^3/18 + 621*x^4/180 + 16526*x^5/2700 +...+ a(n)*x^n/[n!*(n+1)!/2^n] +...
B(x) = 1 + x + 1/3*x^2 + 1/18*x^3 + 1/180*x^4 +...+ x^n/[n!*(n+1)!/2^n] +... where
A(x) = B(x*A(x)) and B(x) = A(x/B(x)) ;
1/B(x) = 1 - x + 2*x^2/3 - 7*x^3/18 + 39*x^4/180 - 321*x^5/2700 +...+ (-1)^n*A103365(n)*x^n/[n!*(n+1)!/2^n] +...
Also, A(x) = C(x*A(x)^2) where:
C(x) = 1 + x - 2*x^2/3 + 19*x^3/18 - 379*x^4/180 + 12726*x^5/2700 +...+ A155927(n)*x^(n+1)/[n!*(n+1)!/2^n] +...
A(x)^2 = 1 + 2*x + 11*x^2/3 + 122*x^3/18 + 2302*x^4/180 + 66482*x^5/2700 +...

Prog

(PARI) {a(n)=local(F=sum(k=0, n, x^k/(k!*(k+1)!/2^k))+x*O(x^n)); polcoeff(serreverse(x/F)/x, n)*n!*(n+1)!/2^n}
(PARI) {a(n)=local(N=matrix(n+1, n+1, m, j, if(m>=j, binomial(m-1, j-1)*binomial(m, j-1)/j))); sum(j=0, n, (N^n)[n+1, j+1])/(n+1)}

Crossrefs

Cf. A105558, A105556, A001263, A103365, A006472, A090181, A155927, A155928 (A^2).

Sequence in context: A097998 A352470 A177385 * A094408 A348860 A352469

Adjacent sequences: A155923 A155924 A155925 * A155927 A155928 A155929

Keyword

nonn

Author

Paul D. Hanna, Jan 30 2009

Status

approved