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.

This page as a plain text file.
%I A364872 #17 Aug 13 2023 03:33:30
%S A364872 2,1,2,2,5,95,137447,19092121105,1456654254113777258001,
%T A364872 8728918703159963392919895262580124849062181
%N A364872 Lexicographically earliest continued fraction which is its own unit fraction series.
%C A364872 Earliest infinite sequence {a0,a1,a2,a3,...} such that: a0+1/(a1+1/(a2+1/(a3+...))) = 1/a0 + 1/a1 + 1/a2 + 1/a3 + ...
%C A364872 There are infinitely many real numbers whose continued fraction is also their unit fraction series - they are dense on the interval (2,oo).
%e A364872 The partial continued fraction must always be strictly larger than the partial unit fractions:
%e A364872 [1]      cannot be  since            1 = 1.
%e A364872 [2]      can be     since            2 > 1/2.
%e A364872 [2,1]    can be     since        2+1/1 > 1/2+1/1.
%e A364872 [2,1,1]  cannot be  since  2+1/(1+1/1) = 1/2+1/1+1/1.
%e A364872 [2,1,2]  can be     since  2+1/(1+1/2) > 1/2+1/1+1/2.
%e A364872 ...
%e A364872 sum(1/a[n]) = 2.71053359137351078733864566... (A364873).
%o A364872 (PARI)
%o A364872 cf(a) = my(m=contfracpnqn(a)); m[1,1]/m[2,1];
%o A364872 uf(a) = sum(i=1, #a, 1/a[i]);
%o A364872 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};
%Y A364872 Cf. A364873.
%K A364872 nonn,cofr
%O A364872 0,1
%A A364872 _Rok Cestnik_, Aug 11 2023