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 21 results. Next

A126473 Number of strings over a 5 symbol alphabet with adjacent symbols differing by three or less.

Original entry on oeis.org

1, 5, 23, 107, 497, 2309, 10727, 49835, 231521, 1075589, 4996919, 23214443, 107848529, 501037445, 2327695367, 10813893803, 50238661313, 233396326661, 1084301290583, 5037394142315, 23402480441009, 108722104190981, 505095858086951, 2346549744920747
Offset: 0

Views

Author

R. H. Hardin, Dec 27 2006

Keywords

Comments

[Empirical] a(base,n) = a(base-1,n) + 7^(n-1) for base >= 3n-2; a(base,n) = a(base-1,n) + 7^(n-1)-2 when base = 3n-3.
From Johannes W. Meijer, Aug 01 2010: (Start)
The a(n) represent the number of n-move routes of a fairy chess piece starting in a given side square (m = 2, 4, 6 or 8) on a 3 X 3 chessboard. This fairy chess piece behaves like a king on the eight side and corner squares but on the central square the king goes crazy and turns into a red king, see A179596.
For the side squares the 512 red kings lead to 47 different red king sequences, see the cross-references for some examples.
The sequence above corresponds to four A[5] vectors with the decimal [binary] values 367 [1,0,1,1,0,1,1,1,1], 463 [1,1,1,0,0,1,1,1,1], 487 [1,1,1,1,0,0,1,1,1] and 493 [1,1,1,1,0,1,1,0,1]. These vectors lead for the corner squares to A179596 and for the central square to A179597.
This sequence belongs to a family of sequences with g.f. (1+x)/(1-4*x-k*x^2). Red king sequences that are members of this family are A003947 (k=0), A015448 (k=1), A123347 (k=2), A126473 (k=3; this sequence) and A086347 (k=4). Other members of this family are A000351 (k=5), A001834 (k=-1), A111567 (k=-2), A048473 (k=-3) and A053220 (k=-4)
Inverse binomial transform of A154244. (End)
Equals the INVERT transform of A055099: (1, 4, 14, 50, 178, ...). - Gary W. Adamson, Aug 14 2010
Number of one-sided n-step walks taking steps from {E, W, N, NE, NW}. - Shanzhen Gao, May 10 2011
For n>=1, a(n) equals the numbers of words of length n-1 on alphabet {0,1,2,3,4} containing no subwords 00 and 11. - Milan Janjic, Jan 31 2015

Crossrefs

Cf. 5 symbol differing by two or less A126392, one or less A057960.
Cf. Red king sequences side squares [numerical value A[5]]: A086347 [495], A179598 [239], A126473 [367], A123347 [335], A179602 [95], A154964 [31], A015448 [327], A152187 [27], A003947 [325], A108981 [11], A007483 [2]. - Johannes W. Meijer, Aug 01 2010
Cf. A055099.

Programs

  • Maple
    with(LinearAlgebra): nmax:=19; m:=2; A[5]:= [1,0,1,1,0,1,1,1,1]: A:=Matrix([[0,1,0,1,1,0,0,0,0],[1,0,1,1,1,1,0,0,0],[0,1,0,0,1,1,0,0,0],[1,1,0,0,1,0,1,1,0],A[5],[0,1,1,0,1,0,0,1,1],[0,0,0,1,1,0,0,1,0],[0,0,0,1,1,1,1,0,1],[0,0,0,0,1,1,0,1,0]]): for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m,k],k=1..9): od: seq(a(n), n=0..nmax); # Johannes W. Meijer, Aug 01 2010
    # second Maple program:
    a:= n-> (M-> M[1,2]+M[2,2])(<<0|1>, <3|4>>^n):
    seq(a(n), n=0..24);  # Alois P. Heinz, Jun 28 2021
  • Mathematica
    LinearRecurrence[{4, 3}, {1, 5}, 24] (* Jean-François Alcover, Dec 10 2024 *)
  • PARI
    a(n)=([0,1; 3,4]^n*[1;5])[1,1] \\ Charles R Greathouse IV, May 10 2016

Formula

From Johannes W. Meijer, Aug 01 2010: (Start)
G.f.: (1+x)/(1-4*x-3*x^2).
a(n) = 4*a(n-1) + 3*a(n-2) with a(0) = 1 and a(1) = 5.
a(n) = ((1+3/sqrt(7))/2)*(A)^(-n) + ((1-3/sqrt(7))/2)*(B)^(-n) with A = (-2 + sqrt(7))/3 and B = (-2-sqrt(7))/3.
Lim_{k->oo} a(n+k)/a(k) = (-1)^(n+1)*A000244(n)/(A015530(n)*sqrt(7)-A108851(n))
(End)
a(n) = A015330(n)+A015330(n+1). - R. J. Mathar, May 09 2023

