%I #96 Sep 18 2024 12:54:15
%S 1,4,20,116,740,5028,35700,261780,1967300,15072836,117297620,
%T 924612532,7367204260,59240277988,480118631220,3917880562644,
%U 32163325863300,265446382860420,2201136740855700,18329850024033012,153225552507991140
%N Expansion of (1-3*x-sqrt(9*x^2-10*x+1))/(2*x).
%C More generally coefficients of (1-m*x-sqrt(m^2*x^2-(2*m+4)*x+1))/(2*x) are given by a(0)=1 and, for n>0, a(n) = (1/n)*Sum_{k=0..n}(m+1)^k*binomial(n,k)*binomial(n,k-1).
%C a(n) = number of lattice paths from (0,0) to (n+1,n+1) that consist of steps (i,0) and (0,j) with i,j>=1 and that stay strictly below the diagonal line y=x except at the endpoints. (See Coker reference.) Equivalently, a(n) = number of marked Dyck (n+1)-paths where the vertices in the middle of each UU and each DD are available to be marked (or not): consider the original path as a Dyck path with a mark at each vertex where two horizontal (or two vertical) steps abut. If only the UU vertices are available for marking, then the counting sequence is the little Schroeder number A001003. - _David Callan_, Jun 07 2006
%C Hankel transform is 4^C(n+1,2). - _Philippe Deléham_, Feb 11 2009
%C a(n) is the number of Schroder paths of semilength n in which the (2,0)-steps come in 3 colors. Example: a(2)=20 because, denoting U=(1,1), H=(2,0), D=(1,-1), we have 3^2=9 paths of shape HH, 3 paths of shape HUD, 3 paths of shape UDH, 3 paths of shape UHD, and 1 path of each of the shapes UDUD, UUDD. - _Emeric Deutsch_, May 02 2011
%C (1 + 4x + 20x^2 + 116x^3 + ...) = (1 + 5x + 29x^2 + 185x^3 + ...) * 1/(1 + x + 5x^2 + 29x^3 +185x^4 + ...); where A059231 = (1, 5, 29, 185, 1257, ...) - _Gary W. Adamson_, Nov 17 2011
%C The first differences between the row sums of the triangle A226392. - _J. M. Bergot_, Jun 21 2013
%H Vincenzo Librandi, <a href="/A082298/b082298.txt">Table of n, a(n) for n = 0..300</a>
%H Paul Barry, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL9/Barry/barry91.html">On Integer-Sequence-Based Constructions of Generalized Pascal Triangles</a>, Journal of Integer Sequences, Vol. 9 (2006), Article 06.2.4.
%H Paul Barry, <a href="https://www.emis.de/journals/JIS/VOL22/Barry3/barry422.html">Generalized Catalan Numbers Associated with a Family of Pascal-like Triangles</a>, J. Int. Seq., Vol. 22 (2019), Article 19.5.8.
%H Paul Barry and Aoife Hennessy, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL14/Barry2/barry126.html">A Note on Narayana Triangles and Related Polynomials, Riordan Arrays, and MIMO Capacity Calculations</a>, J. Int. Seq. 14 (2011), Article 11.3.8.
%H Zhi Chen and Hao Pan, <a href="https://arxiv.org/abs/1608.02448">Identities involving weighted Catalan, Schroder and Motzkin paths</a>, arXiv:1608.02448 [math.CO], 2016. See eq. (1.13), a=4, b=1.
%H Curtis Coker, <a href="http://dx.doi.org/10.1016/S0012-365X(03)00037-2">Enumerating a class of lattice paths</a>, Disc. Math. (2003) Vol. 271, Issues 1-3, 13-28.
%H Aoife Hennessy, <a href="http://repository.wit.ie/1693/1/AoifeThesis.pdf">A Study of Riordan Arrays with Applications to Continued Fractions, Orthogonal Polynomials and Lattice Paths</a>, Ph. D. Thesis, Waterford Institute of Technology, Oct. 2011.
%H John Machacek, <a href="https://arxiv.org/abs/2105.02417">Lattice walks ending on a coordinate hyperplane avoiding backtracking and repeats</a>, arXiv:2105.02417 [math.CO], 2021. See Thm 4.4, G(x,F^1)
%H Greg Morrow, <a href="https://faculty.uccs.edu/gmorrow/wp-content/uploads/sites/19/2024/08/Prob_Distrns__Int_Seqs_related_to_Rook_Paths_Aug_9_2024.pdf">Some probability distributions and integer sequences related to rook paths</a>, Univ. Colorado Springs (2024). See pp. 4, 22
%F a(0)=1, n>0 a(n) = (1/n)*Sum_{k=0..n} 4^k*binomial(n, k)*binomial(n, k-1).
%F a(1)=1, a(n) = 3*a(n-1) + Sum_{i=1..n-1} a(i)*a(n-i). - _Benoit Cloitre_, Mar 16 2004
%F a(n) = Sum_{k=0..n} 1/(n+1) Binomial(n+1,k)Binomial(2n-k,n-k)3^k. - _David Callan_, Jun 07 2006
%F From _Paul Barry_, Feb 01 2009: (Start)
%F G.f.: 1/(1-3x-x/(1-3x-x/(1-3x-x/(1-... (continued fraction);
%F a(n) = Sum_{k=0..n} binomial(n+k,2k)*3^(n-k)*A000108(k). (End)
%F a(n) = Sum_{k=0..n} A060693(n,k)*3^k. - _Philippe Deléham_, Feb 11 2009
%F D-finite with recurrence: (n+1)*a(n) = 5*(2n-1)*a(n-1)-9*(n-2)*a(n-2). - _Paul Barry_, Oct 22 2009
%F G.f.: 1/(1- 4x/(1-x/(1-4x/(1-x/(1-4x/(1-... (continued fraction). - Aoife Hennessy (aoife.hennessy(AT)gmail.com), Dec 02 2009
%F G.f.: (1-3*x-sqrt(9*x^2-10*x+1))/(2*x) = (1-G(0))/x; G(k) = 1+x*3-x*4/G(k+1); (continued fraction, 1-step). - _Sergei N. Gladkovskii_, Jan 05 2012
%F a(n) ~ 3^(2*n+1)/(sqrt(2*Pi)*n^(3/2)). - _Vaclav Kotesovec_, Oct 14 2012
%F a(n) = 4*hypergeom([1 - n, -n], [2], 4)) for n>0. - _Peter Luschny_, May 22 2017
%F G.f. A(x) satisfies: A(x) = (1 + x*A(x)^2) / (1 - 3*x). - _Ilya Gutkovskiy_, Jun 30 2020
%F G.f.: (1+2*x*F(x))^2, where F(x) is the g.f. for A099250. - _Alexander Burstein_, May 11 2021
%p A082298_list := proc(n) local j, a, w; a := array(0..n); a[0] := 1;
%p for w from 1 to n do a[w] := 4*a[w-1]+add(a[j]*a[w-j-1], j=1..w-1) od;convert(a,list)end: A082298_list(20); # _Peter Luschny_, May 19 2011
%p a := n -> `if`(n=0, 1, 4*hypergeom([1 - n, -n], [2], 4)):
%p seq(simplify(a(n)), n=0..20); # _Peter Luschny_, May 22 2017
%t gf[x_] = (1 - 3*x - Sqrt[(9*x^2 - 10*x + 1)])/(2*x); CoefficientList[Series[gf[x], {x, 0, 20}], x] (* _Jean-François Alcover_, Jun 01 2011 *)
%o (PARI) a(n)=if(n<1,1,sum(k=0,n,4^k*binomial(n,k)*binomial(n,k-1))/n)
%o (Magma) Q:=Rationals(); R<x>:=PowerSeriesRing(Q, 40); Coefficients(R!((1-3*x-Sqrt(9*x^2-10*x+1))/(2*x))); // _G. C. Greubel_, Feb 10 2018
%Y Cf. A006318, A047891, A059231, A099250.
%K nonn
%O 0,2
%A _Benoit Cloitre_, May 10 2003