OFFSET
1,3
COMMENTS
a(n) = 1 if and only if n + 1 is prime.
Next term a(13), if it exists, is greater than 1000000. Other unknown terms are at index: 25, 29, 41, 47, 49, 53, 55, 69, 73, 79. (See links). - Jeppe Stig Nielsen, Aug 08 2014
LINKS
Guenter Loeh, Generalized Cullen primes n b^n + 1
Jeppe Stig Nielsen, Table of a(n) for n = 1..100, with unknown entries denoted by question marks
EXAMPLE
1*3^1 + 1 = 4 is not prime. 2*3^2 + 1 = 19 is prime. Thus, a(3) = 2.
PROG
(PARI) a(n) = for(k=1, 4000, if(ispseudoprime(k*n^k+1), return(k)))
n=1; while(n<100, print1(a(n), ", "); n++)
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Derek Orr, Apr 02 2014
STATUS
approved