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

A013662 Decimal expansion of zeta(4).

Original entry on oeis.org

1, 0, 8, 2, 3, 2, 3, 2, 3, 3, 7, 1, 1, 1, 3, 8, 1, 9, 1, 5, 1, 6, 0, 0, 3, 6, 9, 6, 5, 4, 1, 1, 6, 7, 9, 0, 2, 7, 7, 4, 7, 5, 0, 9, 5, 1, 9, 1, 8, 7, 2, 6, 9, 0, 7, 6, 8, 2, 9, 7, 6, 2, 1, 5, 4, 4, 4, 1, 2, 0, 6, 1, 6, 1, 8, 6, 9, 6, 8, 8, 4, 6, 5, 5, 6, 9, 0, 9, 6, 3, 5, 9, 4, 1, 6, 9, 9, 9, 1
Offset: 1

Views

Author

Keywords

Examples

			1.082323233711138191516003696541167...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 811.
  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 89, Exercise.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 262.
  • Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004, Section F17, Series associated with the zeta-function, p. 391.
  • L. D. Landau and E. M. Lifschitz, Band V, Statistische Physik, Akademie Verlag, 1966, pp. 172 and 180-181.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 162.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987, p. 33.

Crossrefs

Programs

  • Magma
    SetDefaultRealField(RealField(110)); L:=RiemannZeta(); Evaluate(L,4); // G. C. Greubel, May 30 2019
    
  • Maple
    evalf(Pi^4/90,120); # Muniru A Asiru, Sep 19 2018
  • Mathematica
    RealDigits[Zeta[4],10,120][[1]] (* Harvey P. Dale, Dec 18 2012 *)
  • Maxima
    ev(zeta(4),numer) ; /* R. J. Mathar, Feb 27 2012 */
    
  • PARI
    default(realprecision, 20080); x=Pi^4/90; for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b013662.txt", n, " ", d)); \\ Harry J. Smith, Apr 29 2009
    
  • Sage
    numerical_approx(zeta(4), digits=100) # G. C. Greubel, May 30 2019

Formula

