OFFSET
1,3
FORMULA
MATHEMATICA
uphi[n_] := Product[{p, e} = pe; p^e - 1, {pe, FactorInteger[n]}];
a[n_] := If[n == 1, 1, Numerator[uphi[n]/n]];
Array[a, 100] (* Jean-François Alcover, Jan 10 2022 *)
PROG
(PARI) a(n)=my(f=factor(n)~); numerator(prod(i=1, #f, f[1, i]^f[2, i]-1)/n);
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Michel Marcus, Sep 26 2018
STATUS
approved