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

A090018 a(n) = 6*a(n-1) + 3*a(n-2) for n > 2, a(0)=1, a(1)=6.

Original entry on oeis.org

1, 6, 39, 252, 1629, 10530, 68067, 439992, 2844153, 18384894, 118841823, 768205620, 4965759189, 32099171994, 207492309531, 1341251373168, 8669985167601, 56043665125110, 362271946253463, 2341762672896108, 15137391876137037, 97849639275510546, 632510011281474387
Offset: 0

Views

Author

Paul Barry, Nov 19 2003

Keywords

Comments

From Johannes W. Meijer, Aug 09 2010: (Start)
a(n) represents the number of n-move routes of a fairy chess piece starting in a given corner or side square on a 3 X 3 chessboard. This fairy chess piece behaves like a white queen on the eight side and corner squares but on the central square the queen explodes with fury and turns into a red queen, see A180032. The central square leads to A180028. (End)

Crossrefs

Sequences with g.f. of the form 1/(1 - 6*x - k*x^2): A106392 (k=-10), A027471 (k=-9), A006516 (k=-8), A081179 (k=-7), A030192 (k=-6), A003463 (k=-5), A084326 (k=-4), A138395 (k=-3), A154244 (k=-2), A001109 (k=-1), A000400 (k=0), A005668 (k=1), A135030 (k=2), this sequence (k=3), A135032 (k=4), A015551 (k=5), A057089 (k=6), A015552 (k=7), A189800 (k=8), A189801 (k=9), A190005 (k=10), A015553 (k=11).

Programs

  • Magma
    [n le 2 select 6^(n-1) else 6*Self(n-1)+3*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 15 2011
    
  • Maple
    a:= n-> (<<0|1>, <3|6>>^n. <<1,6>>)[1,1]:
    seq(a(n), n=0..30);  # Alois P. Heinz, Jan 17 2011
  • Mathematica
    Join[{a=1,b=6},Table[c=6*b+3*a;a=b;b=c,{n,100}]] (* Vladimir Joseph Stephan Orlovsky, Jan 16 2011 *)
    LinearRecurrence[{6,3}, {1,6}, 41] (* G. C. Greubel, Oct 10 2022 *)
  • PARI
    my(x='x+O('x^30)); Vec(1/(1-6*x-3*x^2)) \\ G. C. Greubel, Jan 24 2018
  • Sage
    [lucas_number1(n,6,-3) for n in range(1, 31)] # Zerinvary Lajos, Apr 24 2009
    

Formula

a(n) = (3+2*sqrt(3))^n*(sqrt(3)/4+1/2) + (1/2-sqrt(3)/4)*(3-2*sqrt(3))^n.
a(n) = (-i*sqrt(3))^n * ChebyshevU(n, isqrt(3)), i^2=-1.
From Johannes W. Meijer, Aug 09 2010: (Start)
G.f.: 1/(1 - 6*x - 3*x^2).
Limit_{k->oo} a(n+k)/a(k) = A141041(n) + A090018(n-1)*sqrt(12) for n >= 1.
Limit_{n->oo} A141041(n)/A090018(n-1) = sqrt(12). (End)
a(n) = Sum_{k=0..n} A099089(n,k)*3^k. - Philippe Deléham, Nov 21 2011
E.g.f.: exp(3*x)*(2*cosh(2*sqrt(3)*x) + sqrt(3)*sinh(2*sqrt(3)*x))/2. - Stefano Spezia, Apr 23 2025

Extensions

Typo in Mathematica program corrected by Vincenzo Librandi, Nov 15 2011

A154626 a(n) = 2^n*A001519(n).

Original entry on oeis.org

1, 2, 8, 40, 208, 1088, 5696, 29824, 156160, 817664, 4281344, 22417408, 117379072, 614604800, 3218112512, 16850255872, 88229085184, 461973487616, 2418924584960, 12665653559296, 66318223015936, 347246723858432, 1818207451086848, 9520257811087360
Offset: 0

Views

Author

Paul Barry, Jan 13 2009

Keywords

Comments

