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.

A163521 a(n) = floor(Sum_{k = 2..10^n} k/log(k)).

This page as a plain text file.
%I A163521 #16 Mar 28 2021 19:32:40
%S A163521 30,1255,78698,5762750,455059956,37607986470,3204942375900,
%T A163521 279238346962895,24739954333817884
%N A163521 a(n) = floor(Sum_{k = 2..10^n} k/log(k)).
%C A163521 a(n) = Sum_{x=2..n} x/log(x) closely approximates the number of primes < n^2.
%C A163521 In fact, the sum is as good as Li(n^2) but summing a(n) is rather time consuming for large n.
%C A163521 For n = 10^9,
%C A163521   a(n)    = 24739954333817884,
%C A163521   Pi(n^2) = 24739954287740860,
%C A163521   Li(n^2) = 24739954309690415,
%C A163521   R(n^2)  = 24739954284239494,
%C A163521 where Li = Logarithmic integral approximation of Pi, and R = Riemann's approximation of Pi.
%C A163521 Now x/(log(x)-1) is a much better approximation of Pi(x) than x/log(x):
%C A163521   10^18/(log(10^18)-1) = 24723998785919976,
%C A163521   10^18/log(10^18)     = 24127471216847323.
%C A163521 Ironically, though, a(n) = Sum_{x=2..n} x/(log(x)-1) is far from Pi(n^2).
%e A163521 For n = 9, floor(Sum_{x=2..10^n} x/log(x)) = 24739954333817884, the 9th term.
%t A163521 Table[Floor[Sum[j/Log[j], {j, 2, 10^n}]], {n, 1, 9}] (* _G. C. Greubel_, Jul 27 2017 *)
%o A163521 (PARI) nthsum(n) = for(j=1,n,print1(floor(sum(x=2,10^j,x/log(x)))","));
%K A163521 nonn
%O A163521 1,1
%A A163521 _Cino Hilliard_, Jul 30 2009
%E A163521 Definition clarified by _R. J. Mathar_ and _Omar E. Pol_, Aug 01 2009