OFFSET
1,4
COMMENTS
a(1) = 0 by convention. - Antti Karttunen, Oct 30 2017
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16384
EXAMPLE
For n=4 phi(n')=phi(4')=phi(4)=2.
For n=7 phi(n')=phi(7')=phi(1)=1.
MATHEMATICA
Array[EulerPhi@ If[Abs@ # < 2, 0, # Total[#2/#1 & @@@ FactorInteger[Abs@ #]]] &, 75] (* Michael De Vlieger, Oct 30 2017, after Michael Somos at A003415 *)
PROG
(PARI) rd(n) = {local(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]; )); }
a(n) = eulerphi(rd(n)); \\ Michel Marcus, Sep 24 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Luca Brigada Villa, Sep 24 2013
EXTENSIONS
Description/formula clarified and more terms added by Antti Karttunen, Oct 30 2017
STATUS
approved