zeta(4) = Pi^4/90 = A092425/90. - Harry J. Smith, Apr 29 2009
From Peter Bala, Dec 03 2013: (Start)
Definition: zeta(4) := Sum_{n >= 1} 1/n^4.
zeta(4) = (4/17)*Sum_{n >= 1} ( (1 + 1/2 + ... + 1/n)/n )^2 and
zeta(4) = (16/45)*Sum_{n >= 1} ( (1 + 1/3 + ... + 1/(2*n-1))/n )^2 (see Borwein and Borwein).
zeta(4) = (256/90)*Sum_{n >= 1} n^2*(4*n^2 + 3)*(12*n^2 + 1)/(4*n^2 - 1)^5.
Series acceleration formulas:
zeta(4) = (36/17)*Sum_{n >= 1} 1/( n^4*binomial(2*n,n) ) (Comtet)
= (36/17)*Sum_{n >= 1} P(n)/( (2*n*(2*n - 1))^4*binomial(4*n,2*n) )
= (36/17)*Sum_{n >= 1} Q(n)/( (3*n*(3*n - 1)*(3*n - 2))^4*binomial(6*n,3*n) ),
where P(n) = 80*n^4 - 48*n^3 + 24*n^2 - 8*n + 1 and Q(n) = 137781*n^8 - 275562*n^7 + 240570*n^6 - 122472*n^5 + 41877*n^4 - 10908*n^3 + 2232*n^2 - 288*n + 16 (see section 8 in the Bala link). (End)
zeta(4) = 2/3*2^4/(2^4 - 1)*( Sum_{n even} n^2*p(n)/(n^2 - 1)^5 ), where p(n) = 3*n^4 + 10*n^2 + 3 is a row polynomial of A091043. See A013664, A013666, A013668 and A013670. - Peter Bala, Dec 05 2013
zeta(4) = Sum_{n >= 1} ((floor(sqrt(n))-floor(sqrt(n-1)))/n^2). - Mikael Aaltonen, Jan 18 2015
zeta(4) = Product_{k>=1} 1/(1 - 1/prime(k)^4). - Vaclav Kotesovec, May 02 2020
From Wolfdieter Lang, Sep 16 2020: (Start)
zeta(4) = (1/3!)*Integral_{x=0..oo} x^3/(exp(x) - 1) dx. See Abramowitz-Stegun, 23.2.7., for s=2, p. 807, and Landau-Lifschitz, Band V, p. 172, eq. (4), for x=4. See also A231535.
zeta(4) = (4/21)*Integral_{x=0..oo} x^3/(exp(x) + 1) dx. See Abramowitz-Stegun, 23.2.8., for s=2, p. 807, and Landau-Lifschitz, Band V, p. 172, eq. (1), for x=4. See also A337711. (End)
zeta(4) = (72/17) * Integral_{x=0..Pi/3} x*(log(2*sin(x/2)))^2. See Richard K. Guy reference. - Bernard Schott, Jul 20 2022
From Peter Bala, Nov 12 2023: (Start)
zeta(4) = 1 + (4/3)*Sum_{k >= 1} (1 - 2*(-1)^k)/(k*(k + 1)^4*(k + 2)) = 35053/32400 + 48*Sum_{k >= 1} (1 - 2*(-1)^k)/(k*(k + 1)*(k + 2)*(k + 3)^4*(k + 4)*(k + 5)*(k + 6)).
More generally, it appears that for n >= 0, zeta(4) = c(n) + (4/3)*(2*n + 1)!^2 * Sum_{k >= 1} (1 - 2*(-1)^k)/( (k + 2*n + 1)^3*Product_{i = 0..4*n+2} (k + i) ), where {c(n) : n >= 0} is a sequence of rational approximations to zeta(4) beginning [1, 35053/32400, 2061943067/ 1905120000, 18594731931460103/ 17180389306080000, 257946156103293544441/ 238326360453941760000, ...]. (End)
From Peter Bala, Apr 27 2025: (Start)
zeta(4) = 1/4! * Integral_{x >= 0} x^4 * exp(x)/(exp(x) - 1)^2 dx = 8/7 * 1/4! * Integral_{x >= 0} x^4 * exp(x)/(exp(x) + 1)^2 dx.
zeta(4) = 1/5! * Integral_{x >= 0} x^5 * exp(x)*(exp(x) + 1)/(exp(x) - 1)^3 dx = 1/(3*5*7) * Integral_{x >= 0} x^5 * exp(x)*(exp(x) - 1)/(exp(x) + 1)^3 dx. (End)
10*zeta(4) = Sum_{k>=1} H(k)^3/(k*(k+1)), where H(k) = A001008(k)/A002805(k) is the k-th harmonic number (Ramachandra, 1981). - Amiram Eldar, May 30 2025
zeta(4) = Integral_{x=0..1} Li(3,x)/x dx, where Li(n,x) is the polylogarithm function. - Kritsada Moomuang, Jun 14 2025
zeta(4) = Sum_{i, j >= 1} 1/(i^3*j*binomial(i+j, i)) = 4/3 * Sum_{i, j >= 1} 1/(i^2*j^2*binomial(i+j, i)). - Peter Bala, Aug 03 2025

A013664 Decimal expansion of zeta(6).

Original entry on oeis.org

1, 0, 1, 7, 3, 4, 3, 0, 6, 1, 9, 8, 4, 4, 4, 9, 1, 3, 9, 7, 1, 4, 5, 1, 7, 9, 2, 9, 7, 9, 0, 9, 2, 0, 5, 2, 7, 9, 0, 1, 8, 1, 7, 4, 9, 0, 0, 3, 2, 8, 5, 3, 5, 6, 1, 8, 4, 2, 4, 0, 8, 6, 6, 4, 0, 0, 4, 3, 3, 2, 1, 8, 2, 9, 0, 1, 9, 5, 7, 8, 9, 7, 8, 8, 2, 7, 7, 3, 9, 7, 7, 9, 3, 8, 5, 3, 5, 1, 7
Offset: 1

Views

Author

Keywords

Examples

			1.01734306198444913...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 811.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 262.

Crossrefs

Programs

Formula

