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-6 of 6 results.

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

Original entry on oeis.org

1, 6, 42, 288, 1980, 13608, 93528, 642816, 4418064, 30365280, 208700064, 1434392064, 9858552768, 67757668992, 465697330560, 3200729997312, 21998563967232, 151195763787264, 1039165966526976, 7142170381885440
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->1^6, 1->(1^6)0, starting from 0. The number of 1's and 0's of this word is 6*a(n-1) and 6*a(n-2), resp.

Crossrefs

Programs

Formula

a(n) = 6*a(n-1) + 6*a(n-2); a(0)=1, a(1)=6.
a(n) = S(n, i*sqrt(6))*(-i*sqrt(6))^n with S(n, x) := U(n, x/2), Chebyshev's polynomials of the 2nd kind, A049310.
G.f.: 1/(1-6*x-6*x^2).
a(n) = Sum_{k=0..n} 5^k*A063967(n,k). - Philippe Deléham, Nov 03 2006

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

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

A189800 a(n) = 6*a(n-1) + 8*a(n-2), with a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 6, 44, 312, 2224, 15840, 112832, 803712, 5724928, 40779264, 290475008, 2069084160, 14738305024, 104982503424, 747801460736, 5326668791808, 37942424436736, 270267896954880, 1925146777223168, 13713023838978048, 97679317251653632, 695780094221746176
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    I:=[0,1]; [n le 2 select I[n] else 6*Self(n-1)+8*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 14 2011
    
  • Mathematica
    LinearRecurrence[{6, 8}, {0, 1}, 50]
    CoefficientList[Series[-(x/(-1+6 x+8 x^2)),{x,0,50}],x] (* Harvey P. Dale, Jul 26 2011 *)
  • PARI
    a(n)=([0,1; 8,6]^n*[0;1])[1,1] \\ Charles R Greathouse IV, Oct 03 2016

Formula

G.f.: x/(1 - 2*x*(3+4*x)). - Harvey P. Dale, Jul 26 2011

A368156 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 1 + 2*x, p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where u = p(2,x), v = 1 + x^2.

Original entry on oeis.org

1, 1, 2, 2, 4, 5, 3, 10, 14, 12, 5, 20, 41, 44, 29, 8, 40, 98, 148, 131, 70, 13, 76, 224, 408, 497, 376, 169, 21, 142, 482, 1044, 1542, 1588, 1052, 408, 34, 260, 1003, 2492, 4351, 5456, 4894, 2888, 985, 55, 470, 2026, 5684, 11359, 16790, 18400, 14672, 7813
Offset: 1

Views

Author

Clark Kimberling, Jan 20 2024

Keywords

Comments

Because (p(n,x)) is a strong divisibility sequence, for each integer k, the sequence (p(n,k)) is a strong divisibility sequence of integers.

Examples

			First eight rows:
   1
   1    2
   2    4    5
   3   10   14    12
   5   20   41    44    29
   8   40   98   148   131    70
  13   76  224   408   497   376   169
  21  142  482  1044  1542  1588  1052  408
Row 4 represents the polynomial p(4,x) = 3 + 10*x + 14*x^2 + 12*x^3, so (T(4,k)) = (3,10,14,12), k=0..3.
		

Crossrefs

Cf. A000045 (column 1); A000129, (p(n,n-1)); A007482 (row sums), (p(n,1)); A077925 (alternating row sums), (p(n,-1)); A057088, (p(n,2)); A015523, (p(n,-2)); A015568, (p(n,3)); A180250, (p(n,-3)); A094440, A367208, A367209, A367210, A367211, A367297, A367298, A367299, A367300, A367301, A368150, A368151, A368152, A368153, A368154, A368155.

Programs

  • Mathematica
    p[1, x_] := 1; p[2, x_] := 1 + 2 x; u[x_] := p[2, x]; v[x_] := 1 + 2x^2;
    p[n_, x_] := Expand[u[x]*p[n - 1, x] + v[x]*p[n - 2, x]]
    Grid[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]
    Flatten[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]

Formula

p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where p(1,x) = 1, p(2,x) = 1 + 2*x, u = p(2,x), and v = 1 + x^2.
p(n,x) = k*(b^n - c^n), where k = -1/sqrt(5 + 4*x + 8*x^2), b = (1/2)*(2*x + 1 - 1/k), c = (1/2)*(2*x + 1 + 1/k).

A120717 Expansion of x*(67 + 134*x - 287*x^2 - 378*x^3)/((1+2*x)*(1-3*x)*(1 - 5*x - 10*x^2)).

Original entry on oeis.org

0, 67, 536, 3666, 24834, 163870, 1077594, 7054814, 46137578, 301492462, 1969619930, 12865344702, 84029934282, 548824052494, 3584482608186, 23410842173150, 152899603572266, 998608137861166, 6522041823044762
Offset: 0

Views

Author

Roger L. Bagula, Aug 12 2006

Keywords

Comments

Old name was: Sequence produced by Markov chain based on body-centered pentagonal prism 11 X 11 bond graph.
Characteristic polynomial = -480*x^4 - 1040*x^5 - 752*x^6 - 120*x^7 + 90*x^8 + 35*x^9 - x^11.
The molecule with this structure is known as Ferrocene.

Crossrefs

Cf. A180250.

Programs

  • Magma
    R:=PowerSeriesRing(Integers(), 40); [0] cat Coefficients(R!( x*(67+134*x-287*x^2-378*x^3)/((1+2*x)*(1-3*x)*(1-5*x-10*x^2)) )); // G. C. Greubel, Jul 21 2023
    
  • Mathematica
    M0 = {{0,1,1,1,1,1,0,0,0,0,1}, {1,0,1,1,1,0,1,0,0,0,1}, {1,1,0,1,1,0, 0,1,0,0,1}, {1,1,1,0,1,0,0,0,1,0,1}, {1,1,1,1,0,0,0,0,0,1,1}, {1,0, 0,0,0,0,1,1,1,1,1}, {0,1,0,0,0,1,0,1,1,1,1}, {0,0,1,0,0,1,1,0,1,1, 1}, {0,0,0,1,0,1,1,1,0,1,1}, {0,0,0,0,1,1,1,1,1,0,1}, {1,1,1,1,1,1, 1,1,1,1,0}}; v[1]= Table[Fibonacci[n], {n,0,10}]; v[n_]:= v[n] = M0.v[n-1]; Table[v[n][[1]], {n,50}]
    LinearRecurrence[{6,11,-40,-60},{0,67,536,3666,24834},20] (* Harvey P. Dale, May 25 2023 *)
  • SageMath
    A180250= BinaryRecurrenceSequence(5,10,0,1)
    def A120717(n): return (1/10)*(-63*int(n==0) + 49*(-2)^n - 74*3^n + 22*(4*A180250(n+1) + 25*A180250(n)))
    [A120717(n) for n in range(41)] # G. C. Greubel, Jul 21 2023

Formula

G.f.: x*(67 + 134*x - 287*x^2 - 378*x^3)/((1+2*x)*(1-3*x)*(1 - 5*x - 10*x^2)). - Colin Barker, Nov 29 2012
a(n) = (1/10)*(-63*[n=0] + 49*(-2)^n - 74*3^n + 22*(4*A180250(n) + 25*A180250(n-1))). - G. C. Greubel, Jul 21 2023

Extensions

Edited by N. J. A. Sloane, Jul 13 2007
New name (using g.f. by Colin Barker) by Joerg Arndt, Jan 07 2013
Showing 1-6 of 6 results.