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.

A180048 Coefficient triangle of the denominators of the (n-th convergents to) the continued fraction 1/(w+2/(w+3/(w+4/... . Conjectured to equal unsigned version of A137286.

This page as a plain text file.
%I A180048 #44 Dec 25 2021 21:25:59
%S A180048 1,0,1,2,0,1,0,5,0,1,8,0,9,0,1,0,33,0,14,0,1,48,0,87,0,20,0,1,0,279,0,
%T A180048 185,0,27,0,1,384,0,975,0,345,0,35,0,1,0,2895,0,2640,0,588,0,44,0,1,
%U A180048 3840,0,12645,0,6090,0,938,0,54,0,1,0,35685,0,41685,0,12558,0,1422,0,65
%N A180048 Coefficient triangle of the denominators of the (n-th convergents to) the continued fraction 1/(w+2/(w+3/(w+4/... . Conjectured to equal unsigned version of A137286.
%C A180048 Equivalence to the recurrence formula needs formal proof. This continued fraction converges to 0.525135276160981... for w=1. A conjecture by Ramanujan puts this equal to -1 + 1/(sqrt(e*Pi/2) - Sum_{k>=1} 1/(2k-1)!!).
%C A180048 From _Alexander Kreinin_, Oct 26 2015: (Start)
%C A180048 Let us denote the continued fraction by U(w).
%C A180048 Then it is easy to show that Mill's ratio, R(w) = (1 - Phi(w))/f(w), where Phi is the standard normal distribution function and f is the standard normal density function, satisfies R(w) = 1/(w + U(w)).
%C A180048 Indeed, R(w) = 1/(w+1/(w+2/(w+3/(w+... Then we find U(w) = 1/R(w) - w. It was proved in Alexander Kreinin (arXiv:1405.5852) that R(w+t) + Q(w, t) = exp(w*t + w^2/2)*R(t), where Q(w,t) = Sum_{k>=0} Sum_{m=0..k} q(k,m) * t^m * w^(k+1)/(k+1)!.
%C A180048 Substituting t=0, we obtain R(w) = exp(w^2/2)*sqrt(Pi/2) - Sum_{n>=0} w^(2n+1)/(2n+1)!!. If w=1 we obtain Ramanujan's formula. (End)
%H A180048 G. C. Greubel, <a href="/A180048/b180048.txt">Table of n, a(n) for the first 75 rows, flattened</a>
%H A180048 Authors?, <a href="http://www.komal.hu/forum/forum.cgi?a=to&amp;tid=94&amp;st=50&amp;dr=1&amp;sp=1233">Hungarian discussion forum</a>
%H A180048 Alexander Kreinin, <a href="http://arxiv.org/abs/1405.5852">Combinatorial Properties of Mills' Ratio</a>, arXiv:1405.5852 [math.CO], 2014. See Table 3.
%H A180048 Alexander Kreinin, <a href="https://www.researchgate.net/profile/Alexander_Kreinin/publication/294260037_Integer_Sequences_and_Laplace_Continued_Fraction/links/56bf5f5308ae2f498ef7effa.pdf">Integer Sequences and Laplace Continued Fraction</a>, preprint, 2016.
%H A180048 Alexander Kreinin, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL19/Kreinin/kreinin4.html">Integer Sequences Connected to the Laplace Continued Fraction and Ramanujan's Identity</a>, Journal of Integer Sequences 19 (2016), Article 16.6.2.
%F A180048 p(0)=1; p(1)=w; p(n) = w*p(n-1) + n*p(n-2) (conjecture).
%F A180048 T(n,k) = T(n-1,k-1) + n*T(n-2,k), T(0,0) = 1, T(1,0) = 0, T(1,1) = 1. - _Philippe Deléham_, Oct 28 2013
%F A180048 sum_{k=0..n} T(n,k) = A000932(n). - _Philippe Deléham_, Oct 28 2013
%F A180048 T(2n,0) = A000165(n); T(2n+1,1) = A129890(n); T(2n+2,2) = A035101(n+2). - _Philippe Deléham_, Oct 28 2013
%e A180048 The denominator of 1/(w+2/(w+3/(w+4/(w+5/(w+6/w))))) equals 48 + 87w^2 + 20w^4 + w^6.
%e A180048 From _Joerg Arndt_, Apr 20 2013: (Start)
%e A180048 Triangle begins
%e A180048      1;
%e A180048      0,     1;
%e A180048      2,     0,     1;
%e A180048      0,     5,     0,     1;
%e A180048      8,     0,     9,     0,    1;
%e A180048      0,    33,     0,    14,    0,   1;
%e A180048     48,     0,    87,     0,   20,   0,   1;
%e A180048      0,   279,     0,   185,    0,  27,   0,  1;
%e A180048    384,     0,   975,     0,  345,   0,  35,  0,  1;
%e A180048      0,  2895,     0,  2640,    0, 588,   0, 44,  0, 1;
%e A180048   3840,     0, 12645,     0, 6090,   0, 938,  0, 54, 0, 1;
%e A180048      0, 35685,     0, 41685,    0, ... (End)
%t A180048 Table[ CoefficientList[ Denominator[ Together[ Fold[ #2/(w+#1) &, Infinity, Reverse @ Table[ k, {k, 1, n} ] ] ] ], w ], {n, 16} ] (* or equivalently *) Clear[ p ];p[ 0 ]=1; p[ 1 ]=w; p[ n_ ]:=p[ n ]= w*p[ n-1 ] + n*p[ n-2 ]; Table[ CoefficientList[ p[ k ]//Expand, w ], {k,0,15} ]
%Y A180048 Cf. A137286, A084950, A180047, A180049.
%Y A180048 Cf. A111129.
%K A180048 nonn,tabl
%O A180048 0,4
%A A180048 _Wouter Meeussen_, Aug 08 2010