OFFSET
1,1
COMMENTS
All squares of primes (A001248) are terms.
LINKS
Michel Marcus, Table of n, a(n) for n = 1..1049 (terms up to 10^7)
EXAMPLE
4 is a term since sopf(4)=2 and 4/2 = 2 is a prime.
30 is a term since sopf(30)=10 and 30/10 = 3 is a prime.
PROG
(PARI) sopf(n)=vecsum(factor(n)[, 1]); \\ A008472
isokp(k) = my(q=k/sopf(k)); (denominator(q)==1) && isprime(q);
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Jul 22 2020
STATUS
approved