OFFSET
1,2
COMMENTS
For n > 1: a(n)=n iff n is prime.
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
PROG
(Haskell)
a071889 n = gcd n $ a071888 n -- Reinhard Zumkeller, Jun 10 2012
(PARI) a(n)=if(n<2, return(1)); my(m=n+1); while(gcd(m++, n)==1 || !issquarefree(m), ); gcd(n, m-n) \\ Charles R Greathouse IV, May 29 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jun 10 2002
STATUS
approved