OFFSET
1,3
LINKS
MATHEMATICA
With[{nmax = 100}, Table[Total[PositionIndex[Divisible[n, #[[;; n]]]][True]], {n, nmax}] & [DivisorSigma[1, Range[nmax]]]]
PROG
(PARI) a(n) = vecsum(select(x->(!(n % sigma(x))), [1..n])); \\ Michel Marcus, Dec 11 2024
CROSSREFS
KEYWORD
nonn,new
AUTHOR
Paolo Xausa, Dec 11 2024
STATUS
approved