login
a(n) = Lucas(5^n).
6

%I #28 Apr 01 2023 18:00:23

%S 11,167761,132878596168524201724674011

%N a(n) = Lucas(5^n).

%C Previous name was: a(n) = round(phi^(5^n)) where phi = 1.6180339887... = (sqrt(5) + 1)/2 = A001622.

%C a(4), a 131-digit number, is too large to show here.

%F a(n) = phi^(5^n) + (1-phi)^(5^n) = phi^(5^n) + (-phi)^(-5^n). - _Artur Jasinski_, Oct 05 2008

%F From _Peter Bala_, Nov 14 2022: (Start)

%F a(n) = A000032(5^n).

%F a(n) = a(n-1)^5 + 5*a(n-1)^3 + 5*a(n-1) with a(1) = 11.

%F a(n) == 1 (mod 5).

%F a(n+1) == a(n) (mod 5^(n+1)) for n >= 1 (a particular case of the Gauss congruences for the Lucas numbers).

%F Conjecture: a(n+1) == a(n) (mod 5^(n+r+1)) for n >= r.

%F The smallest positive residue of a(n) mod(5^n) = A268922(n).

%F In the ring of 5-adic integers the limit_{n -> oo} a(n) exists and is equal to A269591. An example is given below. (End)

%e The base 5 representation of a(3) = 132878596168524201724674011 begins 1 + 2*5 + 0*(5^2) + 2*(5^3) + 3*(5^4) + 0*(5^5) + 4*(5^6) + O(5^7) so A269591 begins [1, 2, 0, 2, 3, 0, 4, ...]. - _Peter Bala_, Nov 14 2022

%p a := proc(n) option remember; if n = 1 then 11 else a(n-1)^5 + 5*a(n-1)^3 + 5*a(n-1) end if; end;

%p seq(a(n), n = 1..5); # _Peter Bala_, Nov 14 2022

%t c = N[GoldenRatio, 1000]; Table[Round[c^(5^n)], {n, 1, 5}]

%t c = (1 + Sqrt[5])/2; Table[Expand[c^(5^n) + (1 - c)^(5^n)], {n, 0, 5}] (* _Artur Jasinski_, Oct 05 2008 *)

%t LucasL[5^Range[5]] (* _Harvey P. Dale_, Apr 01 2023 *)

%Y Cf. A000032, A000351, A001622, A001566, A006267, A144836, A144838, A144839, A268922, A269591.

%K nonn,easy

%O 1,1

%A _Artur Jasinski_, Sep 22 2008

%E New name from _Peter Bala_, Nov 10 2022