Extensions

Edited by Johannes W. Meijer, Aug 10 2010

A015523 a(n) = 3*a(n-1) + 5*a(n-2), with a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 3, 14, 57, 241, 1008, 4229, 17727, 74326, 311613, 1306469, 5477472, 22964761, 96281643, 403668734, 1692414417, 7095586921, 29748832848, 124724433149, 522917463687, 2192374556806, 9191710988853, 38537005750589
Offset: 0

Views

Author

Keywords

Comments

From Johannes W. Meijer, Aug 01 2010: (Start)
a(n) represents the number of n-move routes of a fairy chess piece starting in a given corner square (m = 1, 3, 7 and 9) on a 3 X 3 chessboard. This fairy chess piece behaves like a king on the eight side and corner squares but on the central square the king goes crazy and turns into a red king, see A179596.
For n >= 1, the sequence above corresponds to 24 red king vectors, i.e., A[5] vectors, with decimal values 27, 30, 51, 54, 57, 60, 90, 114, 120, 147, 150, 153, 156, 177, 180, 210, 216, 240, 282, 306, 312, 402, 408 and 432. These vectors lead for the side squares to A152187 and for the central square to A179606.
This sequence belongs to a family of sequences with g.f. 1/(1-3*x-k*x^2). Red king sequences that are members of this family are A007482 (k=2), A015521 (k=4), A015523 (k=5; this sequence), A083858 (k=6), A015524 (k=7) and A015525 (k=8). We observe that there is no red king sequence for k=3. Other members of this family are A049072 (k=-4), A057083 (k=-3), A000225 (k=-2), A001906 (k=-1), A000244 (k=0), A006190 (k=1), A030195 (k=3), A099012 (k=9), A015528 (k=10) and A015529 (k=11).
Inverse binomial transform of A052918 (with extra leading 0).
(End)
First differences in A197189. - Bruno Berselli, Oct 11 2011
Pisano period lengths: 1, 3, 4, 6, 4, 12, 3, 12, 12, 12, 120, 12, 12, 3, 4, 24, 288, 12, 72, 12, ... - R. J. Mathar, Aug 10 2012
This is the Lucas U(P=3, Q=-5) sequence, and hence for n >= 0, a(n+2)/a(n+1) equals the continued fraction 3 + 5/(3 + 5/(3 + 5/(3 + ... + 5/3))) with n 5's. - Greg Dresden, Oct 06 2019

Crossrefs

Programs

  • Magma
    [ n eq 1 select 0 else n eq 2 select 1 else 3*Self(n-1)+5*Self(n-2): n in [1..30] ]; // Vincenzo Librandi, Aug 23 2011
    
  • Mathematica
    Join[{a = 0, b = 1}, Table[c = 3 * b + 5 * a; a = b; b = c, {n, 100}]] (* Vladimir Joseph Stephan Orlovsky, Jan 16 2011 *)
    a[0] := 0; a[1] := 1; a[n_] := a[n] = 3a[n - 1] + 5a[n - 2]; Table[a[n], {n, 0, 49}] (* Alonso del Arte, Jan 16 2011 *)
  • PARI
    x='x+O('x^30); concat([0], Vec(x/(1-3*x-5*x^2))) \\ G. C. Greubel, Jan 01 2018
  • Sage
    [lucas_number1(n,3,-5) for n in range(0, 24)] # Zerinvary Lajos, Apr 22 2009
    

Formula

