%I #79 Feb 09 2024 08:05:34
%S 0,1,8,33,98,238,504,966,1716,2871,4576,7007,10374,14924,20944,28764,
%T 38760,51357,67032,86317,109802,138138,172040,212290,259740,315315,
%U 380016,454923,541198,640088,752928,881144,1026256,1189881,1373736,1579641,1809522,2065414
%N a(n) = n*(n+1)*(n+2)*(n+3)*(3*n+2)/120.
%C 5-dimensional version of pentagonal-based pyramidal numbers. - Ben Creech (mathroxmysox(AT)yahoo.com)
%C If Y is a 3-subset of an n-set X then, for n>=7, a(n-6) is the number of 7-subsets of X having at least two elements in common with Y. - _Milan Janjic_, Nov 23 2007
%C Antidiagonal sums of the convolution array A213548. - _Clark Kimberling_, Jun 17 2012
%C After 0, convolution of nonzero triangular numbers (A000217) and nonzero pentagonal numbers (A000326). - _Bruno Berselli_, Jun 27 2013
%C a(n) is also the number of odd chordless cycles in the graph complement of the (n+1)-Andrásfai graph. - _Eric W. Weisstein_, Apr 14 2017
%D Albert H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 194-196.
%D Herbert John Ryser, Combinatorial Mathematics, "The Carus Mathematical Monographs", No. 14, John Wiley and Sons, 1963, pp. 1-8.
%H Vincenzo Librandi, <a href="/A051836/b051836.txt">Table of n, a(n) for n = 0..1000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/AndrasfaiGraph.html">Andrásfai Graph</a>.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/ChordlessCycle.html">Chordless Cycle</a>.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/GraphComplement.html">Graph Complement</a>.
%H <a href="/index/Ps#pyramidal_numbers">Index to sequences related to pyramidal numbers</a>.
%H <a href="/index/Rec#order_06">Index entries for linear recurrences with constant coefficients</a>, signature (6,-15,20,-15,6,-1).
%F a(n) = C(n+4, n)*(3n+5)/5.
%F G.f.: x*(1+2*x)/(1-x)^6. (adapted by _Vincenzo Librandi_, Jul 04 2017)
%F From _Amiram Eldar_, Feb 15 2022: (Start)
%F Sum_{n>=1} 1/a(n) = 135*sqrt(3)*Pi/14 - 1215*log(3)/14 + 925/21.
%F Sum_{n>=1} (-1)^(n+1)/a(n) = 135*sqrt(3)*Pi/7 - 880*log(2)/7 - 355/21. (End)
%F E.g.f.: (1/5!)*x*(120 + 360*x + 240*x^2 + 50*x^3 + 3*x^4)*exp(x). - _G. C. Greubel_, Dec 27 2023
%e By the fourth comment: A000217(1..6) and A000326(1..6) give the term a(6) = 1*21+5*15+12*10+22*6+35*3+51*1 = 504. - _Bruno Berselli_, Jun 27 2013
%p with (combinat):a[0]:=0:for n from 1 to 50 do a[n]:=stirling2(n+2,n)+a[n-1] od: seq(a[n], n=0..34); # _Zerinvary Lajos_, Mar 17 2008
%t Table[n(n + 1)(n + 2)(n + 3)(3n + 2)/120, {n, 0, 60}] (* _Vladimir Joseph Stephan Orlovsky_, Apr 08 2011 *)
%t CoefficientList[Series[x (1 + 2 x) / (1 - x)^6, {x, 0, 33}], x] (* _Vincenzo Librandi_, Jul 04 2017 *)
%t LinearRecurrence[{6,-15,20,-15,6,-1},{0,1,8,33,98,238},40] (* _Harvey P. Dale_, Jun 01 2018 *)
%o (PARI) a(n)=n*(n+1)*(n+2)*(n+3)*(3*n+2)/120 \\ _Charles R Greathouse IV_, Oct 07 2015
%o (Magma) [0] cat [Binomial(n+4, n)*(3*n+5)/5: n in [0..40]]; // _Vincenzo Librandi_, Jul 04 2017
%o (SageMath) [((3*n+2)/(n+4))*binomial(n+4,5) for n in range(41)] # _G. C. Greubel_, Dec 27 2023
%Y Partial sums of A001296.
%Y Cf. A093560 ((3, 1) Pascal, column m=5).
%Y Cf. A000217, A000326, A213548.
%K nonn,easy
%O 0,3
%A _Barry E. Williams_, Dec 12 1999
%E Simpler definition from Ben Creech (mathroxmysox(AT)yahoo.com), Nov 13 2005