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 11-20 of 30 results. Next

A180226 a(n) = 4*a(n-1) + 10*a(n-2), with a(1)=0 and a(2)=1.

Original entry on oeis.org

0, 1, 4, 26, 144, 836, 4784, 27496, 157824, 906256, 5203264, 29875616, 171535104, 984896576, 5654937344, 32468715136, 186424233984, 1070384087296, 6145778689024, 35286955629056, 202605609406464, 1163291993916416, 6679224069730304, 38349816218085376
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    I:=[0,1]; [n le 2 select I[n] else 4*Self(n-1) + 10*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 16 2018
  • Mathematica
    Join[{a=0,b=1},Table[c=4*b+10*a;a=b;b=c,{n,100}]]
    LinearRecurrence[{4,10}, {0,1}, 30] (* G. C. Greubel, Jan 16 2018 *)
  • PARI
    x='x+O('x^30); concat([0], Vec(x^2/(1-4*x-10*x^2))) \\ G. C. Greubel, Jan 16 2018
    

Formula

a(n) = ((2+sqrt(14))^(n-1) - (2-sqrt(14))^(n-1))/(2*sqrt(14)). - Rolf Pleisch, May 14 2011
G.f.: x^2/(1-4*x-10*x^2).

A135030 Generalized Fibonacci numbers: a(n) = 6*a(n-1) + 2*a(n-2).

Original entry on oeis.org

0, 1, 6, 38, 240, 1516, 9576, 60488, 382080, 2413456, 15244896, 96296288, 608267520, 3842197696, 24269721216, 153302722688, 968355778560, 6116740116736, 38637152257536, 244056393778688, 1541612667187200
Offset: 0

Views

Author

Rolf Pleisch, Feb 10 2008, Feb 14 2008

Keywords

Comments

For n>0, a(n) equals the number of words of length n-1 over {0,1,...,7} in which 0 and 1 avoid runs of odd lengths. - Milan Janjic, Jan 08 2017

Crossrefs

Programs

  • Magma
    [n le 2 select n-1 else 6*Self(n-1) + 2*Self(n-2): n in [1..35]]; // Vincenzo Librandi, Sep 18 2016
    
  • Maple
    A:= gfun:-rectoproc({a(0) = 0, a(1) = 1, a(n) = 2*(3*a(n-1) + a(n-2))},a(n),remember):
    seq(A(n),n=1..30); # Robert Israel, Sep 16 2014
  • Mathematica
    Join[{a=0,b=1},Table[c=6*b+2*a;a=b;b=c,{n,100}]] (* Vladimir Joseph Stephan Orlovsky, Jan 16 2011 *)
    LinearRecurrence[{6,2},{0,1},30] (* or *) CoefficientList[Series[ -(x/(2x^2+6x-1)),{x,0,30}],x] (* Harvey P. Dale, Jun 20 2011 *)
  • PARI
    a(n)=([0,1; 2,6]^n*[0;1])[1,1] \\ Charles R Greathouse IV, Oct 03 2016
  • Sage
    [lucas_number1(n,6,-2) for n in range(0, 21)] # Zerinvary Lajos, Apr 24 2009
    

Formula

a(0) = 0; a(1) = 1; a(n) = 2*(3*a(n-1) + a(n-2)).
a(n) = 1/(2*sqrt(11))*( (3 + sqrt(11))^n - (3 - sqrt(11))^n ).
G.f.: x/(1 - 6*x - 2*x^2). - Harvey P. Dale, Jun 20 2011
a(n+1) = Sum_{k=0..n} A099097(n,k)*2^k. - Philippe Deléham, Sep 16 2014
E.g.f.: (1/sqrt(11))*exp(3*x)*sinh(sqrt(11)*x). - G. C. Greubel, Sep 17 2016

Extensions

More terms from Joshua Zucker, Feb 23 2008

A083856 Square array T(n,k) of generalized Fibonacci numbers, read by antidiagonals upwards (n, k >= 0).

