A069637 Number of prime powers <= n with exponents > 1.

0, 0, 0, 1, 1, 1, 1, 2, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

Offset

1,8

Comments

Counts A025475 without 1 = prime^0: a(n) = A085501(n) - 1. - Reinhard Zumkeller, Jul 03 2003

Counts the prime powers (A246655) without the primes. - Peter Luschny, Nov 18 2019

References

H. Sahu, K. Kar and B.S.K.R. Somayajulu, On the average order of pi*(n) - pi(n), Acta Cienc. Indica Math., Vol. 11 (1985), pp. 165-168.

József Sándor, Dragoslav S. Mitrinovic and Borislav Crstici, Handbook of Number Theory I, Springer Science & Business Media, 2005, Chapter VII, p. 237.

Links

Daniel Forgues, Table of n, a(n) for n=1..100000.

Formula

a(n) = A025528(n) - A000720(n) = A000720([n^(1/2)]) + A000720([n^(1/3)]) + ... . - Max Alekseyev, May 11 2009

Sum_{k=1..n} a(k) ~ (4/3) * n^(3/2)/log(n) + O(n^(3/2)/log(n)^2) (Sahu et al., 1985). - Amiram Eldar, Mar 07 2021

Maple

with(numtheory);
A069637 := proc(N) local ct, i; ct:=0;
for i from 1 to N do if not isprime(i) and nops(factorset(i))=1 then ct:=ct+1; fi; od; ct; end; # N. J. A. Sloane, Jun 05 2022

Mathematica

Table[Sum[PrimePi[n^(1/k)], {k, Log[2, n]}]-PrimePi[n], {n, 94}] (* Stefano Spezia, Jun 05 2022 *)

Prog

(SageMath)
[A025528(n) - prime_pi(n)  for n in (1..100)] # Peter Luschny, Nov 18 2019

Crossrefs

Cf. A025528, A000720, A025475, A085501, A246655.

Partial sums of A268340.

Sequence in context: A324608 A237115 A362915 * A072292 A243282 A093390

Adjacent sequences: A069634 A069635 A069636 * A069638 A069639 A069640

Keyword

nonn

Author

Amarnath Murthy, Mar 27 2002

Status

approved