OFFSET
1,1
EXAMPLE
(2^521-1)*6^959+1 is prime
MATHEMATICA
Union[Flatten[Table[Select[p*6^Range[0, 30] + 1, # < 10^20 && PrimeQ[#] &], {p, {3, 7, 31, 127, 8191, 131071, 524287, 2147483647, 2305843009213693951}}]]]
Union[Flatten[Table[Select[(2^MersennePrimeExponent[n]-1)6^Range[ 0, 20]+ 1, PrimeQ], {n, 20}]]] (* Harvey P. Dale, Aug 15 2021 *)
PROG
(Magma) // For the Mersenne prime 2^521-1, for example:
[n: n in [0..1000] | IsPrime((2^521-1)*6^n+1)];
CROSSREFS
KEYWORD
nonn
AUTHOR
Gilbert Mozzo, Feb 18 2011
EXTENSIONS
Partially edited by N. J. A. Sloane, Feb 19 2011
STATUS
approved