a(n) = 3*a(n-1) + 5*a(n-2).
From Paul Barry, Jul 20 2004: (Start)
a(n) = ((3/2 + sqrt(29)/2)^n - (3/2 - sqrt(29)/2)^n)/sqrt(29).
a(n) = Sum_{k=0..floor((n-1)/2)} binomial(n-k-1,k)*5^k*3^(n-2*k-1). (End)
G.f.: x/(1 - 3*x - 5*x^2). - R. J. Mathar, Nov 16 2007
From Johannes W. Meijer, Aug 01 2010: (Start)
Limit_{k->oo} a(n+k)/a(k) = (A072263(n) + a(n)*sqrt(29))/2.
Limit_{n->oo} A072263(n)/a(n) = sqrt(29). (End)
G.f.: G(0)*x/(2-3*x), where G(k) = 1 + 1/(1 - x*(29*k-9)/(x*(29*k+20) - 6/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 17 2013
E.g.f.: 2*exp(3*x/2)*sinh(sqrt(29)*x/2)/sqrt(29). - Stefano Spezia, Oct 06 2019

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

Original entry on oeis.org

0, 1, 4, 19, 88, 409, 1900, 8827, 41008, 190513, 885076, 4111843, 19102600, 88745929, 412291516, 1915403851, 8898489952, 41340171361, 192056155300, 892245135283, 4145149007032, 19257331433977, 89464772757004
Offset: 0

Views

Author

Keywords

Comments

Let b(1)=1, b(k) = floor(b(k-1)) + 3/b(k-1); then for n>1, b(n) = a(n)/a(n-1). - Benoit Cloitre, Sep 09 2002
In general, x/(1 - a*x - b*x^2) has a(n) = Sum_{k=0..floor((n-1)/2)} C(n-k-1,k)*b^k*a^(n-2k-1). - Paul Barry, Apr 23 2005
Pisano period lengths: 1, 2, 1, 4, 24, 2, 21, 4, 3, 24, 40, 4, 84, 42, 24, 8, 288, 6, 18, 24, ... . - R. J. Mathar, Aug 10 2012
This is the Lucas sequence U(4,-3). - Bruno Berselli, Jan 09 2013

Crossrefs

Appears in A179596, A126473 and A179597. - Johannes W. Meijer, Aug 01 2010
Cf. A080042: Lucas sequence V(4,-3).

Programs

  • Magma
    I:=[0, 1]; [n le 2 select I[n] else 4*Self(n-1)+3*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Jun 19 2012
    
  • Mathematica
    LinearRecurrence[{4,3},{0,1},30] (* Vincenzo Librandi, Jun 19 2012 *)
  • PARI
    x='x+O('x^30); concat([0], Vec(x/(1-4*x-3*x^2))) \\ G. C. Greubel, Jan 24 2018
  • Sage
    [lucas_number1(n,4,-3) for n in range(0, 23)]# Zerinvary Lajos, Apr 23 2009
    

Formula

a(n) = 4*a(n-1) + 3*a(n-2).
a(n) = (A086901(n+2) - A086901(n+1))/6. - Ralf Stephan, Feb 01 2004
a(n) = Sum_{k=0..floor((n-1)/2)} C(n-k-1, k)*3^k*4^(n-2k-1). - Paul Barry, Apr 23 2005
a(n) = ((2+sqrt(7))^n - (2-sqrt(7))^n)/sqrt(28). Offset 1. a(3)=19. - Al Hakanson (hawkuu(AT)gmail.com), Jan 05 2009
From Johannes W. Meijer, Aug 01 2010: (Start)
Limit(a(n+k)/a(k), k=infinity) = A108851(n)+a(n)*sqrt(7).
Limit(A108851(n)/a(n), n=infinity) = sqrt(7). (End)
G.f.: x*G(0) where G(k)= 1 + (4*x+3*x^2)/(1 - (4*x+3*x^2)/(4*x + 3*x^2 + 1/G(k+1))); (continued fraction, 3rd kind, 3-step). - Sergei N. Gladkovskii, Jul 28 2012
G.f.: G(0)*x/(2-4*x), where G(k)= 1 + 1/(1 - x*(7*k-4)/(x*(7*k+3) - 2/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 16 2013

A179597 Eight white kings and one red king on a 3 X 3 chessboard. G.f.: (1 + 5*x + 2*x^2)/(1 - 2*x - 11*x^2 - 6*x^3).

Original entry on oeis.org

1, 7, 27, 137, 613, 2895, 13355, 62233, 288741, 1342175, 6233899, 28964169, 134554277, 625117807, 2904117675, 13491856889, 62679715045, 291194561919, 1352817130667, 6284852732713, 29197861274277, 135646005392399
Offset: 0

Views

Author

Johannes W. Meijer, Jul 28 2010, Aug 10 2010

Keywords

Comments

The a(n) represent the number of n-move routes of a fairy chess piece starting in the central square (m = 5) on a 3 X 3 chessboard. This fairy chess piece behaves like a king on the eight side and corner squares but on the central square the king goes crazy and turns into a red king, see A179596.
For the central square the 512 red kings lead to 47 different red king sequences, see the cross-references for some examples.
The sequence above corresponds to four A[5] vectors with decimal [binary] values 367 [1,0,1,1,0,1,1,1,1], 463 [1,1,1,0,0,1,1,1,1], 487 [1,1,1,1,0,0,1,1,1] and 493 [1,1,1,1,0,1,1,0,1]. These vectors lead for the corner squares to A179596 and for the side squares to A126473.
This sequence belongs to a family of sequences with g.f. (1 + (k+2)*x + (2*k-4)*x^2)/(1 - 2*x - (k+8)*x^2 - (2*k)*x^3). Red king sequences that are members of this family are A179607 (k=0), A179605 (k=1), A179601 (k=2), A179597 (k=3; this sequence) and A086348 (k=4). Another member of this family is A179609 (k = -4).

Crossrefs

Red king sequences central square [numerical value A[5]]: A086348 [495], A179599 [239], A179597 [367], A179601 [335], A179603 [95], A154964 [31], A179605 [327], A179606 [27], A179611 [15], A179607 [325], A015521 [11], A007483 [2], A000012 [16], A000007 [0].

Programs

  • Maple
    with(LinearAlgebra): nmax:=21; m:=5; A[1]:= [0,1,0,1,1,0,0,0,0]: A[2]:= [1,0,1,1,1,1,0,0,0]: A[3]:= [0,1,0,0,1,1,0,0,0]: A[4]:=[1,1,0,0,1,0,1,1,0]: A[5]:= [1,0,1,1,0,1,1,1,1]: A[6]:= [0,1,1,0,1,0,0,1,1]: A[7]:= [0,0,0,1,1,0,0,1,0]: A[8]:= [0,0,0,1,1,1,1,0,1]: A[9]:= [0,0,0,0,1,1,0,1,0]: A:=Matrix([A[1],A[2],A[3],A[4],A[5], A[6],A[7],A[8],A[9]]): for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m,k],k=1..9): od: seq(a(n), n=0..nmax);
  • Mathematica
    LinearRecurrence[{2,11,6},{1,7,27},30] (* Harvey P. Dale, Mar 01 2015 *)

Formula

G.f.: (1 + 5*x + 2*x^2)/(1 - 2*x - 11*x^2 - 6*x^3).
a(n) = 2*a(n-1) + 11*a(n-2) + 6*a(n-3) with a(0) = 1, a(1) = 7 and a(2) = 27.
a(n) = 8*(-1/2)^(-n+1)/9 + ((7+11*sqrt(7))*A^(-n-1) + (7-11*sqrt(7))*B^(-n-1))/126 with A = (-2+sqrt(7))/3 and B = (-2-sqrt(7))/3.
Lim_{k->infinity} a(n+k)/a(k) = (-1)^(n+1)*(A000244(n)/(A015530(n)*sqrt(7) - A108851(n))).

A179606 Eight white kings and one red king on a 3 X 3 chessboard. G.f.: (1 + x)/(1 - 3*x - 5*x^2).

Original entry on oeis.org

1, 4, 17, 71, 298, 1249, 5237, 21956, 92053, 385939, 1618082, 6783941, 28442233, 119246404, 499950377, 2096083151, 8788001338, 36844419769, 154473265997, 647641896836, 2715292020493, 11384085545659, 47728716739442
Offset: 0

Views

Author

Johannes W. Meijer, Jul 28 2010

Keywords

Comments

a(n) represents the number of n-move routes of a fairy chess piece starting in the central square (m = 5) on a 3 X 3 chessboard. This fairy chess piece behaves like a king on the eight side and corner squares but on the central square the king goes crazy and turns into a red king, see A179596.
The sequence above corresponds to 24 red king vectors, i.e., A[5] vectors, with decimal values 27, 30, 51, 54, 57, 60, 90, 114, 120, 147, 150, 153, 156, 177, 180, 210, 216, 240, 282, 306, 312, 402, 408 and 432. These vectors lead for the corner squares to A015523 and for the side squares to A152187.
This sequence belongs to a family of sequences with g.f. (1 + (k-4)*x)/(1 - 3*x - k*x^2). Red king sequences that are members of this family are A007483 (k= 2), A015521 (k=4), A179606 (k=5; this sequence), A154964 (k=6), A179603 (k=7) and A179599 (k=8). We observe that there is no red king sequence for k=3. Other members of this family are A006190 (k=1), A133494 (k=0) and A168616 (k=-2).
Inverse binomial transform of A052918.
The sequence b(n+1) = 6*a(n), n >= 0 with b(0)=1, is a berserker sequence, see A180147. The b(n) sequence corresponds to 16 A[5] vectors with decimal values between 111 and 492. These vectors lead for the corner squares to sequence c(n+1)=4*A179606(n), n >= 0 with c(0)=1, and for the side squares to A180140. - Johannes W. Meijer, Aug 14 2010
Equals the INVERT transform of A063782: (1, 3, 10, 32, 104, ...). Example: a(3) = 71 = (1, 1, 4, 7) dot (32, 10, 3, 1) = (32 + 10 + 12 + 17). - Gary W. Adamson, Aug 14 2010

Crossrefs

Cf. A179597 (central square).

Programs

  • Maple
    with(LinearAlgebra): nmax:=22; m:=5; A[1]:= [0,1,0,1,1,0,0,0,0]: A[2]:= [1,0,1,1,1,1,0,0,0]: A[3]:= [0,1,0,0,1,1,0,0,0]: A[4]:= [1,1,0,0,1,0,1,1,0]: A[5]:= [0,0,0,1,1,1,0,0,1]: A[6]:= [0,1,1,0,1,0,0,1,1]: A[7]:= [0,0,0,1,1,0,0,1,0]: A[8]:= [0,0,0,1,1,1,1,0,1]: A[9]:= [0,0,0,0,1,1,0,1,0]: A:=Matrix([A[1],A[2],A[3],A[4],A[5],A[6],A[7],A[8],A[9]]): for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m,k],k=1..9): od: seq(a(n), n=0..nmax);
  • Mathematica
    CoefficientList[Series[(1+x)/(1-3*x-5*x^2), {x, 0, 22}],x] (* or *) LinearRecurrence[{3,5,0},{1,4},23] (* Indranil Ghosh, Mar 05 2017 *)
  • PARI
    print(Vec((1 + x)/(1- 3*x - 5*x^2) + O(x^23))); \\ Indranil Ghosh, Mar 05 2017

