OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
666 = 2*3^2*37, tau(666) = 12, prime(12) = 37, hence 666 is a term.
MATHEMATICA
Select[Range[2200], FactorInteger[#][[-1, 1]]==Prime[DivisorSigma[0, #]]&] (* Harvey P. Dale, Aug 13 2021 *)
PROG
(PARI) for(n=2, 3000, if(component(component(factor(n), 1), omega(n))==prime(numdiv(n)), print1(n, ", ")))
(PARI) is(k) = if(k > 1, my(f = factor(k)); f[#f~, 1] == prime(numdiv(f)), 0); \\ Amiram Eldar, Oct 27 2024
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Jun 08 2002
STATUS
approved
