OFFSET
2,5
COMMENTS
a(p) = 0 for all primes p.
LINKS
Harry J. Smith, Table of n, a(n) for n = 2..1000
PROG
(PARI) lpf(n)= { local(f); f=factor(n); return(f[1, 1]) } { for (n=2, 1000, m=n; a=0; while (!isprime(m), L=lpf(m); m=m / L + L - 1; a++); write("b064918.txt", n, " ", a) ) } \\ Harry J. Smith, Sep 29 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Oct 14 2001
STATUS
approved