cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 15 results. Next

A002294 a(n) = binomial(5*n, n)/(4*n + 1).

Original entry on oeis.org

1, 1, 5, 35, 285, 2530, 23751, 231880, 2330445, 23950355, 250543370, 2658968130, 28558343775, 309831575760, 3390416787880, 37377257159280, 414741863546285, 4628362722856425, 51912988256282175, 584909606696793885, 6617078646960613370
Offset: 0

Views

Author

Keywords

Comments

From Wolfdieter Lang, Sep 14 2007: (Start)
a(n), n >= 1, enumerates quintic trees (rooted, ordered, incomplete) with n vertices (including the root).
This is the Pfaff-Fuss-Catalan sequence C^{m}_n for m = 5. See the Graham et al. reference, p. 347. eq. 7.66. See also the Pólya-Szegő reference.
Also 5-Raney sequence. See the Graham et al. reference, pp. 346-347. (End)
a(n) = A258708(3*n, 2*n) for n > 0. - Reinhard Zumkeller, Jun 23 2015
Conjecturally, a(n) is the number of 4-uniform words on the alphabet [n] that avoid the patterns 231 and 221 (see the Defant and Kravitz link). - Colin Defant, Sep 26 2018
From Stillwell (1995), p. 62: "Eisenstein's Theorem. If y^5 + y = x, then y has a power series expansion y = x - x^5 + 10*x^9/2^1 - 15 * 14 * x^13/3! + 20 * 19 * 18*x^17/4! - ...." - Michael Somos, Sep 19 2019
a(n) is the total number of down steps before the first up step in all 4_1-Dyck paths of length 5*n. A 4_1-Dyck path is a lattice path with steps (1, 4), (1, -1) that starts and ends at y = 0 and stays above the line y = -1. - Sarah Selkirk, May 10 2020
Dropping the first 1 (starting from 1, 5, 35, ... with offset 1), the series reversion gives 1, -5, 15, -35, 70, ... (again offset 1), essentially A000332 and row 5 of A027555. - R. J. Mathar, Aug 17 2023
Number of rooted polyominoes composed of n hexagonal cells of the hyperbolic regular tiling with Schläfli symbol {6,oo}. A rooted polyomino has one external edge identified, and chiral pairs are counted as two. A stereographic projection of the {6,oo} tiling on the Poincaré disk can be obtained via the Christensson link. - Robert A. Russell, Jan 27 2024
This is instance k = 5 of the generalized Catalan family {C(k, n)}_{n>=0} given in a comment of A130564. - Wolfdieter Lang, Feb 05 2024

Examples

			There are a(2) = 5 quintic trees (vertex degree <= 5 and 5 possible branchings) with 2 vertices (one of them the root). Adding one more branch (one more vertex) to these five trees yields 5*5 + binomial(5,2) = 35 = a(3) such trees.
G.f. = 1 + x + 5*x^2 + 35*x^3 + 285*x^4 + 2530*x^5 + 23751*x^6 + 231880*x^7 + ...
G.f. = t + t^5 + 5*t^9 + 35*t^13 + 285*t^17 + 2530*t^21 + 23751*t^25 + 231880*t^29 + ...
		

References

  • Archiv der Mathematik u. Physik, Editor's note: "Über die Bestimmung der Anzahl der verschiedenen Arten, auf welche sich ein n-Eck durch Diagonalen in lauter m-Ecke zerlegen laesst, mit Bezug auf einige Abhandlungen der Herren Lame, Rodrigues, Binet, Catalan und Duhamel in dem Journal de Mathematiques pures et appliquees, publie par Joseph Liouville. T. III. IV.", Archiv der Mathematik u. Physik, 1 (1841), pp. 193ff; see especially p. 198.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 23.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, pp. 200, 347.
  • G. Pólya and G. Szegő, Problems and Theorems in Analysis, Springer-Verlag, Heidelberg, New York, 2 vols., 1972, Vol. 1, problem 211, p. 146 with solution on p. 348.
  • Ulrike Sattler, Decidable classes of formal power series with nice closure properties, Diplomarbeit im Fach Informatik, Univ. Erlangen - Nürnberg, Jul 27 1994.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A001764, A002296, A258708, A346647 (binomial transform), A346665 (inverse binomial transform).
Fourth column of triangle A062993.
Polyominoes: A221184{n-1} (oriented), A004127 (unoriented), A369473 (chiral), A143546 (achiral), A002293 {5,oo}, A002295 {7,oo}.
Cf. A130564.

