OFFSET
1,1
COMMENTS
LINKS
EXAMPLE
MATHEMATICA
Block[{s = Partition[MersennePrimeExponent@ Range@ 45, 2, 1], t}, t = Map[Differences, s][[All, 1]]; Map[s[[FirstPosition[t, #][[1]] ]] &, Union@ FoldList[Max, t]]][[All, 1]] (* Michael De Vlieger, Jan 31 2018 *)
PROG
(PARI) LL(e) = my(n, h); n = 2^e-1; h = Mod(2, n); for (k=1, e-2, h=2*h*h-1); return(0==h) \\ after Joerg Arndt in A000043
my(r=0, p=2); forprime(q=3, , if(LL(q), if(q-p > r, print1(p, ", "); r=q-p); p=q))
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Felix Fröhlich, Jan 30 2018
EXTENSIONS
a(21) from Amiram Eldar, Oct 17 2024
STATUS
approved