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) gpf(n)= { local(f); f=factor(n)~; return(f[1, length(f)]) } { for (n=2, 1000, m=n; a=0; while (!isprime(m), g=gpf(m); m=m / g + g - 1; a++); write("b064922.txt", n, " ", a) ) } \\ Harry J. Smith, Sep 29 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Oct 14 2001
STATUS
approved