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.

Showing 1-10 of 18 results. Next

A091698 Matrix inverse of triangle A063967.

Original entry on oeis.org

1, -1, 1, 1, -3, 1, -1, 8, -5, 1, 1, -23, 19, -7, 1, -1, 74, -69, 34, -9, 1, 1, -262, 256, -147, 53, -11, 1, -1, 993, -986, 615, -265, 76, -13, 1, 1, -3943, 3935, -2571, 1235, -431, 103, -15, 1, -1, 16178, -16169, 10862, -5591, 2216, -653, 134, -17, 1, 1
Offset: 0

Views

Author

Christian G. Bower, Jan 29 2004

Keywords

Comments

Riordan array (1/(1+x), (sqrt(1+6x+5x^2)-x-1)/(2(1+x))). The absolute value array is (1/(1-x),xc(x)/(1-xc(x))) where c(x) is the g.f. of A000108. It factorizes as (1/(1-x),x/(1-x))(1,xc(x)). - Paul Barry, Jun 10 2005

Examples

			From _Paul Barry_, Apr 15 2010: (Start)
Triangle begins
  1,
  -1, 1,
  1, -3, 1,
  -1, 8, -5, 1,
  1, -23, 19, -7, 1,
  -1, 74, -69, 34, -9, 1,
  1, -262, 256, -147, 53, -11, 1,
  -1, 993, -986, 615, -265, 76, -13, 1,
  1, -3943, 3935, -2571, 1235, -431, 103, -15, 1
Production matrix begins
  -1, 1,
  0, -2, 1,
  0, 1, -2, 1,
  0, -1, 1, -2, 1,
  0, 1, -1, 1, -2, 1,
  0, -1, 1, -1, 1, -2, 1,
  0, 1, -1, 1, -1, 1, -2, 1,
  0, -1, 1, -1, 1, -1, 1, -2, 1,
  0, 1, -1, 1, -1, 1, -1, 1, -2, 1,
  0, -1, 1, -1, 1, -1, 1, -1, 1, -2, 1 (End)
		

Crossrefs

Row sums: A091699. Row sums (absolute values): A007317. Column 1: A050511.

Programs

  • Mathematica
    rows = 11; t[n_, k_] := Sum[Binomial[j, n - j]*Binomial[j, k], {j, 0, n}]; T = Table[t[n, k], {n, 0, rows - 1}, {k, 0, rows - 1}] // Inverse; Table[ T[[n, k]], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Oct 11 2017 *)

A091700 Matrix square of triangle A063967.

Original entry on oeis.org

1, 2, 1, 7, 6, 1, 23, 29, 10, 1, 78, 127, 67, 14, 1, 264, 527, 375, 121, 18, 1, 895, 2113, 1892, 831, 191, 22, 1, 3034, 8269, 8922, 4973, 1559, 277, 26, 1, 10286, 31781, 40115, 27139, 10826, 2623, 379, 30, 1, 34872, 120448, 174080, 138617, 67308, 20763, 4087, 497, 34, 1
Offset: 0

Views

Author

Christian G. Bower, Jan 29 2004

Keywords

Examples

			   1;
   2,   1;
   7,   6,  1;
  23,  29, 10,  1;
  78, 127, 67, 14, 1;
  ...
		

Crossrefs

Row sums: A091701.
Column 0: A091702.

A002605 a(n) = 2*(a(n-1) + a(n-2)), a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 2, 6, 16, 44, 120, 328, 896, 2448, 6688, 18272, 49920, 136384, 372608, 1017984, 2781184, 7598336, 20759040, 56714752, 154947584, 423324672, 1156544512, 3159738368, 8632565760, 23584608256, 64434348032, 176037912576, 480944521216, 1313964867584
Offset: 0

Views

Author

Keywords

Comments

Individually, both this sequence and A028859 are convergents to 1 + sqrt(3). Mutually, both sequences are convergents to 2 + sqrt(3) and 1 + sqrt(3)/2. - Klaus E. Kastberg (kastberg(AT)hotkey.net.au), Nov 04 2001
The number of (s(0), s(1), ..., s(n+1)) such that 0 < s(i) < 6 and |s(i) - s(i-1)| <= 1 for i = 1, 2, ..., n + 1, s(0) = 2, s(n+1) = 3. - Herbert Kociemba, Jun 02 2004
The same sequence may be obtained by the following process. Starting a priori with the fraction 1/1, the denominators of fractions built according to the rule: add top and bottom to get the new bottom, add top and 4 times the bottom to get the new top. The limit of the sequence of fractions is sqrt(4). - Cino Hilliard, Sep 25 2005
The Hankel transform of this sequence is [1, 2, 0, 0, 0, 0, 0, 0, 0, ...]. - Philippe Deléham, Nov 21 2007
[1, 3; 1, 1]^n *[1, 0] = [A026150(n), a(n)]. - Gary W. Adamson, Mar 21 2008
(1 + sqrt(3))^n = A026150(n) + a(n)*sqrt(3). - Gary W. Adamson, Mar 21 2008
a(n+1) is the number of ways to tile a board of length n using red and blue tiles of length one and two. - Geoffrey Critzer, Feb 07 2009
Starting with offset 1 = INVERT transform of the Jacobsthal sequence, A001045: (1, 1, 3, 5, 11, 21, ...). - Gary W. Adamson, May 12 2009
Starting with "1" = INVERTi transform of A007482: (1, 3, 11, 39, 139, ...). - Gary W. Adamson, Aug 06 2010
An elephant sequence, see A175654. For the corner squares four A[5] vectors, with decimal values 85, 277, 337 and 340, lead to this sequence (without the leading 0). For the central square these vectors lead to the companion sequence A026150, without the first leading 1. - Johannes W. Meijer, Aug 15 2010
The sequence 0, 1, -2, 6, -16, 44, -120, 328, -896, ... (with alternating signs) is the Lucas U(-2,-2)-sequence. - R. J. Mathar, Jan 08 2013
a(n+1) counts n-walks (closed) on the graph G(1-vertex;1-loop,1-loop,2-loop,2-loop). - David Neil McGrath, Dec 11 2014
Number of binary strings of length 2*n - 2 in the regular language (00+11+0101+1010)*. - Jeffrey Shallit, Dec 14 2015
For n >= 1, a(n) equals the number of words of length n - 1 over {0, 1, 2, 3} in which 0 and 1 avoid runs of odd lengths. - Milan Janjic, Dec 17 2015
a(n+1) is the number of compositions of n into parts 1 and 2, both of two kinds. - Gregory L. Simay, Sep 20 2017
Number of associative, quasitrivial, and order-preserving binary operations on the n-element set {1, ..., n} that have neutral elements. - J. Devillet, Sep 28 2017
(1 + sqrt(3))^n = A026150(n) + a(n)*sqrt(3), for n >= 0; integers in the real quadratic number field Q(sqrt(3)). - Wolfdieter Lang, Feb 10 2018
Starting with 1, 2, 6, 16, ..., number of permutations of length n>0 avoiding the partially ordered pattern (POP) {1>3, 1>4} 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 third and fourth elements. - Sergey Kitaev, Dec 09 2020
a(n) is the number of tilings of a 2 X n board missing one corner cell, with 1 X 1 and L-shaped tiles (where the L-shaped tiles cover 3 squares). Compare to A127864. - Greg Dresden and Yilin Zhu, Jul 17 2025

