A302778 - OEIS

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

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,3

LINKS

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

FORMULA

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

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

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 *)

PROG

(PARI)

A209229(n) = (n && !bitand(n, n-1));

A302777(n) = A209229(isprimepower(n));

s=0; for(n=1, 105, s+=A302777(n); print1(s, ", "));

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.

Sequence in context: A047741 A047784 A047742 * A203967 A365398 A365399

Adjacent sequences: A302775 A302776 A302777 * A302779 A302780 A302781

KEYWORD

nonn

AUTHOR

Antti Karttunen, Apr 16 2018

STATUS

approved