A360164 a(n) is the sum of the square roots of the unitary divisors of n that are odd squares.

1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 8, 6, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 6, 1, 1, 1, 1, 1, 10, 1, 1, 1, 1, 1, 1

Offset

1,9

Comments

First differs from A336649 at n = 27.

The unitary analog of A360163.

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000

Formula

a(n) = Sum_{d|n, gcd(d, n/d)=1, d odd square} sqrt(d).

a(n) = A360162(n) if n is not of the form (2*m - 1)*4^k where m >= 1, k >= 1 (A108269).

Multiplicative with a(2^e) = 1, and for p > 2, a(p^e) = p^(e/2) + 1 if e is even and 1 if e is odd.

Dirichlet g.f.: (zeta(s)*zeta(2*s-1)/zeta(3*s-1))*(2^(3*s)-2^(s+1))/(2^(3*s)-2).

Sum_{k=1..n} a(k) ~ (2*n/Pi^2)*(log(n) + 3*gamma - 1 + log(2) - 3*zeta'(2)/zeta(2)), where gamma is Euler's constant (A001620).

Mathematica

f[p_, e_] := If[OddQ[e], 1, p^(e/2) + 1]; f[2, e_] := 1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

Prog

(PARI) a(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i, 1] == 2, 1, if(f[i, 2]%2, 1, f[i, 1]^(f[i, 2]/2) + 1))); }

Crossrefs

Cf. A001620, A056624, A069290, A108269, A306016, A336649, A358347, A360163.

Sequence in context: A374325 A368169 A367418 * A360163 A336649 A370703

Adjacent sequences: A360161 A360162 A360163 * A360165 A360166 A360167

Keyword

nonn,easy,mult

Author

Amiram Eldar, Jan 29 2023

Status

approved