OFFSET
1,2
COMMENTS
a(11), if it exists, is greater than 10^9. - Lyle Blosser, Jan 18 2026
EXAMPLE
80 is in the sequence since phi(80) = 32 = 22 + 10 = pi(80) + tau(80).
MATHEMATICA
Select[Range[400], EulerPhi[#] == PrimePi[#] + DivisorSigma[0, #] &] (* Amiram Eldar, Oct 14 2025 *)
PROG
(Python)
from math import prod
from itertools import count, islice
from sympy import factorint, nextprime
def agen(): # generator of terms
pk, pik = 2, 0
for k in count(1):
if k == pk:
pk, pik = nextprime(pk), pik + 1
f = factorint(k)
phik = prod((p - 1)*p**(e-1) for p, e in f.items())
tauk = prod(e+1 for e in f.values())
if phik == pik + tauk:
yield k
print(list(islice(agen(), 10))) # Michael S. Branicky, Jan 19 2026
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
S. I. Dimitrov, Oct 14 2025
STATUS
approved
