OFFSET
1,2
FORMULA
Dirichlet g.f.: zeta(s)^2 * Product_{p prime} (1 + p^(1-s) - p^(-s))^2.
Multiplicative with a(p^e) = (e-1)*p^2 + 2*p. - Amiram Eldar, Nov 25 2021
From Vaclav Kotesovec, Nov 26 2021: (Start)
Dirichlet g.f.: zeta(s-1)^2 * zeta(s)^2 * Product_{primes p} (1 + p^(1-2*s) - p^(2-2*s) - p^(-s))^2.
Let f(s) = Product_{primes p} (1 + p^(1-2*s) - p^(2-2*s) - p^(-s)), then
Sum_{k=1..n} a(k) ~ Pi^2 * f(2)^2 * n^2 / 144 * (Pi^2 * (2*log(n) + 4*gamma - 1 + 4*f'(2)/f(2)) + 24*zeta'(2)), where f(2) = Product_{primes p} (1 - 2/p^2 + 1/p^3) = A065464 = 0.428249505677094440218765707581823546121298513355936144..., f'(2) = f(2) * Sum_{primes p} log(p) * (3*p - 2) / (p^3 - 2*p + 1) = 0.6293283828324697510445630056425352981207558777167836747744750359407300858..., zeta'(2) = -A073002 and gamma is the Euler-Mascheroni constant A001620. (End)
MATHEMATICA
Table[Sum[Last[Select[Divisors[d], SquareFreeQ]] Last[Select[Divisors[n/d], SquareFreeQ]], {d, Divisors[n]}], {n, 1, 60}]
f[p_, e_] := (e - 1)*p^2 + 2*p; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 60] (* Amiram Eldar, Nov 25 2021 *)
PROG
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Ilya Gutkovskiy, Nov 25 2021
STATUS
approved