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

%I #15 Nov 27 2020 04:51:42

%S 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,

%T 14,14,15,15,15,15,15,15,16,16,16,16,17,17,18,18,18,18,19,19,20,20,20,

%U 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

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

%H Antti Karttunen, <a href="/A302778/b302778.txt">Table of n, a(n) for n = 1..65537</a>

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

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

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

%o (PARI)

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

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

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

%Y Partial sums of A302777. A left inverse of A050376.

%Y Cf. A302785, A302786.

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

%Y Cf. also A000720, A025528.

%K nonn

%O 1,3

%A _Antti Karttunen_, Apr 16 2018