OFFSET
2,1
COMMENTS
A064920(a(n)) = n. a(n) = n iff n is prime.
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, k=1; until (m==n, k++; g=gpf(k); m=k / g + g - 1); write("b064923.txt", n, " ", k) ) } \\ Harry J. Smith, Sep 29 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Oct 14 2001
STATUS
approved