OFFSET
1,1
COMMENTS
Previous name was: a(n) = round(phi^(7^n)) where phi = 1.6180339887498948482... = (sqrt(5)+1)/2.
FORMULA
a(n) = G^(7^n) + (1 - G)^(7^n) = G^(7^n) + (-G)^(-7^n) where G is the golden ratio A001622. [Artur Jasinski, Oct 05 2008]
From Peter Bala, Nov 28 2022: (Start)
a(n) = Lucas(7^n).
a(n+1) = a(n)^7 + 7*a(n)^5 + 14*a(n)^3 + 7*a(n) with a(0) = 1.
a(n) == 1 (mod 7).
a(n+1) == a(n) (mod 7^(n+1)) for n >= 1 (a particular case of the Gauss congruences for the Lucas numbers).
In the ring of 7-adic integers, the limit_{n -> oo} a(n) exists and is a root of the quartic equation x^4 + 4*x^2 + 2 = 0. (End)
MAPLE
a := proc(n) option remember; if n = 0 then 1 else a(n-1)^7 + 7*a(n-1)^5 + 14*a(n-1)^3 +7*a(n-1) end if; end;
seq(a(n), n = 1..5); # Peter Bala, Nov 28 2022
MATHEMATICA
c = N[GoldenRatio, 1000]; Table[Round[c^(7^n)], {n, 1, 5}]
c = (1 + Sqrt[5])/2; Table[Expand[c^(7^n) + (1 - c)^(7^n)], {n, 0, 5}] (* Artur Jasinski, Oct 05 2008 *)
CROSSREFS
KEYWORD
nonn,bref
AUTHOR
Artur Jasinski, Sep 22 2008
EXTENSIONS
New name from Peter Bala, Nov 28 2022
STATUS
approved