Programs

  • GAP
    List([0..22],n->Binomial(5*n,n)/(4*n+1)); # Muniru A Asiru, Nov 01 2018
  • Haskell
    a002294 n = a002294_list !! n
    a002294_list = [a258708 (3 * n) (2 * n) | n <- [1..]]
    -- Reinhard Zumkeller, Jun 23 2015
    
  • Magma
    [ Binomial(5*n,n)/(4*n+1): n in [0..100]]; // Vincenzo Librandi, Mar 24 2011
    
  • Maple
    seq(binomial(5*k+1,k)/(5*k+1),k=0..30); # Robert FERREOL, Apr 03 2015
    n:=30:G:=series(RootOf(g = 1+x*g^5, g),x=0,n+1):seq(coeff(G,x,k),k=0..n); # Robert FERREOL, Apr 03 2015
  • Mathematica
    CoefficientList[InverseSeries[ Series[ y - y^5, {y, 0, 100}], x], x][[Range[2, 100, 4]]]
    Table[Binomial[5n,n]/(4n+1),{n,0,20}] (* Harvey P. Dale, Dec 30 2011 *)
    a[ n_] := SeriesCoefficient[ HypergeometricPFQ[ {1, 2, 3, 4}/5, {2, 3, 5}/4, x 5^5/4^4], {x, 0, n}]; (* Michael Somos, May 06 2015 *)
    a[ n_] := With[{m = 4 n + 1}, SeriesCoefficient[ InverseSeries @ Series[ x - x^5, {x, 0, m}], {x, 0, m}]]; (* Michael Somos, May 06 2015 *)
  • PARI
    {a(n) = binomial( 5 * n, n) / (4*n + 1)}; /* Michael Somos, Mar 17 2011 */
    
  • PARI
    {a(n) = if( n<0, 0, n = 4*n + 1; polcoeff( serreverse( x - x^5 + x * O(x^n) ), n))}; /* Michael Somos, Mar 17 2011 */
    

Formula

For the connection with the solution of the quintic, hypergeometric series, and Lagrange inversion, see Beukers (2014). - N. J. A. Sloane, Mar 12 2014
G.f.: hypergeometric([1, 2, 3, 4] / 5, [2, 3, 5] / 4, x * 5^5 / 4^4). - Michael Somos, Mar 17 2011
O.g.f. A(x) satisfies A(x) = 1 + x * A(x)^5 = 1 / (1 - x * A(x)^4).
Given g.f. A(x) then z = t * A(t^4) satisfies 0 = z^5 - z + t. - Michael Somos, Mar 17 2011
a(n) = binomial(5*n, n - 1)/n, n >= 1, a(0) = 1. From the Lagrange series of the o.g.f. A(x) with its above given implicit equation.
a(n) = upper left term in M^n, M = the production matrix:
1, 1;
4, 4, 1;
10, 10, 4, 1;
20, 20, 10, 4, 1;
...
where (1, 4, 10, 20, ...) is the tetrahedral sequence, A000292. - Gary W. Adamson, Jul 08 2011
D-finite with recurrence: 8*n*(4*n+1)*(2*n-1)*(4*n-1)*a(n) - 5*(5*n-4)*(5*n-3)*(5*n-2)*(5*n-1)*a(n-1) = 0. - R. J. Mathar, Dec 02 2014
a(n) = binomial(5*n + 1, n)/(5*n + 1) = A062993(n+3,3). - Robert FERREOL, Apr 03 2015
a(0) = 1; a(n) = Sum_{i1 + i2 + ... + i5 = n - 1} a(i1) * a(i2) * ... *a(i5) for n >= 1. - Robert FERREOL, Apr 03 2015
From Ilya Gutkovskiy, Jan 15 2017: (Start)
O.g.f.: 5F4([1/5, 2/5, 3/5, 4/5, 1]; [1/2, 3/4, 1, 5/4]; 3125*x/256).[Cancellation of the 1s, see G.f. the above. - Wolfdieter Lang, Feb 05 2024]
E.g.f.: 4F4([1/5, 2/5, 3/5, 4/5]; [1/2, 3/4, 1, 5/4]; 3125*x/256).
a(n) ~ 5^(5*n + 1/2)/(sqrt(Pi) * 2^(8*n + 7/2) * n^(3/2)). (End)
x*A'(x)/A(x) = (A(x) - 1)/(- 4*A(x) + 5) = x + 9*x^2 + 91*x^3 + 969*x^4 + ... is the o.g.f. of A163456. Cf. A001764 and A002293 - A002296. - Peter Bala, Feb 04 2022
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^9). - Seiichi Manyama, Jun 16 2025

