cp's OEIS Frontend

A022031 Define the sequence T(a(0),a(1)) by a(n+2) is the greatest integer such that a(n+2)/a(n+1) < a(n+1)/a(n) for n >= 0. This is T(4,17).

%I A022031 #18 Feb 16 2016 06:07:38
%S A022031 4,17,72,304,1283,5414,22845,96397,406757,1716352,7242319,30559689,
%T A022031 128949662,544115986,2295951781,9687997993,40879475731,172495033261,
%U A022031 727860031657,3071278144467,12959565068034,54684179957837,230745362360740,973653116715681,4108426630946045
%N A022031 Define the sequence T(a(0),a(1)) by a(n+2) is the greatest integer such that a(n+2)/a(n+1) < a(n+1)/a(n) for n >= 0. This is T(4,17).
%C A022031 The empirical g.f. / recurrence agrees with the original definition for at least 2000 terms (and a(2000) ~ 10^1250). - _M. F. Hasler_, Feb 11 2016
%H A022031 Colin Barker, <a href="/A022031/b022031.txt">Table of n, a(n) for n = 0..1000</a>
%F A022031 Empirical g.f.: -(x^6+x^5+x^4+x^3-x-4) / ((x-1)*(x^6+2*x^5+3*x^4+4*x^3+4*x^2+3*x-1)). - _Colin Barker_, Sep 18 2015
%F A022031 a(n+1) = ceiling(a(n)^2/a(n-1))-1 for all n > 0. a(n+1)/a(n) ~ 4.219599938... as n -> oo. - _M. F. Hasler_, Feb 11 2016
%o A022031 (PARI) a=[4,17];for(n=2,2000,a=concat(a,ceil(a[n]^2/a[n-1])-1));A022031(n)=a[n+1] \\ _M. F. Hasler_, Feb 11 2016
%Y A022031 Cf. A022018 - A022025, A022026 - A022032.
%K A022031 nonn
%O A022031 0,1
%A A022031 _R. K. Guy_
%E A022031 Edited by _M. F. Hasler_, Feb 11 2016