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+3) 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 (32,-400).
FORMULA
a(n) = (20^n * cos(nA), where A is the angle opposite side BC in a triangle ABC having sidelengths |BC|=3, |CA|=4, |AB|=5; ABC is the smallest integer-sided right triangle.
a(n) = 32 a(n-1) - 400 a(n-2), where a(0) = 1, a(1) = 16.
From Stefano Spezia, Oct 03 2024: (Start)
G.f.: (1 - 16*x)/(1 - 32*x + 400*x^2).
E.g.f.: exp(16*x)*cos(12*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[3, 4, 5]]], {n, 0, 30}]
(*Program 2*)
LinearRecurrence[{32, -400}, {1, 16}, 30]
CROSSREFS
KEYWORD
sign,easy
AUTHOR
Clark Kimberling, Oct 03 2024
STATUS
approved