login
Expansion of (1+x+sqrt(1-2x-3x^2))/2.
5

%I #35 Feb 26 2021 07:50:49

%S 1,0,-1,-1,-2,-4,-9,-21,-51,-127,-323,-835,-2188,-5798,-15511,-41835,

%T -113634,-310572,-853467,-2356779,-6536382,-18199284,-50852019,

%U -142547559,-400763223,-1129760415,-3192727797,-9043402501,-25669818476

%N Expansion of (1+x+sqrt(1-2x-3x^2))/2.

%C A signed variant of the Motzkin numbers A001006. Hankel transform is A168052.

%H G. C. Greubel, <a href="/A168051/b168051.txt">Table of n, a(n) for n = 0..1000</a>

%F D-finite with recurrence: n*a(n) -(2n-3)*a(n-1) -3*(n-3)*a(n-2)=0 if n>2. - _R. J. Mathar_, Dec 20 2011 [Edited by _Michael Somos_, Jan 25 2014]

%F 0 = a(n) * (9*a(n+1) + 15*a(n+2) - 12*a(n+3)) + a(n+1) * (-3*a(n+1) + 10*a(n+2) - 5*a(n+3)) + a(n+2) * (a(n+2) + a(n+3)) if n>0. - _Michael Somos_, Jan 25 2014

%F G.f.: 1 + x - (x + x^2) / (1 + x - (x + x^2) / (1 + x - ...)). - _Michael Somos_, Mar 27 2014

%F Convolution inverse of A005043. - _Michael Somos_, Mar 27 2014

%F a(n) ~ -3^(n - 1/2) / (2 * sqrt(Pi) * n^(3/2)). - _Vaclav Kotesovec_, Jun 05 2018

%F From _Gennady Eremin_, Feb 25 2021: (Start)

%F For n > 1, a(n) = A167022(n) / 2.

%F G.f.: (1 + x + A(x)) / 2, where A(x) is the g.f. of A167022. (End)

%e G.f. = 1 - x^2 - x^3 - 2*x^4 - 4*x^5 - 9*x^6 - 21*x^7 - 51*x^8 - 127*x^9 + ...

%t a[ n_] := SeriesCoefficient[ (1 + x + Sqrt[1 - 2 x - 3 x^2]) / 2, {x, 0, n}] (* _Michael Somos_, Jan 25 2014 *)

%o (PARI) {a(n) = polcoeff( (1 + x + sqrt(1 - 2*x - 3*x^2 + x * O(x^n))) / 2, n)}; /* _Michael Somos_, Jan 25 2014 */

%Y Cf. A168049, A166587, A167022.

%Y Cf. A001006, A005043.

%K easy,sign

%O 0,5

%A _Paul Barry_, Nov 17 2009