OFFSET
2,1
COMMENTS
A064916(a(n)) = n.
LINKS
Harry J. Smith, Table of n, a(n) for n = 2..1000
FORMULA
For n > 2: a(n) = if n is prime then n else 2*n - 2.
PROG
(PARI) lpf(n)= { local(f); f=factor(n); return(f[1, 1]) } { for (n=2, 1000, k=1; until (m==n, k++; L=lpf(k); m=k / L + L - 1); write("b064919.txt", n, " ", k) ) } \\ Harry J. Smith, Sep 29 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Oct 14 2001
STATUS
approved