Formula

G.f.: (1+x)/(1 - 3*x - 5*x^2).
a(n) = A015523(n) + A015523(n+1).
a(n) = 3*a(n-1) + 5*a(n-2) with a(0) = 1 and a(1) = 4.
a(n) = ((29 + 7*sqrt(29))*A^(-n-1) + (29-7*sqrt(29))*B^(-n-1))/290 with A = (-3+sqrt(29))/10 and B = (-3-sqrt(29))/10
Limit_{k->oo} a(n+k)/a(k) = (-1)^(n+1)*A000351(n)*A130196(n)/(A015523(n)*sqrt(29) - A072263(n)) for n >= 1.

A086346 On a 3 X 3 board, the number of n-move paths for a chess king ending in a given corner square.

Original entry on oeis.org

1, 3, 18, 80, 400, 1904, 9248, 44544, 215296, 1039104, 5018112, 24227840, 116985856, 564850688, 2727354368, 13168803840, 63584665600, 307013812224, 1482394042368, 7157631156224, 34560101318656, 166870928850944, 805724122775552, 3890380202311680, 18784417308737536, 90699190027419648
Offset: 0

Views

Author

Zak Seidov, Jul 17 2003

Keywords

Comments

From Johannes W. Meijer, Aug 01 2010: (Start)
The a(n) represent the number of n-move paths of a chess king on a 3 X 3 board that end or start in a given corner square m (m = 1, 3, 7, 9). To determine the a(n) we can either sum the components of the column vector A^n[k,m], with A the adjacency matrix of the king's graph, or we can sum the components of the row vector A^n[m,k], see the Maple program.
Inverse binomial transform of A079291 (without the leading 0).
(End)
From R. J. Mathar, Oct 12 2010: (Start)
The row n=3 of an array counting king walks on an n X n board with k steps, starting from a corner:
1, 3, 9, 27, 81, 243, 729, 2187, 6561, 19683, 59049, ...;
1, 3, 18, 80, 400, 1904, 9248, 44544, 215296, 1039104, 5018112, ...;
1, 3, 18, 105, 615, 3600, 21075, 123375, 722250, 4228125, 24751875, ...;
1, 3, 18, 105, 684, 4359, 28278, 182349, 1179792, 7622667, 49283802, ...;
1, 3, 18, 105, 684, 4550, 30807, 209867, 1434279, 9815190, 67209723, ...;
1, 3, 18, 105, 684, 4550, 31340, 218056, 1533712, 10829360, 76720288, ...;
1, 3, 18, 105, 684, 4550, 31340, 219555, 1559835, 11177190, 80573373, ...;
1, 3, 18, 105, 684, 4550, 31340, 219555, 1564080, 11259785, 81765550, ...;
1, 3, 18, 105, 684, 4550, 31340, 219555, 1564080, 11271876, 82025163, ...;
1, 3, 18, 105, 684, 4550, 31340, 219555, 1564080, 11271876, 82059768, ...;
1, 3, 18, 105, 684, 4550, 31340, 219555, 1564080, 11271876, 82059768, ...;
The partial sums along the rows are documented in A123109 (king walks with between 1 and k steps). (End)

