login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A106627 Product L(n)*L_4(n), where L(n) are Lucas numbers and L_4(n) are Lucas 4-step numbers. 1
8, 1, 9, 28, 105, 286, 918, 2871, 8977, 27892, 87084, 271635, 847182, 2641991, 8240325, 25700488, 80156033, 249994997, 779700654, 2431777739, 7584375260 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,1
COMMENTS
a(n) is semiprime iff n is an element of A001606 (an index of a prime Lucas number) and an element of A104577 (an index of a prime Lucas 4-step number). The only known such are n = 2, 8, 16, 19, 71, (through 145858).
LINKS
FORMULA
a(n) = A000032(n) * A073817(n).
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
EXAMPLE
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.
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 *)
PROG
(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
(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
(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
(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
CROSSREFS
Sequence in context: A298519 A238271 A217430 * A302152 A160925 A345059
KEYWORD
easy,nonn
AUTHOR
Jonathan Vos Post, May 11 2005
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 07:33 EDT 2024. Contains 371235 sequences. (Running on oeis4.)