login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A086616 Partial sums of the large Schroeder numbers (A006318). 14

%I #44 Aug 06 2020 11:35:55

%S 1,3,9,31,121,515,2321,10879,52465,258563,1296281,6589727,33887465,

%T 175966211,921353249,4858956287,25786112993,137604139011,737922992937,

%U 3974647310111,21493266631001,116642921832963,635074797251889,3467998148181631,18989465797056721,104239408386028035

%N Partial sums of the large Schroeder numbers (A006318).

%C Row sums of triangle A086614. - _Paul D. Hanna_, Jul 24 2003

%C Hankel transform is A136577(n+1). - _Paul Barry_, Jun 03 2009

%H Vincenzo Librandi, <a href="/A086616/b086616.txt">Table of n, a(n) for n = 0..200</a>

%H Paul Barry, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL12/Barry3/barry93.html">Continued fractions and transformations of integer sequences</a>, JIS 12 (2009), #09.7.6.

%F G.f.: A(x) = 1/(1 - x)^2 + x*A(x)^2.

%F a(1) = 1 and a(n) = n + Sum_{i=1..n-1} a(i)*a(n-i) for n >= 2. - _Benoit Cloitre_, Mar 16 2004

%F G.f.: (1 - x - sqrt(1 - 6*x + x^2))/(2*x*(1 - x)). Cf. A001003. - _Ralf Stephan_, Mar 23 2004

%F a(n) = Sum_{k=0..n} C(n+k+1, 2*k+1) * A000108(k). - _Paul Barry_, Jun 03 2009

%F Recurrence: (n+1)*a(n) = (7*n-2)*a(n-1) - (7*n-5)*a(n-2) + (n-2)*a(n-3). - _Vaclav Kotesovec_, Oct 14 2012

%F a(n) ~ sqrt(24 + 17*sqrt(2))*(3 + 2*sqrt(2))^n/(4*sqrt(Pi)*n^(3/2)). - _Vaclav Kotesovec_, Oct 14 2012

%e a(1) = 2 + 1 = 3;

%e a(2) = 3 + 4 + 2 = 9;

%e a(3) = 4 + 10 + 12 + 5 = 31;

%e a(4) = 5 + 20 + 42 + 40 + 14 = 121.

%t Table[SeriesCoefficient[(1-x-Sqrt[1-6*x+x^2])/(2*x*(1-x)),{x,0,n}],{n,0,20}] (* _Vaclav Kotesovec_, Oct 14 2012 *)

%o (Sage) # Generalized algorithm of L. Seidel

%o def A086616_list(n) :

%o D = [0]*(n+2); D[1] = 1

%o b = True; h = 2; R = []

%o for i in range(2*n) :

%o if b :

%o for k in range(h,0,-1) : D[k] += D[k-1]

%o else :

%o for k in range(1,h, 1) : D[k] += D[k-1]

%o R.append(D[h-1]); h += 1;

%o b = not b

%o return R

%o A086616_list(23) # _Peter Luschny_, Jun 02 2012

%o (PARI) x='x+O('x^66); Vec((1-x-sqrt(1-6*x+x^2))/(2*x*(1-x))) \\ _Joerg Arndt_, May 10 2013

%Y Cf. A086614 (triangle), A086615 (antidiagonal sums).

%Y Cf. A006318.

%K nonn

%O 0,2

%A _Paul D. Hanna_, Jul 24 2003

%E Name changed using a comment of _Emeric Deutsch_ from Dec 20 2004. - _Peter Luschny_, Jun 03 2012

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 24 10:11 EDT 2024. Contains 371935 sequences. (Running on oeis4.)