OFFSET
1,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..200
G. C. Greubel, Generating functions and recurrence
Index entries for linear recurrences with constant coefficients, signature (462, -91531, 10359990, -751926160, 37211988696, -1302889937812, 33041329661136, -616216166660230, 8532037831808700, -88161849837783250, 681071135526667500, -3928678356172712500, 16849648606559250000, -53333569862180312500, 123260625734287500000, -205051485813712890625, 240815404573652343750, -193749765819873046875, 101292345769042968750, -30943180590820312500, 4187798876953125000).
FORMULA
From G. C. Greubel, Sep 10 2022: (Start)
a(n) = Sum_{j=0..5} (-1)^j*binomial(4*n+1, j)*binomial(9-j, 4)^n.
G.f., e.g.f., and recurrence are in the file "Generating functions and recurrence". (End)
MATHEMATICA
Table[Sum[(-1)^j*Binomial[4*n+1, j]*Binomial[9-j, 4]^n, {j, 0, 5}], {n, 30}] (* G. C. Greubel, Sep 10 2022 *)
PROG
(Magma) [(&+[(-1)^j*Binomial(4*n+1, j)*Binomial(9-j, 4)^n: j in [0..5]]): n in [1..30]]; // G. C. Greubel, Sep 10 2022
(SageMath)
def A151643(n): return sum((-1)^j*binomial(4*n+1, j)*binomial(9-j, 4)^n for j in (0..5))
[A151643(n) for n in (1..30)] # G. C. Greubel, Sep 10 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
R. H. Hardin, May 29 2009
EXTENSIONS
Terms a(8) and beyond from Andrew Howroyd, May 06 2020
STATUS
approved