OFFSET
1,1
COMMENTS
The values n with a(n) = 0 listed here are conjectural. There is no prime < 10^30 appearing in the orbits of these values.
EXAMPLE
The orbit of 4 is 4, 11, 23,.... and the first prime in the orbit is 11. So a(4) = 11.
MATHEMATICA
f[ n_ ] := n + DivisorSigma[ 1, n ]; h[ n_ ] := Module[ { a = n}, While[ ! PrimeQ[ a ] && a < 10^30, a = f[ a ] ]; If [ ! PrimeQ[ a ], a = 0 ]; a ]; Table[ h[ i ], {i, 1, 40} ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Joseph L. Pe, Jan 30 2002
STATUS
approved