OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
17-1 = 2^4, 1297-1 = 2^4*3^4, 1336337-1 = 2^4*17^4, 4477457-1 = 2^4*23^4, ...
MATHEMATICA
Prepend[Select[Table[Prime[n], {n, 600000}], Length[Union[Last/@FactorInteger[#-1]]]==1&&Union[Last/@FactorInteger[#-1]]=={4}&], 2]
seq[lim_] := Select[Select[Range[Floor[Surd[lim-1, 2]]], SquareFreeQ]^4 + 1, PrimeQ]; seq[10^6] (* Amiram Eldar, Jan 18 2025 *)
PROG
(PARI) list(lim) = select(isprime, apply(x -> x^4 + 1, select(issquarefree, vector(sqrtnint(lim-1, 4), i, i)))); \\ Amiram Eldar, Jan 18 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Apr 09 2011
EXTENSIONS
a(12)-a(22) from Donovan Johnson, Apr 10 2011
a(1) = 2 inserted and a(23)-a(25) added by Amiram Eldar, Jan 18 2025
STATUS
approved