OFFSET
1,1
COMMENTS
Inspired by A339466. See the references there.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
alias(pf = NumberTheory:-PrimeFactors): gpf := n -> max(pf(n)):
is_a := n -> isprime(n) and n < (gpf((n-1)/gpf(n-1)))^4:
select(is_a, [$5..1150]);
PROG
(PARI) gpf(n) = if (n==1, 1, vecmax(factor(n)[, 1])); \\ A006530
isok(p) = isprime(p) && (p < (gpf((p - 1)/gpf(p - 1)))^4); \\ Michel Marcus, Dec 14 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Dec 13 2020
STATUS
approved