OFFSET
1,1
MATHEMATICA
DivisorSigma[0, #] &@ Prepend[Select[Range[5000]^2, PrimeQ[DivisorSigma[1, #]] &], 2] (* Michael De Vlieger, Aug 14 2016, after Harvey P. Dale at A023194 *)
PROG
(PARI) lista(nn) = for(n=1, nn, if(isprime(sigma(n)), print1(numdiv(n), ", ")));
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, Aug 13 2016
STATUS
approved