Equals Pi^6/945 = A092732/945. - Mohammad K. Azarian, Mar 03 2008
zeta(6) = 8/3*2^6/(2^6 - 1)*( Sum_{n even} n^2*p(n)/(n^2 - 1)^7 ), where p(n) = n^6 + 7*n^4 + 7*n^2 + 1 is a row polynomial of A091043. See A013662, A013666, A013668 and A013670. - Peter Bala, Dec 05 2013
Definition: zeta(6) = Sum_{n >= 1} 1/n^6. - Bruno Berselli, Dec 05 2013
zeta(6) = Sum_{n >= 1} (A010052(n)/n^3). - Mikael Aaltonen, Feb 20 2015
zeta(6) = Sum_{n >= 1} (A010057(n)/n^2). - A.H.M. Smeets, Sep 19 2018
zeta(6) = Product_{k>=1} 1/(1 - 1/prime(k)^6). - Vaclav Kotesovec, May 02 2020
From Wolfdieter Lang, Sep 16 2020: (Start)
zeta(6) = (1/5!)*Integral_{x=0..infinity} x^5/(exp(x) - 1) dx. See Abramowitz-Stegun, 23.2.7., for s=6, p. 807. See also A337710 for the value of the integral.
zeta(6) = (4/465)*Integral_{x=0..infinity} x^5/(exp(x) + 1) dx. See Abramowitz-Stegun, 23.2.8., for s=6, p. 807. The value of the integral is (31/252)*Pi^6 = 118.2661309... . (End)
From Peter Bala, Apr 27 2025: (Start)
zeta(6) = 1/6! * Integral_{x >= 0} x^6 * exp(x)/(exp(x) - 1)^2 dx = 2^5/(2^5 - 1) * 1/6! * Integral_{x >= 0} x^6 * exp(x)/(exp(x) + 1)^2 dx.
zeta(6) = 1/7! * Integral_{x >= 0} x^7 * exp(x)*(exp(x) + 1) /(exp(x) - 1)^3 dx = 2/(3*7*15*31) * Integral_{x >= 0} x^7 * exp(x)*(exp(x) - 1)/(exp(x) + 1)^3 dx. (End)

A013666 Decimal expansion of zeta(8).

Original entry on oeis.org

1, 0, 0, 4, 0, 7, 7, 3, 5, 6, 1, 9, 7, 9, 4, 4, 3, 3, 9, 3, 7, 8, 6, 8, 5, 2, 3, 8, 5, 0, 8, 6, 5, 2, 4, 6, 5, 2, 5, 8, 9, 6, 0, 7, 9, 0, 6, 4, 9, 8, 5, 0, 0, 2, 0, 3, 2, 9, 1, 1, 0, 2, 0, 2, 6, 5, 2, 5, 8, 2, 9, 5, 2, 5, 7, 4, 7, 4, 8, 8, 1, 4, 3, 9, 5, 2, 8, 7, 2, 3, 0, 3, 7, 2, 3, 7, 1, 9, 7
Offset: 1

Views

Author

Keywords

Comments

This sequence is also the decimal expansion of Pi^8/9450. - Mohammad K. Azarian, Mar 03 2008

Examples

			1.00407735619794433937868523850865246525896079064985002032911020265...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 811.

Crossrefs

Programs

  • Maple
    Digits := 100 : evalf(Pi^8/9450) ; # R. J. Mathar, Jan 07 2021
  • Mathematica
    RealDigits[Zeta[8], 10, 100][[1]] (* Vincenzo Librandi, Feb 15 2015 *)

Formula

zeta(8) = 2/3*2^8/(2^8 - 1)*( Sum_{n even} n^2*p(n)/(n^2 - 1)^9 ), where p(n) = 5*n^8 + 60*n^6 + 126*n^4 + 60*n^2 + 5 is a row polynomial of A091043. See A013662, A013664, A013668 and A013670. - Peter Bala, Dec 05 2013
zeta(8) = Sum_{n >= 1} (A010052(n)/n^4). - Mikael Aaltonen, Feb 20 2015
zeta(8) = Product_{k>=1} 1/(1 - 1/prime(k)^8). - Vaclav Kotesovec, May 02 2020
From Wolfdieter Lang, Sep 16 2020 (Start):
zeta(8) = (1/7!)*Integral_{0..infinity} x^7/(exp(x) - 1) dx. See Abramowitz-Stegun, 23.2.7., for s=8, p. 807. The value of the integral is 8*Pi^8/15 = 5060.54987... .
zeta(8) = (2^7/(127*7!))*Integral_{0..infinity} x^7/(exp(x) + 1) dx. See Abramowitz-Stegun, 23.2.8., for s=8, p. 807. The prefactor is 8/40005. The value of the integral is (127/240)*Pi^8 = 5021.014329... .(End)
Equals A092736/9450. - R. J. Mathar, Jan 07 2021
From Peter Bala, Apr 27 2025: (Start)
zeta(8) = 1/8! * Integral_{x >= 0} x^8 * exp(x)/(exp(x) - 1)^2 dx = 2^7/(2^7 - 1) * 1/8! * Integral_{x >= 0} x^8 * exp(x)/(exp(x) + 1)^2 dx.
zeta(8) = 1/9! * Integral_{x >= 0} x^9 * exp(x)*(exp(x) + 1) /(exp(x) - 1)^3 dx = 1/(3*15*63*127) * Integral_{x >= 0} x^9 * exp(x)*(exp(x) - 1)/(exp(x) + 1)^3 dx. (End)

