OFFSET
1,4
COMMENTS
Sum of the divisors d of n such that phi(d) does not divide sigma(d).
Inverse Möbius transform of n * (1 - c(n)), where c = A351114.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
g:= proc(n) option remember; numtheory:-sigma(n) mod numtheory:-phi(n) <> 0 end proc:
f:= n -> convert(select(g, numtheory:-divisors(n)), `+`):
map(f, [$1..100]); # Robert Israel, Aug 26 2025
MATHEMATICA
Table[Sum[d (Ceiling[DivisorSigma[1, d]/EulerPhi[d]] - Floor[DivisorSigma[1, d]/EulerPhi[d]]), {d, Divisors[n]}], {n, 100}]
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Wesley Ivan Hurt, Jul 26 2025
STATUS
approved
