%I #15 May 21 2023 12:05:17
%S 1,3,17,59,145,291,513,827,1249,1795,2481,3323,4337,5539,6945,8571,
%T 10433,12547,14929,17595,20561,23843,27457,31419,35745,40451,45553,
%U 51067,57009,63395,70241,77563,85377,93699,102545,111931,121873,132387,143489,155195
%N a(n) = 1 - 2*n^2 + 4*n*(1 + 2*n^2)/3.
%C Binomial transform of the quasi-finite sequence 1,2,12,16,0,... (0 continued).
%C A bisection of A168582.
%H Vincenzo Librandi, <a href="/A168547/b168547.txt">Table of n, a(n) for n = 0..10000</a>
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (4,-6,4,-1).
%F a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4).
%F G.f.: (1 - x + 11*x^2 + 5*x^3)/(x-1)^4.
%F First differences: a(n+1) - a(n) = 2*A054569(n+1).
%F Second differences: a(n+2) - 2*a(n+1) + a(n) = 4*A004767(n).
%F Third differences: a(n+3) - 3*a(n+2) + 3*a(n+1) - a(n) = 16.
%F a(n+1) = A166464(n) + A035597(n+1).
%F a(n) = 1 - 2*n^2 + 4*A005900(n). - _R. J. Mathar_, Dec 05 2009
%F E.g.f.: (1/3)*(3 + 6*x + 18*x^2 + 8*x^3)*exp(x). - _G. C. Greubel_, Jul 26 2016
%t Table[1-2*n^2+4*n*(1+2*n^2)/3, {n,0,50}] (* _G. C. Greubel_, Jul 26 2016 *)
%t LinearRecurrence[{4,-6,4,-1},{1,3,17,59},60] (* _Harvey P. Dale_, May 21 2023 *)
%o (Magma) [1-2*n^2+4*n*(1+2*n^2)/3: n in [0..50] ]; // _Vincenzo Librandi_, Aug 06 2011
%o (PARI) a(n)=1-2*n^2+4*n*(1+2*n^2)/3 \\ _Charles R Greathouse IV_, Jul 26 2016
%K nonn,easy
%O 0,2
%A _Paul Curtz_, Nov 29 2009
%E Edited and extended by _R. J. Mathar_, Dec 05 2009
|