STATUS
reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
(Python)
from sympy import factorint
def ok(n): return n>1 and (n-1)%sum(p*e for p, e in factorint(n).items())==0
print([k for k in range(10**4) if ok(k)]) # Michael S. Branicky, Sep 30 2022
approved
editing
(MAGMAMagma) [k:k in [2..6200]| IsIntegral((k-1)/( &+[m[1]*m[2]: m in Factorization(k)]))]; // Marius A. Burtea, Sep 16 2019
proposed
approved
editing
proposed
(MAGMA) [k:k in [2..6200]| IsIntegral((k-1)/( &+[m[1]*m[2]: m in Factorization(k)]))]; // Marius A. Burtea, Sep 16 2019
proposed
editing
editing
proposed
For example, 21=7x3, 7+3=10 which divides 21-1=20.
proposed
editing