A013668 Decimal expansion of zeta(10).

Original entry on oeis.org

1, 0, 0, 0, 9, 9, 4, 5, 7, 5, 1, 2, 7, 8, 1, 8, 0, 8, 5, 3, 3, 7, 1, 4, 5, 9, 5, 8, 9, 0, 0, 3, 1, 9, 0, 1, 7, 0, 0, 6, 0, 1, 9, 5, 3, 1, 5, 6, 4, 4, 7, 7, 5, 1, 7, 2, 5, 7, 7, 8, 8, 9, 9, 4, 6, 3, 6, 2, 9, 1, 4, 6, 5, 1, 5, 1, 9, 1, 2, 9, 5, 4, 3, 9, 7, 0, 4, 1, 9, 6, 8, 6, 1, 0, 3, 8, 5, 6, 5
Offset: 1

Views

Author

Keywords

Examples

			1.0009945751278180853371459589003190170060195315644775172577889946362914...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 811.

Crossrefs

Programs

Formula

Equals Pi^10/93555.
zeta(10) = 4/3*2^10/(2^10 - 1)*( Sum_{n even} n^2*p(n)/(n^2 - 1)^11 ), where p(n) = 3*n^10 + 55*n^8 + 198*n^6 + 198*n^4 + 55*n^2 + 3 is a row polynomial of A091043. - Peter Bala, Dec 05 2013
zeta(10) = Sum_{n >= 1} (A010052(n)/n^5) = Sum_{n >= 1} ( (floor(sqrt(n))-floor(sqrt(n-1)))/n^5 ). - Mikael Aaltonen, Feb 20 2015
zeta(10) = Product_{k>=1} 1/(1 - 1/prime(k)^10). - Vaclav Kotesovec, May 02 2020
From Wolfdieter Lang, Sep 16 2020: (Start)
zeta(10) = (1/9!)*Integral_{0..infinity} x^9/(exp(x) - 1). See Abramowitz-Stegun, 23.2.7., for s=10, p. 807. The value of the integral is (128/33)*Pi^10 = (3.6324091...)*10^5.
zeta(10) = (4/1448685)*Integral_{0..infinity} x^9/(exp(x) + 1). See Abramowitz-Stegun, 23.2.8., for s=10, p. 807. The value of the integral is (511/132)*Pi^10 = (3.625314565...)*10^5. (End)

A013670 Decimal expansion of zeta(12).

Original entry on oeis.org

1, 0, 0, 0, 2, 4, 6, 0, 8, 6, 5, 5, 3, 3, 0, 8, 0, 4, 8, 2, 9, 8, 6, 3, 7, 9, 9, 8, 0, 4, 7, 7, 3, 9, 6, 7, 0, 9, 6, 0, 4, 1, 6, 0, 8, 8, 4, 5, 8, 0, 0, 3, 4, 0, 4, 5, 3, 3, 0, 4, 0, 9, 5, 2, 1, 3, 3, 2, 5, 2, 0, 1, 9, 6, 8, 1, 9, 4, 0, 9, 1, 3, 0, 4, 9, 0, 4, 2, 8, 0, 8, 5, 5, 1, 9, 0, 0, 6, 9
Offset: 1

Views

Author

Keywords

Examples

			1.0002460865533080482986379980477396709604160884580034045330409521332520...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 811.

Crossrefs

Programs

Formula

zeta(12) = 2/3*2^12/(2^12 - 1)*( Sum_{n even} n^2*p(n)/(n^2 - 1)^13 ), where p(n) = 7*n^12 + 182*n^10 + 1001*n^8 + 1716*n^6 + 1001*n^4 + 182*n^2 + 7 is a row polynomial of A091043. - Peter Bala, Dec 05 2013
zeta(12) = Sum_{n >= 1} (A010052(n)/n^6) = Sum {n >= 1} ( (floor(sqrt(n))-floor(sqrt(n-1)))/n^6 ). - Mikael Aaltonen, Feb 20 2015
zeta(12) = 691/638512875*Pi^12 (see A002432). - Rick L. Shepherd, May 30 2016
zeta(12) = Product_{k>=1} 1/(1 - 1/prime(k)^12). - Vaclav Kotesovec, May 02 2020

A034870 Even-numbered rows of Pascal's triangle.

