login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A101226
Numbers formed by the fourth nesting of pi(10^n).
1
0, 2, 5, 14, 44, 165, 725, 3607, 19624, 115127, 716177, 4675144, 31743253, 222710685, 1606581461, 11868655658, 89506275569, 687275710608, 5361744667232, 42423406909091, 339921273158836, 2754682455861719, 22553369622588850, 186375668479844528
OFFSET
1,2
FORMULA
a(n) = pi(pi(pi(pi(10^n)))) where pi(x) is the number of primes <= x.
a(n) = A000720(A101225(n)). - Amiram Eldar, Jun 29 2024
EXAMPLE
a(3) = pi(pi(pi(pi(10^3)))) = 5, the third entry in the table.
MATHEMATICA
f[n_] := Nest[PrimePi, 10^n, 4]; Table[ f[n], {n, 13}] (* Robert G. Wilson v, Dec 21 2004 *)
PROG
(PARI) nestpi(n, m) = { my(x, y, z); for(x=1, n, z=10^x; for(y=1, m, z=primepi(z)); print1(z", ")) } \\ try call nestpi(11, 4)
CROSSREFS
KEYWORD
nonn
AUTHOR
Cino Hilliard, Dec 15 2004
EXTENSIONS
a(10)-a(13) from Robert G. Wilson v, Dec 21 2004
a(14)-a(24) using Kim Walisch's primecount added by Amiram Eldar, Jun 29 2024
STATUS
approved