OFFSET
1,2
COMMENTS
a(20) > 10^12. [From Donovan Johnson, Mar 14 2010]
EXAMPLE
p(5) = 11 is congruent to 5 mod sigma(5) = 6; so 5 is a term of the sequence.
MATHEMATICA
Select[ Range[ 1, 10^6 ], Mod[ #, DivisorSigma[ 1, # ] ] == Mod[ Prime[ # ], DivisorSigma[ 1, # ] ] & ]
Select[ Range[ 1, 4 10^7 ], # == Mod[ Prime[ # ], DivisorSigma[ 1, # ]] & ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Joseph L. Pe, Jan 11 2002
EXTENSIONS
More terms from Robert G. Wilson v, Jan 13 2002
a(13)-a(19) from Donovan Johnson, Mar 14 2010
STATUS
approved