cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 61-70 of 142 results. Next

A002571 From a definite integral.

Original entry on oeis.org

1, 5, 10, 30, 74, 199, 515, 1355, 3540, 9276, 24276, 63565, 166405, 435665, 1140574, 2986074, 7817630, 20466835, 53582855, 140281751, 367262376, 961505400, 2517253800, 6590256025, 17253514249, 45170286749, 118257345970
Offset: 1

Views

Author

Keywords

Comments

a(n) are the row sums of the elements of the Golden Triangle (A180662) with alternating signs. - Alexander Adamchuk, Oct 18 2010
Limit_{n->oo} A002570(n)/A002571(n) = 1/sqrt(5). - Sean A. Irvine, Apr 09 2014

Examples

			From _Paul D. Hanna_, Feb 20 2009: (Start)
G.f.: A(x) = x + 5*x^2 + 10*x^3 + 30*x^4 + 74*x^5 + 199*x^6 + ...
log(1+A(x)) = x + 3^2*x^2/2 + 4^2*x^3/3 + 7^2*x^4/4 + 11^2*x^5/5 + ... (End)
G.f.: A(x) = -1 + 1/((1-x-x^2) * (1-3*x^2+x^4) * (1-4*x^3-x^6) * (1-7*x^4+x^8) * (1-11*x^5-x^10)^2 * (1-18*x^6+x^12)^2 * (1-29*x^7-x^14)^4 * (1-47*x^8+x^16)^5 * (1-76*x^9-x^18)^8 * ...* (1 - Lucas(n)*x^n + (-1)^n*x^(2*n))^A006206(n) * ...). - _Paul D. Hanna_, Jan 07 2012
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A001654, A180662 - The Golden Triangle. - Alexander Adamchuk, Oct 18 2010

Programs

  • Maple
    A002571:=-(-1-4*z-z**2+z**3)/(z**2-3*z+1)/(1+z)**2; # conjectured (probably correctly) by Simon Plouffe in his 1992 dissertation
  • PARI
    {a(n)=polcoeff(exp(sum(m=1,n,(fibonacci(m+1)+fibonacci(m-1))^2*x^m/m)+x*O(x^n)),n)} \\ Paul D. Hanna, Feb 20 2009

Formula

Appears to have g.f. x/((1-3x+x^2)*(1+x)^2). - Ralf Stephan, Apr 14 2004
a(n) = (-1)^n*Sum_{i=1..n+1} (-1)^(i+1)*Fibonacci(i)*Fibonacci(i+1). - Alexander Adamchuk, Jun 16 2006
From Paul D. Hanna, Feb 20 2009: (Start)
Given g.f. A(x), then log(1+A(x)) = Sum_{n>=1} A000204(n)^2 * x^n/n where A000204 is the Lucas numbers.
a(n) = (1/n)*(A000204(n)^2 + Sum_{k=1..n-1} A000204(k)^2*a(n-k)) for n>1, with a(1) = 1. (End)
G.f.: -1 + 1/Product_{n>=1} (1 - Lucas(n)*x^n + (-1)^n*x^(2*n))^A006206(n), where A006206(n) is the number of aperiodic binary necklaces of length n with no subsequence 00. - Paul D. Hanna, Jan 07 2012
a(n) = 8*a(n-2) - 8*a(n-4) + a(n-6) + 2(-1)^n, n>6. - Sean A. Irvine, Apr 09 2014
a(n) - a(n-2) = Fibonacci(n+1)^2. - Peter Bala, Aug 30 2015

Extensions

More terms from Max Alekseyev and Alexander Adamchuk, Oct 18 2010

A144677 Related to enumeration of quantum states (see reference for precise definition).

Original entry on oeis.org

1, 2, 3, 6, 9, 12, 18, 24, 30, 40, 50, 60, 75, 90, 105, 126, 147, 168, 196, 224, 252, 288, 324, 360, 405, 450, 495, 550, 605, 660, 726, 792, 858, 936, 1014, 1092, 1183, 1274, 1365, 1470, 1575, 1680, 1800, 1920, 2040, 2176, 2312, 2448, 2601, 2754, 2907, 3078, 3249, 3420
Offset: 0

Views

Author

N. J. A. Sloane, Feb 06 2009

Keywords

Comments

Equals (1, 2, 3, ...) convolved with (1, 0, 0, 2, 0, 0, 3, ...) = (1 + 2*x + 3*x^2 + ...) * (1 + 2*x^3 + 3*x^6 + ...). - Gary W. Adamson, Feb 23 2010
The Ca2 and Ze4 triangle sums, see A180662 for their definitions, of the Connell-Pol triangle A159797 are linear sums of shifted versions of the sequence given above, e.g., Ca2(n) = a(n-1) + 2*a(n-2) + 3*a(n-3) + a(n-4). - Johannes W. Meijer, May 20 2011

Crossrefs

