A137966 G.f. satisfies A(x) = 1+x + x^2*A(x)^6.

1, 1, 1, 6, 21, 86, 396, 1812, 8607, 41958, 207333, 1040234, 5281965, 27078756, 140021248, 729369474, 3823598232, 20158251814, 106809280563, 568471343322, 3037782047947, 16292380484454, 87669285293451, 473172657154822

Offset

0,4

Links

Seiichi Manyama, Table of n, a(n) for n = 0..1000

Formula

a(n) = Sum_{k=0..n-1} C(n-k,k)/(n-k) * C(6*k,n-k-1) for n>0 with a(0)=1. - Paul D. Hanna, Jun 16 2009

Recurrence: 5*(n-1)*n*(5*n - 26)*(5*n - 21)*(5*n - 19)*(5*n - 16)*(5*n - 14)*(5*n - 12)*(5*n - 11)*(5*n - 9)*(5*n - 7)*(5*n - 6)*(5*n - 4)*(5*n - 2)*(5*n + 2)*a(n) = + 576*(n-1)^2*(3*n - 5)*(3*n - 4)*(3*n - 2)*(3*n - 1)*(5*n - 26)*(5*n - 21)*(5*n - 19)*(5*n - 16)*(5*n - 14)*(5*n - 12)*(5*n - 11)*(5*n - 9)*(5*n - 7)*a(n-2) + 288*(3*n - 5)*(3*n - 4)*(5*n - 26)*(5*n - 21)*(5*n - 19)*(5*n - 16)*(5*n - 14)*(5*n - 12)*(11250*n^7 - 111375*n^6 + 440175*n^5 - 888545*n^4 + 975241*n^3 - 574177*n^2 + 165869*n - 18018)*a(n-3) + 144*(5*n - 26)*(5*n - 21)*(5*n - 19)*(5*n - 6)*(10125000*n^11 - 218700000*n^10 + 2075585625*n^9 - 11378954250*n^8 + 39836289925*n^7 - 92894908470*n^6 + 145953551806*n^5 - 152681445300*n^4 + 102505633480*n^3 - 41086190160*n^2 + 8557182144*n - 670602240)*a(n-4) + 72*(5*n - 26)*(5*n - 11)*(5*n - 6)*(5*n - 4)*(20250000*n^11 - 569025000*n^10 + 7025658750*n^9 - 50083579125*n^8 + 227686012400*n^7 - 687547140050*n^6 + 1391232445598*n^5 - 1854143517725*n^4 + 1550931293540*n^3 - 737424345140*n^2 + 162058858752*n - 10360465920)*a(n-5) + 144*(5*n - 16)*(5*n - 11)*(5*n - 9)*(5*n - 6)*(5*n - 4)*(5*n - 2)*(202500*n^9 - 5953500*n^8 + 74924775*n^7 - 526434885*n^6 + 2255339082*n^5 - 6025054075*n^4 + 9796892735*n^3 - 8893818500*n^2 + 3545754268*n - 142331280)*a(n-6) + 72*n*(2*n - 9)*(3*n - 17)*(3*n - 10)*(5*n - 21)*(5*n - 16)*(5*n - 14)*(5*n - 11)*(5*n - 9)*(5*n - 7)*(5*n - 6)*(5*n - 4)*(5*n - 2)*(6*n - 41)*(6*n - 13)*a(n-7). - Vaclav Kotesovec, Nov 18 2017

a(n) ~ sqrt((1 + 2*r*s^6)/(15*Pi)) / (2*s^2 * n^(3/2) * r^(n + 1/2)), where r = 0.1734895129039028676461340698295316044509963479582... and s = 1.408187415484683441175360883795437925341195617549... are roots of the system of equations 1 + r + r^2*s^6 = s, 6*r^2*s^5 = 1. - Vaclav Kotesovec, Nov 18 2017

Mathematica

Flatten[{1, Table[Sum[Binomial[n-k, k]/(n-k) * Binomial[6*k, n-k-1], {k, 0, n-1}], {n, 1, 30}]}] (* Vaclav Kotesovec, Nov 18 2017 *)

Prog

(PARI) {a(n)=local(A=1+x*O(x^n)); for(i=0, n, A=1+x*(1+x*A^6)^1); polcoeff(A, n)}
(PARI) a(n)=if(n==0, 1, sum(k=0, n-1, binomial(n-k, k)/(n-k)*binomial(6*k, n-k-1))) \\ Paul D. Hanna, Jun 16 2009

Crossrefs

Cf. A137967, A137965; A019497, A137954, A137959.

Sequence in context: A191597 A088556 A316105 * A304187 A005498 A002222

Adjacent sequences: A137963 A137964 A137965 * A137967 A137968 A137969

Keyword

nonn

Author

Paul D. Hanna, Feb 26 2008

Status

approved