OFFSET
1,3
COMMENTS
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..10000 (terms 1..160 from Hal M. Switkay)
MATHEMATICA
s = {1}; Do[If[GreaterEqual @@ (f = FactorInteger[n])[[;; , 2]] && PrimePi[f[[-1, 1]]] == Length[f], AppendTo[s, EulerPhi[n]]], {n, 2, 3500}]; s (* Amiram Eldar, Dec 19 2022 *)
PROG
(PARI) is(n) = {if(n==1, return(1)); my(f = factor(n)); f[#f~, 1] == prime(#f~) && vecsort(f[, 2], , 4) == f[, 2]}; \\ A025487
lista(nn) = apply(eulerphi, select(is, [1..nn])); \\ Michel Marcus, Dec 19 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Hal M. Switkay, Dec 18 2022
STATUS
approved