Programs

  • Magma
    I:=[1,2,3,6,9,12,18,24]; [n le 8 select I[n] else 2*Self(n-1)-Self(n-2)+2*Self(n-3)-4*Self(n-4)+2*Self(n-5)-Self(n-6)+2*Self(n-7)-Self(n-8): n in [1..60]]; // Vincenzo Librandi, Mar 28 2015
    
  • Maple
    n:=80; lambda:=3; S10b:=[];
    for ii from 0 to n do
    x:=floor(ii/lambda);
    snc:=1/6*(x+1)*(x+2)*(3*ii-2*x*lambda+3);
    S10b:=[op(S10b),snc];
    od:
    S10b;
    A144677 := proc(n) option remember; local k1; sum(A190717(n-k1),k1=0..2) end: A190717:= proc(n) option remember; A190717(n):= binomial(floor(n/3)+3,3) end: seq(A144677(n), n=0..53); # Johannes W. Meijer, May 20 2011
  • Mathematica
    CoefficientList[Series[1/((x - 1)^4*(x^2 + x + 1)^2), {x, 0, 50}], x] (* Wesley Ivan Hurt, Mar 28 2015 *)
    LinearRecurrence[{2, -1, 2, -4, 2, -1, 2, -1}, {1, 2, 3, 6, 9, 12, 18, 24}, 60 ] (* Vincenzo Librandi, Mar 28 2015 *)
  • Sage
    @CachedFunction
    def a(n): return sum( ((j+3)//3)*((n-j+3)//3) for j in (0..n) )
    [a(n) for n in (0..60)] # G. C. Greubel, Oct 18 2021

Formula

From Johannes W. Meijer, May 20 2011: (Start)
a(n) = A190717(n-2) + A190717(n-1) + A190717(n).
a(n-2) + a(n-1) + a(n) = A014125(n).
G.f.: 1/((1-x)^4*(1+x+x^2)^2). (End)
From Wesley Ivan Hurt, Mar 28 2015: (Start)
a(n) = 2*a(n-1) -a(n-2) +2*a(n-3) -4*a(n-4) +2*a(n-5) -a(n-6) +2*a(n-7) -a(n-8).
a(n) = ((2 + floor(n/3))^3 - floor((n+4)/3) + floor((n+4)/3)^3 - floor((n+5)/3) + floor((n+5)/3)^3 - floor((n+6)/3))/6. (End)
a(n) = Sum_{j=0..n} floor((j+3)/3)*floor((n-j+3)/3). - G. C. Greubel, Oct 18 2021
a(n) = (132+129*n+36*n^2+3*n^3+6*(n+5)*cos(2*n*Pi/3)+2*sqrt(3)*(3*n+11)*sin(2*n*Pi/3))/162. - Wesley Ivan Hurt, Sep 01 2025

A194005 Triangle of the coefficients of an (n+1)-th order differential equation associated with A103631.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 6, 3, 1, 1, 6, 5, 10, 6, 4, 1, 1, 7, 6, 15, 10, 10, 4, 1, 1, 8, 7, 21, 15, 20, 10, 5, 1, 1, 9, 8, 28, 21, 35, 20, 15, 5, 1, 1, 10, 9, 36, 28, 56, 35, 35, 15, 6, 1, 1, 11, 10, 45, 36, 84, 56, 70, 35, 21, 6, 1
Offset: 0

Views

Author

Johannes W. Meijer and A. Hirschberg (a.hirschberg(AT)tue.nl), Aug 11 2011

Keywords

Comments

This triangle is a companion to Parks' triangle A103631.
The coefficients of triangle A103631(n,k) appear in appendix 2 of Park’s remarkable article “A new proof of the Routh-Hurwitz stability criterion using the second method of Liapunov” if we assume that the b(n) coefficients are all equal to 1, see the second Maple program.
The a(n,k) coefficients of the triangle given above are related to the coefficients of a linear (n+1)-th order differential equation for the case b(n)=1, see the examples.
a(n,k) is also the number of symmetric binary strings of odd length n with Hamming weight k>0 and no consecutive 1's. - Christian Barrientos and Sarah Minion, Feb 27 2018

Examples

			For the 5th-order linear differential equation the coefficients a(k) are: a(0) = 1, a(1) = a(4,0) = 1, a(2) = a(4,1) = 4, a(3) = a(4,2) = 3, a(4) = a(4,3) = 3 and a(5) = a(4,4) = 1.
The corresponding Hurwitz matrices A(k) are, see Parks: A(5) = Matrix([[a(1),a(0),0,0,0], [a(3),a(2),a(1),a(0),0], [a(5),a(4),a(3),a(2),a(1)], [0,0,a(5),a(4),a(3)], [0,0,0,0,a(5)]]), A(4) = Matrix([[a(1),a(0),0,0], [a(3),a(2),a(1),a(0)], [a(5),a(4),a(3),a(2)], [0,0,a(5),a(4)]]), A(3) = Matrix([[a(1),a(0),0], [a(3),a(2),a(1)], [a(5),a(4),a(3)]]), A(2) = Matrix([[a(1),a(0)], [a(3),a(2)]]) and A(1) = Matrix([[a(1)]]).
The values of b(k) are, see Parks: b(1) = d(1), b(2) = d(2)/d(1), b(3) = d(3)/(d(1)*d(2)), b(4) = d(1)*d(4)/(d(2)*d(3)) and b(5) = d(2)*d(5)/(d(3)*d(4)).
These a(k) values lead to d(k) = 1 and subsequently to b(k) = 1 and this confirms our initial assumption, see the comments.
'
Triangle starts:
  [0] 1;
  [1] 1, 1;
  [2] 1, 2, 1;
  [3] 1, 3, 2,  1;
  [4] 1, 4, 3,  3,  1;
  [5] 1, 5, 4,  6,  3,  1;
  [6] 1, 6, 5, 10,  6,  4,  1;
  [7] 1, 7, 6, 15, 10, 10,  4,  1;
  [8] 1, 8, 7, 21, 15, 20, 10,  5, 1;
  [9] 1, 9, 8, 28, 21, 35, 20, 15, 5, 1;
		

Crossrefs

Cf. A065941 and A103631.
Triangle sums (see A180662): A000071 (row sums; alt row sums), A075427 (Kn22), A000079 (Kn3), A109222(n+1)-1 (Kn4), A000045 (Fi1), A034943 (Ca3), A001519 (Gi3), A000930 (Ze3)
Interesting diagonals: T(n,n-4) = A189976(n+5) and T(n,n-5) = A189980(n+6)
Cf. A052509.

Programs

  • Haskell
    a194005 n k = a194005_tabl !! n !! k
    a194005_row n = a194005_tabl !! n
    a194005_tabl = [1] : [1,1] : f [1] [1,1] where
       f row' row = rs : f row rs where
         rs = zipWith (+) ([0,1] ++ row') (row ++ [0])
    -- Reinhard Zumkeller, Nov 22 2012
  • Maple
    A194005 := proc(n, k): binomial(floor((2*n+1-k)/2), n-k) end:
    for n from 0 to 11 do seq(A194005(n, k), k=0..n) od;
    seq(seq(A194005(n,k), k=0..n), n=0..11);
    nmax:=11: for n from 0 to nmax+1 do b(n):=1 od:
    A103631 := proc(n,k) option remember: local j: if k=0 and n=0 then b(1)
    elif k=0 and n>=1 then 0 elif k=1 then b(n+1) elif k=2 then b(1)*b(n+1)
    elif k>=3 then expand(b(n+1)*add(procname(j,k-2), j=k-2..n-2)) fi: end:
    for n from 0 to nmax do for k from 0 to n do
    A194005(n,k):= add(A103631(n1,k), n1=k..n) od: od:
    seq(seq(A194005(n,k),k=0..n), n=0..nmax);
  • Mathematica
    Flatten[Table[Binomial[Floor[(2n+1-k)/2],n-k],{n,0,20},{k,0,n}]] (* Harvey P. Dale, Apr 15 2012 *)

Formula

T(n,k) = binomial(floor((2*n+1-k)/2), n-k).
T(n,k) = sum(A103631(n1,k), n1=k..n), 0<=k<=n and n>=0.
T(n,k) = sum(binomial(floor((2*n1-k-1)/2), n1-k), n1=k..n).
T(n,0) = T(n,n) = 1, T(n,k) = T(n-2,k-2) + T(n-1,k), 0 < k < n. - Reinhard Zumkeller, Nov 23 2012

A026003 a(n) = T([n/2],[(n+1)/2]), where T = Delannoy triangle (A008288).

Original entry on oeis.org

1, 1, 3, 5, 13, 25, 63, 129, 321, 681, 1683, 3653, 8989, 19825, 48639, 108545, 265729, 598417, 1462563, 3317445, 8097453, 18474633, 45046719, 103274625, 251595969, 579168825, 1409933619, 3256957317, 7923848253, 18359266785, 44642381823
Offset: 0

Views

Author

Keywords

Comments

Number of lattice paths from (0,0) to the line x=n consisting of U=(1,1), D=(1,-1) and H=(2,0) steps and never going below the x-axis (i.e. left factors of Schroeder paths); for example, a(3)=5, counting the paths UUU,UUD,UDU,HU and UH. - Emeric Deutsch, Oct 27 2002
Transform of A001405 by |A049310(n,k)|, that is, transform of central binomial coefficients C(n,floor(n/2)) by Chebyshev mapping which takes a sequence with g.f. g(x) to the sequence with g.f. (1/(1-x^2))g(x/(1-x^2)). - Paul Barry, Jul 30 2005
The Kn1p sums, p >= 1, see A180662, of the Schroeder triangle A033877 (offset 0) are all related to A026003, e.g. Kn11(n) = A026003(n), Kn12(n) = A026003(n+2) - 1, Kn13(n) = A026003(n+4) - (2*n+7), Kn14(n) = A026003(n+6) - (2*n^2+18*n+41), Kn15(n) = A026003(n+8) - (4*n^3+66*n^2+368*n+693)/3, etc.. - Johannes W. Meijer, Jul 15 2013

References

  • L. Ericksen, Lattice path combinatorics for multiple product identities, J. Stat. Plan. Infer. 140 (2010) 2213-2226 doi:10.1016/j.jspi.2010.01.017
  • Lin Yang and S.-L. Yang, The parametric Pascal rhombus. Fib. Q., 57:4 (2019), 337-346.

Crossrefs

Bisections are the central Delannoy numbers A001850 and A002002 respectively.

Programs

  • Maple
    A026003 :=n -> add(binomial(n-k, k) * binomial(n-2*k, floor((n-2*k)/2)), k=0..floor(n/2)): seq(A026003(n), n=0..30); # Johannes W. Meijer, Jul 15 2013
  • Mathematica
    CoefficientList[Series[(Sqrt[(x^2-2*x-1)/(x^2+2*x-1)]-1)/2/x, {x, 0, 20}], x] (* Vaclav Kotesovec, Feb 13 2014 *)

Formula

G.f.: (sqrt((x^2-2*x-1)/(x^2+2*x-1))-1)/2/x. - Vladeta Jovovic, Apr 27 2003
a(n) = Sum_{k=0..floor(n/2)} C(n-k, k)*C(n-2k, floor((n-2k)/2)). - Paul Barry, Jul 30 2005
From Paul Barry, Mar 01 2010: (Start)
G.f.: 1/(1-x-2x^2/(1-x^2/(1-2x^2/(1-x^2/(1-2x^2/(1-... (continued fraction),
G.f.: 1/(1-x-x^2-x^2/(1-x^2-x^2/(1-x^2-x^2/(1-x^2-x^2/(1-... (continued fraction). (End)
D-finite with recurrence (n+1)*a(n) -2*a(n-1) +6*(-n+1)*a(n-2) -2*a(n-3) +(n-3)*a(n-4)=0. - R. J. Mathar, Nov 30 2012
a(n) ~ (1+sqrt(2))^(n+1) / (2^(3/4) * sqrt(Pi*n)). - Vaclav Kotesovec, Feb 13 2014

A070550 a(n) = a(n-1) + a(n-3) + a(n-4), starting with a(0..3) = 1, 2, 2, 3.

Original entry on oeis.org

1, 2, 2, 3, 6, 10, 15, 24, 40, 65, 104, 168, 273, 442, 714, 1155, 1870, 3026, 4895, 7920, 12816, 20737, 33552, 54288, 87841, 142130, 229970, 372099, 602070, 974170, 1576239, 2550408, 4126648, 6677057, 10803704, 17480760, 28284465, 45765226
Offset: 0

Views

Author

Sreyas Srinivasan (sreyas_srinivasan(AT)hotmail.com), May 02 2002

Keywords

Comments

Shares some properties with Fibonacci sequence.
The sum of any two alternating terms (terms separated by one other term) produces a Fibonacci number (e.g., 2+6=8, 3+10=13, 24+65=89). The product of any two consecutive or alternating Fibonacci terms produces a term from this sequence (e.g., 5*8 = 40, 13*5 = 65, 21*8 = 168).
In Penney's game (see A171861), the number of ways that HTH beats HHH on flip 3,4,5,... - Ed Pegg Jr, Dec 02 2010
The Ca2 sums (see A180662 for the definition of these sums) of triangle A035607 equal the terms of this sequence. - Johannes W. Meijer, Aug 05 2011

Examples

			G.f.: 1 + 2*x + 2*x^2 + 3*x^3 + 6*x^4 + 10*x^5 + 15*x^6 + 24*x^7 + ...
		

Crossrefs

Bisections: A001654, A059929.

Programs

  • Haskell
    a070550 n = a070550_list !! n
    a070550_list = 1 : 2 : 2 : 3 :
       zipWith (+) a070550_list
                   (zipWith (+) (tail a070550_list) (drop 3 a070550_list))
    -- Reinhard Zumkeller, Aug 06 2011
    
  • Maple
    with(combinat): A070550 := proc(n): fibonacci(floor(n/2)+1) * fibonacci(ceil(n/2)+2) end: seq(A070550(n),n=0..37); # Johannes W. Meijer, Aug 05 2011
  • Mathematica
    LinearRecurrence[{1, 0, 1, 1}, {1, 2, 2, 3}, 40] (* Jean-François Alcover, Jan 27 2018 *)
    nxt[{a_,b_,c_,d_}]:={b,c,d,a+b+d}; NestList[nxt,{1,2,2,3},40][[;;,1]] (* Harvey P. Dale, Jul 16 2024 *)
  • PARI
    A070550(n) = fibonacci(n\2+1)*fibonacci((n+5)\2) \\ M. F. Hasler, Aug 06 2011
    
  • PARI
    x='x+O('x^100); Vec((1+x)/(1-x-x^3-x^4)) \\ Altug Alkan, Dec 24 2015

Formula

a(n) = F(floor(n/2)+1)*F(ceiling(n/2)+2), with F(n) = A000045(n). - Ralf Stephan, Apr 14 2004
G.f.: (1+x)/(1-x-x^3-x^4) = (1+x)/((1+x^2)*(1-x-x^2))
a(n) = A126116(n+4) - F(n+3). - Johannes W. Meijer, Aug 05 2011
a(n) = (1+3*i)/10*(-i)^n + (1-3*i)/10*(i)^n + (2+sqrt(5))/5*((1+sqrt(5))/2)^n + (2-sqrt(5))/5*((1-sqrt(5))/2)^n, where i = sqrt(-1). - Sergei N. Gladkovskii, Jul 16 2013
a(n+1)*a(n+3) = a(n)*a(n+2) + a(n+1)*a(n+2) for all n in Z. - Michael Somos, Jan 19 2014
Sum_{n>=1} 1/a(n) = A290565. - Amiram Eldar, Feb 17 2021

A089068 a(n) = a(n-1)+a(n-2)+a(n-3)+2 with a(0)=0, a(1)=0 and a(2)=1.

Original entry on oeis.org

0, 0, 1, 3, 6, 12, 23, 43, 80, 148, 273, 503, 926, 1704, 3135, 5767, 10608, 19512, 35889, 66011, 121414, 223316, 410743, 755475, 1389536, 2555756, 4700769, 8646063, 15902590, 29249424, 53798079, 98950095, 181997600, 334745776, 615693473
Offset: 0

Views

Author

Roger L. Bagula, Dec 03 2003

Keywords

Comments

The a(n+2) represent the Kn12 and Kn22 sums of the square array of Delannoy numbers A008288. See A180662 for the definition of these knight and other chess sums. - Johannes W. Meijer, Sep 21 2010

Crossrefs

Cf. A000073 (Kn11 & Kn21), A089068 (Kn12 & Kn22), A180668 (Kn13 & Kn23), A180669 (Kn14 & Kn24), A180670 (Kn15 & Kn25). - Johannes W. Meijer, Sep 21 2010

Programs

  • Mathematica
    Join[{a=0,b=0,c=1},Table[d=a+b+c+2;a=b;b=c;c=d,{n,50}]] (* Vladimir Joseph Stephan Orlovsky, Apr 19 2011 *)
    RecurrenceTable[{a[0]==a[1]==0,a[2]==1,a[n]==a[n-1]+a[n-2]+a[n-3]+2}, a[n],{n,40}] (* or *) LinearRecurrence[{2,0,0,-1},{0,0,1,3},40] (* Harvey P. Dale, Sep 19 2011 *)

Formula

a(n) = A008937(n-2)+A008937(n-1). - Johannes W. Meijer, Sep 21 2010
a(n) = A018921(n-5)+A018921(n-4), n>4. - Johannes W. Meijer, Sep 21 2010
a(n) = A000073(n+2)-1. - R. J. Mathar, Sep 22 2010
From Johannes W. Meijer, Sep 22 2010: (Start)
a(n) = a(n-1)+A001590(n+1).
a(n) = Sum_{m=0..n} A040000(m)*A000073(n-m).
a(n+2) = Sum_{k=0..floor(n/2)} A008288(n-k+1,k+1).
G.f. = x^2*(1+x)/((1-x)*(1-x-x^2-x^3)). (End)
a(n) = 2*a(n-1)-a(n-4), a(0)=0, a(1)=0, a(2)=1, a(3)=3. - Bruno Berselli, Sep 23 2010

Extensions

Corrected and information added by Johannes W. Meijer, Sep 22 2010, Oct 22 2010
Definition based on arbitrarily set floating-point precision removed by R. J. Mathar, Sep 30 2010

A006603 Generalized Fibonacci numbers.

Original entry on oeis.org

1, 2, 7, 26, 107, 468, 2141, 10124, 49101, 242934, 1221427, 6222838, 32056215, 166690696, 873798681, 4612654808, 24499322137, 130830894666, 702037771647, 3783431872018, 20469182526595, 111133368084892, 605312629105205, 3306633429423460, 18111655081108453
Offset: 0

Views

Author

Keywords

Comments

The Kn21 sums, see A180662, of the Schroeder triangle A033877 equal A006603(n) while the Kn3 sums equal A006603(2*n). The Kn22 sums, see A227504, and the Kn23 sums, see A227505, are also related to the sequence given above. - Johannes W. Meijer, Jul 15 2013
Typo on the right-hand side of Rogers's equation (1-x+x^2+x^3)*R^*(x) = R(x) + x: the sign in front of the x should be switched. - R. J. Mathar, Nov 23 2018

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 50);
    Coefficients(R!( (1-x-2*x^2 -Sqrt(1-6*x+x^2))/(2*x*(1-x+x^2+x^3)) )); // G. C. Greubel, Oct 27 2024
    
  • Maple
    A006603 := n-> add((k*add(binomial(n-k+2, i)*binomial(2*n-3*k-i+3, n-k+1), i= 0.. n-2*k+2))/(n-k+2), k= 1.. n/2+1): seq(A006603(n), n=0..24); # Johannes W. Meijer, Jul 15 2013
  • Mathematica
    CoefficientList[Series[(1-x-2x^2-Sqrt[1-6x+x^2])/(2x(1-x+x^2+x^3)),{x,0,30}],x] (* Harvey P. Dale, Jun 12 2016 *)
  • Maxima
    a(n):=sum((k*sum(binomial(n-k+2,i)*binomial(2*n-3*k-i+3,n-k+1),i,0,n-2*k+2))/(n-k+2),k,1,n/2+1); /* Vladimir Kruchinin, Oct 23 2011 */
    
  • SageMath
    def A006603_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( (1-x-2*x^2 -sqrt(1-6*x+x^2))/(2*x*(1-x+x^2+x^3)) ).list()
    A006603_list(50) # G. C. Greubel, Oct 27 2024

Formula

a(n) = abs(A080244(n-1)).
G.f.: (1 - x - 2*x^2 - sqrt(1 - 6*x + x^2))/(2*x*(1 - x + x^2 + x^3)).
G.f.: (A006318(x) - x)/(1 - x + x^2 + x^3).
a(n) = Sum_{k=1..floor(n/2)+1} k*(1/(n-k+2))*Sum_{i=0..n-2*k+2} C(n-k+2,i)*C(2*n-3*k-i+3,n-k+1). - Vladimir Kruchinin, Oct 23 2011
(n+1)*a(n) -(7*n-2)*a(n-1) +4*(2*n-1)*a(n-2) -6*(n-1)*a(n-3) -(5*n-1)*a(n-4) +(n-2)*a(n-5) = 0. - R. J. Mathar, Nov 23 2018

Extensions

More terms from Emeric Deutsch, Feb 28 2004

A052544 Expansion of (1-x)^2/(1 - 4*x + 3*x^2 - x^3).

Original entry on oeis.org

1, 2, 6, 19, 60, 189, 595, 1873, 5896, 18560, 58425, 183916, 578949, 1822473, 5736961, 18059374, 56849086, 178955183, 563332848, 1773314929, 5582216355, 17572253481, 55315679788, 174128175064, 548137914373, 1725482812088
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

Equals INVERT transform of (1, 1, 3, 8, 21, 55, 144, ...). - Gary W. Adamson, May 01 2009
The Ze2 sums, see A180662, of triangle A065941 equal the terms (doubled) of this sequence. - Johannes W. Meijer, Aug 16 2011
Equals the partial sums of A052529 starting (1, 1, 4, 13, 41, 129, ...). - Gary W. Adamson, Feb 15 2012
First trisection of Narayana's cows sequence A000930. - Oboifeng Dira, Aug 03 2016
From Peter Bala, Nov 03 2017: (Start)
Let f(x) = x/(1 - x^3), the characteristic function of numbers of the form 3*n + 1. Then f(f(x)) = Sum_{n >= 0} a(n)*x^(3*n+1).
a(n) = the number of compositions of 3*n + 1 into parts of the form 3*m + 1. For example, a(2) = 6 and the six compositions of 7 into parts of the form 3*m + 1 are 7, 4 + 1 + 1 + 1, 1 + 4 + 1 + 1, 1 + 1 + 4 + 1, 1 + 1 + 1 + 4 and 1 + 1 + 1 + 1 + 1 + 1 + 1. Cf. A001519, which gives the number of compositions of an odd number into odd parts. (End)
a(n-1) is the number of permutations of length n avoiding the partially ordered pattern (POP) {1>2, 1>3, 4>2} of length 4. That is, number of length n permutations having no subsequences of length 4 in which the first element is larger than the second and third elements, and the fourth element is larger than the second element. - Sergey Kitaev, Dec 09 2020

Examples

			G.f. = 1 + 2*x + 6*x^2 + 19*x^3 + 60*x^4 + 189*x^5 + 595*x^6 + 1873*x^7 + ...
		

Crossrefs

Cf. A124820 (partial sums).

Programs

  • GAP
    a:=[1,2,6];; for n in [4..30] do a[n]:=4*a[n-1]-3*a[n-2]+a[n-3]; od; a; # G. C. Greubel, May 09 2019
  • Magma
    I:=[1, 2, 6]; [n le 3 select I[n] else 4*Self(n-1)-3*Self(n-2) +Self(n-3): n in [1..30]]; // Vincenzo Librandi, Jan 12 2012
    
  • Maple
    spec := [S,{S=Sequence(Union(Z,Prod(Z,Sequence(Z),Sequence(Z))))}, unlabeled]: seq(combstruct[count](spec,size=n), n=0..25);
    A052544 := proc(n): add(binomial(n+2*k, 3*k), k=0...n)  end: seq(A052544(n), n=0..25); # Johannes W. Meijer, Aug 16 2011
  • Mathematica
    LinearRecurrence[{4,-3,1},{1,2,6},30] (* Harvey P. Dale, Jul 13 2011 *)
    Table[Sum[Binomial[n + 2 k, 3 k], {k, 0, n}], {n, 0, 30}] (* or *)
    CoefficientList[Series[(1-x)^2/(1-4x+3x^2-x^3), {x, 0, 30}], x] (* Michael De Vlieger, Aug 03 2016 *)
  • PARI
    {a(n) = sum(k=0, n, binomial(n + 2*k, 3*k))}; /* Michael Somos, Jan 12 2012 */
    
  • PARI
    Vec((1-x)^2/(1-4*x+3*x^2-x^3)+O(x^30)) \\ Charles R Greathouse IV, Jan 12 2012
    
  • Sage
    ((1-x)^2/(1-4*x+3*x^2-x^3)).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, May 09 2019
    

Formula

G.f.: (1-x)^2/(1 -4*x +3*x^2 -x^3).
a(n) = 4*a(n-1) - 3*a(n-2) + a(n-3).
a(n) = Sum(-1/31*(-4-7*_alpha+2*_alpha^2)*_alpha^(-1-n), _alpha=RootOf(-1+4*_Z-3*_Z^2+_Z^3)).
a(n) = Sum_{k=0..n} binomial(n+2*k, 3*k). - Richard L. Ollerton, May 12 2004
G.f.: 1 / (1 - x - x / (1 - x)^2). - Michael Somos, Jan 12 2012
a(n) = hypergeom([(n+1)/2, n/2+1, -n], [1/3, 2/3], -4/27). - Peter Luschny, Nov 03 2017

Extensions

More terms from James Sellers, Jun 06 2000

A080247 Formal inverse of triangle A080246. Unsigned version of A080245.

Original entry on oeis.org

1, 2, 1, 6, 4, 1, 22, 16, 6, 1, 90, 68, 30, 8, 1, 394, 304, 146, 48, 10, 1, 1806, 1412, 714, 264, 70, 12, 1, 8558, 6752, 3534, 1408, 430, 96, 14, 1, 41586, 33028, 17718, 7432, 2490, 652, 126, 16, 1, 206098
Offset: 0

Views

Author

Paul Barry, Feb 15 2003

Keywords

Comments

Row sums are little Schroeder numbers A001003. Diagonal sums are generalized Fibonacci numbers A006603. Columns include A006318, A006319, A006320, A006321.
T(n,k) is the number of dissections of a convex (n+3)-gon by nonintersecting diagonals with exactly k diagonals emanating from a fixed vertex. Example: T(2,1)=4 because the dissections of the convex pentagon ABCDE having exactly one diagonal emanating from the vertex A are: {AC}, {AD}, {AC,EC} and {AD,BD}. - Emeric Deutsch, May 31 2004
For more triangle sums, see A180662, see the Schroeder triangle A033877 which is the mirror of this triangle. - Johannes W. Meijer, Jul 15 2013

Examples

			Triangle starts:
[0]    1
[1]    2,    1
[2]    6,    4,   1
[3]   22,   16,   6,   1
[4]   90,   68,  30,   8,  1
[5]  394,  304, 146,  48, 10,  1
[6] 1806, 1412, 714, 264, 70, 12, 1
...
From _Gary W. Adamson_, Jul 25 2011: (Start)
n-th row = top row of M^n, M = the following infinite square production matrix:
  2, 1, 0, 0, 0, ...
  2, 2, 1, 0, 0, ...
  2, 2, 2, 1, 0, ...
  2, 2, 2, 2, 1, ...
  ... (End)
		

Crossrefs

Cf. A000007, A033877 (mirror), A084938.

Programs

  • Maple
    A080247:=(n,k)->(k+1)*add(binomial(n+1,k+j+1)*binomial(n+j,j),j=0..n-k)/(n+1):
    seq(seq(A080247(n,k),k=0..n),n=0..9);
  • Mathematica
    Clear[w] w[n_, k_] /; k < 0 || k > n := 0 w[0,0]=1 ; w[n_, k_] /; 0 <= k <= n && !n == k == 0 := w[n, k] = w[n-1,k-1] + w[n-1,k] + w[n,k+1] Table[w[n,k],{n,0,10},{k,0,n}] (* David Callan, Jul 03 2006 *)
    T[n_, k_] := Binomial[n, k] Hypergeometric2F1[k - n, n + 1, k + 2, -1];
    Table[T[n, k], {n, 0, 8}, {k, 0, n}] // Flatten (* Peter Luschny, Jan 08 2018 *)
  • Maxima
    T(n,k):=((k+1)*sum(2^m*binomial(n+1,m)*binomial(n-k-1,n-k-m),m,0,n-k))/(n+1); /* Vladimir Kruchinin, Jan 10 2022 */
  • Sage
    def A080247_row(n):
        @cached_function
        def prec(n, k):
            if k==n: return 1
            if k==0: return 0
            return prec(n-1,k-1)-2*sum(prec(n,k+i-1) for i in (2..n-k+1))
        return [(-1)^(n-k)*prec(n, k) for k in (1..n)]
    for n in (1..10): print(A080247_row(n)) # Peter Luschny, Mar 16 2016
    

Formula

G.f.: 2/(2+y*x-y+y*(x^2-6*x+1)^(1/2))/y/x. - Vladeta Jovovic, Feb 16 2003
Essentially same triangle as triangle T(n,k), n > 0 and k > 0, read by rows; given by [0, 2, 1, 2, 1, 2, 1, 2, 1, 2, ...] DELTA A000007 where DELTA is Deléham's operator defined in A084938.
T(n, k) = T(n-1, k-1) + 2*Sum_{j>=0} T(n-1, k+j) with T(0, 0) = 1 and T(n, k)=0 if k < 0. - Philippe Deléham, Jan 19 2004
T(n, k) = (k+1)*Sum_{j=0..n-k} (binomial(n+1, k+j+1)*binomial(n+j, j))/(n+1). - Emeric Deutsch, May 31 2004
Recurrence: T(0,0)=1; T(n,k) = T(n-1,k-1) + T(n-1,k) + T(n,k+1). - David Callan, Jul 03 2006
T(n, k) = binomial(n, k)*hypergeom([k - n, n + 1], [k + 2], -1). - Peter Luschny, Jan 08 2018
T(n,k) = (k+1)/(n+1)*Sum_{m=0..n-k} 2^m*binomial(n+1,m)*binomial(n-k-1,n-k-m). - Vladimir Kruchinin, Jan 10 2022
From Peter Bala, Sep 16 2024: (Start)
Riordan array (S(x), x*S(x)), where S(x) = (1 - x - sqrt(1 - 6*x + x^2))/(2*x) is the g.f. of the large Schröder numbers A006318.
For integer m and n >= 1, (m + 2)*[x^n] S(x)^(m*n) = m*[x^n] (1/S(-x))^((m+2)*n). For cases of this identity see A103885 (m = 1), A333481 (m = 2) and A370102 (m = 3). (End)

A094967 Right-hand neighbors of Fibonacci numbers in Stern's diatomic series.

Original entry on oeis.org

1, 1, 2, 2, 5, 5, 13, 13, 34, 34, 89, 89, 233, 233, 610, 610, 1597, 1597, 4181, 4181, 10946, 10946, 28657, 28657, 75025, 75025, 196418, 196418, 514229, 514229, 1346269, 1346269, 3524578, 3524578, 9227465, 9227465, 24157817, 24157817, 63245986, 63245986, 165580141, 165580141
Offset: 0

Views

Author

Paul Barry, May 26 2004

Keywords

Comments

Fibonacci(2*n+1) repeated. a(n) is the right neighbor of Fibonacci(n+2) in A049456 and A002487 (starts 2,2,5,...). A000045(n+2) = A094966(n) + a(n).
Diagonal sums of A109223. - Paul Barry, Jun 22 2005
The Fi2 sums, see A180662, of triangle A065941 equal the terms of this sequence. - Johannes W. Meijer, Aug 11 2011
a(n) is the last term of (n+1)-th row in Wythoff array A003603. -Reinhard Zumkeller, Jan 26 2012

Crossrefs

Programs

  • GAP
    List([0..50], n -> Fibonacci(n)*(1-(-1)^n)/2 + Fibonacci(n+1)*(1+(-1)^n)/2); # G. C. Greubel, Nov 18 2018
  • Magma
    [IsEven(n) select Fibonacci(n+1) else Fibonacci(n): n in [0..70]]; // Vincenzo Librandi, Nov 18 2018
    
  • Maple
    A094967 := proc(n) combinat[fibonacci](2*floor(n/2)+1) ; end proc: seq(A094967(n), n=0..37);
  • Mathematica
    LinearRecurrence[{0,3,0,-1},{1,1,2,2},40] (* Harvey P. Dale, Apr 05 2015 *)
    f[n_]:=If[OddQ@n, (Fibonacci[n]), Fibonacci[n+1]]; Array[f, 100, 0] (* Vincenzo Librandi, Nov 18 2018 *)
    Table[Fibonacci[n, 0]*Fibonacci[n] + LucasL[n, 0]*Fibonacci[n + 1]/2, {n, 0, 50}] (* G. C. Greubel, Nov 18 2018 *)
  • PARI
    vector(50, n, n--; fibonacci(n)*(1-(-1)^n)/2 + fibonacci(n+1)*(1+(-1)^n)/2) \\ G. C. Greubel, Nov 18 2018
    
  • Sage
    [fibonacci(n)*(1-(-1)^n)/2 + fibonacci(n+1)*(1+(-1)^n)/2 for n in range(50)] # G. C. Greubel, Nov 18 2018
    

Formula

G.f.: (1+x-x^2-x^3)/(1-3*x^2+x^4).
a(n) = Fibonacci(n)*(1-(-1)^n)/2 + Fibonacci(n+1)*(1+(-1)^n)/2.
a(n) = Sum_{k=0..floor(n/2)} binomial(floor(n/2)+k, 2*k). - Paul Barry, Jun 22 2005
Starting (1, 2, 2, 5, 5, 13, 13, ...) = A133080 * A000045, where A000045 starts with "1". - Gary W. Adamson, Sep 08 2007
a(n) = Fibonacci(n+1)^(4*k+3) mod Fibonacci(n+2), for any k>-1, n>0. - Gary Detlefs, Nov 29 2010
Previous Showing 61-70 of 142 results. Next