OFFSET
1,1
PROG
(Python)
from sympy.ntheory.factor_ import antidivisors
def A388007(n):
ad = antidivisors(n)
sad = sum(ad)
return sum(sad%d for d in ad)
def ok(n): return A388007(n) == n
print([k for k in range(1, 2000) if ok(k)]) # Michael S. Branicky, Oct 17 2025
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Michel Marcus, Oct 17 2025
EXTENSIONS
a(7)-a(15) from Michael S. Branicky, Oct 18 2025
STATUS
approved
