A318449 Numerators of the sequence whose Dirichlet convolution with itself yields A001511, the 2-adic valuation of 2n.

1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 3, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 3, 3, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 3, 1, 1, 1, 1, 1, 35, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 3, 3, 3, 1, 1, 1, 1, 1

Offset: 1

Author

Antti Karttunen and Andrew Howroyd, Aug 29 2018

Keywords

• nonn
• frac

Links

• Antti Karttunen, Table of n, a(n) for n = 1..65537

Crossrefs

Cf. A001511, A318450 (denominators).

Programs

Mathematica

a1511[n_] := IntegerExponent[2n, 2];
f[1] = 1; f[n_] := f[n] = 1/2 (a1511[n] - Sum[f[d] f[n/d], {d, Divisors[ n][[2 ;; -2]]}]);
Table[f[n] // Numerator, {n, 1, 105}] (* Jean-François Alcover, Sep 13 2018 *)

PARI

up_to = 65537;
A001511(n) = 1+valuation(n,2);
DirSqrt(v) = {my(n=#v, u=vector(n)); u[1]=1; for(n=2, n, u[n]=(v[n]/v[1] - sumdiv(n, d, if(d>1&&dA317937.
v318449_51 = DirSqrt(vector(up_to, n, A001511(n)));
A318449(n) = numerator(v318449_51[n]);

Formula

a(n) = numerator of f(n), where f(1) = 1, f(n) = (1/2) * (A001511(n) - Sum_{d|n, d>1, d 1.

Sum_{k=1..n} A318449(k) / A318450(k) ~ n * sqrt(2/(Pi*log(n))) * (1 + (1 - gamma/2 + log(2)/2)/(2*log(n))), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, May 10 2025