OFFSET
1,1
COMMENTS
Note that for all primes p > 3, sigma(3p) = sigma(p) + sigma(2p).
MATHEMATICA
lst={}; Do[p=Prime[n]; k=1; While[DivisorSigma[1, p+k]!=DivisorSigma[1, p]+DivisorSigma[1, k], k++ ]; If[k==2p, AppendTo[lst, p]], {n, 3, 200}]; lst
CROSSREFS
KEYWORD
nonn
AUTHOR
T. D. Noe, Jan 20 2004
STATUS
approved