A180033 Eight white queens and one red queen on a 3 X 3 chessboard. G.f.: (1 + x)/(1 - 5*x - 5*x^2).
1, 6, 35, 205, 1200, 7025, 41125, 240750, 1409375, 8250625, 48300000, 282753125, 1655265625, 9690093750, 56726796875, 332084453125, 1944056250000, 11380703515625, 66623798828125, 390022511718750, 2283231552734375
Offset: 0
Links
- Vincenzo Librandi, Table of n, a(n) for n = 0..200
- Jean-Paul Allouche, Jeffrey Shallit, and Manon Stipulanti, Combinatorics on words and generating Dirichlet series of automatic sequences, arXiv:2401.13524 [math.CO], 2025. See p. 14.
- D. Birmajer, J. B. Gil, and M. D. Weiner, On the Enumeration of Restricted Words over a Finite Alphabet, J. Int. Seq. 19 (2016) # 16.1.3, Example 7.
- Milan Janjic, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
- Index entries for linear recurrences with constant coefficients, signature (5,5).
Crossrefs
Programs
-
Magma
I:=[1,6]; [n le 2 select I[n] else 5*Self(n-1)+5*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 15 2011
-
Maple
with(LinearAlgebra): nmax:=20; m:=1; A[5]:= [0,0,0,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,5},{1,6}, 30] (* Vincenzo Librandi, Nov 15 2011 *)
-
PARI
my(x='x+O('x^30)); Vec((1+x)/(1-5*x-5*x^2)) \\ G. C. Greubel, Apr 07 2019
-
Sage
((1+x)/(1-5*x-5*x^2)).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, Apr 07 2019
Formula
G.f.: (1+x)/(1 - 5*x - 5*x^2).
a(n) = 5*a(n-1) + 5*a(n-2) with a(0) = 1 and a(1) = 6.
a(n) = ((7+5*A)*A^(-n-1) + (7+5*B)*B^(-n-1))/45 with A = (-5+3*sqrt(5))/10 and B = (-5-3*sqrt(5))/10.
Limit_{k->oo} a(n+k)/a(k) = 2*5^(n/2)/(L(2*n) - F(2*n)*sqrt(5)) with L(n) = A000032(n) and F(n) = A000045(n).
a(n) = A057088(n+1)/5. a(2*n) = 5^n*F(4*(n+1))/3, a(2*n+1) = 5^n*L(2*(2*n+3))/3. - Ehren Metcalfe, Apr 04 2019
E.g.f.: exp(5*x/2)*(15*cosh(3*sqrt(5)*x/2) + 7*sqrt(5)*sinh(3*sqrt(5)*x/2))/15. - Stefano Spezia, Mar 17 2025
Comments