Original entry on oeis.org

1, 1, 2, 1, 1, 4, 6, 4, 1, 1, 6, 15, 20, 15, 6, 1, 1, 8, 28, 56, 70, 56, 28, 8, 1, 1, 10, 45, 120, 210, 252, 210, 120, 45, 10, 1, 1, 12, 66, 220, 495, 792, 924, 792, 495, 220, 66, 12, 1, 1, 14, 91, 364, 1001, 2002, 3003, 3432, 3003, 2002, 1001, 364, 91, 14, 1
Offset: 0

Views

Author

Keywords

Comments

The sequence of row lengths of this array is [1,3,5,7,9,11,13,...]= A005408(n), n>=0.
Equals X^n * [1,0,0,0,...] where X = an infinite tridiagonal matrix with (1,1,1,...) in the main and subsubdiagonal and (2,2,2,...) in the main diagonal. X also = a triangular matrix with (1,2,1,0,0,0,...) in each column. - Gary W. Adamson, May 26 2008
a(n,m) has the following interesting combinatoric interpretation. Let s(n,m) equal the set of all base-4, n-digit numbers with n-m more 1-digits than 2-digits. For example s(2,1) = {10,01,13,31} (note that numbers like 1 are left-padded with 0's to ensure that they have 2 digits). Notice that #s(2,1) = a(2,1) with # indicating cardinality. This is true in general. a(n,m)=#s(n,m). In words, a(n,m) gives the number of n-digit, base-4 numbers with n-m more 1 digits than 2 digits. A proof is provided in the Links section. - Russell Jay Hendel, Jun 23 2015

Examples

			Triangle begins:
  1;
  1,  2,  1;
  1,  4,  6,   4,   1;
  1,  6, 15,  20,  15,   6,   1;
  1,  8, 28,  56,  70,  56,  28,   8,   1;
  1, 10, 45, 120, 210, 252, 210, 120,  45,  10,  1;
  1, 12, 66, 220, 495, 792, 924, 792, 495, 220, 66, 12, 1;
		

Crossrefs

Cf. A000302 (row sums, powers of 4), alternating row sums are 0, except for n=0 which gives 1.

Programs

  • Haskell
    a034870 n k = a034870_tabf !! n !! k
    a034870_row n = a034870_tabf !! n
    a034870_tabf = map a007318_row [0, 2 ..]
    -- Reinhard Zumkeller, Apr 19 2012, Apr 02 2011
    
  • Magma
    /* As triangle: */ [[Binomial(n,k): k in [0..n]]: n in [0.. 15 by 2]]; // Vincenzo Librandi, Jul 16 2015
    
  • Maple
    T := (n,k) -> simplify(GegenbauerC(`if`(kPeter Luschny, May 08 2016
  • Mathematica
    Flatten[Table[Binomial[n,k],{n,0,20,2},{k,0,n}]] (* Harvey P. Dale, Dec 15 2014 *)
  • Maxima
    taylor(1/(1-x*(y+1)^2),x,0,10,y,0,10); /* Vladimir Kruchinin, Nov 22 2020 */
    
  • Sage
    flatten([[binomial(2*n, k) for k in (0..2*n)] for n in (0..12)]) # G. C. Greubel, Mar 18 2022

Formula

T(n, m) = binomial(2*n, m), 0<= m <= 2*n, 0<=n, else 0.
G.f. for column m=2*k sequence: (x^k)*Pe(k, x)/(1-x)^(2*k+1), k>=0; for column m=2*k-1 sequence (x^k)*Po(k, x)/(1-x)^(2*k), k>=1, with the row polynomials Pe(k, x) := sum(A091042(k, m)*x^m, m=0..k) and Po(k, x) := 2*sum(A091044(k, m)*x^m, m=0..k-1); see also triangle A091043.
From Paul D. Hanna, Apr 18 2012: (Start)
Let A(x) be the g.f. of the flattened sequence, then:
G.f.: A(x) = Sum_{n>=0} x^(n^2) * (1+x)^(2*n).
G.f.: A(x) = Sum_{n>=0} x^n*(1+x)^(2*n) * Product_{k=1..n} (1 - (1+x)^2*x^(4*k-3)) / (1 - (1+x)^2*x^(4*k-1)).
G.f.: A(x) = 1/(1 - x*(1+x)^2/(1 + x*(1-x^2)*(1+x)^2/(1 - x^5*(1+x)^2/(1 + x^3*(1-x^4)*(1+x)^2/(1 - x^9*(1+x)^2/(1 + x^5*(1-x^6)*(1+x)^2/(1 - x^13*(1+x)^2/(1 + x^7*(1-x^8)*(1+x)^2/(1 - ...))))))))), a continued fraction.
(End)
From Peter Bala, Jul 14 2015: (Start)
Denote this array by P. Then P * transpose(P) is the square array ( binomial(2*n + 2*k, 2*k) )n,k>=0, which, read by antidiagonals, is A086645.
Transpose(P) is a generalized Riordan array (1, (1 + x)^2) as defined in the Bala link.
Let p(x) = (1 + x)^2. P^2 gives the coefficients in the expansion of the polynomials ( p(p(x)) )^n, P^3 gives the coefficients in the expansion of the polynomials ( p(p(p(x))) )^n and so on.
Row sums are 2^(2*n); row sums of P^2 are 5^(2*n), row sums of P^3 are 26^(2*n). In general, the row sums of P^k, k = 0,1,2,..., are equal to A003095(k)^(2*n).
The signed version of this array ( (-1)^k*binomial(2*n,k) )n,k>=0 is a left-inverse for A034839.
A034839 * P = A080928. (End)
T(n, k) = GegenbauerC(m, -n, -1) where m = k if kPeter Luschny, May 08 2016
G.f.: 1/(1-x*(y+1)^2). - Vladimir Kruchinin, Nov 22 2020

A091044 One half of odd-numbered entries of even-numbered rows of Pascal's triangle A007318.

Original entry on oeis.org

1, 2, 2, 3, 10, 3, 4, 28, 28, 4, 5, 60, 126, 60, 5, 6, 110, 396, 396, 110, 6, 7, 182, 1001, 1716, 1001, 182, 7, 8, 280, 2184, 5720, 5720, 2184, 280, 8, 9, 408, 4284, 15912, 24310, 15912, 4284, 408, 9, 10, 570, 7752, 38760, 83980, 83980, 38760, 7752, 570, 10, 11
Offset: 1

Views

Author

Wolfdieter Lang, Jan 23 2004

Keywords

Comments

The odd-numbered columns of this triangle can be reduced: see triangle A091043.
The odd-numbered rows coincide with the ones of the reduced triangle A091043.
binomial(2*n,2*m+1) is even for n >= m + 1 >= 1, hence every T(n,m) is a positive integer.
The GCD (greatest common divisor) of the entries of each odd-numbered row n=2*k+1, k>=0, is 1.
The GCD of the entries of the even-numbered row n=2*k, k>=1, is A006519(n) (highest power of 2 in n=2*k).

Examples

			Triangle begins:
  [1];
  [2,2];
  [3,10,3];
  [4,28,28,4];
  [5,60,126,60,5];
  [6,110,396,396,110,6];
  ...
n = 6 = 2*3: gcd(6,110,396) = 2 = A006519(6);
n = 5: gcd(5,60,126) = 1 = A006519(5).
		

Crossrefs

Programs

  • Mathematica
    Flatten[Table[Binomial[2n,2m+1]/2,{n,1,11},{m,0,n-1}]] (* Indranil Ghosh, Feb 22 2017 *)
  • PARI
    {A(i, j) = binomial(2*i + 2*j - 2, 2*i - 1) / 2}; /* Michael Somos, Oct 15 2017 */

Formula

T(n, m)= binomial(2*n, 2*m+1)/2, n >= m + 1 >= 1, else 0.
Put a(n) = n!*(n+1/2)!/(1/2)!. T(n+1,k) = (n+1)*a(n)/(a(k)*a(n-k)).
T(n-1,k-1)*T(n,k+1)*T(n+1,k) = T(n-1,k)*T(n,k-1)*T(n+1,k+1). Cf. A111910. - Peter Bala, Oct 13 2011
From Peter Bala, Jul 29 2013: (Start)
O.g.f.: 1/(1 - 2*t*(x + 1) + t^2*(x - 1)^2)= 1 + (2 + 2*x)*t + (3 + 10*x + 3*x^2)*t^2 + ....
The n-th row polynomial R(n,x) = 1/(4*sqrt(x))*( (1 + sqrt(x))^(2*n) - (sqrt(x) - 1)^(2*n) ) and has n-1 real zeros given by the formula -cot^2(k*Pi/(2*n)) for k = 1,2,...,n-1. Cf A091042.
The row polynomial R(n,x) satisfies (x - 1)^n*R(n,x/(x - 1)) = U(n,2*x - 1), the n-th row polynomial of A053124.
Row sums A000302. Sum {k = 0..n-1} 2^k*T(n,k) = A001109(n). (End)
From Werner Schulte, Jan 13 2017: (Start)
(1) T(n,m) = T(n-1,m) + T(n-1,m-1)*(2*n-1-m)/m for 0 < m < n-1 with T(n,0) = n and T(n,n) = 0;
(2) T(n,m) = 2*T(n-1,m) + 2*T(n-1,m-1) - T(n-2,m) + 2*T(n-2,m-1) - T(n-2,m-2) for 0 < m < n-1 with T(n,0) = T(n,n-1) = n and T(n,m) = 0 if m < 0 or m >= n;
(3) The row polynomials p(n,x) = Sum_{m=0..n-1} T(n,m)*x^m satisfy the recurrence equation p(n+2,x) = (2+2*x)*p(n+1,x) - (x-1)^2*p(n,x) for n >= 1 with initial values p(1,x) = 1 and p(2,x) = 2+2*x.
(End)
G.f.: x*y /(1 - 2*(x+y) + (x-y)^2) with the entries regarded as an infinite square array A(i, j) read by antidiagonals. - Michael Somos, Oct 15 2017

A059978 a(n) = binomial(n+2,n)^6.

Original entry on oeis.org

1, 729, 46656, 1000000, 11390625, 85766121, 481890304, 2176782336, 8303765625, 27680640625, 82653950016, 225199600704, 567869252041, 1340095640625, 2985984000000, 6327518887936, 12827693806929, 25002110044521, 47045881000000, 85766121000000, 151939915084881
Offset: 0

Views

Author

Robert G. Wilson v, Mar 06 2001

Keywords

Comments

Number of 6-dimensional cage assemblies.

References

  • Clifford A. Pickover, "Wonders of Numbers, Adventures in Mathematics, Mind and Meaning," Oxford University Press, 2001, p. 325.

Crossrefs

Programs

  • Maple
    with (combinat):seq(mul(stirling2(n+1,n),k=1..6),n=1..18); # Zerinvary Lajos, Dec 14 2007
  • Mathematica
    m = 6; Table[n^m (n + 1)^m/2^m, {n, 1, 24}]

Formula

G.f.: (x^10 + 716*x^9 + 37257*x^8 + 450048*x^7 + 1822014*x^6 +2864328*x^5 + 1822014*x^4 + 450048*x^3 + 37257*x^2 + 716*x + 1)/(1-x)^13. - Colin Barker, Jul 09 2012
G.f.: 6F5([3,3,3,3,3,3], [1,1,1,1,1], z). - Benedict W. J. Irwin, Mar 14 2016
a(n) = (1/16)*( 3*S(7,n+1) + 10*S(9,n+1) + 3*S(11,n+1) ), where S(r,n) = Sum_{k = 1..n} k^r. Cf. A059977 and A059980. - Peter Bala, Jul 02 2019
From Amiram Eldar, May 15 2022: (Start)
Sum_{n>=0} 1/a(n) = 2688*Pi^2 + 448*Pi^4/15 + 128*Pi^6/945 - 29568.
Sum_{n>=0} (-1)^n/a(n) = 29568 - 32256*log(2) - 5376*zeta(3) - 720*zeta(5). (End)

Extensions

Better definition from Zerinvary Lajos, May 23 2006

A281123 Triangle T read by rows: n-th row (n>=0) gives the non-vanishing coefficients of the polynomial q(n,x) = 2^(-n)*((x+1)^(2^n) - (x-1)^(2^n))/2.

Original entry on oeis.org

1, 1, 1, 1, 1, 7, 7, 1, 1, 35, 273, 715, 715, 273, 35, 1, 1, 155, 6293, 105183, 876525, 4032015, 10855425, 17678835, 17678835, 10855425, 4032015, 876525, 105183, 6293, 155, 1, 1, 651, 119133, 9706503, 430321633, 11618684091, 205263418941, 2492484372855, 21552658988805, 136248095712855
Offset: 0

Views

Author

Martin Renner, Jan 15 2017

Keywords

Comments

Row n has length 1 for n = 1 and 2^(n-1) = A000079(n-1) for n >= 1.
The triangle T gives the non-vanishing coefficients of the polynomial q(0,x) = 1 and q(n,x) = 2^(-n)*Sum_{k=0..2^(n-1)-1} A281122(n,k) * x^(2^n-1-2*k), n >= 1.
The polynomial q(n,x) = product_{k=0..n-1} p(k,x) with polynomial p(n,x) = ((x+1)^(2^n) + (x-1)^(2^n))/2, whose coefficients are tabulated in A201461.
The algorithm r(n) = 1/2*(r(n-1) + A/r(n-1)), starting with r(0) = A, used for approximating sqrt(A), which is known as the Babylonian method or Hero's method after the first-century Greek mathematician Hero of Alexandria and which can be derived from Newton's method, generates fractions beginning with (A+1)/2, (A^2 + 6*A + 1)/(4*(A + 1)), (A^4 + 28*A^3 + 70*A^2 + 28*A + 1)/(8*(A^3 + 7*A^2 + 7*A + 1)), ... This is sqrt(A)*p(n,sqrt(A))/(2^n*q(n,sqrt(A))) with the given polynomials p(n,x) and q(n,x).

Examples

			The triangle T(n, k) starts with:
  1
  1
  1, 1
  1, 7, 7, 1
  1, 35, 273, 715, 715, 273, 35, 1
  1, 155, 6293, 105183, 876525, 4032015, 10855425, 17678835, 17678835, 10855425, 4032015, 876525, 105183, 6293, 155, 1
etc., since the first few polynomials are
q(0,x) = 1,
q(1,x) = x,
q(2,x) = x^3 + x = x*(x^2 + 1),
q(3,x) = x^7 + 7*x^5 + 7*x^3 + x = x*(x^2 + 1)*(x^4 + 6*x^2 + 1),
q(4,x) = x^15 + 35*x^13 + 273*x^11 + 715*x^9 + 715*x^7 + 273*x^5 + 35*x^3 + x = x*(x^2 + 1)*(x^4 + 6*x^2 + 1)*(x^8 + 28*x^6 + 70*x^4 + 28*x^2 + 1),
etc.
		

Crossrefs

Programs

  • Mathematica
    t={1};T[n_,k_]:=Table[2^(-n)Binomial[2^n,2k+1],{n,1,6},{k,0,2^(n-1)-1}];Do[AppendTo[t,T[n,k]]];Flatten[t] (* Indranil Ghosh, Feb 22 2017 *)

Formula

T(n, k) = 1 for n=0, k=0, and T(n, k) = 2^(-n) * binomial(2^n,2*k+1) = A103328(2^(n-1),k) for k = 0..2^(n-1)-1 and n >= 1. - Wolfdieter Lang, Jan 20 2017

Extensions

Edited. - Wolfdieter Lang, Jan 20 2017
More terms from Indranil Ghosh, Feb 22 2017

A059980 Number of 8-dimensional cage assemblies.

Original entry on oeis.org

1, 6561, 1679616, 100000000, 2562890625, 37822859361, 377801998336, 2821109907456, 16815125390625, 83733937890625, 360040606269696, 1370114370683136, 4702525276151521, 14774554437890625, 42998169600000000, 117033789351264256, 300283484326400961
Offset: 1

Views

Author

Robert G. Wilson v, Mar 06 2001

Keywords

References

  • Clifford A. Pickover, "Wonders of Numbers, Adventures in Mathematics, Mind and Meaning," Oxford University Press, 2001, p. 325.

Crossrefs

Programs

  • Mathematica
    m = 8; Table[n^m (n + 1)^m/2^m, {n, 1, 18}]

Formula

G.f.: -x*(x^14 +6544*x^13 +1568215*x^12 +72338144*x^11 +1086859301*x^10 +6727188848*x^9 +19323413187*x^8 +27306899520*x^7 +19323413187*x^6 +6727188848*x^5 +1086859301*x^4 +72338144*x^3 +1568215*x^2 +6544*x +1)/(x-1)^17. - Colin Barker, Jul 09 2012
From Peter Bala, Jul 02 2019 (Start)
a(n) = (n*(n + 1)/2)^8.
a(n) = (1/16)*( S(9,n) + 7*S(11,n) + 7*S(13,n) + S(15,n) ), where S(r,n) = Sum_{k = 1..n} k^r. Cf. A059977 and A059978. (End)
From Amiram Eldar, May 15 2022: (Start)
Sum_{n>=1} 1/a(n) = 146432*Pi^2 + 5632*Pi^4/3 + 2048*Pi^6/105 + 256*Pi^8/4725 - 1647360.
Sum_{n>=1} (-1)^(n+1)/a(n) = 1647360 - 1757184*log(2) - 304128*zeta(3) - 57600*zeta(5) - 4032*zeta(7). (End)
Showing 1-10 of 10 results.