OFFSET
1,1
COMMENTS
Primes p such that the largest prime factor of p+1 is less than p gives A065091, odd primes.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) max(numtheory:-factorset(n^2+1))<n end proc:
select(filter, [seq(ithprime(i), i=1..1000)]); # Robert Israel, Aug 07 2019
MATHEMATICA
<<NumberTheory`NumberTheoryFunctions` Select[Prime[Range[250]], Max[PrimeFactorList[1 + #^2]] < # &] (* Ray Chandler, Jan 08 2005 *)
Select[Prime[Range[212]], Max[First /@ FactorInteger[#^2 + 1]] < # &] (* Jayanta Basu, Jul 01 2013 *)
PROG
(Magma) [p:p in PrimesUpTo(1500)|Max(PrimeDivisors(p^2+1)) lt p]; // Marius A. Burtea, Aug 07 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Aug 27 2002
STATUS
approved