OFFSET
1,18
EXAMPLE
a(18)=3 as sigma(18)+1=40, sigma(40)+1=91, sigma(91)+1=113, which is prime
PROG
(PARI) sr(m, c)=if (c>0, m=sigma(m)+1; c--; sr(m, c), m) for (n=1, 100, x=1; while (!isprime(sr(n, x)) && x<200, x++); if (x==200, x=0); print1(x", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Jon Perry, Jan 09 2003
STATUS
approved