A302089 Primes of form 5*k^k + 4.
139, 15629, 1937102449
Offset: 1
Keywords
Crossrefs
Programs
-
PARI
lista(nn) = forstep(n=1, nn, 2, if(ispseudoprime(p=5*n^n+4), print1(p, ", "))); \\ Altug Alkan, Apr 01 2018
This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
lista(nn) = forstep(n=1, nn, 2, if(ispseudoprime(p=5*n^n+4), print1(p, ", "))); \\ Altug Alkan, Apr 01 2018
a(2) = 9 = 3^2, a(3) = 55 = 5*11, a(4) = 513 = 3 ^ 3 * 19.
with(numtheory):for n from 0 to 50 do: x:=2*n^n + 1 : if type(x,prime)=false then print (x):else fi:od:
Select[Table[2n^n+1,{n,20}],CompositeQ] (* Harvey P. Dale, Jun 21 2015 *)
Comments