A302778 Number of "Fermi-Dirac primes" (A050376) <= n.

0, 1, 2, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8, 8, 8, 9, 10, 10, 11, 11, 11, 11, 12, 12, 13, 13, 13, 13, 14, 14, 15, 15, 15, 15, 15, 15, 16, 16, 16, 16, 17, 17, 18, 18, 18, 18, 19, 19, 20, 20, 20, 20, 21, 21, 21, 21, 21, 21, 22, 22, 23, 23, 23, 23, 23, 23, 24, 24, 24, 24, 25, 25, 26, 26, 26, 26, 26, 26, 27, 27, 28, 28, 29, 29, 29, 29, 29, 29, 30

Offset: 1

Author

Antti Karttunen, Apr 16 2018

Keywords

• nonn

Links

• Antti Karttunen, Table of n, a(n) for n = 1..65537

Crossrefs

Partial sums of A302777. A left inverse of A050376.

Cf. A302785, A302786.

Differs from A203967 for the first time at n=64, where a(64) = 23, while A203967(64) = 24.

Cf. also A000720, A025528.

Programs

Mathematica

s[n_] := Boole[n > 1 && Length[(f = FactorInteger[n])] == 1 && (e = f[[;; , 2]]) == 2^IntegerExponent[e, 2]]; Accumulate @ Array[s, 100] (* Amiram Eldar, Nov 27 2020 *)

PARI

A209229(n) = (n && !bitand(n,n-1));
A302777(n) = A209229(isprimepower(n));
s=0; for(n=1,105,s+=A302777(n); print1(s,", "));

Python

from sympy import primepi, integer_nthroot
def A302778(n): return sum(primepi(integer_nthroot(n,1<Chai Wah Wu, Feb 18-19 2025

Formula

a(1) = 0; for n > 1, a(n) = A302777(n) + a(n-1).

For all n >= 1, a(A050376(n)) = n.