OFFSET
0,3
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..500 (terms 0..200 from Andrew Howroyd)
M. Dziemianczuk, Enumerations of plane trees with multiple edges and Raney lattice paths, Discrete Mathematics 337 (2014): 9-24.
Vaclav Kotesovec, Recurrence (of order 10)
FORMULA
a(n) = Sum_{k=1..n+1} Sum_{i=1..k-1} Sum_{j=0..floor((n-i)/4)} (-1)^j*binomial(k, i)*binomial(i, j)*binomial(n-i, k-i-1)*binomial(n-4*j-1, i-1)/k for n > 0. - Andrew Howroyd, Feb 24 2020
a(n) ~ c * d^n / n^(3/2), where d = 4.60104107391278274868305060007404493889954... is the largest positive real root of the equation 16 - 48*d + 4*d^2 + 28*d^3 - 181*d^4 + 976*d^5 - 666*d^6 - 222*d^7 - 69*d^8 - 94*d^9 + 27*d^10 = 0 and c = 0.91053827246762074415410258833655792134... - Vaclav Kotesovec, Feb 08 2026
MATHEMATICA
Join[{1}, Table[Sum[Sum[Sum[(-1)^j * Binomial[k, i] * Binomial[i, j] * Binomial[n-i, k-i-1] * Binomial[n-4*j-1, i-1]/k, {j, 0, Floor[(n-i)/4]}], {i, 1, k-1} ], {k, 1, n+1}], {n, 1, 30}]] (* Vaclav Kotesovec, Feb 08 2026 *)
PROG
(PARI) a(n)={my(m=4); if(n<1, n==0, sum(k=1, n+1, sum(i=1, k-1, sum(j=0, (n-i)\m, (-1)^j*binomial(k, i)*binomial(i, j)*binomial(n-i, k-i-1)*binomial(n-m*j-1, i-1)))/k))} \\ Andrew Howroyd, Feb 24 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Sep 14 2014
EXTENSIONS
Terms a(11) and beyond from Andrew Howroyd, Feb 24 2020
STATUS
approved