Original entry on oeis.org

0, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 2, 1, 0, 1, 1, 3, 3, 1, 0, 1, 1, 4, 5, 5, 1, 0, 1, 1, 5, 7, 11, 8, 1, 0, 1, 1, 6, 9, 19, 21, 13, 1, 0, 1, 1, 7, 11, 29, 40, 43, 21, 1, 0, 1, 1, 8, 13, 41, 65, 97, 85, 34, 1, 0, 1, 1, 9, 15, 55, 96, 181, 217, 171, 55, 1
Offset: 0

Views

Author

Paul Barry, May 06 2003

Keywords

Comments

Row n >= 0 of the array gives the solution to the recurrence b(k) = b(k-1) + n*b(k-2) for k >= 2 with b(0) = 0 and b(1) = 1.

Examples

			Array T(n,k) (with rows n >= 0 and columns k >= 0) begins as follows:
  0, 1, 1,  1,  1,   1,   1,    1,    1,     1, ... [A057427]
  0, 1, 1,  2,  3,   5,   8,   13,   21,    34, ... [A000045]
  0, 1, 1,  3,  5,  11,  21,   43,   85,   171, ... [A001045]
  0, 1, 1,  4,  7,  19,  40,   97,  217,   508, ... [A006130]
  0, 1, 1,  5,  9,  29,  65,  181,  441,  1165, ... [A006131]
  0, 1, 1,  6, 11,  41,  96,  301,  781,  2286, ... [A015440]
  0, 1, 1,  7, 13,  55, 133,  463, 1261,  4039, ... [A015441]
  0, 1, 1,  8, 15,  71, 176,  673, 1905,  6616, ... [A015442]
  0, 1, 1,  9, 17,  89, 225,  937, 2737, 10233, ... [A015443]
  0, 1, 1, 10, 19, 109, 280, 1261, 3781, 15130, ... [A015445]
  ...
		

Crossrefs

Rows include A057427 (n=0), A000045 (n=1), A001045 (n=2), A006130 (n=3), A006131 (n=4), A015440 (n=5), A015441 (n=6), A015442 (n=7), A015443 (n=8), A015445 (n=9).
Columns include A000012 (k=1,2), A000027 (k=3), A005408 (k=4), A028387 (k=5), A000567 (k=6), A106734 (k=7).
Cf. A083857 (binomial transform), A083859 (main diagonal), A083860 (first subdiagonal), A083861 (second binomial transform), A110112, A110113 (diagonal sums), A193376 (transposed variant), A172237 (transposed variant).

