OFFSET
0,4
FORMULA
EXAMPLE
a(0..2) = 0 since P(0..2) = {1, 2, 6}, and the smallest number in s is 12.
a(3) = 5 since P(3) = 30, and the set s(1..6) = {12, 18, 20, 24, 28} contains k <= 30.
a(4) = 67 since P(4) = 210, and the set s(1..67) = {12, 18, 20, ..., 207, 208} contains k <= 210, etc.
MATHEMATICA
Table[# - Sum[PrimePi@ Floor[#^(1/k)], {k, 2, Floor[Log2[#]]}] - Sum[MoebiusMu[k]*Floor[#/(k^2)], {k, Floor[Sqrt[#]]}] &[Product[Prime[i], {i, n}]], {n, 0, 12}]
PROG
(Python)
from math import isqrt
from sympy import primorial, primepi, integer_nthroot, mobius
def A380403(n):
if n == 0: return 0
m = primorial(n)
return int(-sum(primepi(integer_nthroot(m, k)[0]) for k in range(2, m.bit_length()))-sum(mobius(k)*(m//k**2) for k in range(2, isqrt(m)+1))) # Chai Wah Wu, Jan 24 2025
CROSSREFS
KEYWORD
nonn,hard,more,new
AUTHOR
Michael De Vlieger, Jan 23 2025
EXTENSIONS
Offset changed to 0 by Jinyuan Wang, Jan 24 2025
a(16) from Chai Wah Wu, Jan 24 2025
a(17) from Chai Wah Wu, Jan 25 2025
STATUS
approved