cp's OEIS Frontend

A247698 Brady numbers: B(n) = B(n - 1) + B(n - 2) with B(1) = 2308 and B(2) = 4261.

%I A247698 #67 Jul 08 2022 08:21:18
%S A247698 2308,4261,6569,10830,17399,28229,45628,73857,119485,193342,312827,
%T A247698 506169,818996,1325165,2144161,3469326,5613487,9082813,14696300,
%U A247698 23779113,38475413,62254526,100729939,162984465,263714404,426698869,690413273,1117112142,1807525415,2924637557,4732162972,7656800529
%N A247698 Brady numbers: B(n) = B(n - 1) + B(n - 2) with B(1) = 2308 and  B(2) = 4261.
%C A247698 B(n) / B(n - 1) approaches the golden ratio as n approaches infinity.
%H A247698 Logan Cooper, <a href="/A247698/b247698.txt">Table of n, a(n) for n = 1..1000</a> (truncated from 9966 to 1000 terms by _M. F. Hasler_, May 10 2017)
%H A247698 Brady Haran and Matt Parker, <a href="https://www.youtube.com/watch?v=D8ntDpBm6Ok">Brady Numbers</a>, Numberphile video (2014).
%H A247698 <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (1,1).
%F A247698 a(n) = a(n-1) + a(n-2).
%F A247698 G.f.: x*(2308 + 1953*x) / (1-x-x^2). - _Colin Barker_, Sep 23 2014
%F A247698 a(n) = k*phi^n + o(1), where k = 976.5 + sqrt(354578.45) = 1571.96.... - _Charles R Greathouse IV_, Sep 28 2014
%F A247698 a(n) = 2308*A000045(n-2) + 4261*A000045(n-1) = 1953*A000045(n+1) + 355*A000045(n). - _M. F. Hasler_, May 10 2017
%F A247698 a(n) = F(n+17) - F(n+8) - 9*F(n) - F(n-14) for F(n) = A000045(n). - _Greg Dresden_, Jul 07 2022
%p A247698 Brady1 := proc(n::posint)
%p A247698 option remember, system;
%p A247698 if n = 1 then
%p A247698   2308
%p A247698 elif n = 2 then
%p A247698   4261
%p A247698 else
%p A247698   thisproc( n - 1 ) + thisproc( n - 2 )
%p A247698 end if
%p A247698 end proc:
%p A247698 seq( Brady1( n ), n = 1 .. 100 );
%p A247698 # _James McCarron_, Oct 05 2019
%p A247698 # alternate program
%p A247698 Brady2 := ( n :: posint ) -> coeff( series(x*(2308+1953*x)/(1-x-x^2),x,n+1), x^n ):
%p A247698 seq( Brady2( n ), n = 1 .. 100 );
%p A247698 # _James McCarron_, Oct 05 2019
%t A247698 LinearRecurrence[{1, 1}, {2308, 4261}, n]
%t A247698 Rest[CoefficientList[Series[x*(2308+1953*x)/(1-x-x^2), {x,0,50}], x]] (* _G. C. Greubel_, Sep 07 2018 *)
%o A247698 (Haskell) brady = let makeSeq a b = a : makeSeq b (a + b) in makeSeq 2308 4261
%o A247698 (PARI) Vec(-x*(1953*x+2308)/(x^2+x-1) + O(x^50)) \\ _Colin Barker_, Sep 23 2014
%o A247698 (PARI) a(n)=([1,1;1,0]^n*[1953;355])[1,1] \\ _Charles R Greathouse IV_, Jan 20 2016
%o A247698 (Magma) m:=50; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(x*(2308+1953*x)/(1-x-x^2))); // _G. C. Greubel_, Sep 07 2018
%o A247698 (Python)
%o A247698 def A247698_list(n):
%o A247698     list = [2308, 4261] + [0] * (n - 2)
%o A247698     for i in range(2, n):
%o A247698         list[i] = list[i - 1] + list[i - 2]
%o A247698     return list
%o A247698 print(A247698_list(32)) # _M. Eren Kesim_, Jun 28 2021
%K A247698 nonn,easy
%O A247698 1,1
%A A247698 _Sebastian Zimmer_, Sep 22 2014