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 17 results. Next

A002162 Decimal expansion of the natural logarithm of 2.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Newton calculated the first 16 terms of this sequence.
Area bounded by y = tan x, y = cot x, y = 0. - Clark Kimberling, Jun 26 2020
Choose four values independently and uniformly at random from the unit interval [0,1]. Sort them, and label them a,b,c,d from least to greatest (so that a b^2+c^2. - Akiva Weinberger, Dec 02 2024
Define the trihyperboloid to be the intersection of the three solid hyperboloids x^2+y^2-z^2<1, x^2-y^2+z^2<1, and -x^2+y^2+z^2<1. This fits perfectly within the cube [-1,1]^3. Then this is the ratio of the volume of the trihyperboloid to its bounding cube. - Akiva Weinberger, Dec 02 2024

Examples

			0.693147180559945309417232121458176568075500134360255254120680009493393...
		

References

  • G. Boros and V. H. Moll, Irresistible Integrals: Symbolics, Analysis and Experiments in the Evaluation of Integrals, Cambridge University Press, 2004.
  • Calvin C. Clawson, Mathematical Mysteries: The Beauty and Magic of Numbers, Springer, 2013. See p. 227.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 24, 250.
  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, Sections 1.3.3, 2.21, 6.2, and 7.2.
  • 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).
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 25 and appendix A, equations 25:14:3 and A:7:3 at pages 232, 670.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987, p. 29.

Crossrefs

Cf. A016730 (continued fraction), A002939, A008288, A142979, A142992.

Programs

  • Mathematica
    RealDigits[N[Log[2],200]][[1]] (* Vladimir Joseph Stephan Orlovsky, Feb 21 2011 *)
    RealDigits[Log[2],10,120][[1]] (* Harvey P. Dale, Jan 25 2024 *)
  • PARI
    { default(realprecision, 20080); x=10*log(2); for (n=0, 20000, d=floor(x); x=(x-d)*10; write("b002162.txt", n, " ", d)); } \\ Harry J. Smith, Apr 21 2009

Formula