References

  • Gary Chartrand, Introductory Graph Theory, pp. 217-221, 1984. [From Johannes W. Meijer, Aug 01 2010]

Crossrefs

Programs

  • Magma
    [2^(n-3)*(Evaluate(DicksonFirst(n+2,-1), 2) +2*(-1)^n): n in [0..30]]; // G. C. Greubel, Aug 18 2022
    
  • Maple
    with(LinearAlgebra):
    nmax:=19; m:=1;
    A[5]:= [1, 1, 1, 1, 0, 1, 1, 1, 1]:
    A:=Matrix([[0, 1, 0, 1, 1, 0, 0, 0, 0], [1, 0, 1, 1, 1, 1, 0, 0, 0], [0, 1, 0, 0, 1, 1, 0, 0, 0], [1, 1, 0, 0, 1, 0, 1, 1, 0], A[5], [0, 1, 1, 0, 1, 0, 0, 1, 1], [0, 0, 0, 1, 1, 0, 0, 1, 0], [0, 0, 0, 1, 1, 1, 1, 0, 1], [0, 0, 0, 0, 1, 1, 0, 1, 0]]):
    for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m, k], k=1..9): od: seq(a(n), n=0..nmax); # Johannes W. Meijer, Aug 01 2010
  • Mathematica
    Table[(1/32)(2(-2)^(n+2)+(2+Sqrt[8])^(n+2)+(2-Sqrt[8])^(n+2)), {n, 0, 19}] // FullSimplify
    LinearRecurrence[{2,12,8}, {1,3,18}, 31] (* G. C. Greubel, Aug 18 2022 *)
  • PARI
    Vec((1+x)/((1+2*x)*(1-4*x-4*x^2))+O(x^30)) \\ Joerg Arndt, Jan 29 2024
  • SageMath
    [2^(n-3)*(lucas_number2(n+2,2,-1) +2*(-1)^n) for n in (0..30)] # G. C. Greubel, Aug 18 2022
    

