OFFSET
1,1
MATHEMATICA
us[n_Integer] := (d = Divisors[n]; l = Length[d]; k = 1; s = n; While[k < l, If[ GCD[ d[[k]], n/d[[k]] ] == 1, s = s + d[[k]]]; k++ ]; s); Do[m = n; If[ PrimeQ[ us[ DivisorSigma[1, n]]], Print[n]], {n, 1, 10^7/4} ]
CROSSREFS
KEYWORD
nonn,bref
AUTHOR
Robert G. Wilson v, Aug 21 2001
STATUS
approved