|
%I #10 Sep 08 2022 08:45:39
%S 4,38,136,330,652,1134,1808,2706,3860,5302,7064,9178,11676,14590,
%T 17952,21794,26148,31046,36520,42602,49324,56718,64816,73650,83252,
%U 93654,104888,116986,129980,143902,158784,174658,191556,209510,228552,248714
%N G.f.: -2*(-2 - 11*x - 4*x^2 + x^3)/(x - 1)^4.
%H Vincenzo Librandi, <a href="/A152110/b152110.txt">Table of n, a(n) for n = 0..1000</a>
%F a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). a(n) = 2*(n+1)*(8*n^2 + 8*n + 3)/3 = A152133(2*n-1). - _R. J. Mathar_, Sep 22 2009
%t CoefficientList[Series[-2*(-2 - 11*x - 4*x^2 + x^3)/(x - 1)^4, {x, 0, 40}], x] (* _Vincenzo Librandi_, Dec 11 2012 *)
%o (Magma) I:=[4,38,136,330]; [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..40]]; // _Vincenzo Librandi_, Dec 11 2012
%Y Cf. A006071.
%K nonn,easy
%O 0,1
%A _N. J. A. Sloane_, Sep 21 2009, based on email from _R. J. Mathar_, Mar 22 2009
%E More terms from _R. J. Mathar_, Sep 22 2009
|