Formula

a(n) = (1/32)*(2*(-2)^(n+2) + (2+sqrt(8))^(n+2) + (2-sqrt(8))^(n+2)).
From R. J. Mathar, Jul 22 2010: (Start)
a(n) = 2*a(n-1) + 12*a(n-2) + 8*a(n-3).
G.f.: (1+x) / ( (1+2*x)*(1-4*x-4*x^2) ).
a(n) = (2*A057087(n-1) + 3*A057087(n) + (-2)^n)/4. (End)
Limit_{k->oo} a(n+k)/a(k) = A084128(n) + 2*A057087(n-1)*sqrt(2). - Johannes W. Meijer, Aug 01 2010
a(n) = A110048(n) + A110048(n-1). - R. J. Mathar, Mar 08 2021
a(n) = 2^(n-3)*(A002203(n+2) + 2*(-1)^n). - G. C. Greubel, Aug 18 2022

Extensions

Offset changed and edited by Johannes W. Meijer, Jul 15 2010

A152187 a(n) = 3*a(n-1) + 5*a(n-2), with a(0)=1, a(1)=5.

Original entry on oeis.org

1, 5, 20, 85, 355, 1490, 6245, 26185, 109780, 460265, 1929695, 8090410, 33919705, 142211165, 596232020, 2499751885, 10480415755, 43940006690, 184222098845, 772366329985, 3238209484180, 13576460102465, 56920427728295
Offset: 0

Views

Author

Philippe Deléham, Nov 28 2008

Keywords

Comments

