OFFSET
1,2
COMMENTS
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..16383
MAPLE
a:= n-> Bits[Xor](n, max(1, (numtheory[divisors](n) minus {n})[])):
seq(a(n), n=1..100); # Alois P. Heinz, May 15 2016
MATHEMATICA
a[1] = 0; a[n_] := BitXor[n, Divisors[n][[-2]]]; Array[a, 100] (* Jean-François Alcover, Mar 12 2019 *)
PROG
(PARI) gpd(n) = if(n==1, 1, n/factor(n)[1, 1]); \\ A032742
a(n) = bitxor(n, gpd(n)); \\ Michel Marcus, Mar 12 2019
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Reinhard Zumkeller, May 02 2005
STATUS
approved