OFFSET
1,2
COMMENTS
This corresponds to the numbers such that 3m^m+1 is prime, but these must all be even, m=2k, and therefore it is more natural to record the sequence of k=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.
MAPLE
q:= k-> isprime(3*(2*k)^(2*k)+1):
select(q, [$1..225])[]; # Alois P. Heinz, Aug 04 2025
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
