OFFSET
1,3
COMMENTS
a(p) = pi(p) for p prime.
FORMULA
a(n) = n - Sum_{k=1..n} [u(n/k) = v(k)], where u(n) = 1 - ceiling(n) + floor(n), v is the prime characteristic (A010051), and [ ] is the Iverson bracket.
EXAMPLE
a(8) = 6; 1 is nonprime and divides 8, 3 does not divide 8 and is prime, 4 is not prime and divides 8, 5 is prime and does not divide 8, 7 is prime and does not divide 8, and 8 (nonprime) divides itself. So a(8) = 6.
PROG
(PARI) a(n) = sum(k=1, n, bitxor(isprime(k)&&(n%k), !isprime(k)&&!(n%k))); \\ Michel Marcus, Jan 31 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Jan 30 2022
STATUS
approved