Unsigned version of A152185.
From Johannes W. Meijer, Aug 01 2010: (Start)
The a(n) represent the number of n-move routes of a fairy chess piece starting in a given side square (m = 2, 4, 6 and 8) on a 3 X 3 chessboard. This fairy chess piece behaves like a king on the eight side and corner squares but on the central square the king goes crazy and turns into a red king, see A179596.
The sequence above corresponds to 24 red king vectors, i.e., A[5] vectors, with decimal values 27, 30, 51, 54, 57, 60, 90, 114, 120, 147, 150, 153, 156, 177, 180, 210, 216, 240, 282, 306, 312, 402, 408 and 432. These vectors lead for the corner squares to A015523 and for the central square to A179606.
This sequence belongs to a family of sequences with g.f. (1+2*x)/(1 - 3*x - k*x^2). Red king sequences that are members of this family are A007483 (k=2), A108981 (k=4), A152187 (k=5; this sequence), A154964 (k=6), A179602 (k=7) and A179598 (k=8). We observe that there is no red king sequence for k=3. Other members of this family are A036563 (k=-2), A054486 (k=-1), A084244 (k=0), A108300 (k=1) and A000351 (k=10).
Inverse binomial transform of A015449 (without the first leading 1).
(End)

Crossrefs

Programs

  • Mathematica
    LinearRecurrence[{3,5},{1,5},40] (* Harvey P. Dale, May 03 2013 *)

Formula

