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

A323541 a(n) = Product_{k=0..n} (k^3 + (n-k)^3).

Original entry on oeis.org

0, 1, 128, 59049, 51380224, 80869140625, 207351578198016, 811509810302822449, 4603095542875667038208, 36344623587588604291790241, 386644580358400000000000000000, 5395532942025804980378907333844441, 96578621213529440721046520779140759552
Offset: 0

Views

Author

Vaclav Kotesovec, Jan 17 2019

Keywords

Crossrefs

Cf. 2*A000537 and A163102 (with sum instead of product).

Programs

  • Magma
    m:=3; [(&*[k^m + (n-k)^m: k in [0..n]]): n in [0..15]]; // G. C. Greubel, Jan 18 2019
    
  • Mathematica
    Table[Product[k^3+(n-k)^3, {k, 0, n}], {n, 0, 15}]
  • PARI
    m=3; vector(15, n, n--; prod(k=0,n, k^m + (n-k)^m)) \\ G. C. Greubel, Jan 18 2019
    
  • Sage
    m=3; [product(k^m +(n-k)^m for k in (0..n)) for n in (0..15)] # G. C. Greubel, Jan 18 2019

Formula

a(n) ~ exp(2*(Pi/(3*sqrt(3))-1)*n) * n^(3*n + 3).

A324426 a(n) = Product_{i=1..n, j=1..n} (i^3 + j^3).

Original entry on oeis.org

1, 2, 2592, 134425267200, 3120795915109442519040000, 180825857777547616919759624941086965760000000, 99356698720512072045648926659510730227553351200000695922065408000000000
Offset: 0

Views

Author

Vaclav Kotesovec, Feb 27 2019

Keywords

Crossrefs

Programs

  • Maple
    a:= n-> mul(mul(i^3+j^3, i=1..n), j=1..n):
    seq(a(n), n=0..7);  # Alois P. Heinz, Jun 24 2023
  • Mathematica
    Table[Product[i^3+j^3, {i, 1, n}, {j, 1, n}], {n, 1, 10}]
  • PARI
    a(n) = prod(i=1, n, prod(j=1, n, i^3+j^3)); \\ Michel Marcus, Feb 27 2019
    
  • Python
    from math import prod, factorial
    def A324426(n): return prod(i**3+j**3 for i in range(1,n) for j in range(i+1,n+1))**2*factorial(n)**3<Chai Wah Wu, Nov 26 2023

Formula

a(n) ~ A * 2^(2*n*(n+1) + 1/4) * exp(Pi*(n*(n+1) + 1/6)/sqrt(3) - 9*n^2/2 - 1/12) * n^(3*n^2 - 3/4) / (3^(5/6) * Pi^(1/6) * Gamma(2/3)^2), where A is the Glaisher-Kinkelin constant A074962.
a(n) = A079478(n) * A367543(n). - Vaclav Kotesovec, Nov 22 2023
For n>0, a(n)/a(n-1) = A272246(n)^2 / (2*n^9). - Vaclav Kotesovec, Dec 02 2023

Extensions

a(0)=1 prepended by Alois P. Heinz, Jun 24 2023

A272244 a(n) = Product_{k=0..n} (n^2 + k^2).

Original entry on oeis.org

0, 2, 160, 21060, 4352000, 1313845000, 547573478400, 301758856490000, 212663770808320000, 186659516597629140000, 199722414913149440000000, 255947740845844788169000000, 387074162712817024892928000000, 682170272459193898736228210000000
Offset: 0

Views

Author

Vaclav Kotesovec, Apr 23 2016

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Product[n^2+k^2,{k,0,n}],{n,0,15}]

Formula

a(n) ~ 2^(n + 1/2) * n^(2*(n+1)) / exp((4-Pi)*n/2).

A272247 a(n) = Product_{k=0..n} (n^4 + k^4).

Original entry on oeis.org

0, 2, 8704, 104372388, 3087748038656, 194985808028125000, 23467500289618753093632, 4938279594477505466022892304, 1699491802948673617630927695904768, 907214902722906584628050661111614570016, 719684491044699824651608981274624000000000000
Offset: 0

Views

Author

Vaclav Kotesovec, Apr 23 2016

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Product[n^4+k^4,{k,0,n}],{n,0,10}]

Formula

a(n) ~ 2^(n+1/2) * (1+sqrt(2))^(sqrt(2)*n) * n^(4*n + 4) / exp((4 - Pi/sqrt(2))*n).

A272248 a(n) = Product_{k=0..n} (n^5 + k^5).

Original entry on oeis.org

0, 2, 67584, 7924375800, 2876035930521600, 2693451205324687500000, 5648896640332217707816550400, 23819277009290664033936067933468800, 185754160490281505676324140907107450880000, 2507604631016507710974687639612411760216253760000
Offset: 0

Views

Author

Vaclav Kotesovec, Apr 23 2016

Keywords

Comments

In general, for p>=1, Product_{k=0..n} (n^p + k^p) ~ sqrt(2) * n^(p*(n+1)) * exp(n*Sum_{j>=1} (-1)^(j+1) / (j*(1 + j*p))).

Crossrefs

Programs

  • Mathematica
    Table[Product[n^5+k^5,{k,0,n}],{n,0,10}]

Formula

a(n) ~ 2^(2*n+1/2) * phi^(sqrt(5)*n) * n^(5*n+5) / exp((5-sqrt(phi)*Pi/5^(1/4))*n), where phi = A001622 = (1+sqrt(5))/2 is the golden ratio.

A367823 a(n) = Product_{k=0..n} (n^6 + k^6).

Original entry on oeis.org

0, 2, 532480, 615291760980, 2759348178649088000, 38607000141072009765625000, 1421393314362014000255258433945600, 120987688662917113214734283836109985530000, 21538956798613753302336535534471908360548515840000, 7407676630843157974199152503244776880017496531101658740000
Offset: 0

Views

Author

Vaclav Kotesovec, Dec 02 2023

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Product[n^6 + k^6, {k, 0, n}], {n, 0, 10}]
  • PARI
    a(n) = prod(k=0, n, n^6 + k^6); \\ Michel Marcus, Dec 02 2023

Formula

a(n) ~ 2^(n + 1/2) * (2 + sqrt(3))^(sqrt(3)*n) * n^(6*n + 6) / exp((6 - Pi)*n).

A367833 a(n) = Product_{k=0..n} (n^8 + k^8).

Original entry on oeis.org

0, 2, 33685504, 3851231199376068, 2670340833224951538384896, 8442437541556000328575439453125000, 96982590530284083757788173242580213780447232, 3406021755102121469537208768079471859185253483110025744
Offset: 0

Views

Author

Vaclav Kotesovec, Dec 02 2023

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Product[n^8 + k^8, {k, 0, n}], {n, 0, 10}]

Formula

a(n) ~ 2^(n + 1/2) * ((2 + sqrt(2 + sqrt(2)))/(2 - sqrt(2 + sqrt(2))))^(sqrt(2 + sqrt(2))*n/2) * (sqrt(2) - 1 + sqrt(4 - 2*sqrt(2)))^(sqrt(2 - sqrt(2))*n) * exp((Pi/sqrt(2 - sqrt(2)) - 8)*n) * n^(8*(n + 1)).
Showing 1-7 of 7 results.