OFFSET
1,3
EXAMPLE
a(1) = 1 since 3 = 2^A000043(1) - 1 and 4*(1*3)^2 + 1 = 37 is prime.
MATHEMATICA
f[n_] := Module[{k = 1}, While[!PrimeQ[4*(k*n)^2 + 1], k++]; k]; f /@ (2^MersennePrimeExponent[Range[15]] - 1)(* Amiram Eldar, Jul 17 2021 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Pierre CAMI, Nov 05 2007
EXTENSIONS
Data corrected and a(23)-a(25) added by Amiram Eldar, Jul 17 2021
STATUS
approved