Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A161471
Primes of the form k^k + k + 1.
4
2, 3, 7, 31, 46663, 387420499
OFFSET
1,1
COMMENTS
The associated k are in A058911. - R. J. Mathar, Jun 12 2009
a(7) = 116035988662615798148247830...0301775295923724561430603629007 (1232 digits).
FORMULA
Intersection of A000040 and A066279. - R. J. Mathar, Jun 12 2009
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 *)
PROG
(PARI) lista(nn) = for(k=0, nn, if(ispseudoprime(q=k^k+k+1), print1(q, ", "))); \\ Jinyuan Wang, Mar 01 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Definition simplified by R. J. Mathar, Jun 12 2009
STATUS
approved