OFFSET
0,2
COMMENTS
The Hankel transform (see A001906 for definition) of this sequence is A000302 (powers of 4): 1, 4, 16, 64, 256, 1024, ... - Philippe Deléham, Aug 17 2005
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..200
Vincent Pilaud, V Pons, Permutrees, arXiv preprint arXiv:1606.09643, 2016
FORMULA
Differential equation: (16*x^3 + 12*x^2 - 8*x + 1) * x*(d/dx)A(x) + (8x^3 - 12*x^2 + 6*x - 1) * A(x) + (8x^2 - 6*x + 1) = 0.
G.f.: ((1 - 4*x) + 2*x * sqrt(1 - 4*x)) / (1 - 4*x - 4*x^2). a(n) * (n-1) = a(n-1) * (8*n - 14) - a(n-2) * 12*(n-3) - a(n-3) * 8*(2*n - 5), n > 2. Hankel number wall zig-zag diagonal is A011782. - Michael Somos, Sep 14 2003
INVERT transform of A028329 (offset 1). - Michael Somos, Jan 05 2012
G.f.: (1-2*x*f(x))/(1-2*x*f(x)-2*x) where f(x) is the g.f. of A000108 (Catalan numbers). - Philippe Deléham, Jan 30 2012
a(n) ~ (1-1/sqrt(2))*(2+2*sqrt(2))^n. - Vaclav Kotesovec, Oct 14 2012
From Peter Bala, Feb 05 2017: (Start)
EXAMPLE
1 + 2*x + 8*x^2 + 36*x^3 + 168*x^4 + 796*x^5 + 3800*x^6 + 18216*x^7 + ...
MAPLE
1/(1-x/(sqrt(1/4-x))): series(%, x, 23): seq(coeff(%, x, n), n=0..22); # Peter Luschny, Feb 06 2017
MATHEMATICA
Table[SeriesCoefficient[((1-4*x)+2*x*Sqrt[1-4*x])/(1-4*x-4*x^2), {x, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Oct 14 2012 *)
PROG
(PARI) {a(n) = if( n<0, 0, polcoeff((1 - 4*x + 2*x * sqrt(1 - 4*x + x * O(x^n))) /(1 - 4*x - 4*x^2), n))} /* Michael Somos, Jan 05 2012 */
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paul D. Hanna, Jun 10 2003, Jun 11 2003
STATUS
approved