|
%I #23 Apr 15 2019 12:38:16
%S 1,4,16,67,283,1198,5074,21493,91045,385672,1633732,6920599,29316127,
%T 124185106,526056550,2228411305,9439701769,39987218380,169388575288,
%U 717541519531,3039554653411,12875760133174,54542595186106,231046140877597
%N a(n) = floor(a(n-1)/(sqrt(5) - 2)) for n > 0 and a(0) = 1.
%H Clark Kimberling, <a href="/A024551/b024551.txt">Table of n, a(n) for n = 0..250</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (5,-3,-1).
%F a(n) = 5*a(n-1) - 3*a(n-2) - a(n-3). - _Clark Kimberling_, Aug 16 2012
%F G.f.: (-x^2-x+1)/[(1-x)(1-4x-x^2)].
%F a(n) = (3*Fibonacci(3*n+2) + 1)/4 = 1 + 3*Sum_{k=0..n} A001076(k). - _Ehren Metcalfe_, Apr 15 2019
%t a[0] = 1;
%t a[n_] := Floor[a[n - 1]/FractionalPart[Sqrt[5]]]
%t Table[a[n], {n, 0, 60}]
%t (* _Clark Kimberling_, Aug 16 2012 *)
%t a[0]=1;
%t a[1]=4;
%t a[2]=16;
%t a[n_]:=Floor[a[n-1]^2/a[n-2]]+3
%t Table[a[n],{n,0,60}]
%t With[{c=Sqrt[5]-2},NestList[Floor[#/c]&,1,30]] (* _Harvey P. Dale_, Jul 18 2018 *)
%o (PARI) a(n)=([0,1,0; 0,0,1; -1,-3,5]^n*[1;4;16])[1,1] \\ _Charles R Greathouse IV_, Jan 20 2017
%o (PARI) step(n)=2*n + sqrtint(5*n^2)
%o a(n)=if(n, step(a(n-1)), 1) \\ _Charles R Greathouse IV_, Jan 20 2017
%K nonn,easy
%O 0,2
%A _Clark Kimberling_
|
