OFFSET
0,2
MAPLE
rad := n -> ifelse(n = 0, 1, NumberTheory:-Radical(n)):
a := n -> local k; add(rad(igcd(n, k)), k = 0..n):
seq(a(n), n = 0..60);
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Jul 12 2024
STATUS
approved