OFFSET
1,4
FORMULA
a(n) = prime(n) - A027883(n). - Chai Wah Wu, Dec 08 2024
EXAMPLE
The non prime powers counted under each term:
n=1 n=2 n=3 n=4 n=5 n=6 n=7 n=8 n=9 n=10
-------------------------------------------------
1 1 1 6 10 12 15 18 22 28
1 6 10 14 15 21 26
1 6 12 14 20 24
1 10 12 18 22
6 10 15 21
1 6 14 20
1 12 18
10 15
6 14
1 12
10
6
1
MATHEMATICA
Table[Length[Select[Range[Prime[n]], Not@*PrimePowerQ]], {n, 100}]
PROG
(Python)
from sympy import prime, primepi, integer_nthroot
def A378615(n): return int((p:=prime(n))-n-sum(primepi(integer_nthroot(p, k)[0]) for k in range(2, p.bit_length()))) # Chai Wah Wu, Dec 07 2024
CROSSREFS
KEYWORD
nonn,new
AUTHOR
Gus Wiseman, Dec 06 2024
STATUS
approved