login
Numerator of best approximation to Pi with denominator <= 10^n.
8

%I #31 Jul 17 2023 06:00:55

%S 3,22,22,355,355,312689,1146408,5419351,245850922,2549491779,

%T 21053343141,21053343141,1783366216531,8958937768937,139755218526789,

%U 428224593349304,30246273033735921,66627445592888887

%N Numerator of best approximation to Pi with denominator <= 10^n.

%H Daniel Mondot, <a href="/A072398/b072398.txt">Table of n, a(n) for n = 0..999</a>

%e A072398(5) = 312689 because A072398(5)/A072399(5) = 312689/99532 is the best rational approximation to Pi with positive denominator <= 10^5 = 100000. This approximation is accurate to 0.00000000092766%.

%t nmax = 17; cv = Convergents[Pi, 2*nmax] // Reverse; a[n_] := Select[cv, Denominator[#] <= 10^n &, 1] // Numerator // First; Table[a[n], {n, 0, nmax}] (* _Jean-François Alcover_, Jan 04 2013 *)

%o (PARI) for(n=0,40,print1(numerator(bestappr(Pi,10^n)),",")) \\ Finds these approximations very quickly.

%Y Cf. A072399 (denominators), A000796 (Pi), A068089, A002485/A002486.

%K easy,frac,nice,nonn

%O 0,1

%A _Rick L. Shepherd_, Jun 15 2002