%I #35 Sep 08 2022 08:46:01
%S 0,1,6,20,48,95,166,266,400,573,790,1056,1376,1755,2198,2710,3296,
%T 3961,4710,5548,6480,7511,8646,9890,11248,12725,14326,16056,17920,
%U 19923,22070,24366,26816,29425,32198,35140,38256,41551,45030,48698
%N a(n) = n*(5*n^2 - 3*n + 4) / 6.
%H Vincenzo Librandi, <a href="/A203552/b203552.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) = Sum_{k = 1..n} A(k-1, n-k) where A(i, j) = i^2 + i*j + j^2 + i + j + 1.
%F G.f.: x * (1 + 2*x + 2*x^2) / (1 - x)^4.
%F a( n) = -A203551(-n) for all n in Z.
%F a(n)-a(n-1) = A134238(n). - _Bruno Berselli_, Jan 03 2012
%F a(n) = 4*A000125(n) + 2*A000125(n+1) - A000125(n+3). - _Ivan N. Ianakiev_, Aug 21 2013
%F E.g.f.: x*(5*x^2 + 12*x + 6)*exp(x)/6. - _G. C. Greubel_, Aug 12 2018
%e G.f. = x + 6*x^2 + 20*x^3 + 48*x^4 + 95*x^5 + 166*x^6 + 266*x^7 + 400*x^8 + ...
%t LinearRecurrence[{4,-6,4,-1},{0,1,6,20},40] (* _Vincenzo Librandi_, Jan 07 2012 *)
%o (PARI) {a(n) = n * (5*n^2 - 3*n + 4) / 6};
%o (Magma) I:=[0, 1, 6, 20]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..50]]; // _Vincenzo Librandi_, Jan 07 2012
%Y Cf. A203551.
%K nonn,easy
%O 0,3
%A _Michael Somos_, Jan 02 2012