A269993 Denominators of r-Egyptian fraction expansion for sqrt(1/2), where r = (1,1/2,1/3,1/4,...)
2, 3, 9, 74, 8098, 101114070, 10080916639334518, 234737156891222571756748160861129, 104728182461244680288139397973895577148266725366426255244889745185
Offset: 1
Examples
sqrt(1/2) = 1/2 + 1/(2*3) + 1/(3*9) + ...
Links
- Clark Kimberling, Table of n, a(n) for n = 1..12
- Eric Weisstein's World of Mathematics, Egyptian Fraction
- Index entries for sequences related to Egyptian fractions
Programs
-
Mathematica
r[k_] := 1/k; f[x_, 0] = x; z = 10; n[x_, k_] := n[x, k] = Ceiling[r[k]/f[x, k - 1]] f[x_, k_] := f[x, k] = f[x, k - 1] - r[k]/n[x, k] x = Sqrt[1/2]; Table[n[x, k], {k, 1, z}]
-
PARI
r(k) = 1/k; x = sqrt(1/2); f(x, k) = if(k<1, x, f(x, k - 1) - r(k)/n(x, k)); n(x, k) = ceil(r(k)/f(x, k - 1)); for(k = 1, 10, print1(n(x, k),", ")) \\ Indranil Ghosh, Mar 27 2017, translated from Mathematica code
Comments