OFFSET
1,21
LINKS
FORMULA
a(n) = Sum_{p over distinct primes dividing n} [sigma(n) != 0 mod p].
MATHEMATICA
Table[Length[Select[FactorInteger[n][[All, 1]], Mod[DivisorSigma[ 1, n], #]!= 0&]], {n, 110}] (* Harvey P. Dale, Jan 15 2022 *)
PROG
(PARI) A336387(n) = if(1==n, 0, my(s=sigma(n)); #select(p -> (s%p), factor(n)[, 1]));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jul 25 2020
STATUS
approved
