login
Denominators of signed Egyptian fractions 1/(1+2*a(n)) with sums converging to sqrt(2).
0

%I #4 May 30 2013 16:58:28

%S 17,38,2558,53508058,10183965708276283,

%T 833167602683818992272386593114136,

%U 7008824222617646742567474710166940582408242437487752499090108838014

%N Denominators of signed Egyptian fractions 1/(1+2*a(n)) with sums converging to sqrt(2).

%C See A226049.

%e The algorithm at A226049, with r = sqrt(2), f(n) = 1/(2n+1), gives

%e 1/3 + 1/5 + 1/7 + ... + 1/(1+2*17) - 1/(1+2*38) + 1/(1+2*2558) - ...

%e which converges to sqrt(2). The 17th partial sum differs from sqrt(2) by less than 10^(-500).

%t $MaxExtraPrecision = Infinity; z = 9; f[n_] := 1/(2 n + 1); g[n_] := (1/n - 1)/2; r =

%t Sqrt[2]; s = 0; a[1] = NestWhile[# + 1 &, 1, ! (s += f[#]) > r &]; p = Sum[f[n], {n, 1, a[1]}]; a[2] = Floor[g[p - r]]; a[n_] := Floor[g[((-1)^n) (p - r - Sum[((-1)^k) f[a[k]], {k, 2, n - 1}])]]; Table[a[k], {k, 1, z}]

%Y Cf. A226049.

%K nonn

%O 1,1

%A _Clark Kimberling_, May 27 2013