References

  • John Derbyshire, Prime Obsession, Joseph Henry Press, April 2004, p. 16.

Crossrefs

First differences are given by A026150.
a(n) = A073387(n, 0), n>=0 (first column of triangle).
Equals (1/3) A083337. First differences of A077846. Pairwise sums of A028860 and abs(A077917).
a(n) = A028860(n)/2 apart from the initial terms.
Row sums of A081577 and row sums of triangle A156710.
The following sequences (and others) belong to the same family: A001333, A000129, A026150, A046717, A015518, A084057, A063727, A002533, A002532, A083098, A083099, A083100, A015519.
Cf. A175289 (Pisano periods).
Cf. A002530.
Cf. A127864.

Programs

  • Haskell
    a002605 n = a002605_list !! n
    a002605_list =
       0 : 1 : map (* 2) (zipWith (+) a002605_list (tail a002605_list))
    -- Reinhard Zumkeller, Oct 15 2011
    
  • Magma
    [Floor(((1 + Sqrt(3))^n - (1 - Sqrt(3))^n)/(2*Sqrt(3))): n in [0..30]]; // Vincenzo Librandi, Aug 18 2011
    
  • Magma
    [n le 2 select n-1 else 2*Self(n-1) + 2*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 07 2018
  • Maple
    a[0]:=0:a[1]:=1:for n from 2 to 50 do a[n]:=2*a[n-1]+2*a[n-2]od: seq(a[n], n=0..33); # Zerinvary Lajos, Dec 15 2008
    a := n -> `if`(n<3, n, 2^(n-1)*hypergeom([1-n/2, (1-n)/2], [1-n], -2));
    seq(simplify(a(n)), n=0..29); # Peter Luschny, Dec 16 2015
  • Mathematica
    Expand[Table[((1 + Sqrt[3])^n - (1 - Sqrt[3])^n)/(2Sqrt[3]), {n, 0, 30}]] (* Artur Jasinski, Dec 10 2006 *)
    a[n_]:=(MatrixPower[{{1,3},{1,1}},n].{{1},{1}})[[2,1]]; Table[a[n],{n,-1,40}] (* Vladimir Joseph Stephan Orlovsky, Feb 19 2010 *)
    LinearRecurrence[{2, 2}, {0, 1}, 30] (* Robert G. Wilson v, Apr 13 2013 *)
    Round@Table[Fibonacci[n, Sqrt[2]] 2^((n - 1)/2), {n, 0, 20}] (* Vladimir Reshetnikov, Oct 15 2016 *)
    nxt[{a_,b_}]:={b,2(a+b)}; NestList[nxt,{0,1},30][[All,1]] (* Harvey P. Dale, Sep 17 2022 *)
  • PARI
    Vec(x/(1-2*x-2*x^2)+O(x^99)) \\ Charles R Greathouse IV, Jun 10 2011
    
  • PARI
    A002605(n)=([2,2;1,0]^n)[2,1] \\ M. F. Hasler, Aug 06 2018
    
  • Sage
    [lucas_number1(n,2,-2) for n in range(0, 30)] # Zerinvary Lajos, Apr 22 2009
    
  • Sage
    a = BinaryRecurrenceSequence(2,2)
    print([a(n) for n in (0..29)])  # Peter Luschny, Aug 29 2016
    

Formula

a(n) = (-I*sqrt(2))^(n-1)*U(n-1, I/sqrt(2)) where U(n, x) is the Chebyshev U-polynomial. - Wolfdieter Lang
G.f.: x/(1 - 2*x - 2*x^2).
From Paul Barry, Sep 17 2003: (Start)
E.g.f.: x*exp(x)*(sinh(sqrt(3)*x)/sqrt(3) + cosh(sqrt(3)*x)).
a(n) = (1 + sqrt(3))^(n-1)*(1/2 + sqrt(3)/6) + (1 - sqrt(3))^(n-1)*(1/2 - sqrt(3)/6), for n>0.
Binomial transform of 1, 1, 3, 3, 9, 9, ... Binomial transform is A079935. (End)
a(n) = Sum_{k=0..floor(n/2)} binomial(n - k, k)*2^(n - k). - Paul Barry, Jul 13 2004
a(n) = A080040(n) - A028860(n+1). - Creighton Dement, Jan 19 2005
a(n) = Sum_{k=0..n} A112899(n,k). - Philippe Deléham, Nov 21 2007
a(n) = Sum_{k=0..n} A063967(n,k). - Philippe Deléham, Nov 03 2006
a(n) = ((1 + sqrt(3))^n - (1 - sqrt(3))^n)/(2*sqrt(3)).
a(n) = Sum_{k=0..n} binomial(n, 2*k + 1) * 3^k.
Binomial transform of expansion of sinh(sqrt(3)x)/sqrt(3) (0, 1, 0, 3, 0, 9, ...). E.g.f.: exp(x)*sinh(sqrt(3)*x)/sqrt(3). - Paul Barry, May 09 2003
a(n) = (1/3)*Sum_{k=1..5} sin(Pi*k/2)*sin(2*Pi*k/3)*(1 + 2*cos(Pi*k/6))^n, n >= 1. - Herbert Kociemba, Jun 02 2004
a(n+1) = ((3 + sqrt(3))*(1 + sqrt(3))^n + (3 - sqrt(3))*(1 - sqrt(3))^n)/6. - Al Hakanson (hawkuu(AT)gmail.com), Jun 29 2009
Antidiagonals sums of A081577. - J. M. Bergot, Dec 15 2012
G.f.: Q(0)*x/2, where Q(k) = 1 + 1/(1 - x*(4*k + 2 + 2*x)/(x*(4*k + 4 + 2*x) + 1/Q(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 30 2013
a(n) = 2^(n - 1)*hypergeom([1 - n/2, (1 - n)/2], [1 - n], -2) for n >= 3. - Peter Luschny, Dec 16 2015
Sum_{k=0..n} a(k)*2^(n-k) = a(n+2)/2 - 2^n. - Greg Dresden, Feb 11 2022
a(n) = 2^floor(n/2) * A002530(n). - Gregory L. Simay, Sep 22 2022
From Peter Bala, May 08 2024: (Start)
G.f.: x/(1 - 2*x - 2*x^2) = Sum_{n >= 0} x^(n+1) *( Product_{k = 1..n} (k + 2*x + 1)/(1 + k*x) )
Also x/(1 - 2*x - 2*x^2) = Sum_{n >= 0} (2*x)^n *( x*Product_{k = 1..n} (m*k + 2 - m + x)/(1 + 2*m*k*x) ) for arbitrary m (both series are telescoping). (End)
a(n) = A127864(n-1) + A127864(n-2). - Greg Dresden and Yilin Zhu, Jul 17 2025

Extensions

Edited by N. J. A. Sloane, Apr 15 2009

A002478 Bisection of A000930.

Original entry on oeis.org

1, 1, 3, 6, 13, 28, 60, 129, 277, 595, 1278, 2745, 5896, 12664, 27201, 58425, 125491, 269542, 578949, 1243524, 2670964, 5736961, 12322413, 26467299, 56849086, 122106097, 262271568, 563332848, 1209982081, 2598919345, 5582216355, 11990037126, 25753389181
Offset: 0

Views

Author

Keywords

Comments

Number of ways to tile a 3 X n region with 1 X 1, 2 X 2 and 3 X 3 tiles.
Number of ternary words with subwords (0,0), (0,1) and (1,1) not allowed. - Olivier Gérard, Aug 28 2012
Diagonal sums of A063967. - Paul Barry, Nov 09 2005
Row sums of number triangle A116088. - Paul Barry, Feb 04 2006
Sequence is identical to its second differences negated, minus the first 3 terms. - Paul Curtz, Feb 10 2008
a(n) = term (3,3) in the 3 X 3 matrix [0,1,0; 0,0,1; 1,2,1]^n. - Gary W. Adamson, May 30 2008
a(n)/a(n-1) tends to 2.147899035..., an eigenvalue of the matrix and a root to x^3 - x^2 - 2x - 1 = 0. - Gary W. Adamson, May 30 2008
INVERT transform of (1, 2, 1, 0, 0, 0, ...) = (1, 3, 6, 13, 28, ...); such that (1, 2, 1, 0, 0, 0, ...) convolved with (1, 1, 3, 6, 13, 28, 0, 0, 0, ...) shifts to the left. - Gary W. Adamson, Apr 18 2010
a(n) is the top left entry in the n-th power of the 3 X 3 matrix [1, 1, 1; 1, 0, 1; 1, 0, 0] or of the 3 X 3 matrix [1, 1, 1; 1, 0, 0; 1, 1, 0]. - R. J. Mathar, Feb 03 2014

Examples

			a(3)=6 as there is one tiling of a 3 X 3 region with only 1 X 1 tiles, 4 tilings with exactly one 2 X 2 tile and 1 tiling consisting of the 3 X 3 tile.
		

References

  • Kenneth Edwards, Michael A. Allen, A new combinatorial interpretation of the Fibonacci numbers squared, Part II, Fib. Q., 58:2 (2020), 169-177.
  • L. Euler, (E388) Vollstaendige Anleitung zur Algebra, Zweiter Theil, reprinted in: Opera Omnia. Teubner, Leipzig, 1911, Series (1), Vol. 1, p. 322.
  • S. Heubach, Tiling an m X n Area with Squares of Size up to k X k (m<=5), Congressus Numerantium 140 (1999), pp. 43-64.
  • 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. A000930, A054856, A054857, A025234, A078007, A078039, A226546, A077936 (INVERT transform), A008346 (inverse INVERT transform).

Programs

  • Magma
    I:=[1,1,3]; [n le 3 select I[n] else Self(n-1) +2*Self(n-2) +Self(n-3): n in [1..41]]; // G. C. Greubel, Apr 14 2023
    
  • Mathematica
    f[A_]:= Module[{til = A}, AppendTo[til, A[[-1]] + 2A[[-2]] + A[[-3]]]]; NumOfTilings[ n_ ]:= Nest[ f, {1,1,3}, n-2]; NumOfTilings[30]
    LinearRecurrence[{1,2,1},{1,1,3},40] (* Vladimir Joseph Stephan Orlovsky, Jan 28 2012 *)
    CoefficientList[Series[1/(1-x-2x^2-x^3),{x,0,40}],x] (* Harvey P. Dale, Oct 17 2024 *)
  • PARI
    a(n)=([0,1,0; 0,0,1; 1,2,1]^n*[1;1;3])[1,1] \\ Charles R Greathouse IV, Apr 08 2016
    
  • SageMath
    @CachedFunction
    def a(n): # A002478
        if (n<3): return (1,1,3)[n]
        else: return sum(binomial(2,j)*a(n-j) for j in range(1,4))
    [a(n) for n in (0..40)] # G. C. Greubel, Apr 14 2023

Formula

G.f.: 1 / (1-x-2*x^2-x^3). [Simon Plouffe in his 1992 dissertation.]
a(n) = a(n-1) + 2*a(n-2) + a(n-3).
a(n) = Sum_{k=0..n} binomial(2*n-2*k, k). - Paul Barry, Nov 13 2004
a(n) = Sum_{k=0..floor(n/2)} Sum_{j=0..n-k} C(j, n-k-j)*C(j, k). - Paul Barry, Nov 09 2005
a(n) = Sum_{k=0..n} C(2*k,n-k) = Sum_{k=0..n} C(n,k)*C(3*k,n)/C(3*k,k). - Paul Barry, Feb 04 2006
a(n) = A000930(n) + 2*Sum_{i=0..n-2} a(i)*A000930(n-2-i). - Michael Tulskikh, Jun 07 2020

Extensions

Additional comments from Silvia Heubach (silvi(AT)cine.net), Apr 21 2000

A030195 a(n) = 3*a(n-1) + 3*a(n-2), a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 3, 12, 45, 171, 648, 2457, 9315, 35316, 133893, 507627, 1924560, 7296561, 27663363, 104879772, 397629405, 1507527531, 5715470808, 21668995017, 82153397475, 311467177476, 1180861724853, 4476986706987, 16973545295520
Offset: 0

Views

Author

Keywords

Comments

Scaled Chebyshev U-polynomials evaluated at I*sqrt(3)/2.
Number of zeros in the substitution system {0 -> 1111100, 1 -> 10} at step n from initial string "1" (1 -> 10 -> 101111100 -> ...). - Ilya Gutkovskiy, Apr 10 2017
a(n+1) is the number of compositions of n having parts 1 and 2, both of three kinds. - Gregory L. Simay, Sep 21 2017
More generally, define a(n) = k*a(n-1) + k*a(n-2), a(0) = 0 and a(1) = 1. Then g.f. a(n) = 1/(1 - k*x - k*x^2) and a(n+1) is the number of compositions of n having parts 1 and 2, both of k kinds. - Gregory L. Simay, Sep 22 2017

Examples

			G.f. = x + 3*x^2 + 12*x^3 + 45*x^4 + 171*x^5 + 648*x^6 + 2457*x^7 + ...
		

Crossrefs

Programs

  • Haskell
    a030195 n = a030195_list !! n
    a030195_list =
       0 : 1 : map (* 3) (zipWith (+) a030195_list (tail a030195_list))
    -- Reinhard Zumkeller, Oct 14 2011
    
  • Magma
    I:=[0,1]; [n le 2 select I[n] else 3*Self(n-1) + 3*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 24 2018
  • Mathematica
    CoefficientList[Series[1/(1-3x-3x^2), {x, 0, 25}], x] (* Zerinvary Lajos, Mar 22 2007 *)
    LinearRecurrence[{3, 3}, {0, 1}, 24] (* Or *)
    RecurrenceTable[{a[n] == 3 a[n - 1] + 3 a[n - 2], a[0] == 0, a[1] == 1}, a, {n, 0, 23}] (* Robert G. Wilson v, Aug 18 2012 *)
  • PARI
    {a(n) = n--; polchebyshev(n, 2, I*sqrt(3)/2) * (-I*sqrt(3))^n};
    
  • Sage
    [lucas_number1(n,3,-3) for n in range(0, 25)] # Zerinvary Lajos, Apr 22 2009
    

Formula

a(n+1) = (-I*sqrt(3))^n*U(n, I*sqrt(3)/2).
G.f.: x / (1 - 3*x - 3*x^2).
a(n+1) = Sum_{k=0..floor(n/2)} 3^(n-k)*binomial(n-k, k). - Emeric Deutsch, Nov 14 2001
a(n) = (p^n - q^n)/sqrt(21); p = (3 + sqrt 21)/2, q = (3 - sqrt 21)/2. - Gary W. Adamson, Jul 02 2003
For n > 0, a(n) = Sum_{k=0..n-1} (2^k)*A063967(n-1,k). - Gerald McGarvey, Jul 23 2006
a(n+1) = Sum_{k=0..n} 2^k*A063967(n,k). - Philippe Deléham, Nov 03 2006

Extensions

Edited by Ralf Stephan, Aug 02 2004
I simplified the definition. As a result the offsets in some of the formulas may need to shifted by 1. - N. J. A. Sloane, Apr 01 2006
Formulas shifted to match offset. - Charles R Greathouse IV, Jan 31 2011

A321620 The Riordan square of the Riordan numbers, triangle read by rows, T(n, k) for 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 3, 2, 1, 1, 3, 5, 5, 3, 1, 1, 6, 13, 10, 7, 4, 1, 1, 15, 29, 26, 16, 9, 5, 1, 1, 36, 73, 61, 42, 23, 11, 6, 1, 1, 91, 181, 157, 103, 61, 31, 13, 7, 1, 1, 232, 465, 398, 271, 156, 83, 40, 15, 8, 1, 1, 603, 1205, 1040, 702, 419, 221, 108, 50, 17, 9, 1, 1
Offset: 0

Views

Author

Peter Luschny, Nov 15 2018

Keywords

Comments

If gf is a generating function of the sequence a then by the 'Riordan square of a' we understand the integer triangle given by the Riordan array (gf, gf). This mapping can be seen as an operator RS: Z[[x]] -> Mat[Z].
For instance A039599 is the Riordan square of the Catalan numbers, A172094 is the Riordan square of the little Schröder numbers and A063967 is the Riordan square of the Fibonacci numbers. The Riordan square of the simplest sequence of positive numbers (A000012) is the Pascal triangle.
The generating functions used are listed in the table below. They may differ slightly from those defined elsewhere in order to ensure that RS(a) is invertible (as a matrix). We understand this as a kind of normalization.
---------------------------------------------------------------------------
Sequence a | RS(a) | gf(a)
---------------------------------------------------------------------------
Catalan | A039599 | (1 - sqrt(1 - 4*x))/(2*x).
1, Riordan | A321620 | 1 + 2*x/(1 + x + sqrt(1 - 2*x - 3*x^2)).
Motzkin | A321621 | (1 - x - sqrt(1 - 2*x - 3*x^2))/(2*x^2).
Fine | A321622 | 1 + (1 - sqrt(1 - 4*x))/(3 - sqrt(1 - 4*x)).
large Schröder | A321623 | (1 - x - sqrt(1 - 6*x + x^2))/(2*x).
little Schröder | A172094 | (1 + x - sqrt(1 - 6*x + x^2))/(4*x).
Lucas | A321624 | 1 + x*(1 + 2*x)/(1 - x - x^2).
swinging factorial | A321625 | (1 + x/(1 - 4*x^2))/sqrt(1 - 4*x^2).
ternary trees ||A109956|| u = sqrt(x*3/4); sin(arcsin(3*u)/3)/u.
central trinomial | A116392 | 1/sqrt(1 - 2*x - 3*x^2)
Bell | A154380 | Sum_{k>=0} x^k/Product_{j=1..k}(1 - j*x).
(2*n-1)!! | A321627 | 1/(1-x/(1-2*x/(1-3*x/(1-4*x/(1-5*x/...
powers of 2 | A038208 | 1/(1 - 2*x).
the all 1's seq. | A007318 | 1/(1 - x).
Fibonacci | A063967 | 1/(1 - x - x^2).
tribonacci | A187889 | 1/(1 - x - x^2 - x^3).
tetranacci | A353593 | 1/(1 - x - x^2 - x^3 - x^4).
Jacobsthal | A322942 | (2*x^2-1)/((x + 1)*(2*x - 1))

Examples

			The triangle starts:
   [ 0]   1
   [ 1]   1    1
   [ 2]   0    1    1
   [ 3]   1    1    1   1
   [ 4]   1    3    2   1   1
   [ 5]   3    5    5   3   1   1
   [ 6]   6   13   10   7   4   1   1
   [ 7]  15   29   26  16   9   5   1  1
   [ 8]  36   73   61  42  23  11   6  1  1
   [ 9]  91  181  157 103  61  31  13  7  1 1
   [10] 232  465  398 271 156  83  40 15  8 1 1
		

Crossrefs

Programs

  • Maple
    RiordanSquare := proc(d, n, exp:=false) local td, M, k, m, u, j;
        series(d, x, n+1); td := [seq(coeff(%, x, j), j = 0..n)];
        M := Matrix(n); for k from 1 to n do M[k, 1] := td[k] od;
        for k from 1 to n-1 do for m from k to n-1 do
           M[m+1, k+1] := add(M[j, k]*td[m-j+2], j = k..m) od od;
        if exp then u := 1;
           for k from 1 to n-1 do u := u * k;
              for m from 1 to k do j := `if`(m = 1, u, j/(m-1));
                 M[k+1, m] := M[k+1, m] * j od od fi;
    M end:
    RiordanSquare(1 + 2*x/(1 + x + sqrt(1 - 2*x - 3*x^2)), 8);
  • Mathematica
    RiordanSquare[gf_, len_] := Module[{T}, T[n_, k_] := T[n, k] = If[k == 0, SeriesCoefficient[gf, {x, 0, n}], Sum[T[j, k-1] T[n-j, 0], {j, k-1, n-1}]]; Table[T[n, k], {n, 0, len-1}, {k, 0, n}]];
    M = RiordanSquare[1 + 2x/(1 + x + Sqrt[1 - 2x - 3x^2]), 12];
    M // Flatten (* Jean-François Alcover, Nov 24 2018 *)
  • Sage
    # uses[riordan_array from A256893]
    def riordan_square(gf, len, exp=false):
        return riordan_array(gf, gf, len, exp)
    riordan_square(1 + 2*x/(1 + x + sqrt(1 - 2*x - 3*x^2)), 10)
    # Alternatively, given a list S:
    def riordan_square_array(S):
        N = len(S)
        M = matrix(ZZ, N, N)
        for n in (0..N-1): M[n, 0] = S[n]
        for k in (1..N-1):
            for m in (k..N-1):
                M[m, k] = sum(M[j, k-1]*S[m-j] for j in (k-1..m-1))
        return M
    riordan_square_array([1, 1, 0, 1, 1, 3, 6, 15, 36]) # Peter Luschny, Apr 03 2020

Formula

Given a generating function g and an positive integer N compute the Taylor expansion at the origin t(k) = [x^k] g(x) for k in [0...N-1] and set T(n, 0) = t(n) for n in [0...N-1]. Then compute T(m, k) = Sum_{j in [k-1...m-1]} T(j, k - 1) t(m - j) for k in [1...N-1] and for m in [k...N-1]. The resulting (0, 0)-based lower triangular array is the Riordan square generated by g.

A057087 Scaled Chebyshev U-polynomials evaluated at i. Generalized Fibonacci sequence.

Original entry on oeis.org

1, 4, 20, 96, 464, 2240, 10816, 52224, 252160, 1217536, 5878784, 28385280, 137056256, 661766144, 3195289600, 15428222976, 74494050304, 359689093120, 1736732573696, 8385686667264, 40489676963840, 195501454524416
Offset: 0

Views

Author

Wolfdieter Lang, Aug 11 2000

Keywords

Comments

a(n) gives the length of the word obtained after n steps with the substitution rule 0->1111, 1->11110, starting from 0. The number of 1's and 0's of this word is 4*a(n-1) and 4*a(n-2), respectively.
Inverse binomial transform of odd Pell bisection A001653. With a leading zero, inverse binomial transform of even Pell bisection A001542, divided by 2. - Paul Barry, May 16 2003
For positive n, a(n) equals the permanent of the n X n tridiagonal matrix with 4's along the main diagonal, and 2's along the superdiagonal and the subdiagonal. - John M. Campbell, Jul 19 2011
Pisano period lengths: 1, 1, 8, 1, 3, 8, 6, 1, 24, 3, 120, 8, 21, 6, 24, 1, 16, 24, 360, 3, ... . - R. J. Mathar, Aug 10 2012
Exponential convolution of Pell numbers (A000129) and companion Pell numbers (A002203), divided by 2 and leading zero dropped. - Vladimir Reshetnikov, Oct 07 2016

Crossrefs

Pairwise sums are in A086347.
Appears in A086346, A086347 and A086348. - Johannes W. Meijer, Aug 01 2010

Programs

  • Magma
    I:=[1,4]; [n le 2 select I[n] else 4*Self(n-1) + 4*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 16 2018
  • Maple
    A057087 := n -> `if`(n=0, 1, 4^n*hypergeom([1/2-n/2, -n/2], [-n], -1)):
    seq(simplify(A057087(n)), n=0..21); # Peter Luschny, Dec 17 2015
  • Mathematica
    Table[Fibonacci[n + 1, 2] 2^n, {n, 0, 20}] (* Vladimir Reshetnikov, Oct 08 2016 *)
    LinearRecurrence[{4,4},{1,4},30] (* Harvey P. Dale, Aug 17 2017 *)
  • PARI
    a(n)=if(n<0, 0, (2*I)^n*subst(I*poltchebi(n+1)+poltchebi(n),'x,-I)/2) /* Michael Somos, Sep 16 2005 */
    
  • PARI
    Vec(1/(1-4*x-4*x^2) + O(x^100)) \\ Altug Alkan, Dec 17 2015
    
  • Sage
    [lucas_number1(n,4,-4) for n in range(1, 23)] # Zerinvary Lajos, Apr 23 2009
    

Formula

a(n) = 4*(a(n-1) + a(n-2)), a(-1)=0, a(0)=1.
G.f.: 1/(1 - 4*x - 4*x^2).
a(n) = S(n, 2*i)*(-2*i)^n with S(n, x) := U(n, x/2), Chebyshev's polynomials of the 2nd kind, A049310.
a(n) = Sum_{k=0..n} 3^k*A063967(n,k). - Philippe Deléham, Nov 03 2006
a(n) = A000129(n+1)*A000079(n). - R. J. Mathar, Jul 08 2009
From Johannes W. Meijer, Aug 01 2010: (Start)
Limit_{k->oo} a(n+k)/a(k) = A084128(n) + 2*A057087(n-1)*sqrt(2);
Limit_{n->oo} A084128(n)/A057087(n-1) = sqrt(2). (End)
a(n) = 4^n*hypergeom([1/2-n/2, -n/2], [-n], -1) for n>=1. - Peter Luschny, Dec 17 2015

A026729 Square array of binomial coefficients T(n,k) = binomial(n,k), n >= 0, k >= 0, read by downward antidiagonals.

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, Jan 19 2003

Keywords

Comments

The signed triangular matrix T(n,k)*(-1)^(n-k) is the inverse matrix of the triangular Catalan convolution matrix A106566(n,k), n=k>=0, with A106566(n,k) = 0 if nPhilippe Deléham, Aug 01 2005
As a number triangle: unsigned version of A109466. - Philippe Deléham, Oct 26 2008
A063967*A130595 as infinite lower triangular matrices. - Philippe Deléham, Dec 11 2008
Modulo 2, this sequence becomes A106344. - Philippe Deléham, Dec 18 2008
Let {a_(k,i)}, k>=1, i=0,...,k, be the k-th antidiagonal of the array. Then s_k(n) = Sum_{i=0..k}a_(k,i)* binomial(n,k) is the n-th element of the k-th column of A111808. For example, s_1(n) = binomial(n,1) = n is the first column of A111808 for n>1, s_2(n) = binomial(n,1) + binomial(n,2) is the second column of A111808 for n>1, etc. Therefore, in cases k=3,4,5,6,7,8, s_k(n) is A005581(n), A005712(n), A000574(n), A005714(n), A005715(n), A005716(n), respectively. Besides, s_k(n+5) = A064054(n). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
As a triangle, T(n,k) = binomial(k,n-k). - Peter Bala, Nov 27 2015
For all n >= 0, k >= 0, the k-th homology group of the n-torus H_k(T^n) is the free abelian group of rank T(n,k) = binomial(n,k). See the Math Stack Exchange link below. - Jianing Song, Mar 13 2023

Examples

			Array begins
  1 0 0 0 0 0 ...
  1 1 0 0 0 0 ...
  1 2 1 0 0 0 ...
  1 3 3 1 0 0 ...
  1 4 6 4 1 0 ...
As a triangle, this begins
  1
  0 1
  0 1 1
  0 0 2 1
  0 0 1 3 1
  0 0 0 3 4 1
  0 0 0 1 6 5 1
  ...
Production array is
  0    1
  0    1   1
  0   -1   1   1
  0    2  -1   1  1
  0   -5   2  -1  1  1
  0   14  -5   2 -1  1  1
  0  -42  14  -5  2 -1  1  1
  0  132 -42  14 -5  2 -1  1  1
  0 -429 132 -42 14 -5  2 -1  1  1
  ... (Cf. A000108)
		

Crossrefs

The official entry for Pascal's triangle is A007318. See also A052553 (the same array read by upward antidiagonals).
Cf. A030528 (subtriangle for 1<=k<=n).

Programs

  • GAP
    nmax:=15;; T:=List([0..nmax],n->List([0..nmax],k->Binomial(n,k)));;
    b:=List([2..nmax],n->OrderedPartitions(n,2));;
    a:=Flat(List([1..Length(b)],i->List([1..Length(b[i])],j->T[b[i][j][1]][b[i][j][2]]))); # Muniru A Asiru, Jul 17 2018
  • Magma
    /* As triangle: */ [[Binomial(k, n-k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Nov 29 2015
    
  • Maple
    seq(seq(binomial(k,n-k),k=0..n),n=0..12); # Peter Luschny, May 31 2014
  • Mathematica
    Table[Binomial[k, n - k], {n, 0, 12}, {k, 0, n}] // Flatten (* Michael De Vlieger, Nov 28 2015 *)

Formula

As a number triangle, this is defined by T(n,0) = 0^n, T(0,k) = 0^k, T(n,k) = T(n-1,k-1) + Sum_{j, j>=0} (-1)^j*T(n-1,k+j)*A000108(j) for n>0 and k>0. - Philippe Deléham, Nov 07 2005
As a triangle read by rows, it is [0, 1, -1, 0, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 22 2006
As a number triangle, this is defined by T(n, k) = Sum_{i=0..n} (-1)^(n+i)*binomial(n, i)*binomial(i+k, i-k) and is the Riordan array ( 1, x*(1+x) ). The row sums of this triangle are F(n+1). - Paul Barry, Jun 21 2004
Sum_{k=0..n} x^k*T(n,k) = A000007(n), A000045(n+1), A002605(n), A030195(n+1), A057087(n), A057088(n), A057089(n), A057090(n), A057091(n), A057092(n), A057093(n) for n=0,1,2,3,4,5,6,7,8,9,10. - Philippe Deléham, Oct 16 2006
T(n,k) = A109466(n,k)*(-1)^(n-k). - Philippe Deléham, Dec 11 2008
G.f. for the triangular interpretation: -1/(-1+x*y+x^2*y). - R. J. Mathar, Aug 11 2015
For T(0,0) = 0, the triangle below has the o.g.f. G(x,t) = [t*x(1+x)]/[1-t*x(1+x)]. See A109466 for a signed version and inverse, A030528 for reverse and A102426 for a shifted version. - Tom Copeland, Jan 19 2016

A057088 Scaled Chebyshev U-polynomials evaluated at i*sqrt(5)/2. Generalized Fibonacci sequence.

Original entry on oeis.org

1, 5, 30, 175, 1025, 6000, 35125, 205625, 1203750, 7046875, 41253125, 241500000, 1413765625, 8276328125, 48450468750, 283633984375, 1660422265625, 9720281250000, 56903517578125, 333118994140625, 1950112558593750, 11416157763671875, 66831351611328125, 391237546875000000
Offset: 0

Views

Author

Wolfdieter Lang, Aug 11 2000

Keywords

Comments

a(n) gives the length of the word obtained after n steps with the substitution rule 0->11111, 1->111110, starting from 0. The number of 1's and 0's of this word is 5*a(n-1) and 5*a(n-2), resp.
a(n) / a(n-1) converges to (5 + (3 * sqrt(5))) / 2 as n approaches infinity. (5 + (3 * sqrt(5))) / 2 can also be written as phi^2 + (2 * phi), phi^3 + phi, phi + sqrt(5) + 2, (3 * phi) + 1, (3 * phi^2) - 2, phi^4 - 1 and (5 + (3 * (L(n) / F(n)))) / 2, where L(n) is the n-th Lucas number and F(n) is the n-th Fibonacci number as n approaches infinity. - Ross La Haye, Aug 18 2003, on another version
Pisano period lengths: 1, 3, 3, 6, 1, 3, 24, 12, 9, 3, 10, 6, 56, 24, 3, 24,288, 9, 18, 6, ... - R. J. Mathar, Aug 10 2012

Crossrefs

Programs

  • Magma
    I:=[1, 5]; [n le 2 select I[n] else 5*Self(n-1) + 5*Self(n-2): n in [0..30]]; // G. C. Greubel, Jan 16 2018
  • Maple
    a[0]:=0:a[1]:=1:for n from 2 to 50 do a[n]:=5*a[n-1]+5*a[n-2]od: seq(a[n], n=1..33); # Zerinvary Lajos, Dec 14 2008
  • Mathematica
    LinearRecurrence[{5,5}, {1,5}, 30] (* G. C. Greubel, Jan 16 2018 *)
  • PARI
    x='x+O('x^30); Vec(1/(1 - 5*x - 5*x^2)) \\ G. C. Greubel, Jan 16 2018
    
  • Sage
    [lucas_number1(n,5,-5) for n in range(1, 22)] # Zerinvary Lajos, Apr 24 2009
    

Formula

a(n) = 5*(a(n-1) + a(n-2)), a(-1)=0, a(0)=1.
a(n) = S(n, i*sqrt(5))*(-i*sqrt(5))^n with S(n, x) := U(n, x/2), Chebyshev's polynomials of the 2nd kind, A049310.
G.f.: 1/(1 - 5*x - 5*x^2).
a(n) = (1/3)*Sum_{k=0..n} binomial(n, k)*Fibonacci(k)*3^k. - Benoit Cloitre, Oct 25 2003
a(n) = ((5 + 3*sqrt(5))/2)^n(1/2 + sqrt(5)/6) + (1/2 - sqrt(5)/6)((5 - 3*sqrt(5))/2)^n. - Paul Barry, Sep 22 2004
(a(n)) appears to be given by the floretion - 0.75'i - 0.5'j + 'k - 0.75i' + 0.5j' + 0.5k' + 1.75'ii' - 1.25'jj' + 1.75'kk' - 'ij' - 0.5'ji' - 0.75'jk' - 0.75'kj' - 1.25e ("jes"). - Creighton Dement, Nov 28 2004
a(n) = Sum_{k=0..n} 4^k*A063967(n,k). - Philippe Deléham, Nov 03 2006
G.f.: G(0)/(2-5*x), where G(k)= 1 + 1/(1 - x*(9*k-5)/(x*(9*k+4) - 2/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 17 2013
From Ehren Metcalfe, Nov 18 2017: (Start)
With F(n) = A000045(n), L(n) = A000032(n), beta = (1-sqrt(5))/2:
a(2*n-1) = 5^n*F(4*n)/3 = (5^(n-1/2)*L(4*n) - 2*5^(n-1/2)*beta^(4*n))/3.
a(2*n) = 5^n*L(4*n+2)/3 = (5^(n+1/2)*F(4*n+2) + 2*5^n*beta^(4*n+2))/3.
a(n) = round 5^((n+1)/2)*F(2*(n+1))/3.
a(n) = round 5^(n/2)*L(2*(n+1))/3. (End)

A023610 Convolution of Fibonacci numbers and {F(2), F(3), F(4), ...}.

Original entry on oeis.org

1, 3, 7, 15, 30, 58, 109, 201, 365, 655, 1164, 2052, 3593, 6255, 10835, 18687, 32106, 54974, 93845, 159765, 271321, 459743, 777432, 1312200, 2211025, 3719643, 6248479, 10482351, 17562870, 29391490, 49132669, 82048737, 136884293, 228160495, 379975140, 632293452
Offset: 0

Views

Author

Keywords

Comments

a(n-2) + 1 is the number of (3412,1243)-, (3412,2134)- and (3412,1324)-avoiding involutions in S_n, n>1. - Ralf Stephan, Jul 06 2003
The number of terms in all ordered partitions of (n+1) using only ones and twos. For example, a(3)=15 because there are 15 terms in 1+1+1+1;2+1+1;1+2+1;1+1+2;2+2 - Geoffrey Critzer, Apr 07 2008
a(n) is the number of n-matchings in the graph obtained by a zig-zag triangulation of a convex (2n+1)-gon. Example: a(2)=7 because in the triangulation of the convex pentagon ABCDEA with diagonals AD and AC we have 7 2-matchings: {AB,CD},{AB,DE},{BC,AD},{BC,DE},{BC,EA},{CD,EA} and {DE,AC}. - Emeric Deutsch, Dec 25 2004
Partial sums of A029907. First differences of A002940. - Peter Bala, Oct 24 2007
Equals row sums of triangle A144154. - Gary W. Adamson, Sep 12 2008
Equals the number of 1's in Fibonacci Maximal notation for subsets of
(1, 2, 3, 5, 8, 13, ...) terms. For example (cf. A181630): 4, 5, and 6 are the 3 terms 101, 110, and 111 in Fibonacci Maximal. Total number of 1's for those terms = 7 = a(2). - Gary W. Adamson, Nov 02 2010
a(n) is half the number of strokes needed to draw all the domino tilings of a 2 X (n+2) rectangle. - Roberto Tauraso, Mar 15 2014
a(n) is the total number of 1's in all (n+1)-bit dual Zeckendorf representations of integers (A104326). For example, a(2) = 7 counts the 1's in 101, 110, 111. - Shenghui Yang, Feb 09 2025

Crossrefs

Cf. A000045 (Fibonacci numbers).
Column 1 of triangle A063967.

Programs

  • Haskell
    a023610 n = a023610_list !! n
    a023610_list = f [1] $ drop 3 a000045_list where
       f us (v:vs) = (sum $ zipWith (*) us $ tail a000045_list) : f (v:us) vs
    -- Reinhard Zumkeller, Jan 18 2014
    
  • Mathematica
    Table[Sum[Binomial[n - i, i]*(n - i), {i, 0, n}], {n, 1, 33}] (* Geoffrey Critzer, May 04 2009 *)
  • PARI
    a(n)=(n+2)*fibonacci(n+4)/5+(n-1)*fibonacci(n+2)/5 \\ Charles R Greathouse IV, Jun 11 2015
  • Sage
    def A023610():
        a, b, c, d = 1, 3, 7, 15
        while True:
            yield a
            a, b, c, d = b, c, d, 2*(d-b)+c-a
    a = A023610(); [next(a) for i in range(33)]  # Peter Luschny, Nov 20 2013
    

Formula

O.g.f.: (x+1)/(1-x-x^2)^2. - Len Smiley, Dec 11 2001
a(n) = (1/5)*((n+2)*F(n+4) + (n-1)*F(n+2)), with F(n)=A000045(n). - Ralf Stephan, Jul 06 2003
a(n) = Sum_{k=0..n+1} (n-k+1)*binomial(n-k+1, k). - Paul Barry, Nov 05 2005
Recurrence: a(n+2) = a(n+1) + a(n) + Fib(n+4), n >= 0. For n >= 2, a(n-2) = (-1)^n*((-2n+3)*Fib(-n) - (-n)*Fib(-n-1))/5 = (-1)^n*A010049(-n), the second-order Fibonacci numbers of negative index, where Fib(-n) = (-1)^(n+1)*Fib(n). - Peter Bala, Oct 24 2007
a(n) = (n+1)*F(n+2) - A001629(n+1) where F(n) is the n-th Fibonacci number. - Geoffrey Critzer, Apr 07 2008
a(n) = 2*a(n-1) + a(n-2) - 2*a(n-3) - a(n-4), n >= 4. - L. Edson Jeffery, Mar 29 2013
a(n+1) = A004798(n) + A000045(n+2) for n >= 0. - John Molokach, Jul 04 2013
a(n) = A001629(n+1) + A001629(n+2). - Philippe Deléham, Oct 30 2013
E.g.f.: exp(x/2)*(5*(5 + 7*x)*cosh(sqrt(5)*x/2) + sqrt(5)*(11 + 15*x)*sinh(sqrt(5)*x/2))/25. - Stefano Spezia, Dec 04 2023
Showing 1-10 of 18 results. Next