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.

A056811 Number of primes not exceeding square root of n: primepi(sqrt(n)).

This page as a plain text file.
%I A056811 #49 Aug 15 2024 06:38:39
%S A056811 0,0,0,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,
%T A056811 3,3,3,3,3,3,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,
%U A056811 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4
%N A056811 Number of primes not exceeding square root of n: primepi(sqrt(n)).
%C A056811 Number of primes among factors of LCM(1,...,n) whose exponent is > 1, i.e., number of non-unitary prime factors of LCM(1,...,n).
%C A056811 Number of positive integers <= n with exactly 3 divisors.
%C A056811 Number of squared primes not exceeding n. - _Wesley Ivan Hurt_, May 24 2013
%C A056811 Maximum number of composite numbers not exceeding n that are all coprime to each other. - _Yifan Xie_, Jul 07 2024
%H A056811 T. D. Noe, <a href="/A056811/b056811.txt">Table of n, a(n) for n = 1..10000</a>
%F A056811 a(n) = A056170(A003418(n)) = A000720(A000196(n)).
%F A056811 For k = 1, 2, ..., repeat k A069482(k) (that is, prime(k+1)^2 - prime(k)^2) times, and add 0 three times at the beginning (or begin the preceding by k = 0, with prime(0) set to 1). - _Jean-Christophe Hervé_, Oct 30 2013
%F A056811 G.f.: (1/(1 - x)) * Sum_{k>=1} x^(prime(k)^2). - _Ilya Gutkovskiy_, Sep 14 2019
%F A056811 a(n) ~ 2*n^(1/2)/log(n), by the prime number theorem. - _Harry Richman_, Jan 19 2022
%e A056811 If n=169,...,288 = p()^2,...,p(7)^2-1, then only the first 6 primes have exponents larger than 1, resulting in powers: 128, 81, 125, 49, 121, 169. So a(n)=6 for as much as 288-169+1 = 120 values of n.
%t A056811 Table[PrimePi[Sqrt[n]], {n, 100}] (* _T. D. Noe_, Mar 13 2013 *)
%o A056811 (PARI) a(n) = primepi(sqrt(n)); \\ _Michel Marcus_, Apr 11 2016
%o A056811 (Python)
%o A056811 from math import isqrt
%o A056811 from sympy import primepi
%o A056811 def a(n): return primepi(isqrt(n))
%o A056811 print([a(n) for n in range(1, 88)]) # _Michael S. Branicky_, Jan 19 2022
%Y A056811 Cf. A000196, A000720, A003418, A056170.
%Y A056811 Cf. A069482, A056813.
%K A056811 nonn,easy
%O A056811 1,9
%A A056811 _Labos Elemer_, Aug 28 2000