|
%I #34 Sep 08 2022 08:46:06
%S 0,1,20,74,180,355,616,980,1464,2085,2860,3806,4940,6279,7840,9640,
%T 11696,14025,16644,19570,22820,26411,30360,34684,39400,44525,50076,
%U 56070,62524,69455,76880,84816,93280,102289,111860,122010,132756,144115,156104,168740
%N a(n) = n*(n + 1)*(17*n - 14)/6.
%C Also 19-gonal (or nonadecagonal) pyramidal numbers.
%C This sequence is related to A180232 by 2*a(n) = n*A180232(n) - Sum_{i=0..n-1} A180232(i).
%D E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93 (seventeenth row of the table).
%H Bruno Berselli, <a href="/A237617/b237617.txt">Table of n, a(n) for n = 0..1000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PyramidalNumber.html">Pyramidal Number</a>.
%H <a href="/index/Ps#pyramidal_numbers">Index to sequences related to pyramidal numbers</a>
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (4,-6,4,-1).
%F G.f.: x*(1 + 16*x)/(1 - x)^4.
%F For n>0, a(n) = Sum_{i=0..n-1} (n-i)*(17*i+1); see the generalization in A237616 (Formula field).
%F E.g.f.: x*(6 + 54*x + 17*x^2)*exp(x)/6. - _G. C. Greubel_, Aug 30 2019
%e After 0, the sequence is provided by the row sums of the triangle:
%e 1;
%e 2, 18;
%e 3, 36, 35;
%e 4, 54, 70, 52;
%e 5, 72, 105, 104, 69;
%e 6, 90, 140, 156, 138, 86;
%e 7, 108, 175, 208, 207, 172, 103;
%e 8, 126, 210, 260, 276, 258, 206, 120;
%e 9, 144, 245, 312, 345, 344, 309, 240, 137;
%e 10, 162, 280, 364, 414, 430, 412, 360, 274, 154; etc.,
%e where (r = row index, c = column index):
%e T(r,r) = T(c,c) = 17*r - 16 and T(r,c) = T(r-1,c) + T(r,r) = (r-c+1) * T(r,r), with r>=c>0.
%p seq(n*(n+1)*(17*n-14)/6, n=0..40); # _G. C. Greubel_, Aug 30 2019
%t Table[n(n+1)(17*n-14)/6, {n, 0, 40}]
%t CoefficientList[Series[x(1+16x)/(1-x)^4, {x,0,40}], x] (* _Vincenzo Librandi_, Feb 12 2014 *)
%t LinearRecurrence[{4,-6,4,-1},{0,1,20,74},40] (* _Harvey P. Dale_, Aug 04 2021 *)
%o (Magma) [n*(n+1)*(17*n-14)/6: n in [0..40]];
%o (Magma) I:=[0,1,20,74]; [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_, Feb 12 2014
%o (PARI) vector(40, n, n*(n-1)*(17*n-31)/6) \\ _G. C. Greubel_, Aug 30 2019
%o (Sage) [n*(n+1)*(17*n-14)/6 for n in (0..40)] # _G. C. Greubel_, Aug 30 2019
%o (GAP) List([0..40], n-> n*(n+1)*(17*n-14)/6); # _G. C. Greubel_, Aug 30 2019
%Y Cf. A051871, A180232.
%Y Cf. similar sequences listed in A237616.
%K nonn,easy
%O 0,3
%A _Bruno Berselli_, Feb 11 2014
|
