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-8 of 8 results.

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

Original entry on oeis.org

1, 1, 1, 2, 5, 11, 25, 64, 169, 443, 1181, 3224, 8897, 24701, 69161, 195255, 554577, 1583109, 4541461, 13086574, 37856437, 109892403, 320034309, 934774902, 2737689189, 8037746691, 23652564261, 69749727716, 206091735797, 610061655665, 1808962146529
Offset: 0

Views

Author

Karol A. Penson, Jun 25 2013

Keywords

Crossrefs

Programs

  • Maple
    A226974 := proc(n)
        hypergeom([-n/3,-n/3+2/3,-n/3+1/3,1/4,1/2,3/4],[1/3,2/3,2/3,1,4/3],-256/27) ;
        simplify(%) ;
    end proc:
    seq(A226974(n),n=0..40) ; # R. J. Mathar, Jan 10 2023
  • Mathematica
    Table[Sum[Binomial[n,3*k]*Binomial[4*k,k]/(3*k+1), {k,0,Floor[n/3]}],{n,0,20}] (* Vaclav Kotesovec, Jun 28 2013 *)
  • Maxima
    a(n):=if n<0 then 0 else if n=0 then 1 else sum(sum(sum(a(l)*a(i)*a(j)*a(n-i-j-l-3),l,0,n-3-i-j),j,0,n-3-i),i,0,n-3)+1; /* Vladimir Kruchinin, May 17 2020 */
    
  • PARI
    a(n) = sum(k=0, n\3, binomial(n,3*k)*binomial(4*k,k)/(3*k+1)); \\ Michel Marcus, Sep 16 2021

Formula

a(n) = Sum_{k=0..floor(n/3)} binomial(n,3*k)*A002293(k).
Representation in terms of special values of hypergeometric function of type 6F5: a(n) = hypergeom([1/4, 1/2, 3/4, -(1/3)*n, -(1/3)*n+2/3, -(1/3)*n+1/3], [1/3, 2/3, 2/3, 1, 4/3],-4^4/3^3), n>=0.
Recurrence: 27*n*(n+1)*(n-1)*a(n) = 162*n*(n-1)^2*a(n-1) - 81*(5*n^2-15*n+12)*(n-1)*a(n-2) + 4*(199*n^3 - 1098*n^2 + 2043*n - 1296)*a(n-3) - (n-3)*(1173*n^2 - 5097*n + 5584)*a(n-4) + 6*(n-4)*(n-3)*(155*n-401)*a(n-5) - 283*(n-5)*(n-4)*(n-3)*a(n-6). - Vaclav Kotesovec, Jun 28 2013
a(n) ~ (3+4^(1+1/3))^(n+3/2)/(8*3^(n+1)*sqrt(2*Pi)*n^(3/2)). - Vaclav Kotesovec, Jun 28 2013
G.f. satisfies A(x)=1+x^3*A(x)^4+x/(1-x). - Vladimir Kruchinin, May 17 2020
From Peter Bala, Sep 15 2021: (Start)
O.g.f.: A(x) = (1/x)*series reversion( x*(1 - x^3)/(1 + x*(1 - x^3)) ).
The g.f. of the m-th binomial transform of this sequence is equal to (1/x)*series reversion( x*(1 - x^3)/(1 + (m + 1)*x*(1 - x^3)) ). The case m = -1 gives the sequence [1,0,0,1,0,0,4,0,0,22,0,0,140,...] - an aerated version of A002293. (End)

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

Original entry on oeis.org

1, 1, 1, 1, 2, 6, 16, 36, 76, 172, 436, 1156, 3006, 7606, 19202, 49466, 130156, 345356, 915196, 2421532, 6427001, 17163581, 46087911, 124133531, 334850208, 904691576, 2449891276, 6651540676, 18100561856, 49344295152, 134719523056, 368350942416, 1008680051756
Offset: 0

Views

Author

Vaclav Kotesovec, Jun 28 2013

Keywords

Comments

Generally, Sum(binomial(n,p*k)*binomial((p+1)*k,k)/(p*k+1), k=0..floor(n/p)) is asymptotic to (p+(p+1)^(1+1/p))^(n+3/2)/(p^(n+1)*(p+1)^(1+3/(2*p))*n^(3/2)*sqrt(2*Pi)).

Crossrefs

Cf. A002294, A007317 (p=1), A049130 (p=2), A226974 (p=3), A226910 (p=5).

Programs

  • Mathematica
    Table[Sum[Binomial[n,4*k]*Binomial[5*k,k]/(4*k+1),{k,0,Floor[n/4]}],{n,0,20}]
  • PARI
    a(n)=sum(k=0,n\4,binomial(n,4*k)*binomial(5*k,k)/(4*k+1)) \\ Charles R Greathouse IV, Jun 28 2013

Formula

