OFFSET
1,2
EXAMPLE
MAPLE
f:= proc(n) local s, D, t;
D:= numtheory:-divisors(n);
s:= convert(D, `+`);
add(s mod t, t=D)
end proc:
select(t -> f(t) mod t = 0, [$1..20000]);
PROG
(PARI) f(n) = my(sn = sigma(n)); sumdiv(n, d, sn % d); \\ A318996
isok(k) = !(f(k) % k); \\ Michel Marcus, Feb 03 2021
CROSSREFS
KEYWORD
nonn,more
AUTHOR
J. M. Bergot and Robert Israel, Feb 03 2021
STATUS
approved