cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-1 of 1 results.

A383278 The number of integers k such that A034444(k) * k <= n.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 5, 5, 6, 6, 7, 7, 8, 8, 8, 8, 9, 9, 10, 10, 11, 11, 11, 11, 11, 11, 12, 12, 13, 13, 13, 13, 14, 14, 15, 15, 15, 15, 15, 15, 16, 16, 17, 17, 18, 18, 18, 18, 19, 19, 20, 20, 21, 21, 22, 22, 23, 23, 24, 24, 24, 24, 24, 24, 24, 24
Offset: 1

Views

Author

Amiram Eldar, Apr 21 2025

Keywords

Comments

The number of terms of A383276 not exceeding n.

References

  • József Sándor, Dragoslav S. Mitrinovic, and Borislav Crstici, Handbook of Number Theory I, Springer Science & Business Media, 2005, Chapter II, section 31, page 72.

Crossrefs

Partial sums of A383277.
The unitary analog of A356005.

Programs

  • Mathematica
    Accumulate[Table[DivisorSum[n, 1 &, # * 2^PrimeNu[#] == n &], {n, 1, 100}]]
    (* second program: *)
    f[n_] := Module[{e = IntegerExponent[n, 2], w}, w = PrimeNu[n/2^e]; If[e > w + 1 || e == w, 1, 0]]; Accumulate[Array[f, 100]]
  • PARI
    list(lim) = my(s = 0); for(n = 1, lim, s += sumdiv(n, d, (1 << omega(d)) * d == n); print1(s, ", "));
    
  • PARI
    f(n) = {my(e = valuation(n, 2), w = omega(n >> e)); e > w + 1 || e == w;}
    list(lim) = my(s = 0); for(n = 1, lim, s += f(n); print1(s, ", "));

Formula

a(n) = Sum_{k=1..n} A383277(k).
a(n) = (c + o(1)) * n / sqrt(log(n)), where c = (1/sqrt(Pi)) * Product_{p prime} (p-1/2)/sqrt(p*(p-1)) = A087197 * A345288 = 0.61890644913204789046... (Abbott and Subbarao, 1989).
Showing 1-1 of 1 results.