OFFSET
1,1
COMMENTS
No more terms less than 1.6*10^7.
MATHEMATICA
Select[Range[2, 1000], EulerPhi[#] == Plus @@ Times @@@ FactorInteger[#] &] (* Amiram Eldar, Jun 27 2023 *)
PROG
(Python)
from sympy import factorint, totient
A001414 = lambda k: sum(p*e for p, e in factorint(k).items())
def g():
k = 2
while True:
if A001414(k) == totient(k): yield(k)
k += 1
for a_n in g():
print(a_n)
(PARI) is(k) = my(f=factor(k)); f[, 1]~*f[, 2] == eulerphi(f); \\ Amiram Eldar, Jun 27 2023
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Darío Clavijo, Jun 26 2023
STATUS
approved