Recurrence: -2869*(n-7)*(n-6)*(n-5)*(n-4)*a(n-8) + 2*(n-6)*(n-5)*(n-4)*(5226*n-17267)*a(n-7) - (n-5)*(n-4)*(11582*n^2-55156*n+50139)*a(n-6) - 3*(n-4)*(612*n^3 - 18926*n^2 + 102684*n - 155665)*a(n-5) + 5*(n-4)*(2959*n^3 - 26172*n^2 + 77408*n - 76800)*a(n-4) - 1024*(n-2)*(2*n-5)*(7*n^2-35*n+48)*a(n-3) + 1024*(n-2)*(n-1)*(7*n^2-28*n+30)*a(n-2) - 1024*(n-2)*(n-1)*n*(2*n-3)*a(n-1) + 256*(n-2)*(n-1)*n*(n+1)*a(n) = 0.
a(n) ~ (4+5^(1+1/4))^(n+3/2)/(4^(n+1)*5^(1+3/8)*n^(3/2)*sqrt(2*Pi)).
G.f. A(x) satisfies: A(x) = 1 / (1 - x) + x^4 * A(x)^5. - Ilya Gutkovskiy, Jul 25 2021
From Peter Bala, Sep 15 2021: (Start)
O.g.f.: A(x) = (1/x)*series reversion ( x*(1 - x^4)/(1 + x*(1 - x^4) )).
The g.f. of the m-th binomial transform of this sequence is equal to (1/x)*series reversion ( x*(1 - x^4)/(1 + (m + 1)*x*(1 - x^4)) ). The case m = -1 gives the sequence [1,0,0,0,1,0,0,0,5,0,0,0,35,0,0,0,285,...] - an aerated version of A002294. (End)

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

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 7, 22, 57, 127, 259, 529, 1189, 3004, 8009, 21073, 53233, 129813, 312733, 763573, 1915251, 4914736, 12720841, 32800186, 83869501, 213261712, 542609237, 1388542312, 3579043987, 9273567337, 24075321925, 62475528190, 161969731985, 419914766965
Offset: 0

Views

Author

Karol A. Penson, Jun 22 2013

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Sum[Binomial[n,5*k]*Binomial[6*k,k]/(5*k+1),{k,0,Floor[n/5]}],{n,0,20}] (* Vaclav Kotesovec, Jun 28 2013 *)
  • PARI
    a(n)=sum(k=0,n\5,binomial(n,5*k)*binomial(6*k,k)/(5*k+1)) \\ Charles R Greathouse IV, Jun 24 2013

Formula

Representation in terms of special values of generalized hypergeometric function of type 10F9: a(n) = hypergeom([1/6, 1/3, 1/2, 2/3, 5/6, -(1/5)*n, -(1/5)*n+4/5, -(1/5)*n+3/5, -(1/5)*n+2/5, 1/5-(1/5)*n], [1/5, 2/5, 2/5, 3/5, 3/5, 4/5, 4/5, 1, 6/5], -6^6/5^5), n>=0.
Recurrence: -49781*(n-9)*(n-8)*(n-7)*(n-6)*(n-5)*a(n-10) + 10*(n-8)*(n-7)*(n-6)*(n-5)*(26453*n - 123726)*a(n-9) - 15*(n-7)*(n-6)*(n-5)*(40479*n^2 - 351957*n + 782140)*a(n-8) + 120*(n-6)*(n-5)*(7013*n^3 - 87699*n^2 + 378278*n - 565577)*a(n-7) - 6*(n-5)*(148255*n^4 - 2435310*n^3 + 15491085*n^2 - 45173430*n + 50791476)*a(n-6) + 12*(69513*n^5 - 1361100*n^4 + 10838875*n^3 - 43818750*n^2 + 89776250*n - 74437500)*a(n-5) - 93750*(7*n^4 - 98*n^3 + 525*n^2 - 1274*n + 1180)*(n-3)*a(n-4) + 375000*(n-2)*(n^2-6*n+10)*(n-3)^2*a(n-3) - 46875*(n-2)*(n-1)*(3*n^2-15*n+20)*(n-3)*a(n-2) + 31250*(n-2)^2*(n-1)*n*(n-3)*a(n-1) - 3125*(n-2)*(n-1)*n*(n+1)*(n-3)*a(n) = 0. - Vaclav Kotesovec, Jun 28 2013
a(n) ~ (5+6^(1+1/5))^(n+3/2)/(5^(n+1)*6^(1+3/10)*sqrt(2*Pi)*n^(3/2)). - Vaclav Kotesovec, Jun 28 2013
G.f. A(x) satisfies: A(x) = 1 / (1 - x) + x^5 * A(x)^6. - Ilya Gutkovskiy, Jul 25 2021
From Peter Bala, Sep 15 2021: (Start)
O.g.f.: A(x) = (1/x)*series reversion ( x*(1 - x^5)/(1 + x*(1 - x^5)) ).
The g.f. of the m-th binomial transform of this sequence is equal to (1/x)*series reversion ( x*(1 - x^5)/(1 + (m + 1)*x*(1 - x^5)) ). The case m = -1 gives the sequence [1, 0, 0, 0, 0, 1, 0, 0,0, 0, 6, 0, 0, 0, 0, 51, 0, 0, 0, 0, 506, ...] - an aerated version of A002295. (End)

