OFFSET
1,2
EXAMPLE
30915 is a term since sigma(30915) = 55200 = 49680 + 3333 + 3^7 = psi(30915) + pi(30915) + omega(30915)^7.
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)^7, print1(k, ", "))); q = p); } \\ Amiram Eldar, Feb 01 2026
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
S. I. Dimitrov, Feb 01 2026
EXTENSIONS
a(4)-a(5) from Amiram Eldar, Feb 02 2026
STATUS
approved
