OFFSET
1,1
COMMENTS
From Aloe Poliszuk, Oct 23 2025: (Start)
a(n) is never squarefree.
Also, numbers k such that Sum_{d|k} (k/d * A107078(d) - 1) = pi(k). (End)
EXAMPLE
56 is in the sequence since sigma(56) = 120 = 96 + 8 + 16 = psi(56) + tau(56) + pi(56).
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); forprime(p = 3, lim, m++; for(k = q, p-1, f = factor(k); if(sigma(f) == dedpsi(f) + numdiv(f) + m, print1(k, ", "))); q = p); } \\ Amiram Eldar, Oct 17 2025
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
S. I. Dimitrov, Oct 17 2025
EXTENSIONS
a(7)-a(12) from Amiram Eldar, Oct 17 2025
STATUS
approved
