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.

A038627 Number of solutions x to n * pi(x) = x, where pi(x) = number of primes <= x.

Original entry on oeis.org

0, 4, 3, 3, 6, 7, 6, 6, 3, 9, 1, 18, 11, 12, 21, 3, 10, 33, 31, 32, 24, 8, 13, 32, 35, 4, 15, 9, 15, 26, 22, 24, 9, 3, 14, 55, 36, 3, 65, 52, 33, 139, 42, 2, 85, 25, 7, 96, 16, 33
Offset: 1

Views

Author

Keywords

Comments

Equivalently, a(n) is number of solutions x to the equation pi(n*x) = x. - Farideh Firoozbakht, Jan 09 2005 [For example, a(2) = 4 because 1, 2, 3 & 4 are all solutions of pi(2*x) = x and a(11) = 1 because 15927 is the only solution of the equation pi(11*x) = x.]
S. W. Golomb proved that a(n) > 0 for each integer n > 1. - Carlo Sanna, Nov 09 2015

Examples

			11*pi(x) = x has only 1 solution, so a(11) = 1.
		

Crossrefs

Programs

  • Mathematica
    (* Assumes upper and lower bounds are as defined in A038626. *)
    xmin = .5; xmax = 2;
    Join[{0},Table[c = 0; x = Floor[2.4*xmin]; x1 = 2.7*xmax + 7;
      xmin = Infinity; xmax = 0; While[x <= x1,
       If[x == PrimePi[n x], c++; xmin = Min[x, xmin];
    xmax = Max[x, xmax]]; x++]; c, {n, 2, 15}]] (* Robert Price, Mar 28 2020 *)

Extensions

One more term from Labos Elemer, Sep 05 2003
a(24)-a(26) from Labos Elemer, Sep 12 2003
Edited by N. J. A. Sloane, Sep 06 2008 at the suggestion of R. J. Mathar
a(27)-a(29) from David Radcliffe, Sep 10 2014
a(29) corrected and a(30)-a(50) obtained from the A038625 values computed by Jan Büthe. - Giovanni Resta, Aug 31 2018