OFFSET
1,2
LINKS
N. J. A. Sloane, Table of n, a(n) for n = 1..10000
FORMULA
G.f.: (1/(1 - x))*Sum_{k>=1} (1 - mu(k))*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, Sep 28 2017
a(n) ~ c * n^2 + O(n*log(n)), where c = Pi^2/12 - 3/Pi^2 = 0.518503... . - Amiram Eldar, Oct 27 2022
MATHEMATICA
Accumulate[Table[DivisorSigma[1, n] - EulerPhi[n], {n, 1, 50}]] (* Amiram Eldar, Oct 27 2022 *)
PROG
(PARI) lista(nmax) = {my(s = 0, f); for(n = 1, nmax, f = factor(n); s += sigma(f) - eulerphi(f); print1(s, ", ")); } \\ Amiram Eldar, Oct 27 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Sep 28 2017
STATUS
approved