G.f.: (1+2*x)/(1 - 3*x - 5*x^2).
Lim_{k->infinity} a(n+k)/a(k) = (A072263(n) + A015523(n)*sqrt(29))/2. - Johannes W. Meijer, Aug 01 2010
G.f.: G(0)*(1+2*x)/(2-3*x), where G(k) = 1 + 1/(1 - x*(29*k-9)/(x*(29*k+20) - 6/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 17 2013

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

Original entry on oeis.org

1, 2, 11, 50, 233, 1082, 5027, 23354, 108497, 504050, 2341691, 10878914, 50540729, 234799658, 1090820819, 5067682250, 23543191457, 109375812578, 508132824683, 2360658736466, 10967033419913, 50950109889050
Offset: 0

Views

Author

Philippe Deléham, Jul 11 2005

Keywords

Comments

Binomial transform of A083098, second binomial transform of (1, 0, 7, 0, 49, 0, 243, 0, ...).

Crossrefs

Cf. A080042. - Zerinvary Lajos, May 14 2009
Appears in A179596, A179597 and A126473. - Johannes W. Meijer, Aug 01 2010

Programs

  • Magma
    [Floor(((2 + Sqrt(7))^n + (2 - Sqrt(7))^n) / 2): n in [0..30]]; // Vincenzo Librandi, Jul 18 2011
    
  • Mathematica
    LinearRecurrence[{4,3},{1,2},30] (* Harvey P. Dale, Jan 02 2022 *)
  • PARI
    a(n)=round(((2+sqrt(7))^n+(2-sqrt(7))^n)/2) \\ Charles R Greathouse IV, Dec 06 2011
  • Sage
    [lucas_number2(n,4,-3)/2 for n in range(0, 22)] # Zerinvary Lajos, May 14 2009
    

Formula

a(n) = ((2 + sqrt(7))^n + (2 - sqrt(7))^n) / 2.
G.f.: (1 - 2*x) / (1 - 4*x - 3*x^2).
E.g.f.: exp(2*x)*cosh(sqrt(7)*x).
a(n+1)/a(n) converges to 2 + sqrt(7) = 4.645751311064...
Limit_{k->oo} a(n+k)/a(k) = A108851(n) + A015530(n)*sqrt(7); also lim_{n->oo} A108851(n)/A015530(n) = sqrt(7). - Johannes W. Meijer, Aug 01 2010
a(n) = Sum_{k=0..n} A201730(n,k)*6^k. - Philippe Deléham, Dec 06 2011
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x*(7*k-4)/(x*(7*k+3) - 2/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 27 2013
a(n) = (2 + sqrt(7))^n - A015530(n)*sqrt(7). - Robert FERREOL, Aug 04 2025

A173300 a(n) is the denominator of the fraction f = x^n + y^n given that x + y = 1 and x^2 + y^2 = 2.

Original entry on oeis.org

1, 1, 2, 2, 4, 2, 8, 8, 16, 8, 32, 32, 64, 32, 128, 128, 256, 128, 512, 512, 1024, 512, 2048, 2048, 4096, 2048, 8192, 8192, 16384, 8192, 32768, 32768, 65536, 32768, 131072, 131072, 262144, 131072, 524288, 524288, 1048576, 524288, 2097152, 2097152, 4194304, 2097152
Offset: 1

Views

Author

J. Lowell, Feb 15 2010

Keywords

Comments

The denominators of the coefficients of the Taylor series representation of (1+x)/(1-2*x-11*x^2-6*x^3) around x=-1 lead to this sequence, see the Maple program. - Johannes W. Meijer, Aug 16 2010

Examples

			a(3) = 2 because x^3 + y^3 = 5/2.
		

Crossrefs

Cf. A173989 (2-adic valuations).

Programs

  • Maple
    nmax:=45: f:=n-> coeftayl((1+x)/(1-2*x-11*x^2-6*x^3), x=-1, n): a(1):=1: for n from 0 to nmax do a(n+2):= denom(f(n)) od: seq(a(n),n=1..nmax); # Johannes W. Meijer, Aug 16 2010
  • Mathematica
    Denominator[Map[First, NestList[{Last[#], Last[#] + First[#]/2} &, {1, 2}, 50]]] (* Paolo Xausa, Feb 01 2024, after Nick Hobson *)
  • PARI
    a(n) = denominator(2*polcoeff( lift( Mod((1+x)/2,x^2-3)^n ), 0)) \\ Max Alekseyev, Feb 23 2010
    
  • Python
    from fractions import Fraction
    def a173300_gen(a, b):
        while True:
            yield a.denominator
            b, a = b + Fraction(a, 2), b
    for n, a_n in zip(range(1, 47), a173300_gen(1, 2)):
        print(n, a_n)  # Nick Hobson, Jan 30 2024

Formula

a(n) = denominator of ((1+sqrt(3))/2)^n + ((1-sqrt(3))/2)^n. - Max Alekseyev, Feb 23 2010
Conjecture: a(n) = 4*a(n-4), for n >= 7. - Paolo Xausa, Feb 02 2024

Extensions

More terms from Max Alekseyev, Feb 23 2010

A179598 Eight white kings and one red king on a 3 X 3 chessboard. G.f.: (1 + 2*x)/(1 - 3*x - 8*x^2).

Original entry on oeis.org

1, 5, 23, 109, 511, 2405, 11303, 53149, 249871, 1174805, 5523383, 25968589, 122092831, 574027205, 2698824263, 12688690429, 59656665391, 280479519605, 1318691881943, 6199911802669, 29149270463551, 137047105812005
Offset: 0

Views

Author

Johannes W. Meijer, Jul 28 2010

Keywords

Comments

The a(n) represent the number of n-move routes of a fairy chess piece starting in a given side square (m = 2, 4, 6 or 8) on a 3 X 3 chessboard. This fairy chess piece behaves like a king on the eight side and corner squares but on the central square the king goes crazy and turns into a red king, see A179596.
The sequence above corresponds to 10 red king vectors, i.e., A[5] vectors, with decimal values 239, 351, 375, 381, 431, 471, 477, 491, 494, and 501. These vectors lead for the corner squares to A015525 and for the central square to A179599.
Inverse binomial transform of A126501.

Crossrefs

Cf. A126473 (side squares).

Programs

  • Maple
    with(LinearAlgebra): nmax:=21; m:=2; A[1]:= [0,1,0,1,1,0,0,0,0]: A[2]:= [1,0,1,1,1,1,0,0,0]: A[3]:= [0,1,0,0,1,1,0,0,0]: A[4]:=[1,1,0,0,1,0,1,1,0]: A[5]:= [1,0,1,1,1,1,1,0,1]: A[6]:= [0,1,1,0,1,0,0,1,1]: A[7]:= [0,0,0,1,1,0,0,1,0]: A[8]:= [0,0,0,1,1,1,1,0,1]: A[9]:= [0,0,0,0,1,1,0,1,0]: A:=Matrix([A[1],A[2],A[3],A[4],A[5],A[6],A[7],A[8],A[9]]): for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m,k],k=1..9): od: seq(a(n), n=0..nmax);

Formula

G.f.: (1+2*x)/(1 - 3*x - 8*x^2).
a(n) = 3*a(n-1) + 8*a(n-2) with a(0) = 1 and a(1) = 5.
a(n) = ((41+5*sqrt(41))*A^(-n-1) + (41-5*sqrt(41))*B^(-n-1))/328 with A = (-3+sqrt(41))/16 and B = (-3-sqrt(41))/16.
Showing 1-10 of 21 results. Next