%I #18 Aug 20 2019 01:13:28
%S 1,2,3,6,7,13,21,112,243,275,466,761,1128,4704,9523,10730,17579,28085,
%T 41041,165312,331299,372815,607754,967441,1410360,5648160,11300259,
%U 12713402,20707831,32942845,48005301,192060400,384143763,432165299,703818922,1119543881,1631318640
%N a(n) is the smallest positive integer x such that sqrt(2) + sqrt(x) is closer to an integer than any other value already in the sequence.
%C If b(n) = round(sqrt(2) + sqrt(a(n))), then (b(n)^2 + 2 - a(n))/(2*b(n)) is an approximation for sqrt(2). Conjecture: all convergents of the continued fraction of sqrt(2) except 1 arise in this way. - _Robert Israel_, Aug 18 2019
%e a(6) = 13 because sqrt(2)+sqrt(13) is closer to an integer than any of the previous 5 terms.
%p R:= 1: delta:= sqrt(2)-1:
%p for r from 2 to 10000 do
%p x0:= ceil((r - sqrt(2)-delta)^2);
%p x1:= floor((r-sqrt(2)+delta)^2);
%p for x from x0 to x1 do
%p dx:= abs(sqrt(2)+sqrt(x)-r);
%p if is(dx < delta) then
%p delta:= dx;
%p R:= R, x;
%p fi
%p od
%p od:
%p R; # _Robert Israel_, Aug 18 2019
%t d[x_] := Abs[x - Round[x]]; dm = 1; s = {}; Do[If[(d1 = d[Sqrt[2] + Sqrt[n]]) < dm, dm = d1; AppendTo[s, n]], {n, 1, 10^5}]; s (* _Amiram Eldar_, Aug 18 2019 *)
%o (Python) import math
%o a = 2**(1/2)
%o l = []
%o closest = 1.0
%o for i in range(1, 100000000):
%o b = i**(1/2)
%o c = abs(a+b - round(a+b))
%o if c < closest:
%o print(i, c)
%o closest = c
%o l.append(i)
%o print(l)
%K nonn
%O 1,2
%A _Ben Paul Thurston_, Aug 18 2019
%E More terms from _Giovanni Resta_, Aug 19 2019
|