Hankel transform of 1,1,3,11,45,... (see A026375). Binomial transform of A015448.
From Gary W. Adamson, Jul 22 2016: (Start)
A production matrix for the sequence is M =
1, 1, 0, 0, 0, ...
1, 0, 5, 0, 0, ...
1, 0, 0, 5, 0, ...
1, 0, 0, 0, 5, ...
...
Take powers of M, extracting the upper left terms; getting
the sequence starting (1, 1, 2, 8, 40, 208, ...). (End)
The sequence is N=5 in an infinite set of INVERT transforms of powers of N prefaced with a "1". (1, 2, 8, 40, ...) is the INVERT transform of (1, 1, 5, 25, 125, ...). The first six of such sequences are shown in A006012 (N=3). - Gary W. Adamson, Jul 24 2016
From Gary W. Adamson, Jul 27 2016: (Start)
The sequence is the first in an infinite set in which we perform the operation for matrix M (Cf. Jul 22 2016), but change the left border successively from (1, 1, 1, 1, ...) then to (1, 2, 2, 2, ...), then (1, 3, 3, 3, ...) ...; generally (1, N, N, N, ...). Extracting the upper left terms of each matrix operation, we obtain the infinite set beginning:
N=1 (A154626): 1, 2, 8, 40, 208, 1088, ...
N=2 (A084120): 1, 3, 15, 81, 441, 1403, ...
N=3 (A180034): 1, 4, 22, 124, 700, 3952, ...
N=4 (A001653): 1, 5, 29, 169, 985, 5741, ...
N=5 (A000400): 1, 6, 36, 216, 1296, 7776, ...
N=6 (A015451): 1, 7, 43, 265, 1633, 10063, ...
N=7 (A180029): 1, 8, 50, 316, 1996, 12608, ...
N=8 (A180028): 1, 9, 57, 369, 1285, 15417, ...
N=9 (.......): 1, 10, 64, 424, 2800, 18496, ...
N=10 (A123361): 1, 11, 71, 481, 3241, 21851, ...
N=11 (.......): 1, 12, 78, 540, 3708, 25488, ...
... Each of the sequences begins (1, (N+1), (7*N + 1),
(40*N + (N-1)^2), ... (End)
The set of infinite sequences shown (Cf. comment of Jul 27 2016), can be generated from the matrices P = [(1,N; 1,5]^n, (N=1,2,3,...) by extracting the upper left terms. Example: N=6 sequence (A015451): (1, 7, 43, 265, ...) can be generated from the matrix P = [(1,6); (1,5)]^n. - Gary W. Adamson, Jul 28 2016

Crossrefs

Programs

  • Magma
    [n le 2 select (n) else 6*Self(n-1)-4*Self(n-2): n in [1..25]]; // Vincenzo Librandi, May 15 2015
    
  • Mathematica
    LinearRecurrence[{6, -4}, {1, 2}, 30] (* Vincenzo Librandi, May 15 2015 *)
  • PARI
    Vec((1-4*x) / (1-6*x+4*x^2) + O(x^30)) \\ Colin Barker, Sep 22 2017

Formula

G.f.: (1 - 4*x) / (1 - 6*x + 4*x^2).
a(n) = A084326(n+1) - 4*A084326(n). - R. J. Mathar, Jul 19 2012
From Colin Barker, Sep 22 2017: (Start)
a(n) = (((3-sqrt(5))^n*(1+sqrt(5)) + (-1+sqrt(5))*(3+sqrt(5))^n)) / (2*sqrt(5)).
a(n) = 6*a(n-1) - 4*a(n-2) for n>1. (End)
E.g.f.: exp(3*x)*(5*cosh(sqrt(5)*x) - sqrt(5)*sinh(sqrt(5)*x))/5. - Stefano Spezia, Aug 26 2025

A180032 Eight white queens and one red queen on a 3 X 3 chessboard. G.f.: (1+x)/(1-5*x-7*x^2).

Original entry on oeis.org

1, 6, 37, 227, 1394, 8559, 52553, 322678, 1981261, 12165051, 74694082, 458625767, 2815987409, 17290317414, 106163498933, 651849716563, 4002393075346, 24574913392671, 150891318490777, 926480986202582, 5688644160448349
Offset: 0

Views

Author

Johannes W. Meijer, Aug 09 2010

Keywords

Comments

The a(n) represent the number of n-move routes of a fairy chess piece starting in a given corner or side square (m = 1, 3, 7, 9; 2, 4, 6, 8) on a 3 X 3 chessboard. This fairy chess piece behaves like a white chess queen on the eight side and corner squares but on the central square the queen explodes with fury and turns into a red queen.
On a 3 X 3 chessboard there are 2^9 = 512 ways to explode with fury on the central square (we assume here that a red queen might behave like a white queen). The red queen is represented by the A[5] vector in the fifth row of the adjacency matrix A, see the Maple program. For the corner and side squares the 512 red queens lead to 17 red queen sequences, see the cross-references for the complete set.
The sequence above corresponds to 8 red queen vectors, i.e., A[5] vectors, with decimal values 239, 367, 431, 463, 487, 491, 493 and 494. The central square leads for these vectors to A152240.
This sequence belongs to a family of sequences with g.f. (1+x)/(1 - 5*x - k*x^2). The members of this family that are red queen sequences are A180030 (k=8), A180032 (k=7; this sequence), A000400 (k=6), A180033 (k=5), A126501 (k=4), A180035 (k=3), A180037 (k=2) A015449 (k=1) and A003948 (k=0). Other members of this family are A030221 (k=-1), A109114 (k=-3), A020989 (k=-4), A166060 (k=-6).
Inverse binomial transform of A054413.

Crossrefs

Cf. A180028 (Central square).
Cf. Red queen sequences corner and side squares [decimal value A[5]]: A090018 [511], A135030 [255], A180030 [495], A005668 [127], A180032 [239], A000400 [63], A180033 [47], A001109 [31], A126501 [15], A154244 [23], A180035 [7], A138395 [19], A180037 [3], A084326 [17], A015449 [1], A003463 [16], A003948 [0].

Programs

  • Magma
    I:=[1,6]; [n le 2 select I[n] else 5*Self(n-1)+7*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 15 2011
  • Maple
    with(LinearAlgebra): nmax:=20; m:=1; A[5]:= [1,1,1,1,0,1,1,1,0]: A:=Matrix([[0,1,1,1,1,0,1,0,1], [1,0,1,1,1,1,0,1,0], [1,1,0,0,1,1,1,0,1], [1,1,0,0,1,1,1,1,0], A[5], [0,1,1,1,1,0,0,1,1], [1,0,1,1,1,0,0,1,1], [0,1,0,1,1,1,1,0,1], [1,0,1,0,1,1,1,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);
  • Mathematica
    LinearRecurrence[{5,7},{1,6},40] (* Vincenzo Librandi, Nov 15 2011 *)
    CoefficientList[Series[(1+x)/(1-5x-7x^2),{x,0,30}],x] (* Harvey P. Dale, Apr 04 2024 *)

Formula

G.f.: (1+x)/(1 - 5*x - 7*x^2).
a(n) = 5*a(n-1) + 7*a(n-2) with a(0) = 1 and a(1) = 6.
a(n) = ((7+9*A)*A^(-n-1) + (7+9*B)*B^(-n-1))/53 with A = (-5+sqrt(53))/14 and B = (-5-sqrt(53))/14.

A180036 Eight white queens and one red queen on a 3 X 3 chessboard. G.f.: (1 - 2*x)/(1 - 5*x - 3*x^2).

Original entry on oeis.org

1, 3, 18, 99, 549, 3042, 16857, 93411, 517626, 2868363, 15894693, 88078554, 488076849, 2704619907, 14987330082, 83050510131, 460214540901, 2550224234898, 14131764797193, 78309496690659, 433942777844874
Offset: 0

Views

Author

Johannes W. Meijer, Aug 09 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 white queen on the eight side and corner squares but on the central square the queen explodes with fury and turns into a red queen, see A180028.
The sequence above corresponds to 56 red queen vectors, i.e., A[5] vector, with decimal values varying between 7 and 448. The corner and side squares lead for these vectors to A180035.

Programs

  • Magma
    I:=[1,3]; [n le 2 select I[n] else 5*Self(n-1)+3*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 15 2011
  • Maple
    with(LinearAlgebra): nmax:=21; m:=5; A[5]:= [0,0,0,0,0,0,1,1,1]: A:=Matrix([[0,1,1,1,1,0,1,0,1], [1,0,1,1,1,1,0,1,0], [1,1,0,0,1,1,1,0,1], [1,1,0,0,1,1,1,1,0], A[5], [0,1,1,1,1,0,0,1,1], [1,0,1,1,1,0,0,1,1], [0,1,0,1,1,1,1,0,1], [1,0,1,0,1,1,1,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);
  • Mathematica
    LinearRecurrence[{5,3},{1,3},201] (* Vincenzo Librandi, Nov 15 2011 *)

Formula

G.f.: (1-2*x)/(1 - 5*x - 3*x^2).
a(n) = 5*a(n-1) + 3*a(n-2) with a(0) = 1 and a(1) = 3.
a(n) = ((1+16*A)*A^(-n-1) + (1+16*B)*B^(-n-1))/37 with A = (-5+sqrt(37))/6 and B = (-5-sqrt(37))/6.
a(n) = Sum_{k=0..n} A202395(n,k)*2^k. - Philippe Deléham, Dec 21 2011

Extensions

Second formula corrected by Vincenzo Librandi, Nov 15 2011

A180038 Eight white queens and one red queen on a 3 X 3 chessboard. G.f.: (1 - 3*x)/(1 - 5*x - 2*x^2).

Original entry on oeis.org

1, 2, 12, 64, 344, 1848, 9928, 53336, 286536, 1539352, 8269832, 44427864, 238678984, 1282250648, 6888611208, 37007557336, 198815009096, 1068090160152, 5738080818952, 30826584415064, 165609083713224, 889698587396248
Offset: 0

Views

Author

Johannes W. Meijer, Aug 09 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 white queen on the eight side and corner squares but on the central square the queen explodes with fury and turns into a red queen, see A180028.
The sequence above corresponds to 28 red queen vectors, i.e., A[5] vector, with decimal values varying between 3 and 384. The corner and side squares lead for these vectors to A180037.

Programs

  • Magma
    I:=[1,2]; [n le 2 select I[n] else 5*Self(n-1)+2*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 15 2011
  • Maple
    with(LinearAlgebra): nmax:=21; m:=5; A[5]:= [0,0,0,0,0,0,0,1,1]: A:=Matrix([[0,1,1,1,1,0,1,0,1], [1,0,1,1,1,1,0,1,0], [1,1,0,0,1,1,1,0,1], [1,1,0,0,1,1,1,1,0], A[5], [0,1,1,1,1,0,0,1,1], [1,0,1,1,1,0,0,1,1], [0,1,0,1,1,1,1,0,1], [1,0,1,0,1,1,1,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);
  • Mathematica
    LinearRecurrence[{5,2},{1,2},50] (* Vincenzo Librandi, Nov 15 2011 *)

Formula

G.f.: (1-3*x)/(1 - 5*x - 2*x^2).
a(n) = 5*a(n-1) + 2*a(n-2) with a(0) = 1 and a(1) = 2.
a(n) = ((19*A-1)*A^(-n-1) + (19*B-1)*B^(-n-1))/33 with A = (-5+sqrt(33))/4 and B = (-5-sqrt(33))/4.

A180031 Number of n-move paths on a 3 X 3 chessboard of a queen starting or ending in the central square.

Original entry on oeis.org

1, 8, 48, 304, 1904, 11952, 74992, 470576, 2952816, 18528688, 116265968, 729559344, 4577924464, 28726097072, 180253881072, 1131078181936, 7097421958256, 44535735246768, 279458051899888, 1753576141473584
Offset: 0

Views

Author

Johannes W. Meijer, Aug 09 2010

Keywords

Comments

The a(n) represent the number of n-move paths of a chess queen starting or ending in the central square (m = 5) on a 3 X 3 chessboard. The other squares lead to A180030.
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 queen's graph, or we can sum the components of the row vector A^n[m,k], see the Maple program.
Closely related with this sequence are the red queen sequences, see A180028 and A180032.
This sequence belongs to a family of sequences with g.f. (1+k*x)/(1 - 5*x - (k+5)*x^2). The members of this family that are red queen sequences are A180031 (k=3; this sequence), A152240 (k=2), A000400 (k=1), A057088 (k=0), A122690 (k=-1), A180036 (k=-2), A180038 (k=-3), A015449 (k=-4) and A000007 (k=-5). Other members of this family are A030221 (k= -6), 3*A109114 (k=-8), 4*A020989 (k=-9), 6*A166060 (k=-11).

Programs

  • Magma
    I:=[1,8]; [n le 2 select I[n] else 5*Self(n-1)+8*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 15 2011
  • Maple
    with(LinearAlgebra): nmax:=19; m:=5; A[5]:= [1,1,1,1,0,1,1,1,1]: A:=Matrix([[0,1,1,1,1,0,1,0,1], [1,0,1,1,1,1,0,1,0], [1,1,0,0,1,1,1,0,1], [1,1,0,0,1,1,1,1,0], A[5], [0,1,1,1,1,0,0,1,1], [1,0,1,1,1,0,0,1,1], [0,1,0,1,1,1,1,0,1], [1,0,1,0,1,1,1,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);
  • Mathematica
    LinearRecurrence[{5,8},{1,8},50] (* Vincenzo Librandi, Nov 15 2011 *)

Formula

G.f.: (1+3*x)/(1 - 5*x - 8*x^2).
a(n) = 5*a(n-1) + 8*a(n-2) with a(0) = 1 and a(1) = 8.
a(n) = ((A+11)*A^(-n-1) + (B+11)*B^(-n-1))/57 with A = (-5+sqrt(57))/16 and B = (-5-sqrt(57))/16.

A180034 Eight white queens and one red queen on a 3 X 3 chessboard. G.f.: (1 - 2*x)/(1 - 6*x + 2*x^2).

Original entry on oeis.org

1, 4, 22, 124, 700, 3952, 22312, 125968, 711184, 4015168, 22668640, 127981504, 722551744, 4079347456, 23030981248, 130027192576, 734101192960, 4144552772608, 23399114249728, 132105579953152, 745835251219456
Offset: 0

Views

Author

Johannes W. Meijer, Aug 09 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 white queen on the eight side and corner squares but on the central square the queen explodes with fury and turns into a red queen, see A180028.
The sequence above corresponds to 56 red queen vectors, i.e., A[5] vector, with decimal values varying between 23 and 464. The corner and side squares lead for these vectors to A154244.

Crossrefs

Programs

  • Magma
    I:=[1,4]; [n le 2 select I[n] else 6*Self(n-1)-2*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 15 2011
  • Maple
    with(LinearAlgebra): nmax:=21; m:=5; A[5]:= [0,0,0,0,1,0,1,1,1]: A:=Matrix([[0,1,1,1,1,0,1,0,1], [1,0,1,1,1,1,0,1,0], [1,1,0,0,1,1,1,0,1], [1,1,0,0,1,1,1,1,0], A[5], [0,1,1,1,1,0,0,1,1], [1,0,1,1,1,0,0,1,1], [0,1,0,1,1,1,1,0,1], [1,0,1,0,1,1,1,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);
  • Mathematica
    LinearRecurrence[{6,-2},{1,4},50] (* Vincenzo Librandi, Nov 15 2011 *)

Formula

G.f.: (1-2*x)/(1 - 6*x + 2*x^2).
a(n) = 6*a(n-1) - 2*a(n-2) with a(0) = 1 and a(1) = 4.
a(n) = ((1+4*A)*A^(-n-1) + (1+4*B)*B^(-n-1))/14 with A = (3+sqrt(7))/2 and B = (3-sqrt(7))/2.
a(n) = A154244(n) - 2*A154244(n-1). - R. J. Mathar, Aug 14 2012

A180029 Eight white queens and one red queen on a 3 X 3 chessboard. G.f.: (1 + 2*x)/(1 - 6*x - 2*x^2).

Original entry on oeis.org

1, 8, 50, 316, 1996, 12608, 79640, 503056, 3177616, 20071808, 126786080, 800860096, 5058732736, 31954116608, 201842165120, 1274961223936, 8053451673856, 50870632491008, 321330698293760, 2029725454744576
Offset: 0

Views

Author

Johannes W. Meijer, Aug 09 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 white queen on the eight side and corner squares but on the central square the queen explodes with fury and turns into a red queen, see A180028.
The sequence above corresponds to 8 red queen vectors, i.e., A[5] vector, with decimal values 255, 383, 447, 479, 503, 507, 509 and 510. The other squares lead for these vectors to A135030.

Programs

  • Magma
    I:=[1,8]; [n le 2 select I[n] else 6*Self(n-1)+2*Self(n-2): n in [1..20]]; // Vincenzo Librandi, Nov 15 2011
  • Maple
    with(LinearAlgebra): nmax:=19; m:=5; A[5]:= [0,1,1,1,1,1,1,1,1]: A:=Matrix([[0,1,1,1,1,0,1,0,1], [1,0,1,1,1,1,0,1,0], [1,1,0,0,1,1,1,0,1], [1,1,0,0,1,1,1,1,0], A[5], [0,1,1,1,1,0,0,1,1], [1,0,1,1,1,0,0,1,1], [0,1,0,1,1,1,1,0,1], [1,0,1,0,1,1,1,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);
  • Mathematica
    LinearRecurrence[{6,2},{1,8},50 ] (* Vincenzo Librandi, Nov 15 2011 *)

Formula

G.f.: (1+2*x)/(1 - 6*x - 2*x^2).
a(n) = 6*a(n-1) + 2*a(n-2) with a(0) = 1 and a(1) = 8.
a(n) = ((5-4*A)*A^(-n-1) + (5-4*B)*B^(-n-1))/22 with A = (-3+sqrt(11))/2 and B = (-3-sqrt(11))/2.
Lim_{k->infinity} a(n+k)/a(k) = (-1)^(n-1)*A016116(n+1)/(A041015(n-1)*sqrt(11) - A041014(n-1)) for n >= 1.

A180030 Number of n-move paths on a 3 X 3 chessboard of a queen starting or ending in a corner or side square.

Original entry on oeis.org

1, 6, 38, 238, 1494, 9374, 58822, 369102, 2316086, 14533246, 91194918, 572240558, 3590762134, 22531735134, 141384772742, 887177744782, 5566966905846, 34932256487486, 219197017684198, 1375443140320878, 8630791843077974
Offset: 0

Views

Author

Johannes W. Meijer, Aug 09 2010

Keywords

Comments

The a(n) represent the number of n-move paths of a chess queen starting or ending in a given corner or side square (m = 1, 3, 7, 9; 2, 4, 6, 8) on a 3 X 3 chessboard. The central square leads to A180031.
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 queen's graph, or we can sum the components of the row vector A^n[m,k], see the Maple program.
Closely related with this sequence are the red queen sequences, see A180028 and A180032.
Inverse binomial transform of A015555 (without the leading 0).

Programs

  • Magma
    I:=[1,6]; [n le 2 select I[n] else 5*Self(n-1)+8*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 15 2011
  • Maple
    with(LinearAlgebra): nmax:=20; m:=1; A[5]:= [1,1,1,1,0,1,1,1,1]: A:=Matrix([[0,1,1,1,1,0,1,0,1], [1,0,1,1,1,1,0,1,0], [1,1,0,0,1,1,1,0,1], [1,1,0,0,1,1,1,1,0], A[5], [0,1,1,1,1,0,0,1,1], [1,0,1,1,1,0,0,1,1], [0,1,0,1,1,1,1,0,1], [1,0,1,0,1,1,1,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);
  • Mathematica
    LinearRecurrence[{5,8},{1,6},201] (* Vincenzo Librandi, Nov 15 2011 *)

Formula

G.f.: (1+x)/(1 - 5*x - 8*x^2).
a(n) = 5*a(n-1) + 8*a(n-2) with a(0) = 1 and a(1) = 6.
a(n) = ((7+11*A)*A^(-n-1) + (7+11*B)*B^(-n-1))/57 with A = (-5+sqrt(57))/16 and B = (-5-sqrt(57))/16.
Showing 1-9 of 9 results.