OFFSET
1,3
COMMENTS
a(n) is the sum of the (n - 1)-th antidiagonal in A074712. - Ctibor O. Zizka, Mar 14 2025
Möbius transform of A189835(n). - Wesley Ivan Hurt, Jul 16 2025
FORMULA
From Wesley Ivan Hurt, Jul 16 2025: (Start)
a(n) = Sum_{d|n} A189835(d) * mu(n/d).
a(p^k) = p^(2*k)-p^k-k*p^k+k*p^(k-1) for p prime and k>=1. (End)
EXAMPLE
a(6) = Sum_{d|6} (6-d) * phi(6/d) = 5*phi(6) + 4*phi(3) + 3*phi(2) + 0*phi(1) = 5*2 + 4*2 + 3*1 + 0*1 = 21.
MAPLE
with(numtheory): seq(add((n-d)*phi(n/d), d in divisors(n)), n=1..80); # Ridouane Oudra, Jan 21 2024
MATHEMATICA
Table[Sum[(n - k)*EulerPhi[n/k^(1 - Ceiling[n/k] + Floor[n/k])] (1 - Ceiling[n/k] + Floor[n/k]), {k, n}], {n, 80}]
PROG
(PARI) a(n) = sumdiv(n, d, (n-d) * eulerphi(n/d)); \\ Michel Marcus, May 21 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, May 20 2021
STATUS
approved
