A161471 Primes of the form k^k + k + 1.
2, 3, 7, 31, 46663, 387420499
Offset: 1
Keywords
Programs
-
Mathematica
Unprotect[Power]; Power[0, 0] = 1; Protect[Power]; lst={}; Do[p=n^n+n+1; If[PrimeQ[p], AppendTo[lst,p]], {n,0,100}]; lst Join[{2},Select[Table[k^k+k+1,{k,1000}],PrimeQ]] (* Harvey P. Dale, Oct 09 2022 *)
-
PARI
lista(nn) = for(k=0, nn, if(ispseudoprime(q=k^k+k+1), print1(q, ", "))); \\ Jinyuan Wang, Mar 01 2020
Formula
Extensions
Definition simplified by R. J. Mathar, Jun 12 2009
Comments