cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A364872 Lexicographically earliest continued fraction which is its own unit fraction series.

Original entry on oeis.org

2, 1, 2, 2, 5, 95, 137447, 19092121105, 1456654254113777258001, 8728918703159963392919895262580124849062181
Offset: 0

Views

Author

Rok Cestnik, Aug 11 2023

Keywords

Comments

Earliest infinite sequence {a0,a1,a2,a3,...} such that: a0+1/(a1+1/(a2+1/(a3+...))) = 1/a0 + 1/a1 + 1/a2 + 1/a3 + ...
There are infinitely many real numbers whose continued fraction is also their unit fraction series - they are dense on the interval (2,oo).

Examples

			The partial continued fraction must always be strictly larger than the partial unit fractions:
[1]      cannot be  since            1 = 1.
[2]      can be     since            2 > 1/2.
[2,1]    can be     since        2+1/1 > 1/2+1/1.
[2,1,1]  cannot be  since  2+1/(1+1/1) = 1/2+1/1+1/1.
[2,1,2]  can be     since  2+1/(1+1/2) > 1/2+1/1+1/2.
...
sum(1/a[n]) = 2.71053359137351078733864566... (A364873).
		

Crossrefs

Cf. A364873.

Programs

  • PARI
    cf(a) = my(m=contfracpnqn(a)); m[1,1]/m[2,1];
    uf(a) = sum(i=1, #a, 1/a[i]);
    A364872(N) = {a=[2]; for(i=2, N, a=concat(a, if(cf(a)==uf(a), a[i-1], ceil(1/(cf(a)-uf(a))))); while(cf(a)<=uf(a), a[i]++)); a};