OFFSET
1,1
COMMENTS
Sequence consists of 2 and the union of the Mersenne primes (A000668) and the Fermat primes (A019434).
a(18) has 157 digits and is too large to include. - Ray Chandler, Jun 22 2009
From Wolfdieter Lang, Mar 28 2012: (Start)
The sequence of exponents k of 2 for these Mersenne or Fermat primes is k=[0, 1, 2, 3, 4, 5, 7, 8, 13, 16, 17, 19, 31, 61, 89, 107, 127,...], n>=1, if 3 is considered as Fermat prime.
The second exponent is 2 if 3 is considered as Mersenne prime. The next exponent k(18)=521 for the Mersenne prime a(18).
r(a(n)) := sqrt(a(n)*2^(k(n)+2) + 1) is a solution of the congruence x^2 == 1 (mod a(n)*2^(k(n)+1)), n>=1. The sequence k is given in the preceding comment. If n=2 then, besides r(3) = 5, also sqrt(3*2^(2+2) + 1) = 7 is an incongruent solution mod 12 if the exponent 2 is chosen. For n=2 there are altogether four incongruent solutions of x^2 == 1 (mod 12), namely 1, 12-1 = 11, r(3)=5 and 12-5 = 7. If n=1 there are altogether two incongruent solutions, namely 1 and r(2) = 3 = 4-1 (a trivial solution == -1 (mod 4)). For n>=3 there are eight incongruent solutions, and besides the trivial (positive) ones, 1 and a(n)*2^(k(n)+1) - 1, one has a nontrivial pair r(a(n)) and a(n)*2^(k(n)+1) - r(a(n)). For the number of incongruent solutions of the congruence x^2 == 1 (mod n) see A060594. For the r(a(n)) values see A210844.
r(a(n)), n>=1, also solves the congruence x^2 == 1 (Modd a(n)*2^(k(n)+1)), because floor((r(a(n))^2)/(a(n)*2^(k(n)+1)) is even. For Modd n (not to be confused with mod n) see a comment on A203571.
(End)
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..22
EXAMPLE
From Wolfdieter Lang, Mar 28 2012: (Start)
Solutions to the congruence x^2 == 1 (mod a(n)*2^(k(n)+1):
n=3: r(5) = sqrt(5*2^(2+2) + 1) = 9. 9^2 = 81 == 1 (mod 5*8).
The companion solution is 40-9 = 31. Because floor(81/40)=2 is even, 81 == 1 (Modd 40) also.
n=4: r(7) = sqrt(7*2^(3+2) + 1) = 15. 15^2 = 225 == 1 (mod 7*16). The companion solution is 112-15 = 97. Because floor(225/112)=2 is even, 225 == 1 (Modd 112) also.
n=7: r(127) = sqrt(127*2^(7+2) + 1) = 255. 255^2 == 1 (mod 127*2^8). The companion solution is 32512-255 = 32257.
(End)
MATHEMATICA
Select[Prime[Range[30000]], Length[FactorInteger[#-1]]==1 || Length[FactorInteger[#+1]]==1&] (* Vladimir Joseph Stephan Orlovsky, Feb 03 2012 *)
Select[Union[Join @@ Array[2^# + {-1, +1} &, 140, 0]], PrimeQ] (* Michael De Vlieger, Oct 23 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Aug 07 2008
EXTENSIONS
More terms from R. J. Mathar, Jan 23 2009
a(17) from Ray Chandler, Jun 22 2009
STATUS
approved