OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)
MAPLE
with(numtheory); A048968:=n->`if`(type(sigma(n)/tau(n), prime), n, NULL); seq(A048968(n), n=1..2400); # Wesley Ivan Hurt, Feb 04 2014
MATHEMATICA
Select[ Range[2400], PrimeQ[ DivisorSigma[1, #] / DivisorSigma[0, #] ]& ] (* Jean-François Alcover, Sep 24 2012 *)
PROG
(PARI) isok(k) = {my(f = factor(k), d = numdiv(f), s = sigma(f)); !(s % d) && isprime(s / d); } \\ Amiram Eldar, Jan 13 2025
CROSSREFS
KEYWORD
nonn,nice,changed
AUTHOR
EXTENSIONS
More terms from Jud McCranie
STATUS
approved