OFFSET
2,1
COMMENTS
Well-defined since A064916(n) < n for nonprimes.
a(p) = p for all primes p.
LINKS
Harry J. Smith, Table of n, a(n) for n = 2..1000
Michael Gilleland, Some Self-Similar Integer Sequences
PROG
(PARI) lpf(n)= { local(f); f=factor(n); return(f[1, 1]) } { for (n=2, 1000, m=n; while (!isprime(m), L=lpf(m); m=m / L + L - 1); write("b064917.txt", n, " ", m) ) } \\ Harry J. Smith, Sep 29 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Oct 14 2001
STATUS
approved