OFFSET
0,2
COMMENTS
If a prime p divides a term, then the indices n such that p divides a(n) comprise an arithmetic sequence; e.g., 7 divides a(4n+2) for n >= 0; 17 divides a(8n+4) for n>= 0. See the Renault paper in References.
REFERENCES
Marc Renault, "The Period, Rank, and Order of the (a,b)-Fibonacci Sequence mod m", Math. Mag. 86 (2013) 372 - 380.
LINKS
Index entries for linear recurrences with constant coefficients, signature (48,-1600).
FORMULA
a(n) = (40)^n * cos(nB), where B is the angle opposite side CA in a triangle ABC having sidelengths |BC|=3, |CA|=4, |AB|=5; ABC is the smallest integer-sided right triangle.
a(n) = 48 a(n-1) - 1600 a(n-2), where a(0) = 1, a(1) = 24.
From Stefano Spezia, Oct 03 2024: (Start)
G.f.: (1 - 24*x)/(1 - 48*x + 1600*x^2).
E.g.f.: exp(24*x)*cos(32*x). (End)
MATHEMATICA
(*Program 1*)
A[a_, b_, c_] := ArcCos[(b^2 + c^2 - a^2)/(2 b c)];
Table[TrigExpand[(20)^n Cos[n A[4, 5, 3]]], {n, 0, 30}]
(*Program 2*)
LinearRecurrence[{48, -1600}, {1, 24}, 30]
CROSSREFS
KEYWORD
sign,easy
AUTHOR
Clark Kimberling, Oct 02 2024
STATUS
approved