OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (4,2,-20,-1,40,8,-32,-16).
FORMULA
a(n) = Sum_{k=0..floor(n/2)} binomial(n-k+3, 3) * binomial(n-k, k) * 2^k.
a(n) = ((350+177*n+21*n^2)*(n+1)*U(n+1) + 2*(277+132*n+15*n^2)*(n+2)*U(n))/ (2*9^3) with U(n) = A001045(n+1), n>=0.
G.f.: 1/(1-(1+2*x)*x)^4 = 1/ ( (1+x)^4*(1-2*x)^4 ).
E.g.f.: (1/4374)*(32*(113 + 366*x + 234*x^2 + 36*x^3)*exp(2*x) - (-758 + 690*x - 153*x^2 + 9*x^3)*exp(-x)). - G. C. Greubel, Sep 29 2022
MATHEMATICA
Table[(1/4374)*(2^(n+4)*(226 +267*n +90*n^2 +9*n^3) +(-1)^n*(758 +555*n +126*n^2 +9*n^3)), {n, 0, 40}] (* G. C. Greubel, Sep 29 2022 *)
PROG
(Magma) [(1/4374)*(2^(n+4)*(226 +267*n +90*n^2 +9*n^3) +(-1)^n*(758 +555*n +126*n^2 +9*n^3)): n in [0..40]]; // G. C. Greubel, Sep 29 2022
(SageMath)
def A073373(n): return (1/4374)*(2^(n+4)*(226+267*n+90*n^2+9*n^3) +(-1)^n*(758 +555*n+126*n^2+9*n^3))
[A073373(n) for n in range(40)] # G. C. Greubel, Sep 29 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Aug 02 2002
STATUS
approved