cp's OEIS Frontend

A213358 a(n) = maximum(continuedfraction(F(n+1)^n/F(n)^n)) - L(n) + (1-(-1)^n)/2, where F(n) is Fibonacci(n) and L(n) is the n-th Lucas number.

%I A213358 #22 Mar 05 2017 00:57:06
%S A213358 1,1,0,4,7,16,0,0,397,17,0,8047,0,0,0,0,0,0,0,6633,0,0,0,0,0,0,0,0,0,
%T A213358 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
%U A213358 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
%N A213358 a(n) = maximum(continuedfraction(F(n+1)^n/F(n)^n)) - L(n) + (1-(-1)^n)/2, where F(n) is Fibonacci(n) and L(n) is the n-th Lucas number.
%C A213358 This is a variation on A113506 and A113500.
%C A213358 If we consider the sequence of quotients [1,1,1,1] of a simple continued fraction, we can express this by the continued fraction of F(4+1)/F(4), where we have chosen not to write the last quotient as 2 but as 1,1.
%C A213358 It is a very interesting fact that if we take the fifth power of F(n+1)/F(n), for varying n, we very soon observe extremely large values in its continued fraction expansion. That's why Benoit Cloitre, author of A113560, takes vecmax, the max of all these quotients, in order to capture the large rate of growth of these quotients.  It would be better, in order to emphasize the largeness of these big quotients, to somehow quantify the ratio between the max and the average of all the quotients, or the average of all the others.
%C A213358 Since Cloitre found a Fibonacci-like behavior of the vecmax, I thought to replace the power 5 by the power n, thus getting the sequence vecmax(contfrac(F(n+1)^n/F(n)^n)) instead of vecmax(contfrac(F(n+1)^5/F(n)^5)). After I noticed its behavior was almost that of the Lucas sequence, I took the sequence which subtracted the almost-Lucas sequence from it, thus eventually getting 0 after the 20th term.
%H A213358 G. C. Greubel, <a href="/A213358/b213358.txt">Table of n, a(n) for n = 1..1000</a>
%t A213358 Table[Max[ContinuedFraction[Fibonacci[n+1]^n/Fibonacci[n]^n]] - LucasL[n] + (1 - (-1)^n)/2, {n, 200}] (* _T. D. Noe_, Jun 28 2012 *)
%o A213358 (PARI)
%o A213358 F(n) = fibonacci(n);
%o A213358 L(n) = F(n-1) + F(n+1);
%o A213358 a(n) = vecmax(contfrac(F(n+1)^n/F(n)^n)) - L(n) + (1-(-1)^n)/2;
%Y A213358 Cf. A113506, A113500.
%K A213358 nonn
%O A213358 1,4
%A A213358 _Art DuPre_, Jun 10 2012