|
%I #15 Oct 21 2022 22:10:07
%S 357,14450,221650,3508050,55975250,894989650,14317376850,229068199250,
%T 3665051866450,58640672576850,938250132084050,15011999596762450,
%U 240191983481869650,3843071695444596050,61489146966052263250
%N a(n) = 20*a(n-1) - 64*a(n-2) - 150 for n > 2; a(0) = 357, a(1) = 14450, a(2) = 221650.
%C Related to Reverse and Add trajectory of 318 in base 4: A075153(6*n+1) = 3*a(n).
%C lim_{n -> infinity} a(n)/a(n-1) = 16.
%H G. C. Greubel, <a href="/A166913/b166913.txt">Table of n, a(n) for n = 0..500</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (21, -84, 64).
%F a(n) = (2560*16^n + 600*4^n - 10)/3 for n > 0.
%F G.f.: (357 + 6953*x - 51812*x^2 + 44352*x^3)/((1-x)*(1-4*x)*(1-16*x)).
%F a(0)=357, a(1)=14450, a(2)=221650, a(3)=3508050, a(n)=21*a(n-1)- 84*a(n-2)+ 64*a(n-3). - _Harvey P. Dale_, Jun 18 2014
%F E.g.f.: (1/3)*(-10*exp(x) + 600*exp(4*x) + 2560*exp(16*x)) - 693. - _G. C. Greubel_, May 28 2016
%t Join[{357},LinearRecurrence[{21,-84,64},{14450,221650,3508050},20]] (* _Harvey P. Dale_, Jun 18 2014 *)
%o (PARI) m=15; v=concat([357, 14450, 221650], vector(m-3)); for(n=4, m, v[n]=20*v[n-1]-64*v[n-2]-150); v
%Y Cf. A075153, A166912, A166914, A166915, A166916, A166917, A167120, A167121, A167122.
%K nonn,easy
%O 0,1
%A _Klaus Brockhaus_, Oct 27 2009
|
