OFFSET
1,4
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = A001615(n) mod n.
EXAMPLE
MAPLE
MATHEMATICA
a253628[n_] :=
Mod[DirichletConvolve[j, MoebiusMu[j]^2, j, #], #] & /@ Range@n; a253628[75] (* Michael De Vlieger, Jan 07 2015, after Jan Mangaldan at A001615 *)
PROG
(Sage)
[(n*mul(1+1/p for p in prime_divisors(n)))%n for n in [1..100]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Tom Edgar, Jan 06 2015
STATUS
approved