%I #39 Jan 10 2023 01:50:05
%S 2,14,45,104,200,342,539,800,1134,1550,2057,2664,3380,4214,5175,6272,
%T 7514,8910,10469,12200,14112,16214,18515,21024,23750,26702,29889,
%U 33320,37004,40950,45167,49664,54450,59534,64925,70632,76664,83030,89739,96800,104222,112014
%N Product of n with sum of next n consecutive integers.
%H Vincenzo Librandi, <a href="/A036659/b036659.txt">Table of n, a(n) for n = 1..1000</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) = 2 + (3*n*n + 10*n + 11)*n/2. - _Frank Ellermann_, Mar 16 2002
%F a(n) = n^2(3n + 1)/2. - Ronnie B Kon (ronniekon(AT)yahoo.com), Jun 09 2006
%F G.f.: x*(2 + 6*x + x^2)/(1 - x)^4. - _Colin Barker_, Mar 27 2012
%F From _Amiram Eldar_, Jan 10 2023: (Start)
%F Sum_{n>=1} 1/a(n) = Pi^2/3 + sqrt(3)*Pi + 9*log(3) - 18.
%F Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/6 - 2*sqrt(3)*Pi - 12*log(2) + 18. (End)
%e n=3 -> 3*(4+5+6)=45
%t CoefficientList[Series[(2+6*x+x^2)/(1-x)^4,{x,0,40}],x] (* _Vincenzo Librandi_, Apr 25 2012 *)
%t Table[(n^2*(1+3n))/2,{n,40}] (* or *) LinearRecurrence[{4,-6,4,-1},{2,14,45,104},40] (* _Harvey P. Dale_, Jan 04 2015 *)
%o (PARI) a(n)=n^2*(3*n+1)/2 \\ _Charles R Greathouse IV_, Dec 08 2011
%o (Magma) I:=[2, 14, 45, 104]; [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_, Apr 25 2012
%Y Cf. A036660.
%K nonn,easy
%O 1,1
%A _G. L. Honaker, Jr._, Dec 15 1998
|