cp's OEIS Frontend

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.

A292785 Number of odd squarefree semiprimes < 10^n.

This page as a plain text file.
%I A292785 #22 Oct 17 2024 15:47:19
%S A292785 0,16,194,1932,18181,168330,1555366,14424896,134429269,1258812629,
%T A292785 11840308472,111817802539,1059796159358,10076977878935,96091981692305,
%U A292785 918679869869451,8803388128870716,84537081067757934,813340036390023775,7838825925395981969,75669246174605279757
%N A292785 Number of odd squarefree semiprimes < 10^n.
%F A292785 a(n) = A066265(n) - A122121(n) - A220262(n) + 1 for n > 1.
%e A292785 a(2)=16 because there are 16 squarefree odd semiprimes < 10^2: 15=3*5, 21=3*7, 33=3*11, 35=5*7, 39=3*13, 51=3*17, 55=5*11, 57=3*19, 65=5*13, 69=3*23, 77=7*11, 85=5*17, 87=3*29, 91=7*13, 93=3*31, 95=5*19.
%t A292785 -1 + Accumulate@ Table[Count[Range[10^n + 1, 10^(n + 1) - 1, 2], _?(And[ SquareFreeQ@ #, PrimeNu[#] == 2] &)], {n, 0, 5}] (* _Michael De Vlieger_, Oct 10 2017 *)
%o A292785 (Python)
%o A292785 from math import isqrt
%o A292785 from sympy import primepi, primerange
%o A292785 def A292785(n): return int((-(t:=primepi(s:=isqrt(m:=10**n)))*(t+1)>>1)+1+sum(primepi(m//k) for k in primerange(3, s+1))) if n>1 else 0 # _Chai Wah Wu_, Oct 17 2024
%Y A292785 Cf. A046388, A066265, A122121, A220262.
%K A292785 nonn
%O A292785 1,2
%A A292785 _Hugo Pfoertner_, Oct 10 2017
%E A292785 a(21) from _Jinyuan Wang_, Jul 30 2021