Programs

  • Julia
    function generalized_fibonacci(r, n)
       F = BigInt[1 r; 1 0]
       Fn = F^n
       Fn[2, 1]
    end
    for r in 0:6 println([generalized_fibonacci(r, n) for n in 0:9]) end # Peter Luschny, Mar 06 2017
  • Maple
    A083856_row := proc(r, n) local R; R := proc(n) option remember;
    if n<=1 then n else R(n-1)+r*R(n-2) fi end: R(n) end:
    for r from 0 to 9 do seq(A083856_row(r, n), n=0..9) od; # Peter Luschny, Mar 06 2017
  • Mathematica
    T[, 0] = 0; T[, 1|2] = 1; T[n_, k_] := T[n, k] = T[n, k-1] + n T[n, k-2];
    Table[T[n-k, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 22 2018 *)

Formula

T(n, k) = (((1 + sqrt(4*n + 1))/2)^k - ((1 - sqrt(4*n + 1))/2)^k)/sqrt(4*n + 1). [corrected by Michel Marcus, Jun 25 2018]
From Thomas Baruchel, Jun 25 2018: (Start)
The g.f. for row n >= 0 is x/(1 - x - n*x^2).
The g.f. for column k >= 1 is g(k,x) = 1/(1-x) + Sum_{m = 1..floor((k-1)/2)} (1 - x)^(-1 - m) * binomial(k - 1 - m, m) * Sum_{i = 0..m} x^i * Sum_{j = 0..i} (-1)^j * (i - j)^m * binomial(1 + m, j).
The g.f. for column k >= 1 is also g(k,x) = 1 + Sum_{m = 1..floor((k+1)/2)} ((1 - x)^(-m) * binomial(k-m, m-1) * Sum_{j = 0..m} (-1)^j * binomial(m, j) * x^m * Phi(x, -m+1, -j+m)) + Sum_{s = 0..floor((k-1)/2)} binomial(k-s-1, s) * PolyLog(-s, x), where Phi is the Lerch transcendent function. (End)
T(n,k) = Sum_{i = 0..k} (-1)^(k+i) * binomial(k,i) * A083857(n,i). - Petros Hadjicostas, Dec 24 2019

Extensions

Various sections edited by Petros Hadjicostas, Dec 24 2019

A180250 a(n) = 5*a(n-1) + 10*a(n-2), with a(1)=0 and a(2)=1.

Original entry on oeis.org

0, 1, 5, 35, 225, 1475, 9625, 62875, 410625, 2681875, 17515625, 114396875, 747140625, 4879671875, 31869765625, 208145546875, 1359425390625, 8878582421875, 57987166015625, 378721654296875, 2473479931640625, 16154616201171875, 105507880322265625
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [n le 2 select n-1 else 5*Self(n-1) +10*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 16 2018
    
  • Mathematica
    Join[{a=0,b=1},Table[c=5*b+10*a;a=b;b=c,{n,100}]]
    LinearRecurrence[{5,10}, {0,1}, 30] (* G. C. Greubel, Jan 16 2018 *)
  • PARI
    a(n)=([0,1;10,5]^(n-1))[1,2] \\ Charles R Greathouse IV, Oct 03 2016
    
  • PARI
    my(x='x+O('x^30)); concat([0], Vec(x^2/(1-5*x-10*x^2))) \\ G. C. Greubel, Jan 16 2018
    
  • SageMath
    A180250= BinaryRecurrenceSequence(5,10,0,1)
    [A180250(n-1) for n in range(1,41)] # G. C. Greubel, Jul 21 2023

Formula

a(n) = ((5+sqrt(65))^(n-1) - (5-sqrt(65))^(n-1))/(2^(n-1)*sqrt(65)). - Rolf Pleisch, May 14 2011
G.f.: x^2/(1-5*x-10*x^2).
a(n) = (i*sqrt(10))^(n-1) * ChebyshevU(n-1, -i*sqrt(5/8)). - G. C. Greubel, Jul 21 2023

A015551 Expansion of x/(1 - 6*x - 5*x^2).

Original entry on oeis.org

0, 1, 6, 41, 276, 1861, 12546, 84581, 570216, 3844201, 25916286, 174718721, 1177893756, 7940956141, 53535205626, 360916014461, 2433172114896, 16403612761681, 110587537144566, 745543286675801, 5026197405777636
Offset: 0

Views

Author

Keywords

Comments

Let the generator matrix for the ternary Golay G_12 code be [I|B], where the elements of B are taken from the set {0,1,2}. Then a(n)=(B^n)1,2 for instance. - _Paul Barry, Feb 13 2004
Pisano period lengths: 1, 2, 4, 4, 1, 4, 42, 8, 12, 2, 10, 4, 12, 42, 4, 16, 96, 12, 360, 4, ... - R. J. Mathar, Aug 10 2012

Crossrefs

Programs

  • Magma
    I:=[0,1]; [n le 2 select I[n] else 6*Self(n-1)+5*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 14 2011
    
  • Mathematica
    Join[{a=0,b=1},Table[c=6*b+5*a;a=b;b=c,{n,100}]] (* Vladimir Joseph Stephan Orlovsky, Jan 16 2011 *)
    CoefficientList[Series[x/(1-6x-5x^2),{x,0,20}],x] (* or *) LinearRecurrence[ {6,5},{0,1},30] (* Harvey P. Dale, Oct 30 2017 *)
  • PARI
    a(n)=([0,1; 5,6]^n*[0;1])[1,1] \\ Charles R Greathouse IV, Oct 03 2016
  • Sage
    [lucas_number1(n,6,-5) for n in range(0, 21)] # Zerinvary Lajos, Apr 24 2009
    

Formula

a(n) = 6*a(n-1) + 5*a(n-2).
a(n) = sqrt(14)*(3+sqrt(14))^n/28 - sqrt(14)*(3-sqrt(14))^n/28. - Paul Barry, Feb 13 2004

A193376 T(n,k) = number of ways to place any number of 2 X 1 tiles of k distinguishable colors into an n X 1 grid; array read by descending antidiagonals, with n, k >= 1.

Original entry on oeis.org

1, 1, 2, 1, 3, 3, 1, 4, 5, 5, 1, 5, 7, 11, 8, 1, 6, 9, 19, 21, 13, 1, 7, 11, 29, 40, 43, 21, 1, 8, 13, 41, 65, 97, 85, 34, 1, 9, 15, 55, 96, 181, 217, 171, 55, 1, 10, 17, 71, 133, 301, 441, 508, 341, 89, 1, 11, 19, 89, 176, 463, 781, 1165, 1159, 683, 144, 1, 12, 21, 109, 225, 673
Offset: 1

Views

Author

R. H. Hardin, Jul 24 2011

Keywords

Comments

Transposed variant of A083856. - R. J. Mathar, Aug 23 2011
As to the sequences by columns beginning (1, N, ...), let m = (N-1). The g.f. for the sequence (1, N, ...) is 1/(1 - x - m*x^2). Alternatively, the corresponding matrix generator is [[1,1], [m,0]]. Another equivalency is simply: The sequence beginning (1, N, ...) is the INVERT transform of (1, m, 0, 0, 0, ...). Convergents to the sequences a(n)/a(n-1) are (1 + sqrt(4*m+1))/2. - Gary W. Adamson, Feb 25 2014

Examples

			Array T(n,k) (with rows n >= 1 and column k >= 1) begins as follows:
  ..1...1....1....1.....1.....1.....1......1......1......1......1......1...
  ..2...3....4....5.....6.....7.....8......9.....10.....11.....12.....13...
  ..3...5....7....9....11....13....15.....17.....19.....21.....23.....25...
  ..5..11...19...29....41....55....71.....89....109....131....155....181...
  ..8..21...40...65....96...133...176....225....280....341....408....481...
  .13..43...97..181...301...463...673....937...1261...1651...2113...2653...
  .21..85..217..441...781..1261..1905...2737...3781...5061...6601...8425...
  .34.171..508.1165..2286..4039..6616..10233..15130..21571..29844..40261...
  .55.341.1159.2929..6191.11605.19951..32129..49159..72181.102455.141361...
  .89.683.2683.7589.17621.35839.66263.113993.185329.287891.430739.624493...
  ...
Some solutions for n = 5 and k = 3 with colors = 1, 2, 3 and empty = 0:
..0....2....3....2....0....1....0....0....2....0....0....2....3....0....0....0
..0....2....3....2....2....1....2....3....2....1....0....2....3....1....1....1
..1....0....0....0....2....0....2....3....2....1....0....1....0....1....1....1
..1....2....2....0....3....2....2....3....2....0....3....1....3....3....2....1
..0....2....2....0....3....2....2....3....0....0....3....0....3....3....2....1
		

Crossrefs

Column 1 is A000045(n+1), column 2 is A001045(n+1), column 3 is A006130, column 4 is A006131, column 5 is A015440, column 6 is A015441(n+1), column 7 is A015442(n+1), column 8 is A015443, column 9 is A015445, column 10 is A015446, column 11 is A015447, and column 12 is A053404,
Row 2 is A000027(n+1), row 3 is A004273(n+1), row 4 is A028387, row 5 is A000567(n+1), and row 6 is A106734(n+2).
Diagonal is A171180, superdiagonal 1 is A083859(n+1), and superdiagonal 2 is A083860(n+1).

Programs

  • Maple
    T:= proc(n,k) option remember; `if`(n<0, 0,
          `if`(n<2 or k=0, 1, k*T(n-2, k) +T(n-1, k)))
        end;
    seq(seq(T(n, d+1-n), n=1..d), d=1..12); # Alois P. Heinz, Jul 29 2011
  • Mathematica
    T[n_, k_] := T[n, k] = If[n < 0, 0, If[n < 2 || k == 0, 1, k*T[n-2, k]+T[n-1, k]]]; Table[Table[T[n, d+1-n], {n, 1, d}], {d, 1, 12}] // Flatten (* Jean-François Alcover, Mar 04 2014, after Alois P. Heinz *)

Formula

With z X 1 tiles of k colors on an n X 1 grid (with n >= z), either there is a tile (of any of the k colors) on the first spot, followed by any configuration on the remaining (n-z) X 1 grid, or the first spot is vacant, followed by any configuration on the remaining (n-1) X 1. Thus, T(n,k) = T(n-1,k) + k*T(n-z,k), with T(n,k) = 1 for n = 0, 1, ..., z-1.
The solution is T(n,k) = Sum_r r^(-n-1)/(1 + z*k*r^(z-1)), where the sum is over the roots r of the polynomial k*x^z + x - 1.
For z = 2, T(n,k) = ((2*k / (sqrt(1 + 4*k) - 1))^(n+1) - (-2*k/(sqrt(1 + 4*k) + 1))^(n+1)) / sqrt(1 + 4*k).
T(n,k) = Sum_{s=0..[n/2]} binomial(n-s,s) * k^s.
For z X 1 tiles, T(n,k,z) = Sum_{s = 0..[n/z]} binomial(n-(z-1)*s, s) * k^s. - R. H. Hardin, Jul 31 2011

Extensions

Formula and proof from Robert Israel in the Sequence Fans mailing list.

A100303 Expansion of (1 - x - 4*x^2)/(1 - x - 8*x^2).

Original entry on oeis.org

1, 0, 4, 4, 36, 68, 356, 900, 3748, 10948, 40932, 128516, 455972, 1484100, 5131876, 17004676, 58059684, 194097092, 658574564, 2211351300, 7479947812, 25170758212, 85010340708, 286376406404, 966459132068, 3257470383300, 10989143439844, 37048906506244, 124962054024996
Offset: 0

Views

Author

Paul Barry, Nov 12 2004

Keywords

Comments

Construct a graph as follows: form the graph whose adjacency matrix is the tensor product of that of P_3 and [1,1;1,1], then add a loop at each of the extremity nodes. (Spectrum : [0; 1^3; (1-sqrt(33))/2;(1+sqrt(33))/2]). a(n) counts closed walks of length n at each of the internal nodes.

Crossrefs

Programs

  • Magma
    [1] cat [n le 2 select 4*(n-1) else Self(n-1) +8*Self(n-2): n in [1..30]]; // G. C. Greubel, Feb 04 2023
    
  • Mathematica
    CoefficientList[Series[(1-x-4x^2)/(1-x-8x^2),{x,0,30}],x] (* Harvey P. Dale, Dec 01 2013 *)
    LinearRecurrence[{1,8}, {1,0,4}, 31] (* G. C. Greubel, Feb 04 2023 *)
  • SageMath
    def A100303(n): return (1/2)*int(n==0) + 4*lucas_number1(n-1, 1, -8)
    [A100303(n) for n in range(31)] # G. C. Greubel, Feb 04 2023

Formula

a(n) = 0^n/2 + (4/sqrt(33))*( ((1 + sqrt(33))/2)^(n-1) - ((1 - sqrt(33))/2)^(n-1) ).
a(n) = (1/2)*(A015443(n) - A015443(n-1)), n > 0. - Ralf Stephan, Jul 21 2013
E.g.f.: (33 + exp(x/2)*(33*cosh(sqrt(33)*x/2) - sqrt(33)*sinh(sqrt(33)*x/2)))/66. - Stefano Spezia, Sep 08 2022
a(n) = (1/2)*[n=0] + 4*A015443(n). - G. C. Greubel, Feb 04 2023

A350470 Array read by ascending antidiagonals. T(n, k) = J(k, n) where J are the Jacobsthal polynomials.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 5, 5, 1, 1, 1, 7, 9, 11, 1, 1, 1, 9, 13, 29, 21, 1, 1, 1, 11, 17, 55, 65, 43, 1, 1, 1, 13, 21, 89, 133, 181, 85, 1, 1, 1, 15, 25, 131, 225, 463, 441, 171, 1, 1, 1, 17, 29, 181, 341, 937, 1261, 1165, 341, 1
Offset: 0

Views

Author

Peter Luschny, Mar 19 2022

Keywords

Examples

			Array starts:
n\k 0, 1,  2,  3,   4,    5,    6,     7,      8,      9, ...
---------------------------------------------------------------------
[0] 1, 1,  1,  1,   1,    1,    1,     1,      1,      1, ... A000012
[1] 1, 1,  3,  5,  11,   21,   43,    85,    171,    341, ... A001045
[2] 1, 1,  5,  9,  29,   65,  181,   441,   1165,   2929, ... A006131
[3] 1, 1,  7, 13,  55,  133,  463,  1261,   4039,  11605, ... A015441
[4] 1, 1,  9, 17,  89,  225,  937,  2737,  10233,  32129, ... A015443
[5] 1, 1, 11, 21, 131,  341, 1651,  5061,  21571,  72181, ... A015446
[6] 1, 1, 13, 25, 181,  481, 2653,  8425,  40261, 141361, ... A053404
[7] 1, 1, 15, 29, 239,  645, 3991, 13021,  68895, 251189, ... A350468
[8] 1, 1, 17, 33, 305,  833, 5713, 19041, 110449, 415105, ... A168579
[9] 1, 1, 19, 37, 379, 1045, 7867, 26677, 168283, 648469, ... A350469
      A005408 | A082108 |
           A016813   A014641
		

Crossrefs

Cf. A350467 (main diagonal), A352361 (Fibonacci polynomials), A352362 (Lucas polynomials).

Programs

  • Maple
    J := (n, x) -> add(2^k*binomial(n - k, k)*x^k, k = 0..n):
    seq(seq(J(k, n-k), k = 0..n), n = 0..10);
  • Mathematica
    T[n_, k_] := Hypergeometric2F1[(1 - k)/2, -k/2, -k, -8 n];
    Table[T[n, k], {n, 0, 9}, {k, 0, 9}] // TableForm
    (* or *)
    T[n_, k_] := With[{s = Sqrt[8*n+1]}, ((1+s)^(k+1) - (1-s)^(k+1)) / (2^(k+1)*s)];
    Table[Simplify[T[n, k]], {n, 0, 9}, {k, 0, 9}] // TableForm
  • PARI
    T(n, k) = ([1, 2; k, 0]^n)[1, 1] ;
    export(T)
    for(k = 0, 9, print(parvector(10, n, T(n - 1, k))))

Formula

T(n, k) = Sum_{j=0..k} binomial(k - j, j)*(2*n)^j.
T(n, k) = ((1+s)^(k+1) - (1-s)^(k+1)) / (2^(k+1)*s) where s = sqrt(8*n + 1).
T(n, k) = [x^k] (1 / (1 - x - 2*n*x^2)).
T(n, k) = hypergeom([1/2 - k/2, -k/2], [-k], -8*n).

A100302 Expansion of (1 - x - 6*x^2)/((1 - x)*(1 - x - 8*x^2)).

Original entry on oeis.org

1, 1, 3, 5, 23, 57, 235, 685, 2559, 8033, 28499, 92757, 320743, 1062793, 3628731, 12131069, 41160911, 138209457, 467496739, 1573172389, 5313146295, 17898525401, 60403695755, 203591898957, 686821464991, 2315556656641, 7810128376563, 26334581629685, 88815608642183
Offset: 0

Views

Author

Paul Barry, Nov 12 2004

Keywords

Comments

Construct a graph as follows: form the graph whose adjacency matrix is the tensor product of that of P_3 and [1,1;1,1], then add a loop at each of the extremity nodes. a(n) counts closed walks of length n at each of the extremity nodes.

Crossrefs

Programs

  • Magma
    I:=[1,1,3]; [n le 3 select I[n] else 2*Self(n-1) +7*Self(n-2) -8*Self(n-3): n in [1..31]]; // G. C. Greubel, Feb 04 2023
    
  • Mathematica
    LinearRecurrence[{2,7,-8},{1,1,3},29] (* Stefano Spezia, Sep 08 2022 *)
  • SageMath
    def A100302(n): return (3 + lucas_number1(n+1, 1, -8))/4
    [A100302(n) for n in range(31)] # G. C. Greubel, Feb 04 2023

Formula

a(n) = 2*a(n-1) + 7*a(n-2) - 8*a(n-3).
a(n) = 3/4 + (1/(4*sqrt(33)))*(((1 + sqrt(33))/2)^(n+1) - ((1 - sqrt(33))/2)^(n+1)).
E.g.f.: 3*exp(x)/4 + exp(x/2)*(33*cosh(sqrt(33)*x/2) + sqrt(33)*sinh(sqrt(33)*x/2))/132. - Stefano Spezia, Sep 08 2022
a(n) = (1/4)*(3 + A015443(n+1)). - G. C. Greubel, Feb 04 2023

A100304 Expansion of (1 - x - 6*x^2)/(1 - x - 8*x^2).

Original entry on oeis.org

1, 0, 2, 2, 18, 34, 178, 450, 1874, 5474, 20466, 64258, 227986, 742050, 2565938, 8502338, 29029842, 97048546, 329287282, 1105675650, 3739973906, 12585379106, 42505170354, 143188203202, 483229566034, 1628735191650, 5494571719922, 18524453253122, 62481027012498
Offset: 0

Views

Author

Paul Barry, Nov 12 2004

Keywords

Comments

Construct a graph as follows: form the graph whose adjacency matrix is the tensor product of that of P_3 and [1,1;1,1], then add a loop at each of the 'internal' nodes. (Spectrum : [0^3; 1; (1-sqrt(33))/2;(1+sqrt(33))/2]). a(n) counts closed walks of length n at each of the extremity nodes.

Crossrefs

Essentially half A100303.
Cf. A015443, A100302 (partial sums), A100305.

Programs

  • Magma
    [1] cat [n le 2 select 2*(n-1) else Self(n-1) +8*Self(n-2): n in [1..30]]; // G. C. Greubel, Feb 04 2023
    
  • Mathematica
    LinearRecurrence[{1,8},{1,0,2},27] (* Stefano Spezia, Sep 08 2022 *)
  • SageMath
    def A100304(n): return (3/4)*int(n==0) + 2*lucas_number1(n-1, 1, -8)
    [A100304(n) for n in range(31)] # G. C. Greubel, Feb 04 2023

Formula

a(n) = 3*0^n/4 + (2/sqrt(33))*( ((1 + sqrt(33))/2)^(n-1) - ((1 - sqrt(33))/2)^(n-1) ).
E.g.f.: (99 + exp(x/2)*(33*cosh(sqrt(33)*x/2) - sqrt(33)*sinh(sqrt(33)*x/2)))/132. - Stefano Spezia, Sep 08 2022
a(n) = (3/4)*[n=0] + 2*A015443(n-2). - G. C. Greubel, Feb 04 2023
Previous Showing 11-20 of 30 results. Next