|
|
A247698
|
|
Brady numbers: B(n) = B(n - 1) + B(n - 2) with B(1) = 2308 and B(2) = 4261.
|
|
3
|
|
|
2308, 4261, 6569, 10830, 17399, 28229, 45628, 73857, 119485, 193342, 312827, 506169, 818996, 1325165, 2144161, 3469326, 5613487, 9082813, 14696300, 23779113, 38475413, 62254526, 100729939, 162984465, 263714404, 426698869, 690413273, 1117112142, 1807525415, 2924637557, 4732162972, 7656800529
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
B(n) / B(n - 1) approaches the golden ratio as n approaches infinity.
|
|
LINKS
|
Brady Haran and Matt Parker, Brady Numbers, Numberphile video (2014).
|
|
FORMULA
|
a(n) = a(n-1) + a(n-2).
G.f.: x*(2308 + 1953*x) / (1-x-x^2). - Colin Barker, Sep 23 2014
|
|
MAPLE
|
Brady1 := proc(n::posint)
option remember, system;
if n = 1 then
2308
elif n = 2 then
4261
else
thisproc( n - 1 ) + thisproc( n - 2 )
end if
end proc:
seq( Brady1( n ), n = 1 .. 100 );
# alternate program
Brady2 := ( n :: posint ) -> coeff( series(x*(2308+1953*x)/(1-x-x^2), x, n+1), x^n ):
seq( Brady2( n ), n = 1 .. 100 );
|
|
MATHEMATICA
|
LinearRecurrence[{1, 1}, {2308, 4261}, n]
Rest[CoefficientList[Series[x*(2308+1953*x)/(1-x-x^2), {x, 0, 50}], x]] (* G. C. Greubel, Sep 07 2018 *)
|
|
PROG
|
(Haskell) brady = let makeSeq a b = a : makeSeq b (a + b) in makeSeq 2308 4261
(PARI) Vec(-x*(1953*x+2308)/(x^2+x-1) + O(x^50)) \\ Colin Barker, Sep 23 2014
(Magma) m:=50; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(x*(2308+1953*x)/(1-x-x^2))); // G. C. Greubel, Sep 07 2018
(Python)
list = [2308, 4261] + [0] * (n - 2)
for i in range(2, n):
list[i] = list[i - 1] + list[i - 2]
return list
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|