%I #15 Oct 15 2024 10:31:24
%S 1,24,-448,-59904,-2158592,-7766016,3080978432,160312590336,
%T 2765438844928,-123759079981056,-10365137990975488,
%U -299512095597133824,2207640196898357248,585186082406535266304,24556707640476321185792,242424234892406990831616
%N 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.
%C 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.
%D Marc Renault, "The Period, Rank, and Order of the (a,b)-Fibonacci Sequence mod m", Math. Mag. 86 (2013) 372 - 380.
%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (48,-1600).
%F 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.
%F a(n) = 48 a(n-1) - 1600 a(n-2), where a(0) = 1, a(1) = 24.
%F From _Stefano Spezia_, Oct 03 2024: (Start)
%F G.f.: (1 - 24*x)/(1 - 48*x + 1600*x^2).
%F E.g.f.: exp(24*x)*cos(32*x). (End)
%t (*Program 1*)
%t A[a_, b_, c_] := ArcCos[(b^2 + c^2 - a^2)/(2 b c)];
%t Table[TrigExpand[(20)^n Cos[n A[4, 5, 3]]], {n, 0, 30}]
%t (*Program 2*)
%t LinearRecurrence[{48, -1600}, {1, 24}, 30]
%Y Cf. A374880, A139030, A376285, A376455.
%K sign,easy,changed
%O 0,2
%A _Clark Kimberling_, Oct 02 2024