OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
László Tóth, On the Bi-Unitary Analogues of Euler's Arithmetical Function and the Gcd-Sum Function, JIS 12 (2009), Article 09.5.2, function S**(n).
FORMULA
Sum_{k=1..n} a(k) = c * n^2 / 2 + O(n * log(n)^2), where c = A065486. - Amiram Eldar, Dec 22 2023
MAPLE
MATHEMATICA
beta[n_] := Times @@ Transpose[FactorInteger[n]][[2]]; phi[x_, n_] := Sum[Boole[ GCD[k, n] == 1 ], {k, 1, x}]; a[n_] := DivisorSum[n, beta[#] * phi[n/#, #] &, GCD[#, n/#] == 1 &]; Array[a, 100] (* Amiram Eldar, Sep 22 2019 *)
CROSSREFS
KEYWORD
AUTHOR
R. J. Mathar, Jul 21 2016
STATUS
approved