Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #5 Dec 29 2023 12:59:55
%S 1,1,2,3,1,6,3,11,10,23,26,2,47,70,10,102,176,45,221,449,160,5,493,
%T 1121,539,35,1105,2817,1680,196,2516,7031,5082,868,14,5763,17604,
%U 14856,3486,126,13328,43996,42660,12810,840,30995,110147,120338,44640,4410,42
%N Triangle read by rows: T(n,k) is number of Motzkin paths of length n having k UH's, where U=(1,1), H=(1,0) (0<=k<=floor(n/3)).
%C Row n contains 1+floor(n/3) terms. Row sums are the Motzkin numbers (A001006). Column 0 yields A090344. Sum(k*T(n,k),k=0..floor(n/3))=A014531(n-2).
%F G.f.=[1-z-sqrt(1-2z-3z^2-4tz^3+4z^3)]/[2(1-z+tz)z^2].
%e T(4,1)=3 because we have H(UH)D, (UH)DH and (UH)HD, where U=(1,1), H=(1,0), D=(1,-1) (the UH's are shown between parentheses).
%e Triangle begins:
%e 1;
%e 1;
%e 2;
%e 3,1;
%e 6,3;
%e 11,10;
%e 23,26,2;
%e 47,70,10;
%p G:=(1-z-sqrt(1-2*z-3*z^2-4*z^3*t+4*z^3))/2/z^2/(1-z+t*z): Gser:=simplify(series(G,z=0,20)): P[0]:=1: for n from 1 to 16 do P[n]:=coeff(Gser,z^n) od: for n from 0 to 16 do seq(coeff(t*P[n],t^j),j=1..1+floor(n/3)) od; # yields sequence in triangular form
%Y Cf. A001006, A090344, A014531.
%K nonn,tabf
%O 0,3
%A _Emeric Deutsch_, Dec 09 2005