|
PROG
|
(PARI) isok(m) = my(sn = sigma(m)); sumdiv(m, d, sn % d) == m;
(Python)
from itertools import count, islice
from sympy import divisor_sigma, divisors
def A349866gen(): # generator of terms
return filter(lambda m: sum(divisor_sigma(m) % d for d in divisors(m, generator=True)) == m, count(1))
|