A106627 Product L(n)*L_4(n), where L(n) are Lucas numbers and L_4(n) are Lucas 4-step numbers.
8, 1, 9, 28, 105, 286, 918, 2871, 8977, 27892, 87084, 271635, 847182, 2641991, 8240325, 25700488, 80156033, 249994997, 779700654, 2431777739, 7584375260
Offset: 0
Examples
a(0) = 8 because L(0) * L_4(0) = 2 * 4. a(1) = 1 because L(1) * L_4(1) = 1 * 1. a(2) = 9 because L(2) * L_4(2) = 3 * 3. a(3) = 28 because L(3) * L_4(3) = 4 * 7. a(4) = 105 because L(4) * L_4(4) = 7 * 15. a(5) = 286 because L(5) * L_4(5) = 11 * 26. a(6) = 918 because L(6) * L_4(6) = 18 * 51.
Links
- G. C. Greubel, Table of n, a(n) for n = 0..1000
- Index entries for linear recurrences with constant coefficients, signature (1,4,5,9,3,-2,1,-1).
Programs
-
GAP
a:=[8,1,9,28,105,286,918,2871];; for n in [9..30] do a[n]:=a[n-1] +4*a[n-2]+5*a[n-3]+9*a[n-4]+3*a[n-5]-2*a[n-6]+a[n-7]-a[n-8]; od; a; # G. C. Greubel, Feb 19 2019
-
Magma
I:=[8,1,9,28,105,286,918,2871]; [n le 8 select I[n] else Self(n-1)+4*Self(n-2)+5*Self(n-3)+9*Self(n-4)+3*Self(n-5)-2*Self(n-6) + Self(n-7)-Self(n-8): n in [1..30]]; // G. C. Greubel, Feb 19 2019
-
Mathematica
LinearRecurrence[{1,4,5,9,3,-2,1,-1}, {8,1,9,28,105,286,918,2871}, 40] (* G. C. Greubel, Feb 19 2019 *) -
PARI
my(x='x+O('x^40)); Vec((8-7*x-24*x^2-25*x^3-36*x^4-9*x^5+4*x^6 -x^7)/(1-x-4*x^2-5*x^3-9*x^4-3*x^5+2*x^6-x^7+x^8)) \\ G. C. Greubel, Feb 19 2019 -
Sage
((8-7*x-24*x^2-25*x^3-36*x^4-9*x^5+4*x^6-x^7)/(1-x-4*x^2-5*x^3 -9*x^4-3*x^5+2*x^6-x^7+x^8)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Feb 19 2019
Formula
a(n) = +a(n-1) +4*a(n-2) +5*a(n-3) +9*a(n-4) +3*a(n-5) -2*a(n-6) +a(n-7) -a(n-8). - R. J. Mathar, Dec 22 2010
G.f.: (8 -7*x -24*x^2 -25*x^3 -36*x^4 -9*x^5 +4*x^6 -x^7) / (1 -x -4*x^2 -5*x^3 -9*x^4 -3*x^5 +2*x^6 -x^7 +x^8). - Colin Barker, Jun 17 2012
Comments