%I #11 Oct 15 2024 10:31:05
%S 1,16,112,-2816,-134912,-3190784,-48140288,-264175616,10802495488,
%T 451350102016,10122205069312,143370521411584,538974657445888,
%U -40101019526365184,-1498822487822041088,-31921911799759241216,-421972182463479283712,-734345118927640592384
%N 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.
%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+3) 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 (32,-400).
%F 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.
%F a(n) = 32 a(n-1) - 400 a(n-2), where a(0) = 1, a(1) = 16.
%F From _Stefano Spezia_, Oct 03 2024: (Start)
%F G.f.: (1 - 16*x)/(1 - 32*x + 400*x^2).
%F E.g.f.: exp(16*x)*cos(12*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[3, 4, 5]]], {n, 0, 30}]
%t (*Program 2*)
%t LinearRecurrence[{32, -400}, {1, 16}, 30]
%Y Cf. A066771, A374880, A139030, A376455.
%K sign,easy,changed
%O 0,2
%A _Clark Kimberling_, Oct 03 2024