OFFSET
1,1
COMMENTS
Primes in A309562. - Robert Israel, Aug 09 2019
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000
EXAMPLE
p = 10181, 1+p^4 = 10743894862923122 = 2*17*1657*4657*5113*8009, so the largest prime factor is 8009 < p = 10181.
MAPLE
filter:= proc(p) max(numtheory:-factorset(p^4+1)) < p end proc:
select(filter, [seq(ithprime(i), i=1..20000)]); # Robert Israel, Aug 09 2019
MATHEMATICA
<<NumberTheory`NumberTheoryFunctions` Select[Prime[Range[15000]], Max[PrimeFactorList[1 + #^4]] < # &] (* Ray Chandler, Jan 08 2005 *)
Select[Prime[Range[15000]], FactorInteger[#^4+1][[-1, 1]]<#&] (* Harvey P. Dale, Feb 27 2017 *)
PROG
(PARI) isok(p) = isprime(p) && (vecmax(factor(p^4+1)[, 1]) < p); \\ Michel Marcus, Jul 09 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Jan 05 2005
EXTENSIONS
Extended by Ray Chandler, Jan 08 2005
STATUS
approved