login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A234466
a(n) = 7*binomial(8*n+7,n)/(8*n+7).
13
1, 7, 77, 1015, 14763, 228459, 3689595, 61474519, 1048927880, 18236463245, 321899509386, 5753527081211, 103922382296180, 1893943017506925, 34783258504651434, 643111366544129175, 11960812088346090200, 223614812152492437432, 4200107505573406222425
OFFSET
0,2
COMMENTS
Fuss-Catalan sequence is a(n,p,r) = r*binomial(np+r,n)/(np+r), this is the case p=8, r=7.
LINKS
J-C. Aval, Multivariate Fuss-Catalan Numbers, arXiv:0711.0906v1, Discrete Math., 308 (2008), 4660-4669.
Thomas A. Dowling, Catalan Numbers Chapter 7.
Elżbieta Liszewska, Wojciech Młotkowski, Some relatives of the Catalan sequence, arXiv:1907.10725 [math.CO], 2019.
Wojciech Mlotkowski, Fuss-Catalan Numbers in Noncommutative Probability, Docum. Mathm. 15: 939-955.
FORMULA
G.f. satisfies: B(x) = {1 + x*B(x)^(p/r)}^r, where p=8, r=7.
E.g.f.: hypergeom([7, 9, 10, 11, 12, 13, 14]/8, [8, 9, 10, 11, 12, 13, 14]/7, (8^8/7^7)*x). Cf.: Ilya Gutkovskiy in A118971. - Wolfdieter Lang, Feb 06 2020
D-finite with recurrence: +7*(7*n+4)*(7*n+1)*(7*n+5)*(7*n+2)*(7*n+6)*(7*n+3)*(n+1)*a(n) -128*(8*n+3)*(4*n+3)*(8*n+1)*(2*n+1)*(8*n-1)*(4*n+1)*(8*n+5)*a(n-1)=0. - R. J. Mathar, Feb 21 2020
From Wolfdieter Lang, Feb 15 2024: (Start)
a(n) = binomial(8*n + 6, n+1)/(7*n + 6). This is instance k = 7 of c(k, n+1) given in a comment in A130564.
The compositional inverse of y*(1 - y)^7 is x*G(x), where G is the o.g.f.. That is, G(x)*(1 - x*G(x))^7 = 1. This is equivalent to the formula of the first line above with B = G. Take A = B^(1/7) then A*(1 - x*B) = 1 or B*(1 - x*B)^7 = 1.
The o.g.f is G(x) = 8F7([7..14]/8, [8..14]/7; (8^8/7^7)*x) = (7/(8*x))*(1 - 7F6([-1,1,2,3,4,5,6]/8, [1,2,3,4,5,6]/7; (8^8/7^7)*x)). See the e.g.f. above.(End)
MATHEMATICA
Table[7 Binomial[8 n + 7, n]/(8 n + 7), {n, 0, 40}] (* Vincenzo Librandi, Dec 26 2013 *)
PROG
(PARI) a(n) = 7*binomial(8*n+7, n)/(8*n+7);
(PARI) {a(n)=local(B=1); for(i=0, n, B=(1+x*B^(8/7))^7+x*O(x^n)); polcoeff(B, n)}
(Magma) [7*Binomial(8*n+7, n)/(8*n+7): n in [0..30]]; // Vincenzo Librandi, Dec 26 2013
KEYWORD
nonn,easy
AUTHOR
Tim Fulford, Dec 26 2013
STATUS
approved