OFFSET
1,5
COMMENTS
This is the phi- function in Sandor and Atanassof.
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..10000
József Sándor and Krassimir Atanassov, Some new arithmetic functions, Notes on Number Theory and Discrete Mathematics, Volume 30, 2024, Number 4, Pages 851-856.
FORMULA
From Amiram Eldar, Aug 19 2025: (Start)
Multiplicative with a(p^e) = (p-1)*p^(e-1) - 1.
Dirichlet g.f.: zeta(s-1) * zeta(s) * Product_{p prime} (1 - 3/p^s + 1/p^(2*s-1) + 1/p^(2*s)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{p prime} (1 + 1/(p+1) - (p+1)/p^2) = 0.39439177573628632634... . (End)
MATHEMATICA
A383954[n_] := If[n == 1, 1, Times @@ (EulerPhi[Power @@@ FactorInteger[n]] - 1)];
Array[A383954, 100] (* Paolo Xausa, Aug 19 2025 *)
PROG
(PARI) a(n) = my(f=factor(n)); prod(k=1, #f~, p=f[k, 1]; eulerphi(f[k, 1]^f[k, 2])-1);
CROSSREFS
KEYWORD
nonn,mult,easy
AUTHOR
Michel Marcus, Aug 19 2025
STATUS
approved
