A383657 Numerator of Dirichlet g.f.: Sum_{n>=1} a(n)/n^s = zeta(s)^(3/2).

1, 3, 3, 15, 3, 9, 3, 35, 15, 9, 3, 45, 3, 9, 9, 315, 3, 45, 3, 45, 9, 9, 3, 105, 15, 9, 35, 45, 3, 27, 3, 693, 9, 9, 9, 225, 3, 9, 9, 105, 3, 27, 3, 45, 45, 9, 3, 945, 15, 45, 9, 45, 3, 105, 9, 105, 9, 9, 3, 135, 3, 9, 45, 3003, 9, 27, 3, 45, 9, 27, 3, 525, 3

Offset

1,2

Comments

In general, for m > 0, if Dirichlet g.f. is zeta(s)^m, then Sum_{j=1..n} a(j) ~ n*log(n)^(m-1)/Gamma(m) * (1 + (m-1)*(m*gamma - 1)/log(n)), where gamma is the Euler-Mascheroni constant A001620 and Gamma() is the gamma function.

Links

Vaclav Kotesovec, Table of n, a(n) for n = 1..10000

Formula

Sum_{k=1..n} A383657(k)/A383658(k) ~ 2*n*sqrt(log(n)/Pi) * (1 - (1 - 3*gamma/2)/(2*log(n))), where gamma is the Euler-Mascheroni constant A001620.

Mathematica

coeff=CoefficientList[Series[1/(1-x)^(3/2), {x, 0, 20}]//Normal, x]; dptTerm[n_]:=Module[{flist=FactorInteger[n]}, If[n==1, coeff[[1]], Numerator[Times@@(coeff[[flist[[All, 2]]+1]])]]]; Array[dptTerm, 73] (* Shenghui Yang, May 04 2025 *)

Prog

(PARI) for(n=1, 100, print1(numerator(direuler(p=2, n, 1/(1-X)^(3/2))[n]), ", "))

Crossrefs

Cf. A383658.

Cf. A046643, A046644, A256688, A256689, A256690, A256691, A256692, A256693.

Cf. A000005, A007425, A007426, A061200, A034695.

Sequence in context: A281759 A280467 A282264 * A281213 A282388 A131943

Adjacent sequences: A383654 A383655 A383656 * A383658 A383659 A383660

Keyword

nonn,frac,mult

Author

Vaclav Kotesovec, May 04 2025

Status

approved