login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A179191 Expansion of o.g.f. (1/2)*(-1 + 1/sqrt(1 - 4*x - 4*x^2)). 2
0, 1, 4, 16, 68, 296, 1312, 5888, 26672, 121696, 558464, 2574848, 11917952, 55345408, 257741824, 1203224576, 5629027072, 26383656448, 123868321792, 582414688256, 2742116907008, 12926036258816, 60998951747584 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

COMMENTS

G.f. A(x) satisfies A(x)^2 + A(x) = (x^2 + x)/(1 - 4*x - 4*x^2). - Michael Somos, Jan 28 2019

LINKS

G. C. Greubel, Table of n, a(n) for n = 0..1000

M. Dziemianczuk, On Directed Lattice Paths With Additional Vertical Steps, arXiv:1410.5747 [math.CO], 2014.

FORMULA

G.f.: (1/2)*(-1 + 1/sqrt(1 - 4*x - 4*x^2)).

G.f.: A(x) = x*A001006(A000045(x)/x-1)/(1-x*A001006(A000045(x)/x-1)).

a(n) = Sum_{m=1..n} m*Sum_{k=m..n} Sum_{i=k..n} binomial(i-1,k-1)* binomial(i,n-i))*Sum_{j=0..k} binomial(j,2*j-m-k)*binomial(k,j))/k)). - Vladimir Kruchinin, Mar 11 2011

a(n) = Sum_{k=0..n} 2^(n-k-1)*binomial(n,k)*binomial(n-k,k). - Vladimir Kruchinin, Mar 12 2015

From Vaclav Kotesovec, Jan 26 2019: (Start)

Recurrence: n*a(n) = 2*(2*n - 1)*a(n-1) + 4*(n-1)*a(n-2).

a(n) ~ 2^(n - 7/4) * (1 + sqrt(2))^(n + 1/2) / sqrt(Pi*n). (End)

0 = a(n)*(16*a(n+1) +24*a(n+2) -8*a(n+3)) + a(n+1)*(+8*a(n+1) +16*a(n+2) -6*a(n+3)) + a(n+2)*(-2*a(n+2) +a(n+3)) for all n in Z except n=-1. - Michael Somos, Jan 27 2019

EXAMPLE

G.f. = x + 4*x^2 + 16*x^3 + 68*x^4 + 296*x^5 + 1312*x^6 + 5888*x^7 + ....

MATHEMATICA

CoefficientList[1/2 (-1 + (1-4x-4x^2)^(-1/2)) + O[x]^23, x] (* Jean-Fran├žois Alcover, Jul 27 2018 *)

PROG

(Maxima)

a(n):=sum(m*sum(sum(binomial(i-1, k-1)*binomial(i, n-i), i, k, n)*sum(binomial(j, 2*j-m-k)*binomial(k, j), j, 0, k)/k, k, m, n), m, 1, n); /* Vladimir Kruchinin, Mar 11 2011 */

a(n):=sum(2^(n-k-1)*binomial(n, k)*binomial(n-k, k), k, 0, n); /* Vladimir Kruchinin, Mar 12 2015 */

(PARI) my(x='x+O('x^30)); concat([0], Vec((-1 +1/sqrt(1-4*x-4*x^2))/2)) \\ G. C. Greubel, Jan 25 2019

(MAGMA) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); Coefficients(R!( (-1 + 1/Sqrt(1-4*x-4*x^2))/2 )); // G. C. Greubel, Jan 25 2019

(Sage) ((-1 + 1/sqrt(1-4*x-4*x^2))/2).series(x, 20).coefficients(x, sparse=False) # G. C. Greubel, Jan 25 2019

CROSSREFS

Cf. A178693, A178694, A179190.

Sequence in context: A179611 A290912 A089979 * A128730 A151243 A006319

Adjacent sequences:  A179188 A179189 A179190 * A179192 A179193 A179194

KEYWORD

nonn

AUTHOR

Clark Kimberling, Jul 01 2010

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 26 19:29 EDT 2019. Contains 323597 sequences. (Running on oeis4.)