%I #25 Feb 19 2024 01:53:25
%S 1,9,57,313,1593,7737,36409,167481,757305,3378745,14913081,65244729,
%T 283348537,1222872633,5249404473,22429273657,95443717689,404681363001,
%U 1710351420985,7207909559865,30297653743161,127054676987449
%N a(1)=1, a(n) = n*4^(n-1) + a(n-1).
%H Vincenzo Librandi, <a href="/A014916/b014916.txt">Table of n, a(n) for n = 1..1000</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (9,-24,16).
%F a(1)=1, a(2)=9, a(n) = 8*a(n-1) - 16*a(n-2) + 1. - _Vincenzo Librandi_, Oct 23 2012
%F G.f.: x/((1-x)*(1-4*x)^2). - _Vincenzo Librandi_, Oct 23 2012
%p a:=n->sum (4^n-4^j, j=0..n): seq(a(n)/3, n=1..31); # _Zerinvary Lajos_, Dec 14 2008
%t Join[{a=1,b=9},Table[c=8*b-16*a+1;a=b;b=c,{n,60}]] (* _Vladimir Joseph Stephan Orlovsky_, Feb 07 2011 *)
%t CoefficientList[Series[1/((1 - x)(1 - 4*x)^2), {x, 0, 40}], x] (* _Vincenzo Librandi_, Oct 23 2012 *)
%t LinearRecurrence[{9,-24,16},{1,9,57},30] (* _Harvey P. Dale_, Jul 25 2015 *)
%o (Magma) I:=[1, 9]; [n le 2 select I[n] else 8*Self(n-1)-16*Self(n-2)+ 1: n in [1..30]]; // _Vincenzo Librandi_, Oct 23 2012
%K nonn,easy
%O 1,2
%A _Olivier GĂ©rard_
|