Extensions

More terms from Olivier Gérard, Jul 05 2001

A349332 G.f. A(x) satisfies A(x) = 1 + x * A(x)^5 / (1 - x).

Original entry on oeis.org

1, 1, 6, 46, 406, 3901, 39627, 418592, 4551672, 50610692, 572807157, 6577068383, 76426719408, 897078662538, 10620634999318, 126676885170703, 1520759193166329, 18361269213121164, 222814883564042704, 2716125963857227904, 33244557641365865109
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 15 2021

Keywords

Crossrefs

Programs

  • Maple
    a:= n-> coeff(series(RootOf(1+x*A^5/(1-x)-A, A), x, n+1), x, n):
    seq(a(n), n=0..20);  # Alois P. Heinz, Nov 15 2021
  • Mathematica
    nmax = 20; A[] = 0; Do[A[x] = 1 + x A[x]^5/(1 - x) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    Table[Sum[Binomial[n - 1, k - 1] Binomial[5 k, k]/(4 k + 1), {k, 0, n}], {n, 0, 20}]
  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0);
    A[#A] = 1 + sum(k=1, m-1, (polcoeff(Ser(A)^5, k)) )); A[n+1]}
    for(n=0, 30, print1(a(n), ", ")) \\ Vaclav Kotesovec, Nov 23 2024, after Paul D. Hanna

Formula

a(n) = Sum_{k=0..n} binomial(n-1,k-1) * binomial(5*k,k) / (4*k+1).
a(n) ~ 3381^(n + 1/2) / (25 * sqrt(Pi) * n^(3/2) * 2^(8*n + 7/2)). - Vaclav Kotesovec, Nov 15 2021
Recurrence: 8*n*(2*n - 1)*(4*n - 1)*(4*n + 1)*a(n) = (4405*n^4 - 10346*n^3 + 9575*n^2 - 4354*n + 840)*a(n-1) - 12*(n-2)*(1255*n^3 - 3957*n^2 + 4492*n - 1820)*a(n-2) + 2*(n-3)*(n-2)*(10655*n^2 - 32733*n + 26908)*a(n-3) - 4*(n-4)*(n-3)*(n-2)*(3445*n - 6986)*a(n-4) + 3381*(n-5)*(n-4)*(n-3)*(n-2)*a(n-5). - Vaclav Kotesovec, Nov 17 2021

A346646 a(n) = Sum_{k=0..n} binomial(n,k) * binomial(4*k,k) / (3*k + 1).

Original entry on oeis.org

1, 2, 7, 38, 257, 1935, 15505, 129519, 1115061, 9823160, 88121887, 802227794, 7392428009, 68819554003, 646276497617, 6114880542117, 58237420303109, 557850829527246, 5370956411708779, 51947475492561014, 504492516832543885, 4917564488572565160
Offset: 0

Views

Author

Ilya Gutkovskiy, Jul 26 2021

Keywords

Comments

Binomial transform of A002293.

Crossrefs

Programs

  • Maple
    A346646 := proc(n)
        hypergeom([-n,1/4,1/2,3/4],[2/3,1,4/3],-256/27) ;
        simplify(%) ;
    end proc:
    seq(A346646(n),n=0..40) ; # R. J. Mathar, Jan 10 2023
  • Mathematica
    Table[Sum[Binomial[n, k] Binomial[4 k, k]/(3 k + 1), {k, 0, n}], {n, 0, 21}]
    nmax = 21; A[] = 0; Do[A[x] = 1/(1 - x) + x (1 - x)^2 A[x]^4 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    nmax = 21; CoefficientList[Series[Sum[(Binomial[4 k, k]/(3 k + 1)) x^k/(1 - x)^(k + 1), {k, 0, nmax}], {x, 0, nmax}], x]
    Table[HypergeometricPFQ[{1/4, 1/2, 3/4, -n}, {2/3, 1, 4/3}, -256/27], {n, 0, 21}]
  • PARI
    a(n) = sum(k=0, n, binomial(n,k)*binomial(4*k,k)/(3*k + 1)); \\ Michel Marcus, Jul 26 2021

Formula

G.f. A(x) satisfies: A(x) = 1 / (1 - x) + x * (1 - x)^2 * A(x)^4.
G.f.: Sum_{k>=0} ( binomial(4*k,k) / (3*k + 1) ) * x^k / (1 - x)^(k+1).
a(n) ~ 283^(n + 3/2) / (2048 * sqrt(2*Pi) * n^(3/2) * 3^(3*n + 3/2)). - Vaclav Kotesovec, Jul 30 2021
D-finite with recurrence 3*n*(3*n-1)*(3*n+1)*a(n) -2*(2*n-1) *(91*n^2 -91*n +24)*a(n-1) +6*(n-1) *(155*n^2 -310*n +167)*a(n-2) -438*(n-1) *(n-2)*(2*n-3) *a(n-3) +283*(n-1)*(n-2) *(n-3)*a(n-4)=0. - R. J. Mathar, Aug 17 2023

A349311 G.f. A(x) satisfies: A(x) = (1 + x * A(x)^5) / (1 - x).

Original entry on oeis.org

1, 2, 12, 112, 1232, 14832, 189184, 2512064, 34358784, 480745984, 6848734464, 99003237376, 1448575666176, 21411827808256, 319255531155456, 4796005997940736, 72520546008219648, 1102912584949792768, 16859182461720526848, 258886644574700699648
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 14 2021

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 19; A[] = 0; Do[A[x] = (1 + x A[x]^5)/(1 - x) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    Table[Sum[Binomial[n + 4 k, 5 k] Binomial[5 k, k]/(4 k + 1), {k, 0, n}], {n, 0, 19}]

Formula

a(n) = Sum_{k=0..n} binomial(n+4*k,5*k) * binomial(5*k,k) / (4*k+1).
a(n) = F([(1+n)/4, (2+n)/4, (3+n)/4, (4+n)/4, -n], [1/2, 3/4, 1, 5/4], -1), where F is the generalized hypergeometric function. - Stefano Spezia, Nov 14 2021
a(n) ~ sqrt(1 + 4*r) / (2 * 5^(3/4) * sqrt(2*Pi) * (1-r)^(1/4) * n^(3/2) * r^(n + 1/4)), where r = 0.0600920016324256496641829206872407657377702010870270617... is the real root of the equation 4^4 * (1-r)^5 = 5^5 * r. - Vaclav Kotesovec, Nov 15 2021

A346648 a(n) = Sum_{k=0..n} binomial(n,k) * binomial(6*k,k) / (5*k + 1).

Original entry on oeis.org

1, 2, 9, 73, 751, 8587, 104425, 1323952, 17303503, 231455104, 3153167249, 43597546197, 610232050453, 8629733401556, 123114479858631, 1769728635257503, 25607523627970183, 372688563309335806, 5451995469296025115, 80122698147986922194, 1182341393088427774071
Offset: 0

Views

Author

Ilya Gutkovskiy, Jul 26 2021

Keywords

Comments

Binomial transform of A002295.

Crossrefs

Programs

  • Mathematica
    Table[Sum[Binomial[n, k] Binomial[6 k, k]/(5 k + 1), {k, 0, n}], {n, 0, 20}]
    nmax = 20; A[] = 0; Do[A[x] = 1/(1 - x) + x (1 - x)^4 A[x]^6 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    nmax = 20; CoefficientList[Series[Sum[(Binomial[6 k, k]/(5 k + 1)) x^k/(1 - x)^(k + 1), {k, 0, nmax}], {x, 0, nmax}], x]
    Table[HypergeometricPFQ[{1/6, 1/3, 1/2, 2/3, 5/6, -n}, {2/5, 3/5, 4/5, 1, 6/5}, -46656/3125], {n, 0, 20}]
  • PARI
    a(n) = sum(k=0, n, binomial(n,k)*binomial(6*k,k)/(5*k + 1)); \\ Michel Marcus, Jul 26 2021

Formula

G.f. A(x) satisfies: A(x) = 1 / (1 - x) + x * (1 - x)^4 * A(x)^6.
G.f.: Sum_{k>=0} ( binomial(6*k,k) / (5*k + 1) ) * x^k / (1 - x)^(k+1).
a(n) ~ 49781^(n + 3/2) / (3359232 * sqrt(3*Pi) * n^(3/2) * 5^(5*n + 3/2)). - Vaclav Kotesovec, Jul 30 2021

A346649 a(n) = Sum_{k=0..n} binomial(n,k) * binomial(7*k,k) / (6*k + 1).

Original entry on oeis.org

1, 2, 10, 95, 1146, 15343, 218407, 3241316, 49588850, 776483636, 12383420161, 200444399493, 3284531747403, 54378741581471, 908238222519566, 15284835717461020, 258933935458506210, 4412025177612412048, 75564998345532498844, 1300158755391113561288
Offset: 0

Views

Author

Ilya Gutkovskiy, Jul 26 2021

Keywords

Comments

Binomial transform of A002296.

Crossrefs

Programs

  • Mathematica
    Table[Sum[Binomial[n, k] Binomial[7 k, k]/(6 k + 1), {k, 0, n}], {n, 0, 19}]
    nmax = 19; A[] = 0; Do[A[x] = 1/(1 - x) + x (1 - x)^5 A[x]^7 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    nmax = 19; CoefficientList[Series[Sum[(Binomial[7 k, k]/(6 k + 1)) x^k/(1 - x)^(k + 1), {k, 0, nmax}], {x, 0, nmax}], x]
    Table[HypergeometricPFQ[{1/7, 2/7, 3/7, 4/7, 5/7, 6/7, -n}, {1/3, 1/2, 2/3, 5/6, 1, 7/6}, -823543/46656], {n, 0, 19}]
  • PARI
    a(n) = sum(k=0, n, binomial(n,k)*binomial(7*k,k)/(6*k + 1)); \\ Michel Marcus, Jul 26 2021

Formula

G.f. A(x) satisfies: A(x) = 1 / (1 - x) + x * (1 - x)^5 * A(x)^7.
G.f.: Sum_{k>=0} ( binomial(7*k,k) / (6*k + 1) ) * x^k / (1 - x)^(k+1).
a(n) ~ 870199^(n + 3/2) / (282475249 * sqrt(Pi) * n^(3/2) * 2^(6*n + 2) * 3^(6*n + 3/2)). - Vaclav Kotesovec, Jul 30 2021

A346650 a(n) = Sum_{k=0..n} binomial(n,k) * binomial(8*k,k) / (7*k + 1).

Original entry on oeis.org

1, 2, 11, 120, 1661, 25484, 415619, 7066670, 123865313, 2222178999, 40604688117, 753051711707, 14138552326609, 268204210248763, 5132686807360949, 98973130183436759, 1921142366704203305, 37508070639707177792, 736080632477073862271, 14511777729474947626918
Offset: 0

Views

Author

Ilya Gutkovskiy, Jul 26 2021

Keywords

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

Crossrefs

Programs

  • 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}]
  • PARI
    a(n) = sum(k=0, n, binomial(n,k)*binomial(8*k,k)/(7*k + 1)); \\ Michel Marcus, Jul 26 2021

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

A349290 G.f. A(x) satisfies A(x) = 1 / ((1 - x) * (1 - x * A(x)^4)).

Original entry on oeis.org

1, 2, 11, 96, 1001, 11456, 139013, 1756596, 22867421, 304560171, 4130200726, 56836946342, 791689962811, 11140615233281, 158140107648676, 2261708608884896, 32559326010349817, 471428798399646336, 6860801662510005266, 100302910051255600486
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 13 2021

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 19; A[] = 0; Do[A[x] = 1/((1 - x) (1 - x A[x]^4)) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    Table[Sum[Binomial[n + 3 k, 4 k] Binomial[5 k, k]/(4 k + 1), {k, 0, n}], {n, 0, 19}]
  • PARI
    a(n) = sum(k=0, n, binomial(n+3*k,4*k) * binomial(5*k,k) / (4*k+1)); \\ Michel Marcus, Nov 14 2021

Formula

a(n) = Sum_{k=0..n} binomial(n+3*k,4*k) * binomial(5*k,k) / (4*k+1).
a(n) = F([1/5, 2/5, 3/5, 4/5, (1+n)/3, (2+n)/3, (3+n)/3, -n], [1/4, 1/2, 1/2, 3/4, 3/4, 1, 5/4], -3^3*5^5/2^16), where F is the generalized hypergeometric function. - Stefano Spezia, Nov 13 2021
a(n) ~ sqrt(1 + 3*r) / (2 * 5^(3/4) * sqrt(2*Pi*(1-r)) * n^(3/2) * r^(n + 1/4)), where r = 0.0631152861998150860738633360987635931... is the root of the equation 5^5 * r = 4^4 * (1-r)^4. - Vaclav Kotesovec, Nov 14 2021
a(n) = 1 + Sum_{i, j, k, l, m>=0 and i+j+k+l+m=n-1} a(i) * a(j) * a(k) * a(l) * a(m). - Seiichi Manyama, Jul 10 2025

A346665 a(n) = Sum_{k=0..n} (-1)^(n-k) * binomial(n,k) * binomial(5*k,k) / (4*k + 1).

Original entry on oeis.org

1, 0, 4, 22, 172, 1409, 12216, 109904, 1016876, 9614584, 92490261, 902364918, 8907507708, 88802649446, 892833960460, 9042639746819, 92171773008828, 944819352291920, 9733592874215112, 100725697334689896, 1046535959932600141, 10913073121311627481, 114175868855824821752
Offset: 0

Views

Author

Ilya Gutkovskiy, Jul 27 2021

Keywords

Comments

Inverse binomial transform of A002294.

Crossrefs

Programs

  • Maple
    A346665 := proc(n)
        add((-1)^(n-k)*binomial(n,k)*binomial(5*k,k)/(4*k+1),k=0..n) ;
    end proc:
    seq(A346665(n),n=0..80); # R. J. Mathar, Aug 17 2023
  • Mathematica
    Table[Sum[(-1)^(n - k) Binomial[n, k] Binomial[5 k, k]/(4 k + 1), {k, 0, n}], {n, 0, 22}]
    nmax = 22; A[] = 0; Do[A[x] = 1/(1 + x) + x (1 + x)^3 A[x]^5 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    nmax = 22; CoefficientList[Series[Sum[(Binomial[5 k, k]/(4 k + 1)) x^k/(1 + x)^(k + 1), {k, 0, nmax}], {x, 0, nmax}], x]
    Table[(-1)^n HypergeometricPFQ[{1/5, 2/5, 3/5, 4/5, -n}, {1/2, 3/4, 1, 5/4}, 3125/256], {n, 0, 22}]
  • PARI
    a(n) = sum(k=0, n, (-1)^(n-k)*binomial(n,k)*binomial(5*k,k)/(4*k + 1)); \\ Michel Marcus, Jul 28 2021

Formula

G.f. A(x) satisfies: A(x) = 1 / (1 + x) + x * (1 + x)^3 * A(x)^5.
G.f.: Sum_{k>=0} ( binomial(5*k,k) / (4*k + 1) ) * x^k / (1 + x)^(k+1).
a(n) ~ 2869^(n + 3/2) / (78125 * sqrt(Pi) * n^(3/2) * 2^(8*n + 7/2)). - Vaclav Kotesovec, Jul 30 2021
D-finite with recurrence 8*n*(4*n+1)*(2*n-1)*(4*n-1)*a(n) -(n-1) *(1845*n^3 -1333*n^2 -238*n +240)*a(n-1) -4*(n-1) *(2485*n^3 -7263*n^2 +7388*n -2580) *a(n-2) -2*(n-1) *(n-2) *(8095*n^2 -24029*n +18924) *a(n-3) -4*(n-1) *(n-2) *(n-3) *(2805*n -5578) *a(n-4) -2869*(n-1) *(n-2) *(n-3) *(n-4) *a(n-5)=0. - R. J. Mathar, Aug 17 2023

A345368 a(n) = Sum_{k=0..n} binomial(5*k,k) / (4*k + 1).

Original entry on oeis.org

1, 2, 7, 42, 327, 2857, 26608, 258488, 2588933, 26539288, 277082658, 2936050788, 31494394563, 341325970323, 3731742758203, 41108999917483, 455850863463768, 5084213586320193, 56997201842602368, 641906808539396253, 7258985455500009623, 82393287049581399283
Offset: 0

Views

Author

Ilya Gutkovskiy, Jul 28 2021

Keywords

Comments

Partial sums of A002294.

Crossrefs

Programs

  • Mathematica
    Table[Sum[Binomial[5 k, k]/(4 k + 1), {k, 0, n}], {n, 0, 21}]
    nmax = 21; A[] = 0; Do[A[x] = 1/(1 - x) + x (1 - x)^4 A[x]^5 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
  • PARI
    a(n) = sum(k=0, n, binomial(5*k, k)/(4*k+1)); \\ Michel Marcus, Jul 28 2021

Formula

G.f. A(x) satisfies: A(x) = 1 / (1 - x) + x * (1 - x)^4 * A(x)^5.
a(n) ~ 5^(5*n + 11/2) / (2869 * sqrt(Pi) * n^(3/2) * 2^(8*n + 7/2)). - Vaclav Kotesovec, Jul 28 2021
Showing 1-10 of 15 results. Next