OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
Select[Range[1800], FactorInteger[#][[-1, 1]]FactorInteger[EulerPhi[#]][[-1, 1]] == #&] (* Harvey P. Dale, Mar 25 2023 *)
PROG
(PARI) for(n=1, 3000, if(vecmax(component(factor(n), 1))*vecmax(component(factor(eulerphi(n)), 1))==n, print1(n, ", ")))
(PARI) is(k) = if(k > 2, my(f = factor(k)); k == f[#f~, 1] * vecmax(factor(eulerphi(f))[, 1]), 0); \\ Amiram Eldar, Oct 28 2024
CROSSREFS
KEYWORD
easy,nonn,changed
AUTHOR
Benoit Cloitre, Jun 16 2002
STATUS
approved