A098819 Least k such that k*(Mersenne_prime(n)^2) + 1 is prime.

2, 4, 16, 12, 22, 28, 30, 22, 138, 12, 102, 180, 426, 1444, 9502, 6568, 2212, 5380, 9942, 3034, 24222, 31530, 2962, 7848, 93280, 3718

Offset: 1

Author

Pierre CAMI, Oct 08 2004

Keywords

• more
• nonn

Comments

Subsequence of A035092. - Amiram Eldar, Mar 01 2019

Examples

			2*((2^2-1)^2)+1=19 prime so a(1)=2.
		

Crossrefs

Cf. A000043, A000668, A035092.

Programs

Mathematica

a[n_] := Module[{k = 1, m = (2^MersennePrimeExponent[n] - 1)^2}, While[! PrimeQ[k*m + 1], k++]; k]; Array[a, 20] (* Amiram Eldar, Mar 01 2019 *)

Extensions

a(19)-a(22) from Donovan Johnson, Mar 23 2008

a(23)-a(26) from Amiram Eldar, Mar 01 2019