This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A120033 #20 Mar 05 2019 01:56:04 %S A120033 0,1,1,4,4,12,20,40,75,147,285,535,1062,2006,3918,7548,14595,28293, %T A120033 54761,106452,206421,401522,780966,1520543,2962226,5777162,11272279, %U A120033 22009839,43006972,84077384,164482781,321944211,630487562,1235382703 %N A120033 Number of semiprimes s such that 2^n < s <= 2^(n+1). %C A120033 The partial sum equals the number of Pi_2(2^n) = A125527(n). %H A120033 Dana Jacobsen, <a href="/A120033/b120033.txt">Table of n, a(n) for n = 0..62</a> (first 48 terms from Charles R Greathouse IV, corrected a(47)-a(48)) %e A120033 (2^2, 2^3] there is one semiprime, namely 6. 4 was counted in the previous entry. %t A120033 SemiPrimePi[n_] := Sum[PrimePi[n/Prime[i]] - i + 1, {i, PrimePi[Sqrt[n]]}]; t = Table[SemiPrimePi[2^n], {n, 0, 35}]; Rest@t - Most@t %o A120033 (PARI) pi2(n)=my(s,i); forprime(p=2, sqrt(n), s+=primepi(n\p); i++); s - i * (i-1)/2 %o A120033 a(n)=pi2(2^(n+1))-pi2(2^n) \\ _Charles R Greathouse IV_, May 15 2012 %o A120033 (Perl) use ntheory ":all"; print "$_ ",semiprime_count(1+(1<<$_), 1<<($_+1)),"\n" for 0..48; # _Dana Jacobsen_, Mar 04 2019 %o A120033 (Perl) use ntheory ":all"; my $l=0; for (0..48) { my $c=semiprime_count(1<<($_+1)); print "$_ ",$c-$l,"\n"; $l=$c; } # _Dana Jacobsen_, Mar 04 2019 %Y A120033 Cf. A001358, A072000, A066265, A036378, A120033-A120043. %K A120033 nonn %O A120033 0,4 %A A120033 _Jonathan Vos Post_ and _Robert G. Wilson v_, Mar 20 2006