OFFSET
1,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..990
Index entries for linear recurrences with constant coefficients, signature (40,-675,6294,-35679,127548,-289173,409062,-347112,161056,-31360).
FORMULA
G.f.: x*(1 +126*x -483*x^2 -3884*x^3 +15300*x^4 -10848*x^5 -8960*x^6)/ ( (1-10*x) *(1-7*x)^2 *(1-4*x)^3 *(1-x)^4 ). - R. J. Mathar, Sep 14 2013
a(n) = (1/162)*( 8*10^(n+3) - 30*(3*n +8)*7^(n+2) + 6*(9*n^2 +39*n +40)*4^(n+2) - (27*n^3 +135*n^2 +198*n +80)). - G. C. Greubel, Mar 15 2022
MATHEMATICA
T[n_, k_, m_]:= T[n, k, m]= If[k==1 || k==n, 1, (m*n-m*k+1)*T[n-1, k-1, m] + (m*k - m+1)*T[n-1, k, m]];
A144380[n_]:= T[n+3, n, 3];
Table[A144380[n], {n, 30}] (* modified by G. C. Greubel, Mar 15 2022 *)
PROG
(Magma) [(1/162)*( 8*10^(n+3) - 30*(3*n +8)*7^(n+2) + 6*(9*n^2 +39*n +40)*4^(n+2) - (27*n^3 +135*n^2 +198*n +80)): n in [1..30]]; // G. C. Greubel, Mar 15 2022
(Sage)
@CachedFunction
def T(n, k, m):
if (k==1 or k==n): return 1
else: return (m*(n-k)+1)*T(n-1, k-1, m) + (m*k-m+1)*T(n-1, k, m)
def A144380(n): return T(n+3, n, 3)
[A144380(n) for n in (1..30)] # G. C. Greubel, Mar 15 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula and Gary W. Adamson, Oct 01 2008
STATUS
approved