A333844 G.f.: Sum_{k>=1} k * x^(k^4) / (1 - x^(k^4)).

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3

Offset: 1

Author

Ilya Gutkovskiy, Apr 07 2020

Keywords

• nonn
• mult

Comments

Sum of 4th roots of 4th powers dividing n.

Links

• Amiram Eldar, Table of n, a(n) for n = 1..10000
• A. Dixit, B. Maji, and A. Vatwani, Voronoi summation formula for the generalized divisor function sigma_z^k(n), arXiv:2303.09937 [math.NT], 2023, sigma(z=1,k=4,n).

Crossrefs

Cf. A000203, A053164, A063775, A069290, A300909, A333843.

Programs

Mathematica

nmax = 112; CoefficientList[Series[Sum[k x^(k^4)/(1 - x^(k^4)), {k, 1, Floor[nmax^(1/4)] + 1}], {x, 0, nmax}], x] // Rest
Table[DivisorSum[n, #^(1/4) &, IntegerQ[#^(1/4)] &], {n, 112}]
f[p_, e_] := (p^(Floor[e/4] + 1) - 1)/(p - 1); a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Dec 01 2020 *)

Formula

Dirichlet g.f.: zeta(s) * zeta(4*s-1).

If n = Product (p_j^k_j) then a(n) = Product ((p_j^(floor(k_j/4) + 1) - 1)/(p_j - 1)).

Sum_{k=1..n} a(k) ~ zeta(3)*n + zeta(1/2)*sqrt(n)/2. - Vaclav Kotesovec, Dec 01 2020