OFFSET
0,5
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..400
FORMULA
a(n) = Sum_{k=0..floor(n/4)} binomial(n,k) * binomial(k,n-4*k).
The g.f. exp( Sum_{k>=1} a(k) * x^k/k ) has integer coefficients and equals (1/x) * Series_Reversion( x / (1+x^4+x^5) ). See A365730.
MATHEMATICA
Table[SeriesCoefficient[Series[(1+t^4+t^5)^n, {t, 0, n}], n], {n, 0, 40}] (* Vincenzo Librandi, Sep 20 2025 *)
PROG
(PARI) a(n) = sum(k=0, n\4, binomial(n, k)*binomial(k, n-4*k));
(Magma) R<t> := PolynomialRing(Integers()); seq := [ MonomialCoefficient((1 + t^4 + t^5)^n, t^n) : n in [0..40] ]; seq; // Vincenzo Librandi, Sep 20 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Sep 19 2025
STATUS
approved
