A058911 Numbers k such that k^k + k + 1 is prime.
0, 1, 2, 3, 6, 9, 462
Offset: 1
Examples
a(2) = 2 because 2^2 + 2 + 1 = 7.
Programs
-
Mathematica
Do[ If[ PrimeQ[ n^n + n + 1], Print[n]], {n, 1, 700} ] Join[{0},Select[Range[470],PrimeQ[#^#+#+1]&]] (* Harvey P. Dale, Dec 11 2022 *)
-
PARI
f2(n,k) = for(x=1,n,y=x^x+x+k;if(ispseudoprime(y),print1(x","))) \\ Cino Hilliard, Jan 07 2005
-
PFGW
ABC2 $a^$a + $a + 1 a: from 0 to 1000 // Jinyuan Wang, Mar 01 2020
Comments