A364589 G.f. satisfies A(x) = 1/(1-x) + x^3*A(x)^3.

Original entry on oeis.org

1, 1, 1, 2, 4, 7, 14, 31, 67, 146, 331, 760, 1749, 4072, 9583, 22673, 53929, 129055, 310328, 749152, 1815481, 4415313, 10771564, 26352955, 64644926, 158963191, 391767016, 967523138, 2394060433, 5934576763, 14735792889, 36647185192, 91274339014, 227645446307
Offset: 0

Views

Author

Seiichi Manyama, Jul 29 2023

Keywords

Crossrefs

Programs

  • PARI
    a(n) = sum(k=0, n\3, binomial(n-k, 2*k)*binomial(3*k,k)/(2*k+1));

Formula

a(n) = Sum_{k=0..floor(n/3)} binomial(n-k,2*k) * binomial(3*k,k) / (2*k+1).

A364590 G.f. satisfies A(x) = 1/(1-x) + x^4*A(x)^3.

Original entry on oeis.org

1, 1, 1, 1, 2, 4, 7, 11, 19, 37, 74, 142, 268, 518, 1033, 2077, 4152, 8290, 16687, 33899, 69148, 141160, 288650, 592354, 1220086, 2519226, 5210164, 10794088, 22408556, 46613554, 97125751, 202662419, 423459427, 886048249, 1856448852, 3894362560, 8178530890
Offset: 0

Views

Author

Seiichi Manyama, Jul 29 2023

Keywords

Crossrefs

Programs

  • PARI
    a(n) = sum(k=0, n\4, binomial(n-2*k, 2*k)*binomial(3*k, k)/(2*k+1));

Formula

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

A370836 Expansion of (1/x) * Series_Reversion( x/(x+1/(1+x^2)) ).

Original entry on oeis.org

1, 1, 0, -2, -2, 6, 19, 1, -98, -170, 268, 1464, 967, -7253, -19035, 11497, 142894, 186814, -592148, -2327480, -371472, 14922592, 30367918, -44517534, -291059645, -242260229, 1550840094, 4611423196, -2050694753, -36095033685, -54276040088, 150373292998
Offset: 0

Views

Author

Seiichi Manyama, Mar 03 2024

Keywords

Crossrefs

Programs

  • PARI
    my(N=40, x='x+O('x^N)); Vec(serreverse(x/(x+1/(1+x^2)))/x)
    
  • PARI
    a(n) = sum(k=0, n\2, (-1)^k*binomial(n, 2*k)*binomial(3*k, k)/(2*k+1));

Formula

a(n) = Sum_{k=0..floor(n/2)} (-1)^k * binomial(n,2*k) * binomial(3*k,k)/(2*k+1).

A364588 G.f. satisfies A(x) = 1/(1-x) + x^2*A(x)^4.

Original entry on oeis.org

1, 1, 2, 5, 15, 49, 170, 613, 2275, 8629, 33301, 130333, 516077, 2063685, 8321892, 33803161, 138181521, 568031297, 2346668400, 9737766513, 40569611691, 169632827345, 711611670532, 2994165070045, 12632782541053, 53433933353885, 226540298098019
Offset: 0

Views

Author

Seiichi Manyama, Jul 29 2023

Keywords

Crossrefs

Programs

  • PARI
    a(n) = sum(k=0, n\2, binomial(n+k, 3*k)*binomial(4*k,k)/(3*k+1));

Formula

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

A369688 G.f. satisfies A(x) = 1 + x*A(x) + x^2*(1-x)^3*A(x)^5.

Original entry on oeis.org

1, 1, 2, 4, 12, 36, 126, 442, 1644, 6172, 23801, 92731, 366688, 1462852, 5891808, 23898576, 97600556, 400844140, 1654818768, 6862550360, 28576414261, 119434041561, 500849380048, 2106740001442, 8886482895068, 37580609774876, 159303913630686
Offset: 0

Views

Author

Seiichi Manyama, Jan 28 2024

Keywords

Crossrefs

Programs

  • PARI
    a(n) = sum(k=0, n\2, binomial(n, 2*k)*binomial(5*k, k)/(4*k+1));

Formula

a(n) = Sum_{k=0..floor(n/2)} binomial(n,2*k) * binomial(5*k,k) / (4*k+1).
Showing 1-8 of 8 results.