Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #9 Nov 28 2019 07:59:12
%S 0,1,456,14976,181780,1273970,6293756,24395756,79119496,223821235,
%T 568280240,1321714636,2858876956,5817509516,11237224740,20751835560,
%U 36849296016,63215722181,105182448536,170297734920,269047574180,415753060646,629674964556,936359517556
%N Doubly heptagonal pyramidal numbers.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/HeptagonalPyramidalNumber.html">Heptagonal Pyramidal Number</a>
%H <a href="/index/Ps#pyramidal_numbers">Index to sequences related to pyramidal numbers</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 G.f.: x*(1 + 446*x + 10461*x^2 + 52420*x^3 + 75580*x^4 + 32544*x^5 + 3504*x^6 + 44*x^7)/(1 - x)^10.
%F a(n) = A002413(A002413(n)).
%F a(n) = Sum_{k=0..A002413(n)} A000566(k).
%F a(n) = n *(5*n-2) *(n+1) *(5*n^3+3*n^2-2*n+6) *(25*n^3+15*n^2-10*n-12)/1296. - _R. J. Mathar_, Nov 28 2019
%t A002413[n_] := n (n + 1) (5 n - 2)/6; a[n_] := A002413[A002413[n]]; Table[a[n], {n, 0, 25}]
%t Table[Sum[k (5 k - 3)/2, {k, 0, n (n + 1) (5 n - 2)/6}], {n, 0, 25}]
%t nmax = 25; CoefficientList[Series[x (1 + 446 x + 10461 x^2 + 52420 x^3 + 75580 x^4 + 32544 x^5 + 3504 x^6 + 44 x^7)/(1 - x)^10, {x, 0, nmax}], x]
%t LinearRecurrence[{10, -45, 120, -210, 252, -210, 120, -45, 10, -1}, {0, 1, 456, 14976, 181780, 1273970, 6293756, 24395756, 79119496, 223821235}, 26]
%Y Cf. A000566, A002413, A140236, A264891, A329753, A329754, A329755, A329757.
%K nonn,easy
%O 0,3
%A _Ilya Gutkovskiy_, Nov 20 2019