OFFSET
1,2
COMMENTS
See Bacher link for the definition of Euclid-reduced.
LINKS
Roland Bacher, Euclid meets Popeye: The Euclidean Algorithm for 2X2 matrices, arXiv:2209.09529 [math.NT], 2022.
FORMULA
a(n) = Sum_{d^2|n} moebius(d)*A357259(n/d^2).
MATHEMATICA
f[n_] := DivisorSum[n, # + 1 - n/# &, #^2 >= n &]; a[n_] := DivisorSum[n, MoebiusMu[Sqrt[#]] * f[n/#] &, IntegerQ[Sqrt[#]] &]; Array[a, 100] (* Amiram Eldar, Sep 21 2022 *)
PROG
(PARI) f(n) = sumdiv(n, d, if (d^2 >= n, d + 1 -n/d)); \\ A357259
a(n) = sumdiv(n, d, if (issquare(d), moebius(sqrtint(d))*f(n/d)));
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Sep 21 2022
STATUS
approved