|
%I #38 Mar 02 2023 02:53:56
%S 0,8,35,99,224,440,783,1295,2024,3024,4355,6083,8280,11024,14399,
%T 18495,23408,29240,36099,44099,53360,64008,76175,89999,105624,123200,
%U 142883,164835,189224,216224,246015,278783,314720,354024,396899,443555,494208
%N a(n) = (n-1)*(n+2)*(n^2 + n + 2)/4.
%C The products of two consecutive numbers in this sequence may be evaluated in terms of the Frobenius numbers for 5 consecutive integers, A138985(n) = F(n): for n>0, a(2n-1)*a(2n) = F(4n^2-2)^2 - (2n)^2; a(2n)*a(2n+1) = F(4n^2+4n)^2 - (2n+1)^2. - _Charlie Marion_, Jan 23 2012
%H Vincenzo Librandi, <a href="/A168566/b168566.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (5,-10,10,-5,1).
%F G.f.: x^2*(-8 + 5*x - 4*x^2 + x^3)/(x-1)^5. - _R. J. Mathar_, Jan 04 2011
%F a(n) = A000217(n)^2 - 1 = A000537(n)-1. - _Charlie Marion_, Sep 27 2011
%F From _G. C. Greubel_, Jul 26 2016: (Start)
%F a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5).
%F E.g.f.: (1/4)*(4 - (4 - 4*x + 14*x^2 + 8*x^3)*exp(x)). (End)
%F Sum_{n>=2} 1/a(n) = 49/36 - tanh(sqrt(7)*Pi/2)*Pi/sqrt(7). - _Amiram Eldar_, Mar 02 2023
%t s=0;lst={s};Do[s+=n^3;AppendTo[lst,s],{n,2,5!}];lst (* _Vladimir Joseph Stephan Orlovsky_, Apr 27 2010 *)
%t Table[(n-1)*(n+2)*(n^2 + n + 2)/4, {n,1,25}] (* _G. C. Greubel_, Jul 26 2016 *)
%t LinearRecurrence[{5,-10,10,-5,1},{0,8,35,99,224},40] (* _Harvey P. Dale_, Jan 21 2019 *)
%o (Magma) [(n-1)*(n+2)*(n^2+n+2)/4: n in [1..50]]
%o (PARI) a(n)=(n-1)*(n+2)*(n^2+n+2)/4 \\ _Charles R Greathouse IV_, Jul 26 2016
%Y Cf. A000578 (first differences), A000217, A000537, A138985.
%K nonn,easy
%O 1,2
%A _Vincenzo Librandi_, Nov 30 2009
|
