login
A346650
a(n) = Sum_{k=0..n} binomial(n,k) * binomial(8*k,k) / (7*k + 1).
11
1, 2, 11, 120, 1661, 25484, 415619, 7066670, 123865313, 2222178999, 40604688117, 753051711707, 14138552326609, 268204210248763, 5132686807360949, 98973130183436759, 1921142366704203305, 37508070639707177792, 736080632477073862271, 14511777729474947626918
OFFSET
0,2
COMMENTS
Binomial transform of A007556.
In general, for m > 1, Sum_{k=0..n} binomial(n,k) * binomial(m*k,k) / ((m-1)*k + 1) ~ (m^m + (m-1)^(m-1))^(n + 3/2) / (sqrt(2*Pi) * m^((3*m-1)/2) * n^(3/2) * (m-1)^((m-1)*n + 3/2)). - Vaclav Kotesovec, Jul 30 2021
LINKS
FORMULA
G.f. A(x) satisfies: A(x) = 1 / (1 - x) + x * (1 - x)^6 * A(x)^8.
G.f.: Sum_{k>=0} ( binomial(8*k,k) / (7*k + 1) ) * x^k / (1 - x)^(k+1).
a(n) ~ 17600759^(n + 3/2) / (34359738368 * sqrt(Pi) * n^(3/2) * 7^(7*n + 3/2)). - Vaclav Kotesovec, Jul 30 2021
MATHEMATICA
Table[Sum[Binomial[n, k] Binomial[8 k, k]/(7 k + 1), {k, 0, n}], {n, 0, 19}]
nmax = 19; A[_] = 0; Do[A[x_] = 1/(1 - x) + x (1 - x)^6 A[x]^8 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
nmax = 19; CoefficientList[Series[Sum[(Binomial[8 k, k]/(7 k + 1)) x^k/(1 - x)^(k + 1), {k, 0, nmax}], {x, 0, nmax}], x]
Table[HypergeometricPFQ[{1/8, 1/4, 3/8, 1/2, 5/8, 3/4, 7/8, -n}, {2/7, 3/7, 4/7, 5/7, 6/7, 1, 8/7}, -16777216/823543], {n, 0, 19}]
PROG
(PARI) a(n) = sum(k=0, n, binomial(n, k)*binomial(8*k, k)/(7*k + 1)); \\ Michel Marcus, Jul 26 2021
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jul 26 2021
STATUS
approved