OFFSET
0,2
REFERENCES
Hugh C. Williams, Édouard Lucas and Primality Testing, John Wiley and Sons, 1998, p. 75.
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (8,-8,1).
FORMULA
a(n) = 8*a(n-1) - 8*a(n-2) + a(n-3).
From R. J. Mathar, Sep 03 2010: (Start)
G.f.: 15*x/((1-x)*(1-7*x+x^2)). (End)
From G. C. Greubel, Jun 12 2025: (Start)
a(n) = 15*(-1)^n*A156088(n).
E.g.f.: exp(7*x/2)*( 3*cosh(3*sqrt(5)*x/2) - sqrt(5)*sinh(3*sqrt(5)*x/2) ) - 3*exp(x). (End)
Sum_{n>=1} 1/a(n) = (3 - sqrt(5))/10 = A187799 / 100. - Amiram Eldar, Dec 25 2025
MAPLE
luc := proc(n) option remember: if n=0 then RETURN(2) fi: if n=1 then RETURN(1) fi: luc(n-1)+luc(n-2): end: for n from 0 to 40 do printf(`%d, `, luc(4*n+2)-3) od: # James Sellers, Mar 05 2003
MATHEMATICA
LinearRecurrence[{8, -8, 1}, {0, 15, 120}, 20] (* Jean-François Alcover, Nov 29 2023 *)
PROG
(Magma)
A081079:= func< n | Lucas(4*n+2) -3 >; // G. C. Greubel, Jun 12 2025
(SageMath)
def A081079(n): return lucas_number2(4*n+2, 1, -1) -3 # G. C. Greubel, Jun 12 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
R. K. Guy, Mar 04 2003
STATUS
approved
