

A153197


Binomial transform of Bessel numbers A006789.


5



1, 2, 5, 15, 51, 189, 748, 3128, 13731, 62969, 300552, 1488704, 7634723, 40464741, 221311617, 1247444859, 7238458309, 43196661875, 264878725516, 1667564565616, 10770316016557, 71314258947903, 483765644021787, 3359905164274725
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

0,2


COMMENTS

Prefaced with a 1: (1, 1, 2, 5, 15, 51, ...), convolved with A006789 = A006789; identical to reversing k terms of one sequence then taking the dot product of k terms of the other: e.g., A006789(6) = 143 = (51, 15, 5, 2, 1, 1) dot (1, 1, 2, 5, 14, 43) = (51 + 15 + 10 + 10 + 14 + 43). Equals row sums of triangle A153199. A153197 can be generated from the Hankel transform [1,1,1,...] by taking successive iterates of the operations: (binomial transform of [1,1,1,...]) followed by INVERT transform of the result, then binomial transform of the result, (repeat cycle)...; until the operations converge upon a two sequence fixed limit cycle of A006789 and A153197. Or, the infinite set of operations Q may begin: INVERT transform of [1,1,1,..] followed by binomial transform of the result, INVERT transform of the result, etc; until the operations again converge upon A006789 and A153197. The two sequences A006789 and A153197 have the mutual relationships that binomial transform of A006789 = A153197; while the INVERT transform of A153197 prefaced with a 1 and then prefaced with a 1 afterward = A006789.
Product of the two sequences (A153197 prefaced with a 1) and A006789 = A006789 with offset 1. Or, (1,1,2,5,15,51,...) * (1,1,2,5,14,43,...) = (1,2,5,14,43,...).
Conjecture: Given any sequence with Hankel transform of [1,1,1,...], performing alternate operations: binomial transform followed by INVERT transform, then binomial transform of the last result (repeat); or INVERT transform starting first, will converge upon A006789 and A153197 as a two sequence limit cycle. The conjecture can be extended to any Hankel transform (and their accompanying sequence set): analogous operations will converge upon a Besseltype sequence and its binomial transform.


LINKS

Table of n, a(n) for n=0..23.
Aoife Hennessy, A Study of Riordan Arrays with Applications to Continued Fractions, Orthogonal Polynomials and Lattice Paths, Ph. D. Thesis, Waterford Institute of Technology, Oct. 2011.


FORMULA

INVERT transform of A153197 with 1 prepended is A006789 with initial 1 removed. Bessel numbers with offset 1: (1, 2, 5, 14, 43, 143, ...).
If the offset is 1, the g.f. A(x) satisfies A(x) = x / (1  2*x  x * A( x / (1  x))).  Michael Somos, Mar 06 2011
G.f.: 1/U(0) where U(k)= 1  x*(k+2)  x^2/U(k+1) ; (continued fraction, 1step).  Sergei N. Gladkovskii, Oct 11 2012
G.f.: 1/x^21/xU(0)/x^2 where U(k)= 1  x*(k+1)  x^2/U(k+1); (continued fraction, 1step).  Sergei N. Gladkovskii, Oct 24 2012


EXAMPLE

a(3) = 15 = (1, 3, 3, 1) dot (1, 1, 2, 5) = (1 + 3 + 6 + 5).
G.f. = 1 + 2*x + 5*x^2 + 15*x^3 + 51*x^4 + 189*x^5 + 748*x^6 + 3128*x^7 + ...


MATHEMATICA

a[ n_] := If[ n < 0, 0, With[{m = n + 1}, SeriesCoefficient[ Nest[ Series[ x / (1  2 x  x (Normal[#] /. x > x / (1  x))), {x, 0, m}] &, 0, Ceiling[m/2]], {x, 0, m}]]]; (* Michael Somos, Aug 05 2014 *)


PROG

(PARI) {a(n) = local( A = O(x) ); if( n<0, 0, n++; for( k=1, ceil(n/2), A = x / (1  2*x  x * subst( A, x, x / (1  x)))); polcoeff( A, n))}; /* Michael Somos, Mar 06 2011 */


CROSSREFS

Cf. A153198, A153199, A006789.
Sequence in context: A181768 A279554 A279555 * A299968 A279556 A108307
Adjacent sequences: A153194 A153195 A153196 * A153198 A153199 A153200


KEYWORD

nonn


AUTHOR

Gary W. Adamson, Dec 20 2008


EXTENSIONS

Edited by N. J. A. Sloane, Mar 06 2011


STATUS

approved



