OFFSET
1,2
MATHEMATICA
Do[s=(Prime[n]+1)/DivisorSigma[1, n]; If[IntegerQ[s], Print[n]], {n, 1000000, 10000000}
Select[Range[81*10^5], Mod[Prime[#]+1, DivisorSigma[1, #]]==0&] (* Harvey P. Dale, May 03 2018 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Feb 05 2002
EXTENSIONS
More terms from Labos Elemer, Jun 12 2003
STATUS
approved