log(2) = Sum_{k>=1} 1/(k*2^k) = Sum_{j>=1} (-1)^(j+1)/j.
log(2) = Integral_{t=0..1} dt/(1+t).
log(2) = (2/3) * (1 + Sum_{k>=1} 2/((4*k)^3-4*k)) (Ramanujan).
log(2) = 4*Sum_{k>=0} (3-2*sqrt(2))^(2*k+1)/(2*k+1) (Y. Luke). - R. J. Mathar, Jul 13 2006
log(2) = 1 - (1/2)*Sum_{k>=1} 1/(k*(2*k+1)). - Jaume Oliver Lafont, Jan 06 2009, Jan 08 2009
log(2) = 4*Sum_{k>=0} 1/((4*k+1)*(4*k+2)*(4*k+3)). - Jaume Oliver Lafont, Jan 08 2009
log(2) = 7/12 + 24*Sum_{k>=1} 1/(A052787(k+4)*A000079(k)). - R. J. Mathar, Jan 23 2009
From Alexander R. Povolotsky, Jul 04 2009: (Start)
log(2) = (1/4)*(3 - Sum_{n>=1} 1/(n*(n+1)*(2*n+1))).
log(2) = (230166911/9240 - Sum_{k>=1} (1/2)^k*(11/k + 10/(k+1) + 9/(k+2) + 8/(k+3) + 7/(k+4) + 6/(k+5) - 6/(k+7) - 7/(k+8) - 8/(k+9) - 9/(k+10) - 10/(k+11)))/35917. (End)
log(2) = A052882/A000670. - Mats Granvik, Aug 10 2009
From log(1-x-x^2) at x=1/2, log(2) = (1/2)*Sum_{k>=1} L(k)/(k*2^k), where L(n) is the n-th Lucas number (A000032). - Jaume Oliver Lafont, Oct 24 2009
log(2) = Sum_{k>=1} 1/(cos(k*Pi/3)*k*2^k) (cf. A176900). - Jaume Oliver Lafont, Apr 29 2010
log(2) = (Sum_{n>=1} 1/(n^2*(n+1)^2*(2*n+1)) + 11)/16. - Alexander R. Povolotsky, Jan 13 2011
log(2) = ((Sum_{n>=1} (2*n+1)/(Sum_{k=1..n} k^2)^2)+396)/576. - Alexander R. Povolotsky, Jan 14 2011
From Alexander R. Povolotsky, Dec 16 2008: (Start)
log(2) = 105*(319/44100 - Sum_{n>=1} 1/(2*n*(2*n+1)*(2*n+3)*(2*n+5)*(2*n+7))).
log(2) = 319/420 - (3/2)*Sum_{n>=1} 1/(6*n^2+39*n+63). (End)
log(2) = Sum_{k>=1} A191907(2,k)/k. - Mats Granvik, Jun 19 2011
log(2) = Integral_{x=0..oo} 1/(1 + e^x) dx. - Jean-François Alcover, Mar 21 2013
log(2) = lim_{s->1} zeta(s)*(1-1/2^(s-1)). - Mats Granvik, Jun 18 2013
From Peter Bala, Dec 10 2013: (Start)
log(2) = 2*Sum_{n>=1} 1/( n*A008288(n-1,n-1)*A008288(n,n) ), a result due to Burnside.
log(2) = (1/3)*Sum_{n >= 0} (5*n+4)/( (3*n+1)*(3*n+2)*C(3*n,n) )*(1/2)^n = (1/12)*Sum_{n >= 0} (28*n+17)/( (3*n+1)*(3*n+2)*C(3*n,n) )*(-1/4)^n.
log(2) = (3/16)*Sum_{n >= 0} (14*n+11)/( (4*n+1)*(4*n+3)*C(4*n,2*n) )*(1/4)^n = (1/12)*Sum_{n >= 0} (34*n+25)/( (4*n+1)*(4*n+3)*C(4*n,2*n) )*(-1/18)^n. For more series of this type see the Bala link.
See A142979 for series acceleration formulas for log(2) obtained from the Mercator series log(2) = Sum_{n >= 1} (-1)^(n+1)/n. See A142992 for series for log(2) related to the root lattice C_n. (End)
log(2) = lim_{n->oo} Sum_{k=2^n..2^(n+1)-1} 1/k. - Richard R. Forberg, Aug 16 2014
From Peter Bala, Feb 03 2015: (Start)
log(2) = (2/3)*Sum_{k >= 0} 1/((2*k + 1)*9^k).
Define a pair of integer sequences A(n) = 9^n*(2*n + 1)!/n! and B(n) = A(n)*Sum_{k = 0..n} 1/((2*k + 1)*9^k). Both satisfy the same second-order recurrence equation u(n) = (40*n + 16)*u(n-1) - 36*(2*n - 1)^2*u(n-2). From this observation we obtain the continued fraction expansion log(2) = (2/3)*(1 + 2/(54 - 36*3^2/(96 - 36*5^2/(136 - ... - 36*(2*n - 1)^2/((40*n + 16) - ... ))))). Cf. A002391, A073000 and A105531 for similar expansions. (End)
log(2) = Sum_{n>=1} (Zeta(2*n)-1)/n. - Vaclav Kotesovec, Dec 11 2015
From Peter Bala, Oct 30 2016: (Start)
Asymptotic expansions:
for N even, log(2) - Sum_{k = 1..N/2} (-1)^(k-1)/k ~ (-1)^(N/2)*(1/N - 1/N^2 + 2/N^4 - 16/N^6 + 272/N^8 - ...), where the sequence of unsigned coefficients [1, 1, 2, 16, 272, ...] is A000182 with an extra initial term of 1. See Borwein et al., Theorem 1 (b);
for N odd, log(2) - Sum_{k = 1..(N-1)/2} (-1)^(k-1)/k ~ (-1)^((N-1)/2)*(1/N - 1/N^3 + 5/N^5 - 61/N^7 + 1385/N^9 - ...), by Borwein et al., Lemma 2 with f(x) := 1/(x + 1/2), h := 1/2 and then set x = (N - 1)/2, where the sequence of unsigned coefficients [1, 1, 5, 61, 1385, ...] is A000364. (End)
log(2) = lim_{n->oo} Sum_{k=1..n} sin(1/(n+k)). See Mathematical Reflections link. - Michel Marcus, Jan 07 2017
log(2) = Sum_{n>=1} A006519(n) / ((1 + 2^A006519(n)) * A000265(n) * (1 + A000265(n))). - Nicolas Nagel, Mar 19 2018
From Amiram Eldar, Jul 02 2020: (Start)
Equals Sum_{k>=2} zeta(k)/2^k.
Equals -Sum_{k>=2} log(1 - 1/k^2).
Equals Sum_{k>=1} 1/A002939(k).
Equals Integral_{x=0..Pi/3} tan(x) dx. (End)
log(2) = Integral_{x=0..Pi/2} (sec(x) - tan(x)) dx. - Clark Kimberling, Jul 08 2020
From Peter Bala, Nov 14 2020: (Start)
log(2) = Integral_{x = 0..1} (x - 1)/log(x) dx (Boros and Moll, p. 97).
log(2) = (1/2)*Integral_{x = 0..1} (x + 2)*(x - 1)^2/log(x)^2 dx.
log(2) = (1/4)*Integral_{x = 0..1} (x^2 + 3*x + 4)*(x - 1)^3/log(x)^3 dx. (End)
log(2) = 2*arcsinh(sqrt(2)/4) = 2*sqrt(2)*Sum_{n >= 0} (-1)^n*C(2*n,n)/ ((8*n+4)*32^n) = 3*Sum_{n >= 0} (-1)^n/((8*n+4)*(2^n)*C(2*n,n)). - Peter Bala, Jan 14 2022
log(2) = Integral_{x=0..oo} ( e^(-x) * (1-e^(-2x)) * (1-e^(-4x)) * (1-e^(-6x)) ) / ( x * (1-e^(-14x)) ) dx (see Crux Mathematicorum link). - Bernard Schott, Jul 11 2022
From Peter Bala, Oct 22 2023: (Start)
log(2) = 23/32 + 2!^3/16 * Sum_{n >= 1} (-1)^n * (n + 1)/(n*(n + 1)*(n + 2))^2 = 707/1024 - 4!^3/(16^2 * 2!^2) * Sum_{n >= 1} (-1)^n * (n + 2)/(n*(n + 1)*(n + 2)*(n + 3)*(n + 4))^2 = 42611/61440 + 6!^3/(16^3 * 3!^2) * Sum_{n >= 1} (-1)^n * (n + 3)/(n*(n + 1)*(n + 2)*(n + 3)*(n + 4)*(n + 5)*(n + 6))^2.
More generally, it appears that for k >= 0, log(2) = c(k) + (2*k)!^3/(16^k * k!^2) * Sum_{n >= 1} (-1)^(n+k+1) * (n + k)/(n*(n + 1)*...*(n + 2*k))^2 , where c(k) is a rational approximation to log(2). The first few values of c(k) are [0, 23/32, 707/1024, 42611/61440, 38154331/55050240, 76317139/110100480, 26863086823/38755368960, ...].
Let P(n,k) = n*(n + 1)*...*(n + k).
Conjecture: for k >= 0 and r even with r - 1 <= k, the series Sum_{n >= 1} (-1)^n * (d/dn)^r (P(n,k)) / (P(n,k)^2 = A(r,k)*log(2) + B(r,k), where A(r,k) and B(r,k) are both rational numbers. (End)
From Peter Bala, Nov 13 2023: (Start)
log(2) = 5/8 + (1/8)*Sum_{k >= 1} (-1)^(k+1) * (2*k + 1)^2 / ( k*(k + 1) )^4
= 257/384 + (3!^5/2^9)*Sum_{k >= 1} (-1)^(k+1) * (2*k + 1)*(2*k + 3)^2*(2*k + 5) / ( k*(k + 1)*(k + 2)*(k + 3) )^4
= 267515/393216 + (5!^5/2^19)*Sum_{k >= 1} (-1)^(k+1) * (2*k + 1)*(2*k + 3)*(2*k + 5)^2*(2*k + 7)*(2*k + 9) / ( k*(k + 1)*(k + 2)*(k + 3)*(k + 4)*(k + 5) )^4
log(2) = 3/4 - 1/128 * Sum_{k >= 0} (-1/16)^k * (10*k + 12)*binomial(2*k+2,k+1)/ ((k + 1)*(2*k + 3)). The terms of the series are O(1/(k^(3/2)*4^n)). (End)
log(2) = eta(1) is a period, where eta(x) is the Dirichlet eta function. - Andrea Pinos, Mar 19 2024
log(2) = K_{n>=0} (n^2 + [n=0])/1, where K is the Gauss notation for an infinite continued fraction. In the expanded form, log(2) = 1/(1 + 1/(1 + 4/(1 + 9/1 + 16/(1 + 25/(1 + ... (see Clawson at p. 227). - Stefano Spezia, Jul 01 2024
log(2) = lim_{n->oo} Sum_{k=1..n} 1/(n + k) = lim_{x->0} (2^x - 1)/x = lim_{x->0} (2^x - 2^(-x))/(2*x) (see Finch). - Stefano Spezia, Oct 19 2024
From Colin Linzer, Nov 08 2024: (Start)
log(2) = Integral_{t=0...oo} (1 - tanh(t)) dt.
log(2) = Integral_{t=0...1} arctanh(t) dt.
log(2) = (1/2) * Integral_{t=-1...1} |arctanh(t)| dt. (End)
log(2) = 1 + Sum_{n >= 1} (-1)^n/(n*(4*n^2 - 1)) = 1/2 + (1/2)*Sum_{n >= 1} 1/(n*(4*n^2 - 1)). - Peter Bala, Jan 07 2025
log(2) = Integral_{x=0..1} Integral_{y=0..1} 1/((1 - x*y)*(1 + x)*(1 + y)) dy dx. - Kritsada Moomuang, May 22 2025

A000389 Binomial coefficients C(n,5).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 6, 21, 56, 126, 252, 462, 792, 1287, 2002, 3003, 4368, 6188, 8568, 11628, 15504, 20349, 26334, 33649, 42504, 53130, 65780, 80730, 98280, 118755, 142506, 169911, 201376, 237336, 278256, 324632, 376992, 435897, 501942, 575757, 658008, 749398
Offset: 0

Views

Author

Keywords

Comments

a(n+4) is the number of inequivalent ways of coloring the vertices of a regular 4-dimensional simplex with n colors, under the full symmetric group S_5 of order 120, with cycle index (x1^5 + 10*x1^3*x2 + 20*x1^2*x3 + 15*x1*x2^2 + 30*x1*x4 + 20*x2*x3 + 24*x5)/120.
Figurate numbers based on 5-dimensional regular simplex. According to Hyun Kwang Kim, it appears that every nonnegative integer can be represented as the sum of g = 10 of these 5-simplex(n) numbers (compared with g=3 for triangular numbers, g=5 for tetrahedral numbers and g=8 for pentatope numbers). - Jonathan Vos Post, Nov 28 2004
The convolution of the nonnegative integers (A001477) with the tetrahedral numbers (A000292), which are the convolution of the nonnegative integers with themselves (making appropriate allowances for offsets of all sequences). - Graeme McRae, Jun 07 2006
a(n) is the number of terms in the expansion of (a_1 + a_2 + a_3 + a_4 + a_5 + a_6)^n. - Sergio Falcon, Feb 12 2007
Product of five consecutive numbers divided by 120. - Artur Jasinski, Dec 02 2007
Equals binomial transform of [1, 5, 10, 10, 5, 1, 0, 0, 0, ...]. - Gary W. Adamson, Feb 02 2009
Equals INVERTi transform of A099242 (1, 7, 34, 153, 686, 3088, ...). - Gary W. Adamson, Feb 02 2009
For a team with n basketball players (n>=5), this sequence is the number of possible starting lineups of 5 players, without regard to the positions (center, forward, guard) of the players. - Mohammad K. Azarian, Sep 10 2009
a(n) is the number of different patterns, regardless of order, when throwing (n-5) 6-sided dice. For example, one die can display the 6 numbers 1, 2, ..., 6; two dice can display the 21 digit-pairs 11, 12, ..., 56, 66. - Ian Duff, Nov 16 2009
Sum of the first n pentatope numbers (1, 5, 15, 35, 70, 126, 210, ...), see A000332. - Paul Muljadi, Dec 16 2009
Sum_{n>=0} a(n)/n! = e/120. Sum_{n>=4} a(n)/(n-4)! = 501*e/120. See A067764 regarding the second ratio. - Richard R. Forberg, Dec 26 2013
For a set of integers {1,2,...,n}, a(n) is the sum of the 2 smallest elements of each subset with 4 elements, which is 3*C(n+1,5) (for n>=4), hence a(n) = 3*C(n+1,5) = 3*A000389(n+1). - Serhat Bulut, Mar 11 2015
a(n) = fallfac(n,5)/5! is also the number of independent components of an antisymmetric tensor of rank 5 and dimension n >= 1. Here fallfac is the falling factorial. - Wolfdieter Lang, Dec 10 2015
Number of compositions (ordered partitions) of n+1 into exactly 6 parts. - Juergen Will, Jan 02 2016
Number of weak compositions (ordered weak partitions) of n-5 into exactly 6 parts. - Juergen Will, Jan 02 2016
a(n+3) could be the general number of all geodetic graphs of diameter n>=2 homeomorphic to the Petersen Graph. - Carlos Enrique Frasser, May 24 2018
From Robert A. Russell, Dec 24 2020: (Start)
a(n) is the number of chiral pairs of colorings of the 5 tetrahedral facets (or vertices) of the regular 4-D simplex (5-cell, pentachoron, Schläfli symbol {3,3,3}) using subsets of a set of n colors. Each member of a chiral pair is a reflection but not a rotation of the other.
a(n+4) is the number of unoriented colorings of the 5 tetrahedral facets of the regular 4-D simplex (5-cell, pentachoron) using subsets of a set of n colors. Each chiral pair is counted as one when enumerating unoriented arrangements. (End)
For integer m and positive integer r >= 4, the polynomial a(n) + a(n + m) + a(n + 2*m) + ... + a(n + r*m) in n has its zeros on the vertical line Re(n) = (4 - r*m)/2 in the complex plane. - Peter Bala, Jun 02 2024

Examples

			G.f. = x^5 + 6*x^6 + 21*x^7 + 56*x^8 + 126*x^9 + 252*x^10 + 462*x^11 + ...
For A={1,2,3,4}, the only subset with 4 elements is {1,2,3,4}; sum of 2 minimum elements of this subset: a(4) = 1+2 = 3 = 3*C(4+1,5).
For A={1,2,3,4,5}, the subsets with 4 elements are {1,2,3,4}, {1,2,3,5}, {1,2,4,5}, {1,3,4,5}, {2,3,4,5}; sum of 2 smallest elements of each subset: a(5) = (1+2)+(1+2)+(1+2)+(1+3)+(2+3) = 18 = 3*C(5+1,5). - _Serhat Bulut_, Mar 11 2015
a(6) = 6 from the six independent components of an antisymmetric tensor A of rank 5 and dimension 6: A(1,2,3,4,5), A(1,2,3,4,6), A(1,2,3,5,6), A(1,2,4,5,6), A(1,3,4,5,6), A(2,3,4,5,6). See the Dec 10 2015 comment. - _Wolfdieter Lang_, Dec 10 2015
		

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. 828.
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 196.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 7.
  • Gupta, Hansraj; Partitions of j-partite numbers into twelve or a smaller number of parts. Collection of articles dedicated to Professor P. L. Bhatnagar on his sixtieth birthday. Math. Student 40 (1972), 401-441 (1974).
  • J. C. P. Miller, editor, Table of Binomial Coefficients. Royal Society Mathematical Tables, Vol. 3, Cambridge Univ. Press, 1954.
  • 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. A099242. - Gary W. Adamson, Feb 02 2009
Cf. A242023. A104712 (fourth column, k=5).
5-cell colorings: A337895 (oriented), A132366(n-1) (achiral).
Unoriented colorings: A063843 (5-cell edges, faces), A128767 (8-cell vertices, 16-cell facets), A337957 (16-cell vertices, 8-cell facets), A338949 (24-cell), A338965 (600-cell vertices, 120-cell facets).
Chiral colorings: A331352 (5-cell edges, faces), A337954 (8-cell vertices, 16-cell facets), A234249 (16-cell vertices, 8-cell facets), A338950 (24-cell), A338966 (600-cell vertices, 120-cell facets).

Programs

  • Haskell
    a000389 n = a000389_list !! n
    a000389_list = 0 : 0 : f [] a000217_list where
       f xs (t:ts) = (sum $ zipWith (*) xs a000217_list) : f (t:xs) ts
    -- Reinhard Zumkeller, Mar 03 2015, Apr 13 2012
    
  • Magma
    [Binomial(n, 5): n in [0..40]]; // Vincenzo Librandi, Mar 12 2015
  • Maple
    f:=n->(1/120)*(n^5-10*n^4+35*n^3-50*n^2+24*n): seq(f(n), n=0..60);
    ZL := [S, {S=Prod(B,B,B,B,B,B), B=Set(Z, 1 <= card)}, unlabeled]: seq(combstruct[count](ZL, size=n+1), n=0..42); # Zerinvary Lajos, Mar 13 2007
    A000389:=1/(z-1)**6; # Simon Plouffe, 1992 dissertation
  • Mathematica
    Table[Binomial[n, 5], {n, 5, 50}] (* Stefan Steinerberger, Apr 02 2006 *)
    CoefficientList[Series[x^5 / (1 - x)^6, {x, 0, 40}], x] (* Vincenzo Librandi, Mar 12 2015 *)
    LinearRecurrence[{6,-15,20,-15,6,-1},{0,0,0,0,0,1},50] (* Harvey P. Dale, Jul 17 2016 *)
  • PARI
    (conv(u,v)=local(w); w=vector(length(u),i,sum(j=1,i,u[j]*v[i+1-j])); w);
    (t(n)=n*(n+1)/2); u=vector(10,i,t(i)); conv(u,u)
    

Formula

G.f.: x^5/(1-x)^6.
a(n) = n*(n-1)*(n-2)*(n-3)*(n-4)/120.
a(n) = (n^5-10*n^4+35*n^3-50*n^2+24*n)/120. (Replace all x_i's in the cycle index with n.)
a(n+2) = Sum_{i+j+k=n} i*j*k. - Benoit Cloitre, Nov 01 2002
Convolution of triangular numbers (A000217) with themselves.
Partial sums of A000332. - Alexander Adamchuk, Dec 19 2004
a(n) = -A110555(n+1,5). - Reinhard Zumkeller, Jul 27 2005
a(n+3) = (1/2!)*(d^2/dx^2)S(n,x)|A049310.%20-%20_Wolfdieter%20Lang">{x=2}, n>=2, one half of second derivative of Chebyshev S-polynomials evaluated at x=2. See A049310. - _Wolfdieter Lang, Apr 04 2007
a(n) = A052787(n+5)/120. - Zerinvary Lajos, Apr 26 2007
Sum_{n>=5} 1/a(n) = 5/4. - R. J. Mathar, Jan 27 2009
For n>4, a(n) = 1/(Integral_{x=0..Pi/2} 10*(sin(x))^(2*n-9)*(cos(x))^9). - Francesco Daddi, Aug 02 2011
Sum_{n>=5} (-1)^(n + 1)/a(n) = 80*log(2) - 655/12 = 0.8684411114... - Richard R. Forberg, Aug 11 2014
a(n) = -a(4-n) for all n in Z. - Michael Somos, Oct 07 2014
0 = a(n)*(+a(n+1) + 4*a(n+2)) + a(n+1)*(-6*a(n+1) + a(n+2)) for all n in Z. - Michael Somos, Oct 07 2014
a(n) = 3*C(n+1, 5) = 3*A000389(n+1). - Serhat Bulut, Mar 11 2015
From Ilya Gutkovskiy, Jul 23 2016: (Start)
E.g.f.: x^5*exp(x)/120.
Inverse binomial transform of A054849. (End)
From Robert A. Russell, Dec 24 2020: (Start)
a(n) = A337895(n) - a(n+4) = (A337895(n) - A132366(n-1)) / 2 = a(n+4) - A132366(n-1).
a(n+4) = A337895(n) - a(n) = (A337895(n) + A132366(n-1)) / 2 = a(n) + A132366(n-1).
a(n+4) = 1*C(n,1) + 4*C(n,2) + 6*C(n,3) + 4*C(n,4) + 1*C(n,5), where the coefficient of C(n,k) is the number of unoriented pentachoron colorings using exactly k colors. (End)

Extensions

Corrected formulas that had been based on other offsets. - R. J. Mathar, Jun 16 2009
I changed the offset to 0. This will require some further adjustments to the formulas. - N. J. A. Sloane, Aug 01 2010

A052762 Products of 4 consecutive integers: a(n) = n*(n-1)*(n-2)*(n-3).

Original entry on oeis.org

0, 0, 0, 0, 24, 120, 360, 840, 1680, 3024, 5040, 7920, 11880, 17160, 24024, 32760, 43680, 57120, 73440, 93024, 116280, 143640, 175560, 212520, 255024, 303600, 358800, 421200, 491400, 570024, 657720, 755160, 863040, 982080, 1113024
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

Also, starting with n=4, the square of area of cyclic quadrilateral with sides n, n-1, n-2, n-3. - Zak Seidov, Jun 20 2003
Number of n-colorings of the complete graph on 4 vertices, which is also the tetrahedral graph. - Eric M. Schmidt, Oct 31 2012
Cf. A130534 for relations to colored forests and disposition of flags on flagpoles. - Tom Copeland, Apr 05 2014
Number of 4-permutations of the set {1, 2, ..., n}. - Joerg Arndt, Apr 05 2014

Crossrefs

Programs

  • Magma
    [n*(n-1)*(n-2)*(n-3): n in [0..30]]; // G. C. Greubel, Nov 19 2017
  • Maple
    spec := [S,{B=Set(Z),S=Prod(Z,Z,Z,Z,B)},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
    seq(numbperm (n,4), n=0..34); # Zerinvary Lajos, Apr 26 2007
    G(x):=x^4*exp(x): f[0]:=G(x): for n from 1 to 34 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n],n=0..34); # Zerinvary Lajos, Apr 05 2009
  • Mathematica
    Table[n*(n+1)*(n+2)*(n+3), {n,-3,60}] (* Vladimir Joseph Stephan Orlovsky, Apr 21 2010 *)
    Times@@@Partition[Range[-3,60], 4, 1] (* Harvey P. Dale, May 09 2012 *)
    LinearRecurrence[ {5,-10,10,-5,1}, {0,0,0,0,24}, 60] (* Harvey P. Dale, May 09 2012 *)
  • Maxima
    A052762(n):=n*(n-1)*(n-2)*(n-3)$
    makelist(A052762(n),n,0,30); /* Martin Ettl, Nov 03 2012 */
    
  • PARI
    a(n)=24*binomial(n,4) \\ Charles R Greathouse IV, Nov 20 2011
    

Formula

a(n) = n*(n-1)*(n-2)*(n-3) = n!/(n-4)! (for n >= 4).
a(n) = A001094(n) - n.
E.g.f.: x^4*exp(x).
Recurrence: {a(1)=0, a(2)=0, a(3)=0, a(4)=24, (-1-n)*a(n) + (n-3)*a(n+1)}.
a(n) + 1 = A062938(n-4) for n > 4. - Amarnath Murthy, Dec 13 2003
a(n) = numbperm(n, 4). - Zerinvary Lajos, Apr 26 2007
O.g.f.: -24*x^4/(-1+x)^5. - R. J. Mathar, Nov 23 2007
For n > 4: a(n) = A173333(n, n-4). - Reinhard Zumkeller, Feb 19 2010
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5), with a(0)=0, a(1)=0, a(2)=0, a(3)=0, a(4)=24. - Harvey P. Dale, May 09 2012
a(n) = a(n-1) + 4*A007531(n). - J. M. Bergot, May 30 2012
a(n) - 1 = A069756(n-2) for n >= 4. - Jean-Christophe Hervé, Nov 01 2015
a(n) = 24 * A000332(n). - Bruce J. Nicholson, Apr 03 2017
From R. J. Mathar, Jun 30 2021: (Start)
Sum_{n>=4} 24*(-1)^n/a(n) = A242023.
Sum_{n>=4} 1/a(n) = 1/18. (End)

Extensions

More terms from Henry Bottomley, Mar 20 2000
Formula corrected by Philippe Deléham, Dec 12 2003

A173333 Triangle read by rows: T(n, k) = n! / k!, 1 <= k <= n.

Original entry on oeis.org

1, 2, 1, 6, 3, 1, 24, 12, 4, 1, 120, 60, 20, 5, 1, 720, 360, 120, 30, 6, 1, 5040, 2520, 840, 210, 42, 7, 1, 40320, 20160, 6720, 1680, 336, 56, 8, 1, 362880, 181440, 60480, 15120, 3024, 504, 72, 9, 1, 3628800, 1814400, 604800, 151200, 30240, 5040, 720, 90, 10, 1
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 19 2010

Keywords

Comments

From Wolfdieter Lang, Jun 27 2012: (Start)
T(n-1,k), k=1,...,n-1, gives the number of representative necklaces with n beads (C_N symmetry) of n+1-k distinct colors, say c[1],c[2],...,c[n-k+1], corresponding to the color signature determined by the partition k,1^(n-k) of n. The representative necklaces have k beads of color c[1]. E.g., n=4, k=2: partition 2,1,1, color signature (parts as exponents) c[1]c[1]c[2]c[3], 3=T(3,2) necklaces (write j for color c[j]): cyclic(1123), cyclic(1132) and cyclic(1213). See A212359 for the numbers for general partitions or color signatures. (End)

Examples

			Triangle starts:
n\k      1       2      3      4     5    6   7  8  9 10 ...
1        1
2        2       1
3        6       3      1
4       24      12      4      1
5      120      60     20      5     1
6      720     360    120     30     6    1
7     5040    2520    840    210    42    7   1
8    40320   20160   6720   1680   336   56   8  1
9   362880  181440  60480  15120  3024  504  72  9  1
10 3628800 1814400 604800 151200 30240 5040 720 90 10  1
... - _Wolfdieter Lang_, Jun 27 2012
		

Crossrefs

Row sums give A002627.
Central terms give A006963:
T(2*n-1,n) = A006963(n+1).
T(2*n,n) = A001813(n).
T(2*n,n+1) = A001761(n).
1 < k <= n: T(n,k) = T(n,k-1) / k.
1 <= k <= n: T(n+1,k) = A119741(n,n-k+1).
1 <= k <= n: T(n+1,k+1) = A162995(n,k).
T(n,1) = A000142(n).
T(n,2) = A001710(n) for n>1.
T(n,3) = A001715(n) for n>2.
T(n,4) = A001720(n) for n>3.
T(n,5) = A001725(n) for n>4.
T(n,6) = A001730(n) for n>5.
T(n,7) = A049388(n-7) for n>6.
T(n,8) = A049389(n-8) for n>7.
T(n,9) = A049398(n-9) for n>8.
T(n,10) = A051431(n) for n>9.
T(n,n-7) = A159083(n+1) for n>7.
T(n,n-6) = A053625(n+1) for n>6.
T(n,n-5) = A052787(n) for n>5.
T(n,n-4) = A052762(n) for n>4.
T(n,n-3) = A007531(n) for n>3.
T(n,n-2) = A002378(n-1) for n>2.
T(n,n-1) = A000027(n) for n>1.
T(n,n) = A000012(n).

Programs

  • Haskell
    a173333 n k = a173333_tabl !! (n-1) !! (k-1)
    a173333_row n = a173333_tabl !! (n-1)
    a173333_tabl = map fst $ iterate f ([1], 2)
       where f (row, i) = (map (* i) row ++ [1], i + 1)
    -- Reinhard Zumkeller, Jul 04 2012
  • Mathematica
    Table[n!/k!, {n, 1, 10}, {k, 1, n}] // Flatten (* Jean-François Alcover, Mar 01 2019 *)

Formula

E.g.f.: (exp(x*y) - 1)/(x*(1 - y)). - Olivier Gérard, Jul 07 2011
T(n,k) = A094587(n,k), 1 <= k <= n. - Reinhard Zumkeller, Jul 05 2012

A045619 Numbers that are the products of 2 or more consecutive integers.

Original entry on oeis.org

0, 2, 6, 12, 20, 24, 30, 42, 56, 60, 72, 90, 110, 120, 132, 156, 182, 210, 240, 272, 306, 336, 342, 360, 380, 420, 462, 504, 506, 552, 600, 650, 702, 720, 756, 812, 840, 870, 930, 990, 992, 1056, 1122, 1190, 1260, 1320, 1332, 1406, 1482, 1560, 1640, 1680
Offset: 1

Views

Author

Keywords

Comments

Erdős and Selfridge proved that, apart from the first term, these are never perfect powers (A001597). - T. D. Noe, Oct 13 2002
Numbers of the form x!/y! with y+1 < x. - Reinhard Zumkeller, Feb 20 2008

Examples

			30 is in the sequence as 30 = 5*6 = 5*(5+1). - _David A. Corneth_, Oct 19 2021
		

Crossrefs

Programs

  • Mathematica
    maxNum = 1700; lst = {}; For[i = 1, i <= Sqrt[maxNum], i++, j = i + 1; prod = i*j; While[prod < maxNum, AppendTo[lst, prod]; j++; prod *= j]]; lst = Union[lst]
  • PARI
    list(lim)=my(v=List([0]),P,k=1,t); while(1, k++; P=binomial('n+k-1,k)*k!; if(subst(P,'n,1)>lim, break); for(n=1,lim, t=eval(P); if(t>lim, next(2)); listput(v,t))); Set(v) \\ Charles R Greathouse IV, Nov 16 2021
  • Python
    import heapq
    from sympy import sieve
    def aupton(terms, verbose=False):
        p = 6; h = [(p, 2, 3)]; nextcount = 4; aset = {0, 2}
        while len(aset) < terms:
            (v, s, l) = heapq.heappop(h)
            aset.add(v)
            if verbose: print(f"{v}, [= Prod_{{i = {s}..{l}}} i]")
            if v >= p:
                p *= nextcount
                heapq.heappush(h, (p, 2, nextcount))
                nextcount += 1
            v //= s; s += 1; l += 1; v *= l
            heapq.heappush(h, (v, s, l))
        return sorted(aset)
    print(aupton(52)) # Michael S. Branicky, Oct 19 2021
    

Formula

a(n) = A000142(A137911(n))/A000142(A137912(n)-1) for n>1. - Reinhard Zumkeller, Feb 27 2008
Since the oblong numbers (A002378) have relative density of 100%, we have a(n) ~ (n-1) n ~ n^2. - Daniel Forgues, Mar 26 2012
a(n) = n^2 - 2*n^(5/3) + O(n^(4/3)). - Charles R Greathouse IV, Aug 27 2013

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Jul 20 2000
More terms from Reinhard Zumkeller, Feb 27 2008
Incorrect program removed by David A. Corneth, Oct 19 2021

A265609 Array read by ascending antidiagonals: A(n,k) the rising factorial, also known as Pochhammer symbol, for n >= 0 and k >= 0.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 2, 0, 1, 3, 6, 6, 0, 1, 4, 12, 24, 24, 0, 1, 5, 20, 60, 120, 120, 0, 1, 6, 30, 120, 360, 720, 720, 0, 1, 7, 42, 210, 840, 2520, 5040, 5040, 0, 1, 8, 56, 336, 1680, 6720, 20160, 40320, 40320, 0
Offset: 0

Views

Author

Peter Luschny, Dec 19 2015

Keywords

Comments

The Pochhammer function is defined P(x,n) = x*(x+1)*...*(x+n-1). By convention P(0,0) = 1.
From Antti Karttunen, Dec 19 2015: (Start)
Apart from the initial row of zeros, if we discard the leftmost column and divide the rest of terms A(n,k) with (n+k) [where k is now the once-decremented column index of the new, shifted position] we get the same array back. See the given recursive formula.
When the numbers in array are viewed in factorial base (A007623), certain repeating patterns can be discerned, at least in a few of the topmost rows. See comment in A001710 and arrays A265890, A265892. (End)
A(n,k) is the k-th moment (about 0) of a gamma (Erlang) distribution with shape parameter n and rate parameter 1. - Geoffrey Critzer, Dec 24 2018

Examples

			Square array A(n,k) [where n=row, k=column] is read by ascending antidiagonals as:
A(0,0), A(1,0), A(0,1), A(2,0), A(1,1), A(0,2), A(3,0), A(2,1), A(1,2), A(0,3), ...
Array starts:
n\k [0  1   2    3     4      5        6         7          8]
--------------------------------------------------------------
[0] [1, 0,  0,   0,    0,     0,       0,        0,         0]
[1] [1, 1,  2,   6,   24,   120,     720,     5040,     40320]
[2] [1, 2,  6,  24,  120,   720,    5040,    40320,    362880]
[3] [1, 3, 12,  60,  360,  2520,   20160,   181440,   1814400]
[4] [1, 4, 20, 120,  840,  6720,   60480,   604800,   6652800]
[5] [1, 5, 30, 210, 1680, 15120,  151200,  1663200,  19958400]
[6] [1, 6, 42, 336, 3024, 30240,  332640,  3991680,  51891840]
[7] [1, 7, 56, 504, 5040, 55440,  665280,  8648640, 121080960]
[8] [1, 8, 72, 720, 7920, 95040, 1235520, 17297280, 259459200]
.
Seen as a triangle, T(n, k) = Pochhammer(n - k, k), the first few rows are:
   [0] 1;
   [1] 1, 0;
   [2] 1, 1,  0;
   [3] 1, 2,  2,   0;
   [4] 1, 3,  6,   6,    0;
   [5] 1, 4, 12,  24,   24,    0;
   [6] 1, 5, 20,  60,  120,  120,     0;
   [7] 1, 6, 30, 120,  360,  720,   720,     0;
   [8] 1, 7, 42, 210,  840, 2520,  5040,  5040,     0;
   [9] 1, 8, 56, 336, 1680, 6720, 20160, 40320, 40320, 0.
		

References

  • Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics, Addison-Wesley, 1994.
  • H. S. Wall, Analytic Theory of Continued Fractions, Chelsea 1973, p. 355.

Crossrefs

Triangle giving terms only up to column k=n: A124320.
Row 0: A000007, row 1: A000142, row 3: A001710 (from k=1 onward, shifted two terms left).
Column 0: A000012, column 1: A001477, column 2: A002378, columns 3-7: A007531, A052762, A052787, A053625, A159083 (shifted 2 .. 6 terms left respectively, i.e. without the extra initial zeros), column 8: A239035.
Row sums of the triangle: A000522.
A(n, n) = A000407(n-1) for n>0.
2^n*A(1/2,n) = A001147(n).
Cf. also A007623, A008279 (falling factorial), A173333, A257505, A265890, A265892.

Programs

  • Maple
    for n from 0 to 8 do seq(pochhammer(n,k), k=0..8) od;
  • Mathematica
    Table[Pochhammer[n, k], {n, 0, 8}, {k, 0, 8}]
  • Sage
    for n in (0..8): print([rising_factorial(n,k) for k in (0..8)])
    
  • Scheme
    (define (A265609 n) (A265609bi (A025581 n) (A002262 n)))
    (define (A265609bi row col) (if (zero? col) 1 (* (+ row col -1) (A265609bi row (- col 1)))))
    ;; Antti Karttunen, Dec 19 2015

Formula

A(n,k) = Gamma(n+k)/Gamma(n) for n > 0 and n^k for n=0.
A(n,k) = Sum_{j=0..k} n^j*S1(k,j), S1(n,k) the Stirling cycle numbers A132393(n,k).
A(n,k) = (k-1)!/(Sum_{j=0..k-1} (-1)^j*binomial(k-1, j)/(j+n)) for n >= 1, k >= 1.
A(n,k) = (n+k-1)*A(n,k-1) for k >= 1, A(n,0) = 1. - Antti Karttunen, Dec 19 2015
E.g.f. for row k: 1/(1-x)^k. - Geoffrey Critzer, Dec 24 2018
A(n, k) = FallingFactorial(n + k - 1, k). - Peter Luschny, Mar 22 2022
G.f. for row n as a continued fraction of Stieltjes type: 1/(1 - n*x/(1 - x/(1 - (n+1)*x/(1 - 2*x/(1 - (n+2)*x/(1 - 3*x/(1 - ... ))))))). See Wall, Chapter XVIII, equation 92.5. Cf. A226513. - Peter Bala, Aug 27 2023

A053625 Product of 6 consecutive integers.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 720, 5040, 20160, 60480, 151200, 332640, 665280, 1235520, 2162160, 3603600, 5765760, 8910720, 13366080, 19535040, 27907200, 39070080, 53721360, 72681840, 96909120, 127512000, 165765600, 213127200, 271252800, 342014400, 427518000, 530122320
Offset: 0

Views

Author

Henry Bottomley, Mar 20 2000

Keywords

Crossrefs

Programs

  • GAP
    F:=Factorial;; Concatenation([0,0,0,0,0,0], List([6..30], n-> F(n)/F(n-5) )); # G. C. Greubel, Aug 27 2019
  • Magma
    I:=[0,0,0,0,0,0,720]; [n le 7 select I[n] else 7*Self(n-1) -21*Self(n-2)+35*Self(n-3)-35*Self(n-4)+21*Self(n-5)-7*Self(n-6) +Self(n-7): n in [1..30]]; // Vincenzo Librandi, Apr 28 2012
    
  • Maple
    seq(combinat[numbperm](n, 6), n=0..31); # Zerinvary Lajos, Apr 26 2007
  • Mathematica
    CoefficientList[Series[720*x^6/(1-x)^7,{x,0,30}],x] (* Vincenzo Librandi, Apr 28 2012 *)
    Times@@@Partition[Range[-5,30],6,1] (* or *) LinearRecurrence[{7,-21,35,-35,21,-7,1},{0,0,0,0,0,0,720},30] (* Harvey P. Dale, Nov 13 2015 *)
    Pochhammer[Range[30]-6, 6] (* G. C. Greubel, Aug 27 2019 *)
  • PARI
    a(n)=factorback([n-5..n]) \\ Charles R Greathouse IV, Oct 07 2015
    
  • Sage
    [rising_factorial(n-5,6) for n in (0..30)] # G. C. Greubel, Aug 27 2019
    

Formula

a(n) = n*(n-1)*(n-2)*(n-3)*(n-4)*(n-5) = n!/(n-6)! = A052787(n)*(n-6) = a(n-1)*n/(n-6).
E.g.f.: x^6*exp(x).
a(n) = 720 * A000579(n). - Zerinvary Lajos, Apr 26 2007
For n > 5: a(n) = A173333(n, n-6). - Reinhard Zumkeller, Feb 19 2010
G.f.: 720*x^6/(1-x)^7. - Colin Barker, Mar 27 2012
a(n) = 7*a(n-1) - 21*a(n-2) + 35*a(n-3) - 35*a(n-4) + 21*a(n-5) - 7*a(n-6) + a(n-7). - Vincenzo Librandi, Apr 28 2012
From Amiram Eldar, Mar 08 2022: (Start)
Sum_{n>=6} 1/a(n) = 1/600.
Sum_{n>=6} (-1)^n/a(n) = 4*log(2)/15 - 661/3600. (End)

A054559 Number of labeled pure 2-complexes on n nodes (0-simplexes) with 5 2-simplexes and 8 1-simplexes.

Original entry on oeis.org

30, 180, 630, 1680, 3780, 7560, 13860, 23760, 38610, 60060, 90090, 131040, 185640, 257040, 348840, 465120, 610470, 790020, 1009470, 1275120, 1593900, 1973400, 2421900, 2948400, 3562650, 4275180, 5097330, 6041280, 7120080, 8347680, 9738960, 11309760, 13076910
Offset: 5

Views

Author

Vladeta Jovovic, Apr 10 2000

Keywords

Comments

Number of {T_1,T_2,...,T_k} where T_i,i=1..k are 3-subsets of an n-set such that {D | D is 2-subset of T_i for some i=1..k} has l elements; k=5,l=8.
Let H be the n X n Hilbert matrix H(i,j) = 1/(i+j-1) for 1 <= i,j <= n. Let B be the inverse matrix of H. The sum of the elements in row 3 of B equals -a(n+2). - T. D. Noe, May 01 2011

References

  • V. Jovovic, On the number of two-dimensional simplicial complexes (in Russian), Metody i sistemy tekhnicheskoy diagnostiki, Vypusk 16, Mezhvuzovskiy zbornik nauchnykh trudov, Izdatelstvo Saratovskogo universiteta, 1991.

Crossrefs

Programs

  • Magma
    I:=[30, 180, 630, 1680, 3780, 7560]; [n le 6 select I[n] else 6*Self(n-1)-15*Self(n-2)+20*Self(n-3)-15*Self(n-4)+6*Self(n-5)-Self(n-6): n in [1..30]]; // Vincenzo Librandi, Apr 29 2012
    
  • Mathematica
    Table[n*(n+1)*(n+2)*(n+3)*(n+4)/4, {n,1,100}] (* Vladimir Joseph Stephan Orlovsky, Jul 21 2009 *)
    CoefficientList[Series[30/(1-x)^6,{x,0,30}],x] (* Vincenzo Librandi, Apr 29 2012 *)
  • PARI
    x='x+O('x^30); Vec(serlaplace(x^5*exp(x)/4)) \\ G. C. Greubel, Nov 23 2017

Formula

a(n) = 30*C(n,5) = 30*A000389(n) = n*(n-1)*(n-2)*(n-3)*(n-4)/4.
G.f.: 30*x^5/(1-x)^6. - Colin Barker, Jan 19 2012
a(n) = 6*a(n-1) - 15*a(n-2) + 20*a(n-3) - 15*a(n-4) + 6*a(n-5) - a(n-6). - Vincenzo Librandi, Apr 29 2012
E.g.f.: x^5*exp(x)/4. - G. C. Greubel, Nov 23 2017
From Amiram Eldar, Mar 08 2022: (Start)
Sum_{n>=5} 1/a(n) = 1/24.
Sum_{n>=5} (-1)^(n+1)/a(n) = 8*log(2)/3 - 131/72. (End)

A215228 T(n,k) = number of length-n 0..k arrays connected end-around, with no sequence of L

Original entry on oeis.org

2, 3, 2, 4, 6, 0, 5, 12, 6, 0, 6, 20, 24, 12, 0, 7, 30, 60, 72, 0, 0, 8, 42, 120, 240, 120, 18, 0, 9, 56, 210, 600, 720, 408, 0, 0, 10, 72, 336, 1260, 2520, 2940, 840, 24, 0, 11, 90, 504, 2352, 6720, 12600, 10080, 2448, 0, 0, 12, 110, 720, 4032, 15120, 40110, 57960, 38640
Offset: 1

Views

Author

R. H. Hardin, Aug 06 2012

Keywords

Comments

Table starts
2 3 4 5 6 7 8 9 10
2 6 12 20 30 42 56 72 90
0 6 24 60 120 210 336 504 720
0 12 72 240 600 1260 2352 4032 6480
0 0 120 720 2520 6720 15120 30240 55440
0 18 408 2940 12600 40110 105168 240408 496080
0 0 840 10080 57960 228480 710640 1874880 4379760
0 24 2448 38640 280560 1338120 4883424 14783328 38962080
0 0 5760 140400 1330560 7761600 33384960 116212320 345945600
0 0 15960 529440 6394680 45291120 228945360 915183360 3075040080
0 66 39864 1956900 30548760 263674950 1568401296 7203324744
0 72 108024 7335840 146516040 1537291560 10751253072
Empirical: row n is a polynomial of degree n.
Coefficients for rows 1-10, highest power first:
1 1
1 1 0
1 0 -1 0
1 0 -1 0 0
1 0 -5 0 4 0
1 0 -6 5 5 -5 0
1 0 -7 0 14 0 -8 0
1 0 -8 0 27 -12 -20 12 0
1 0 -9 0 27 0 -31 0 12 0
1 0 -10 0 35 9 -60 -25 34 16 0
Row n is divisible by n.
Column k is divisible by k+1.
From Robert Israel, Nov 23 2017: (Start)
Row n is a monic polynomial of degree n.
Proof: Let b(j,n,k) be the number of such arrays taking exactly j different values.
Then T(n,k) = Sum_{j <= n} b(j,n,k). But since the j values may be any combination of 0..k taken j at a time, b(j,n,k) = binomial(k+1,j)* b(j,n,j-1) which (if nonzero) is a polynomial in k of degree j.
In particular, b(n,n,n-1) = n!, so b(n,n,k) has degree n and leading coefficient 1. (End)

Examples

			Some solutions for n=5, k=4:
  3  0  1  1  1  0  4  4  0  1  3  2  2  3  1  0
  2  4  0  3  0  4  3  2  2  2  4  0  4  4  4  1
  0  2  2  2  2  3  0  3  1  4  0  4  3  1  0  0
  3  0  3  0  3  1  3  4  4  0  3  0  0  3  4  2
  1  3  2  4  0  2  1  0  1  4  2  1  4  0  2  3
		

Crossrefs

Column 2 is A066297.
Row 2 is A002378.
Row 3 is A007531(n+1).
Row 4 is A047928(n+1).
Row 5 is A052787(n+2).

A239035 Product of 8 consecutive integers. a(n) = RisingFactorial(n, 8).

Original entry on oeis.org

0, 40320, 362880, 1814400, 6652800, 19958400, 51891840, 121080960, 259459200, 518918400, 980179200, 1764322560, 3047466240, 5079110400, 8204716800, 12893126400, 19769460480, 29654190720, 43609104000, 62990928000, 89513424000, 125318793600, 173059286400
Offset: 0

Views

Author

Michel Marcus, Mar 09 2014

Keywords

Crossrefs

Cf. Product of n consecutive integers: A002378 (n=2), A007531 (n=3), A052762 (n=4), A052787 (n=5), A053625 (n=6), A159083 (n=7).

Programs

  • Magma
    [(n^4+14*n^3+63*n^2+98*n+28)^2-(28+8*n)^2: n in [0..30]]; // Vincenzo Librandi, Mar 11 2014
    
  • Mathematica
    CoefficientList[Series[40320 x/(1 - x)^9, {x, 0, 30}], x] (* Vincenzo Librandi, Mar 11 2014 *)
  • PARI
    a(n) = prod(k=0, 7, n+k);
    
  • PARI
    concat([0], Vec(40320*x/(1-x)^9  + O(x^100))) \\ Colin Barker, Mar 09 2014
    
  • Sage
    print([rising_factorial(n, 8) for n in range(23)]) # Peter Luschny, Mar 22 2022

Formula

a(n) = n*(n+1)*(n+2)*(n+3)*(n+4)*(n+5)*(n+6)*(n+7).
a(n) = (n^4+14*n^3+63*n^2+98*n+28)^2 - (28+8*n)^2. (see Thill link).
G.f.: 40320*x / (1-x)^9. - Colin Barker, Mar 09 2014
From Amiram Eldar, Mar 08 2022: (Start)
Sum_{n>=1} 1/a(n) = 1/35280.
Sum_{n>=1} (-1)^(n+1)/a(n) = 8*log(2)/315 - 1163/66150. (End)
Showing 1-10 of 17 results. Next