login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 56th year, we are closing in on 350,000 sequences, and we’ve crossed 9,700 citations (which often say “discovered thanks to the OEIS”).

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A270523 Denominators of r-Egyptian fraction expansion for Pi - 3, where r(k) = 1/k!. 1
8, 31, 360, 63288, 3000329177, 2267607071582813683, 7548646359131509583693406626221228733, 22436552662647350051378366551573442407224062622229053640998338187956658409 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

Suppose that r is a sequence of rational numbers r(k) <= 1 for k >= 1, and that x is an irrational number in (0,1).  Let f(0) = x, n(k) = floor(r(k)/f(k-1)), and f(k) = f(k-1) - r(k)/n(k).  Then x = r(1)/n(1) + r(2)/n(2) + r(3)/n(3) + ... , the r-Egyptian fraction for x.

See A269993 for a guide to related sequences.

LINKS

Clark Kimberling, Table of n, a(n) for n = 1..11

Eric Weisstein's World of Mathematics, Egyptian Fraction

Index entries for sequences related to Egyptian fractions

EXAMPLE

Pi - 3 = 1/(1*8) + 1/(2*31) + 1/(6*360) + 1/(24*63288) + ...

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 = Pi - 3; Table[n[x, k], {k, 1, z}]

PROG

(PARI) r(k) = 1/k!;

f(k, x) = if (k==0, x, f(k-1, x) - r(k)/a(k, x); );

a(k, x=Pi-3) = ceil(r(k)/f(k-1, x)); \\ Michel Marcus, Mar 31 2016

CROSSREFS

Cf. A269993, A000142, A000796.

Sequence in context: A303176 A209484 A209343 * A332916 A270353 A270400

Adjacent sequences:  A270520 A270521 A270522 * A270524 A270525 A270526

KEYWORD

nonn,frac,easy

AUTHOR

Clark Kimberling, Mar 30 2016

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified November 30 17:12 EST 2021. Contains 349424 sequences. (Running on oeis4.)