OFFSET
1,2
EXAMPLE
16 is a term since sigma(16) = 31 = 24 + 6 + 1 = psi(16) + pi(16) + omega(16).
PROG
(PARI) dedpsi(f) = prod(k = 1, #f~, (f[k, 1]+1) * f[k, 1]^(f[k, 2]-1));
list(lim) = {my(q = 2, m = 0, f); print1(1, ", "); forprime(p = 3, lim, m++; for(k = q, p-1, f = factor(k); if(sigma(f) == dedpsi(f) + m + omega(f), print1(k, ", "))); q = p); } \\ Amiram Eldar, Jan 20 2026
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
S. I. Dimitrov, Jan 20 2026
EXTENSIONS
a(7)-a(8) from Amiram Eldar, Jan 20 2026
STATUS
approved
