OFFSET

1,2

COMMENTS

This corresponds to the numbers such that 3m^m+1 is prime, but these must all be even, m=2n, and therefore it is more natural to record the sequence of n=m/2.

Next term > 15000. - Matevz Markovic, Oct 09 2012

EXAMPLE

a(1) = 1, because 2^2*3+1 = 13 is the smallest prime of this form.

a(2) = 2, because 4^4*3+1 = 769 is the next smallest prime of this form. a(3) = 3, because 6^6*3+1 = 139969 is again a prime.

PROG

(PARI) for(i=1, 9999, ispseudoprime(i^i*3+1)&print1(i/2, ", "))

CROSSREFS

KEYWORD

hard,more,nonn

AUTHOR

M. F. Hasler, Jul 10 2009

STATUS

approved