%I #38 May 04 2023 12:35:08
%S 1,67,927,6677,32942,126378,404634,1129854,2833479,6515509,13947505,
%T 28115451,53846156,98669156,173975076,296541132,490504893,789878583,
%U 1241708083,1909993393,2880500634,4266609710,6216356510,8920844010,12624212835,17635378761
%N Third partial sums of sixth powers (A001014).
%C This is one of a sequence of arrays that are the convolutions of the zero-padded sequences binomial(2n-1+k,k) with the Eulerian polynomials E(n,x) of A008292, represented by E(n,x) (1-x)^(-2n), which generate increasing partial sums of powers of integers:
%C n= 2) (1 + 4*x + x^2)/(1-x)^4 is the o.g.f. of A000578, the convolution of (1,4,1) with A000292, giving the powers of m^3.
%C n= 3) (1 + 11*x + 11*x^2 + x^3)/(1-x)^6 is the o.g.f. of A000538, convolution of (1,11,11,1) with A000389, giving the partial sums of m^4.
%C n= 4) (1 + 26*x + 66*x^2 + 26*x^3 + x^4)/(1-x)^8, the o.g.f. of A101092, convolution of (1,26,66,26,1) with A000580, the second partial sums of m^5
%C n= 5) (1 + 57*x + 302*x^2 + 302*x^3 + 57*x^4 + x^5)/(1-x)^10, the o.g.f. of A254460, convolution of (1,57,302,302,57,1) with A000582, giving the third partial sums of m^6. - _Tom Copeland_, Dec 07 2015
%H Colin Barker, <a href="/A254640/b254640.txt">Table of n, a(n) for n = 1..1000</a>
%H Luciano Ancora, <a href="/A254640/a254640_1.pdf">Partial sums of m-th powers with Faulhaber polynomials</a>
%H <a href="/index/Rec#order_10">Index entries for linear recurrences with constant coefficients</a>, signature (10,-45,120,-210,252,-210,120,-45,10,-1).
%F a(n) = n*(1+n)*(2+n)*(3+n)*(3+2*n)*(2 -30*n +35*n^2 +30*n^3 +5*n^4)/5040.
%F G.f.: x*(1+x)*(1 +56*x +246*x^2 +56*x^3 +x^4)/(1-x)^10. - _Colin Barker_, Feb 04 2015
%p seq(binomial(n+3,4)*(2*n+3)*(5*n^4 +30*n^3 +35*n^2 -30*n +2)/210, n=1..30); # _G. C. Greubel_, Aug 28 2019
%t Table[n(1+n)(2+n)(3+n)(3+2n)(2 -30n +35n^2 +30n^3 +5n^4)/5040, {n, 30}] (* or *) CoefficientList[Series[(x+1)(x^4 +56x^3 +246x^2 +56x +1)/(x - 1)^10, {x, 0, 30}], x] (* _Vincenzo Librandi_, Feb 05 2015 *)
%o (PARI) vector(30, n, n*(1+n)*(2+n)*(3+n)*(3+2*n)*(2-30*n+35*n^2+30*n^3+5*n^4)/5040) \\ _Colin Barker_, Feb 04 2015
%o (Magma) [n*(1+n)*(2+n)*(3+n)*(3+2*n)*(2-30*n+35*n^2+30*n^3+ 5*n^4)/5040: n in [1..30]]; // _Vincenzo Librandi_, Feb 05 2015
%o (Sage) [binomial(n+3,4)*(2*n+3)*(5*n^4 +30*n^3 +35*n^2 -30*n +2)/210 for n in (1..30)] # _G. C. Greubel_, Aug 28 2019
%o (GAP) List([1..30], n-> Binomial(n+3,4)*(2*n+3)*(5*n^4 +30*n^3 +35*n^2 -30*n +2)/210); # _G. C. Greubel_, Aug 28 2019
%o (Python)
%o def A254640(n): return n*(n*(n*(n*(n*(n*(n*(n*(10*n + 135) + 720) + 1890) + 2394) + 945) - 640) - 450) + 36)//5040 # _Chai Wah Wu_, Dec 07 2021
%Y Cf. A008292, A000578, A000292, A000538, A000389, A101092, A000580, A254460, A000582.
%K nonn,easy
%O 1,2
%A _Luciano Ancora_, Feb 04 2015