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.

A128910 Similar to A057835 except using K * X / log(X), K=1.022.

This page as a plain text file.
%I A128910 #20 Aug 18 2022 21:11:29
%S A128910 0,3,20,119,715,4523,30509,213343,1530983,11203550,83064263,620498643,
%T A128910 4643259527,34592032908,254639722327,1832740718223,12680919388801,
%U A128910 81678704122892,452951221016511,1574800035301944,8395299939524712,282240813012897282,4457697545906326118,58106920364272792945,693274802905577732102,7864635685729658131835
%N A128910 Similar to A057835 except using K * X / log(X), K=1.022.
%C A128910 This is an improvement over the classic X / log(X) approximation in the range many people work with.
%C A128910 pi(x), R(x), and li(x) are all asymptotically x/log x + x/log^2 x + O(x/log^3 x), so this approximation is good around exp(1/.022) ≈ 5 * 10^19. Asymptotically the best value for K would be 1. - _Charles R Greathouse IV_, Aug 18 2022
%F A128910 a(n) = abs(round(1.022*10^n/log(10^n)) - pi(10^n)). - _Charles R Greathouse IV_, Mar 22 2015
%F A128910 a(n) ~ 10^n/kn with k = 104.6629.... - _Charles R Greathouse IV_, Mar 22 2015
%e A128910 a(10)=11203550 via abs (455,052,511 - 443,848,961).
%t A128910 Table[ PrimePi[10^n] - Round[N[1.022*10^n/Log[10^n]]], {n, 23}] (* and absolute value thereof (orig entries 21-23 <0); courtesy of _Robert G. Wilson v_ *)
%o A128910 (PARI) a(n) = abs(round(1.022*10^n/log(10^n)) - primepi(10^n)) \\ _Charles R Greathouse IV_, Mar 22 2015
%Y A128910 Cf. A057835, A006880.
%K A128910 nonn
%O A128910 1,2
%A A128910 _Bill McEachen_, Apr 23 2007