# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a126931
Showing 1-1 of 1

%I A126931 #22 Sep 08 2022 08:45:29
%S A126931 1,3,10,33,110,366,1220,4065,13550,45162,150540,501786,1672620,
%T A126931 5575356,18584520,61948257,206494190,688313490,2294378300,7647926046,
%U A126931 25493086820,84976950468,283256501560,944188318938,3147294396460
%N A126931 a(n) = A127359(n+1)/2 - A127359(n).
%C A126931 Hankel transform is A000012=[1,1,1,1,1,1,1,...].
%C A126931 a(n) is the number of Motzkin paths of length n in which the (1,0)-steps at level 0 come in 3 colors and there are no (1,0)-steps at a higher level. Example: a(3)=33 because, denoting  U=(1,1), H=(1,0), and D=(1,-1), we have 3^3 = 27 paths of shape HHH, 3 paths of shape HUD, and 3 paths of shape UDH. - _Emeric Deutsch_, May 02 2011
%H A126931 Michael De Vlieger, <a href="/A126931/b126931.txt">Table of n, a(n) for n = 0..1912</a>
%H A126931 Isaac DeJager, Madeleine Naquin, Frank Seidl, <a href="https://www.valpo.edu/mathematics-statistics/files/2019/08/Drube2019.pdf">Colored Motzkin Paths of Higher Order</a>, VERUM 2019.
%F A126931 G.f.: 1/(1-3x-x^2/(1-x^2/(1-x^2/(1-x^2/(1-... (continued fraction). - _Paul Barry_, Mar 10 2009
%F A126931 G.f.: 2/(1 - 6*x + sqrt(1-4*x^2)). - _Emeric Deutsch_, May 02 2011
%F A126931 Conjecture: 3*(n+1)*a(n) +10*(-n-1)*a(n-1) +12*(-n+2)*a(n-2) +40*(n-2)*a(n-3)=0. - _R. J. Mathar_, Nov 26 2012
%p A126931 A127359 := proc(n) add(binomial(n,floor(k/2))*3^(n-k), k=0..n) ; end proc: A126931 := proc(n) A127359(n+1)/2-A127359(n) ; end proc: seq(A126931(n),n=0..50) ; # _R. J. Mathar_, Mar 25 2010
%t A126931 With[{s = Partition[#, 2, 1] &@ Array[Sum[Binomial[#, Floor[k/2]]*3^(# - k), {k, 0, #}] &, 26, 0]}, Map[#2/2 - #1 & @@ # &, s]] (* _Michael De Vlieger_, Dec 15 2019 *)
%t A126931 CoefficientList[Series[2/(1-6*x+Sqrt[1-4*x^2]), {x,0,30}], x] (* _G. C. Greubel_, Jan 29 2020 *)
%o A126931 (PARI) my(x='x+O('x^30)); Vec( 2/(1 - 6*x + sqrt(1-4*x^2)) ) \\ _G. C. Greubel_, Jan 29 2020
%o A126931 (Magma) R<x>:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( 2/(1 - 6*x + Sqrt(1-4*x^2)) )); // _G. C. Greubel_, Jan 29 2020
%o A126931 (Sage)
%o A126931 def A126931_list(prec):
%o A126931     P.<x> = PowerSeriesRing(ZZ, prec)
%o A126931     return P( 2/(1 - 6*x + sqrt(1-4*x^2)) ).list()
%o A126931 A126931_list(30) # _G. C. Greubel_, Jan 29 2020
%Y A126931 Cf. A127359.
%K A126931 nonn
%O A126931 0,2
%A A126931 _Philippe Deléham_, Mar 17 2007
%E A126931 More terms from _R. J. Mathar_, Mar 25 2010

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE