%I #29 Sep 25 2024 23:13:04
%S 10,98,9602,92198402,8500545331353602,
%T 72259270930397519221389558374402,
%U 5221402235392591963136699520829303150191924374488750728808857602
%N a(n) = a(n-1)^2 - 2 with a(1) = 10.
%C This is the Lucas-Lehmer sequence with starting value u(1) = 10 and the position of the zeros when it is reduced mod(2^p - 1) also gives the position of the Mersenne primes. As we have started with n = 1, these will occupy the (p - 1)th positions in the sequence. For example, the first 12 terms mod(2^13 - 1) are 10, 98, 1411, 506, 2113, 672, 1077, 4996, 2037, 4721, 128, 0 and hence 8191 is a Mersenne prime. The radicals in the above closed forms are the solutions to x^2 - 10x + 1 = 0.
%H Gabriel Klambauer, <a href="https://doi.org/10.2307/2321992">Summation of Series</a>, Amer. Math. Monthly, Vol. 87, No. 2 (Feb., 1980), pp. 128-130.
%H Raphael M. Robinson, <a href="http://dx.doi.org/10.1090/S0002-9939-1954-0064787-4">Mersenne and Fermat Numbers</a>, Proceedings of the American Mathematical Society, Vol. 5, No. 5. (October 1954), pp. 842-846.
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Engel_expansion">Engel expansion</a>
%F a(n) = 2*cosh(2^(n-1)*log(5 + 2*sqrt(6))) = exp(2^(n-1)*log(5 + 2*sqrt(6))) + exp(2^(n-1)*log(5 - 2*sqrt(6))) = (5 + 2*sqrt(6))^(2^(n-1)) + (5 - 2*sqrt(6))^(2^(n-1)) = ceiling(exp(2^(n-1)*log(5 + 2*sqrt(6)))) = ceiling((5 + 2*sqrt(6))^(2^(n-1))).
%F From _Peter Bala_, Feb 01 2022: (Start)
%F Product_{n >= 1} (1 + 2/a(n)) = (1/2)*sqrt(6); Product_{n >= 1} (1 - 1/a(n)) = (4/11)*sqrt(6).
%F Engel expansion of 5 - sqrt(24) = 1/a(1) + 1/(a(1)*a(2)) + 1/(a(1)*a(2)*a(3)) + .... See Klambauer, p. 130. (End)
%e a(4) = 2*cosh(2^3*log(5 + 2*sqrt(6))) = 92198402.
%t a[1] = 10; a[n_] := a[n] = a[n - 1]^2 - 2; a[#]&/@Range[7]
%o (Python)
%o A135927 = [10]
%o for n in range(1, 8): A135927.append(A135927[-1]**2-2)
%o print(A135927) # _Karl-Heinz Hofmann_, Feb 01 2022
%Y Cf. A000668, A000043, A003010, A095847, A001566, A135928, A246723.
%K nonn,easy
%O 1,1
%A _Ant King_, Dec 07 2007