Binomial coefficient

From testwiki
Jump to navigation Jump to search

Template:Short description Template:Redirect

File:Pascal's triangle 5.svg
The binomial coefficients can be arranged to form Pascal's triangle, in which each entry is the sum of the two immediately above.
File:Binomial theorem visualisation.svg
Visualisation of binomial expansion up to the 4th power

In mathematics, the binomial coefficients are the positive integers that occur as coefficients in the binomial theorem. Commonly, a binomial coefficient is indexed by a pair of integers Template:Math and is written (nk). It is the coefficient of the Template:Math term in the polynomial expansion of the binomial power Template:Math; this coefficient can be computed by the multiplicative formula

(nk)=n×(n1)××(nk+1)k×(k1)××1,

which using factorial notation can be compactly expressed as

(nk)=n!k!(nk)!.

For example, the fourth power of Template:Math is

(1+x)4=(40)x0+(41)x1+(42)x2+(43)x3+(44)x4=1+4x+6x2+4x3+x4,

and the binomial coefficient (42)=4×32×1=4!2!2!=6 is the coefficient of the Template:Math term.

Arranging the numbers (n0),(n1),,(nn) in successive rows for Template:Math gives a triangular array called Pascal's triangle, satisfying the recurrence relation

(nk)=(n1k1)+(n1k).

The binomial coefficients occur in many areas of mathematics, and especially in combinatorics. In combinatorics the symbol (nk) is usually read as "Template:Mvar choose Template:Mvar" because there are (nk) ways to choose an (unordered) subset of Template:Mvar elements from a fixed set of Template:Mvar elements. For example, there are (42)=6 ways to choose Template:Math elements from Template:Math, namely Template:Math, Template:Math, Template:Math, Template:Math, Template:Math and Template:Math.

The first form of the binomial coefficients can be generalized to (zk) for any complex number Template:Mvar and integer Template:Math, and many of their properties continue to hold in this more general form.

History and notation

Andreas von Ettingshausen introduced the notation (nk) in 1826,[1] although the numbers were known centuries earlier (see Pascal's triangle). In about 1150, the Indian mathematician Bhaskaracharya gave an exposition of binomial coefficients in his book Līlāvatī.[2]

Alternative notations include Template:Math, Template:Math, Template:Math, Template:Math,[3] Template:Math, and Template:Math, in all of which the Template:Mvar stands for combinations or choices; the Template:Mvar notation means the number of ways to choose k out of n objects. Many calculators use variants of the Template:Nowrap because they can represent it on a single-line display. In this form the binomial coefficients are easily compared to the numbers of [[permutation#k-permutations of n|Template:Mvar-permutations of Template:Mvar]], written as Template:Math, etc.

Definition and interpretations

Template:Diagonal split header 0 1 2 3 4
0 1 Template:SilverC Template:SilverC Template:SilverC Template:SilverC
1 1 1 Template:SilverC Template:SilverC Template:SilverC
2 1 2 1 Template:SilverC Template:SilverC
3 1 3 3 1 Template:SilverC
4 1 4 6 4 1
The first few binomial coefficients
on a left-aligned Pascal's triangle

For natural numbers (taken to include 0) Template:Mvar and Template:Mvar, the binomial coefficient (nk) can be defined as the coefficient of the monomial Template:Math in the expansion of Template:Math. The same coefficient also occurs (if Template:Math) in the binomial formula Template:NumBlk2 (valid for any elements Template:Mvar, Template:Mvar of a commutative ring), which explains the name "binomial coefficient".

Another occurrence of this number is in combinatorics, where it gives the number of ways, disregarding order, that Template:Mvar objects can be chosen from among Template:Mvar objects; more formally, the number of Template:Mvar-element subsets (or Template:Mvar-combinations) of an Template:Mvar-element set. This number can be seen as equal to the one of the first definition, independently of any of the formulas below to compute it: if in each of the Template:Mvar factors of the power Template:Math one temporarily labels the term Template:Mvar with an index Template:Mvar (running from Template:Math to Template:Mvar), then each subset of Template:Mvar indices gives after expansion a contribution Template:Math, and the coefficient of that monomial in the result will be the number of such subsets. This shows in particular that (nk) is a natural number for any natural numbers Template:Mvar and Template:Mvar. There are many other combinatorial interpretations of binomial coefficients (counting problems for which the answer is given by a binomial coefficient expression), for instance the number of words formed of Template:Mvar bits (digits 0 or 1) whose sum is Template:Mvar is given by (nk), while the number of ways to write k=a1+a2++an where every Template:Math is a nonnegative integer is given by Template:Tmath. Most of these interpretations can be shown to be equivalent to counting Template:Mvar-combinations.

Computing the value of binomial coefficients

Several methods exist to compute the value of (nk) without actually expanding a binomial power or counting Template:Mvar-combinations.

Recursive formula

One method uses the recursive, purely additive formula (nk)=(n1k1)+(n1k) for all integers n,k such that 1k<n, with boundary values (n0)=(nn)=1 for all integers Template:Math.

The formula follows from considering the set Template:Math and counting separately (a) the Template:Mvar-element groupings that include a particular set element, say "Template:Mvar", in every group (since "Template:Mvar" is already chosen to fill one spot in every group, we need only choose Template:Math from the remaining Template:Math) and (b) all the k-groupings that don't include "Template:Mvar"; this enumerates all the possible Template:Mvar-combinations of Template:Mvar elements. It also follows from tracing the contributions to Xk in Template:Math. As there is zero Template:Math or Template:Math in Template:Math, one might extend the definition beyond the above boundaries to include (nk)=0 when either Template:Math or Template:Math. This recursive formula then allows the construction of Pascal's triangle, surrounded by white spaces where the zeros, or the trivial coefficients, would be.

Multiplicative formula

A more efficient method to compute individual binomial coefficients is given by the formula (nk)=nk_k!=n(n1)(n2)(n(k1))k(k1)(k2)1=i=1kn+1ii, where the numerator of the first fraction, nk_, is a falling factorial. This formula is easiest to understand for the combinatorial interpretation of binomial coefficients. The numerator gives the number of ways to select a sequence of Template:Mvar distinct objects, retaining the order of selection, from a set of Template:Mvar objects. The denominator counts the number of distinct sequences that define the same Template:Mvar-combination when order is disregarded. This formula can also be stated in a recursive form. Using the "C" notation from above, Cn,k=Cn,k1(nk+1)/k, where Cn,0=1. It is readily derived by evaluating Cn,k/Cn,k1 and can intuitively be understood as starting at the leftmost coefficient of the n-th row of Pascal's triangle, whose value is always 1, and recursively computing the next coefficient to its right until the k-th one is reached.

Due to the symmetry of the binomial coefficients with regard to Template:Mvar and Template:Math, calculation of the above product, as well as the recursive relation, may be optimised by setting its upper limit to the smaller of Template:Mvar and Template:Math.

Factorial formula

Finally, though computationally unsuitable, there is the compact form, often used in proofs and derivations, which makes repeated use of the familiar factorial function: (nk)=n!k!(nk)!for  0kn, where Template:Math denotes the factorial of Template:Mvar. This formula follows from the multiplicative formula above by multiplying numerator and denominator by Template:Math; as a consequence it involves many factors common to numerator and denominator. It is less practical for explicit computation (in the case that Template:Mvar is small and Template:Mvar is large) unless common factors are first cancelled (in particular since factorial values grow very rapidly). The formula does exhibit a symmetry that is less evident from the multiplicative formula (though it is from the definitions) Template:NumBlk2 which leads to a more efficient multiplicative computational routine. Using the falling factorial notation, (nk)={nk_/k!if  kn2nnk_/(nk)!if  k>n2.

Generalization and connection to the binomial series

Template:Main The multiplicative formula allows the definition of binomial coefficients to be extended[4] by replacing n by an arbitrary number α (negative, real, complex) or even an element of any commutative ring in which all positive integers are invertible: (αk)=αk_k!=α(α1)(α2)(αk+1)k(k1)(k2)1for k and arbitrary α.

With this definition one has a generalization of the binomial formula (with one of the variables set to 1), which justifies still calling the (αk) binomial coefficients: Template:NumBlk2

This formula is valid for all complex numbers α and X with |X| < 1. It can also be interpreted as an identity of formal power series in X, where it actually can serve as definition of arbitrary powers of power series with constant coefficient equal to 1; the point is that with this definition all identities hold that one expects for exponentiation, notably (1+X)α(1+X)β=(1+X)α+βand((1+X)α)β=(1+X)αβ.

If α is a nonnegative integer n, then all terms with Template:Math are zero,[5] and the infinite series becomes a finite sum, thereby recovering the binomial formula. However, for other values of α, including negative integers and rational numbers, the series is really infinite.

Pascal's triangle

File:Pascal's triangle - 1000th row.png
1000th row of Pascal's triangle, arranged vertically, with grey-scale representations of decimal digits of the coefficients, right-aligned. The left boundary of the image corresponds roughly to the graph of the logarithm of the binomial coefficients, and illustrates that they form a log-concave sequence.

Template:Main

Pascal's rule is the important recurrence relation Template:NumBlk2 which can be used to prove by mathematical induction that (nk) is a natural number for all integer n ≥ 0 and all integer k, a fact that is not immediately obvious from formula (1). To the left and right of Pascal's triangle, the entries (shown as blanks) are all zero.

Pascal's rule also gives rise to Pascal's triangle:

0: 1
1: 1 1
2: 1 2 1
3: 1 3 3 1
4: 1 4 6 4 1
5: 1 5 10 10 5 1
6: 1 6 15 20 15 6 1
7: 21 35 35 21
8: 28 56 70 56 28

Row number Template:Mvar contains the numbers (nk) for Template:Math. It is constructed by first placing 1s in the outermost positions, and then filling each inner position with the sum of the two numbers directly above. This method allows the quick calculation of binomial coefficients without the need for fractions or multiplications. For instance, by looking at row number 5 of the triangle, one can quickly read off that

(x+y)5=1_x5+5_x4y+10_x3y2+10_x2y3+5_xy4+1_y5.

Combinatorics and statistics

Binomial coefficients are of importance in combinatorics because they provide ready formulas for certain frequent counting problems:

  • There are (nk) ways to choose k elements from a set of n elements. See Combination.
  • There are (n+k1k) ways to choose k elements from a set of n elements if repetitions are allowed. See Multiset.
  • There are (n+kk) strings containing k ones and n zeros.
  • There are (n+1k) strings consisting of k ones and n zeros such that no two ones are adjacent.[6]
  • The Catalan numbers are 1n+1(2nn).
  • The binomial distribution in statistics is (nk)pk(1p)nk.

Binomial coefficients as polynomials

For any nonnegative integer k, the expression (tk) can be written as a polynomial with denominator Template:Math:

(tk)=tk_k!=t(t1)(t2)(tk+1)k(k1)(k2)21;

this presents a polynomial in t with rational coefficients.

As such, it can be evaluated at any real or complex number t to define binomial coefficients with such first arguments. These "generalized binomial coefficients" appear in Newton's generalized binomial theorem.

For each k, the polynomial (tk) can be characterized as the unique degree k polynomial Template:Math satisfying Template:Math and Template:Math.

Its coefficients are expressible in terms of Stirling numbers of the first kind:

(tk)=i=0ks(k,i)tik!.

The derivative of (tk) can be calculated by logarithmic differentiation:

ddt(tk)=(tk)i=0k11ti.

This can cause a problem when evaluated at integers from 0 to t1, but using identities below we can compute the derivative as:

ddt(tk)=i=0k1(1)ki1ki(ti).

Binomial coefficients as a basis for the space of polynomials

Over any field of characteristic 0 (that is, any field that contains the rational numbers), each polynomial p(t) of degree at most d is uniquely expressible as a linear combination k=0dak(tk) of binomial coefficients, because the binomial coefficients consist of one polynomial of each degree. The coefficient ak is the kth difference of the sequence p(0), p(1), ..., p(k). Explicitly,[7] Template:NumBlk2

Integer-valued polynomials

Template:Main Each polynomial (tk) is integer-valued: it has an integer value at all integer inputs t. (One way to prove this is by induction on k using Pascal's identity.) Therefore, any integer linear combination of binomial coefficient polynomials is integer-valued too. Conversely, (Template:EquationNote) shows that any integer-valued polynomial is an integer linear combination of these binomial coefficient polynomials. More generally, for any subring R of a characteristic 0 field K, a polynomial in K[t] takes values in R at all integers if and only if it is an R-linear combination of binomial coefficient polynomials.

Example

The integer-valued polynomial Template:Math can be rewritten as

9(t2)+6(t1)+0(t0).

Identities involving binomial coefficients

The factorial formula facilitates relating nearby binomial coefficients. For instance, if k is a positive integer and n is arbitrary, then Template:NumBlk2 and, with a little more work,

(n1k)(n1k1)=n2kn(nk).

We can also get

(n1k)=nkn(nk).

Moreover, the following may be useful:

(nk)(kj)=(nj)(njkj)=(nkj)(nk+jj).

For constant n, we have the following recurrence:

(nk)=nk+1k(nk1).

To sum up, we have

(nk)=(nnk)=nk+1k(nk1)=nnk(n1k)
=nk(n1k1)=nn2k((n1k)(n1k1))=(n1k)+(n1k1).

Sums of the binomial coefficients

The formula Template:NumBlk2 says that the elements in the Template:Mvarth row of Pascal's triangle always add up to 2 raised to the Template:Mvarth power. This is obtained from the binomial theorem (Template:EquationNote) by setting Template:Math and Template:Math. The formula also has a natural combinatorial interpretation: the left side sums the number of subsets of {1, ..., n} of sizes k = 0, 1, ..., n, giving the total number of subsets. (That is, the left side counts the power set of {1, ..., n}.) However, these subsets can also be generated by successively choosing or excluding each element 1, ..., n; the n independent binary choices (bit-strings) allow a total of 2n choices. The left and right sides are two ways to count the same collection of subsets, so they are equal.

The formulas Template:NumBlk2 and

k=0nk2(nk)=(n+n2)2n2

follow from the binomial theorem after differentiating with respect to Template:Mvar (twice for the latter) and then substituting Template:Math.

The Chu–Vandermonde identity, which holds for any complex values m and n and any non-negative integer k, is Template:NumBlk2 and can be found by examination of the coefficient of xk in the expansion of Template:Math using equation (Template:EquationNote). When Template:Math, equation (Template:EquationNote) reduces to equation (Template:EquationNote). In the special case Template:Math, using (Template:EquationNote), the expansion (Template:EquationNote) becomes (as seen in Pascal's triangle at right) Template:Image frame Template:NumBlk2 where the term on the right side is a central binomial coefficient.

Another form of the Chu–Vandermonde identity, which applies for any integers j, k, and n satisfying Template:Math, is Template:NumBlk2 The proof is similar, but uses the binomial series expansion (Template:EquationNote) with negative integer exponents. When Template:Math, equation (Template:EquationNote) gives the hockey-stick identity

m=kn(mk)=(n+1k+1)

and its relative

r=0m(n+rr)=(n+m+1m).

Let F(n) denote the n-th Fibonacci number. Then

k=0n/2(nkk)=F(n+1).

This can be proved by induction using (Template:EquationNote) or by Zeckendorf's representation. A combinatorial proof is given below.

Multisections of sums

For integers s and t such that 0t<s, series multisection gives the following identity for the sum of binomial coefficients:

(nt)+(nt+s)+(nt+2s)+=1sj=0s1(2cosπjs)ncosπ(n2t)js.

For small Template:Mvar, these series have particularly nice forms; for example,[8]

(n0)+(n3)+(n6)+=13(2n+2cosnπ3)
(n1)+(n4)+(n7)+=13(2n+2cos(n2)π3)
(n2)+(n5)+(n8)+=13(2n+2cos(n4)π3)
(n0)+(n4)+(n8)+=12(2n1+2n2cosnπ4)
(n1)+(n5)+(n9)+=12(2n1+2n2sinnπ4)
(n2)+(n6)+(n10)+=12(2n12n2cosnπ4)
(n3)+(n7)+(n11)+=12(2n12n2sinnπ4)

Partial sums

Although there is no closed formula for partial sums

j=0k(nj)

of binomial coefficients,[9] one can again use (Template:EquationNote) and induction to show that for Template:Math,

j=0k(1)j(nj)=(1)k(n1k),

with special case[10]

j=0n(1)j(nj)=0

for Template:Math. This latter result is also a special case of the result from the theory of finite differences that for any polynomial P(x) of degree less than n,[11]

j=0n(1)j(nj)P(j)=0.

Differentiating (Template:EquationNote) k times and setting x = −1 yields this for P(x)=x(x1)(xk+1), when 0 ≤ k < n, and the general case follows by taking linear combinations of these.

When P(x) is of degree less than or equal to n, Template:NumBlk2 where an is the coefficient of degree n in P(x).

More generally for (Template:EquationNote),

j=0n(1)j(nj)P(m+(nj)d)=dnn!an

where m and d are complex numbers. This follows immediately applying (Template:EquationNote) to the polynomial Template:Tmath instead of Template:Tmath, and observing that Template:Tmath still has degree less than or equal to n, and that its coefficient of degree n is dnan.

The series k1kj=01(j+xk)=1(x1k1) is convergent for k ≥ 2. This formula is used in the analysis of the German tank problem. It follows from k1kj=0M1(j+xk)=1(x1k1)1(M+xk1) which is proved by induction on M.

Identities with combinatorial proofs

Many identities involving binomial coefficients can be proved by combinatorial means. For example, for nonnegative integers nq, the identity

k=qn(nk)(kq)=2nq(nq)

(which reduces to (Template:EquationNote) when q = 1) can be given a double counting proof, as follows. The left side counts the number of ways of selecting a subset of [n] = {1, 2, ..., n} with at least q elements, and marking q elements among those selected. The right side counts the same thing, because there are (nq) ways of choosing a set of q elements to mark, and 2nq to choose which of the remaining elements of [n] also belong to the subset.

In Pascal's identity

(nk)=(n1k1)+(n1k),

both sides count the number of k-element subsets of [n]: the two terms on the right side group them into those that contain element n and those that do not.

The identity (Template:EquationNote) also has a combinatorial proof. The identity reads

k=0n(nk)2=(2nn).

Suppose you have 2n empty squares arranged in a row and you want to mark (select) n of them. There are (2nn) ways to do this. On the other hand, you may select your n squares by selecting k squares from among the first n and nk squares from the remaining n squares; any k from 0 to n will work. This gives

k=0n(nk)(nnk)=(2nn).

Now apply (Template:EquationNote) to get the result.

If one denotes by Template:Math the sequence of Fibonacci numbers, indexed so that Template:Math, then the identity k=0n2(nkk)=F(n) has the following combinatorial proof.[12] One may show by induction that Template:Math counts the number of ways that a Template:Math strip of squares may be covered by Template:Math and Template:Math tiles. On the other hand, if such a tiling uses exactly Template:Mvar of the Template:Math tiles, then it uses Template:Math of the Template:Math tiles, and so uses Template:Math tiles total. There are (nkk) ways to order these tiles, and so summing this coefficient over all possible values of Template:Mvar gives the identity.

Sum of coefficients row

Template:See also

The number of k-combinations for all k, 0kn(nk)=2n, is the sum of the nth row (counting from 0) of the binomial coefficients. These combinations are enumerated by the 1 digits of the set of base 2 numbers counting from 0 to 2n1, where each digit position is an item from the set of n.

Dixon's identity

Dixon's identity is

k=aa(1)k(2ak+a)3=(3a)!(a!)3

or, more generally,

k=aa(1)k(a+ba+k)(b+cb+k)(c+ac+k)=(a+b+c)!a!b!c!,

where a, b, and c are non-negative integers.

Continuous identities

Certain trigonometric integrals have values expressible in terms of binomial coefficients: For any m,n,

ππcos((2mn)x)cosn(x) dx=π2n1(nm)
ππsin((2mn)x)sinn(x) dx={(1)m+(n+1)/2π2n1(nm),n odd0,otherwise
ππcos((2mn)x)sinn(x) dx={(1)m+(n/2)π2n1(nm),n even0,otherwise

These can be proved by using Euler's formula to convert trigonometric functions to complex exponentials, expanding using the binomial theorem, and integrating term by term.

Congruences

If n is prime, then (n1k)(1)kmodn for every k with 0kn1. More generally, this remains true if n is any number and k is such that all the numbers between 1 and k are coprime to n.

Indeed, we have

(n1k)=(n1)(n2)(nk)12k=i=1kniii=1kii=(1)kmodn.

Generating functions

Ordinary generating functions

For a fixed Template:Mvar, the ordinary generating function of the sequence (n0),(n1),(n2), is

k=0(nk)xk=(1+x)n.

For a fixed Template:Mvar, the ordinary generating function of the sequence (0k),(1k),(2k),, is

n=0(nk)yn=yk(1y)k+1.

The bivariate generating function of the binomial coefficients is

n=0k=0n(nk)xkyn=11yxy.

A symmetric bivariate generating function of the binomial coefficients is

n=0k=0(n+kk)xkyn=11xy.

which is the same as the previous generating function after the substitution xxy.

Exponential generating function

A symmetric exponential bivariate generating function of the binomial coefficients is:

n=0k=0(n+kk)xkyn(n+k)!=ex+y.

Divisibility properties

Template:Main In 1852, Kummer proved that if m and n are nonnegative integers and p is a prime number, then the largest power of p dividing (m+nm) equals pc, where c is the number of carries when m and n are added in base p. Equivalently, the exponent of a prime p in (nk) equals the number of nonnegative integers j such that the fractional part of k/pj is greater than the fractional part of n/pj. It can be deduced from this that (nk) is divisible by n/gcd(n,k). In particular therefore it follows that p divides (prs) for all positive integers r and s such that Template:Math. However this is not true of higher powers of p: for example 9 does not divide (96).

A somewhat surprising result by David Singmaster (1974) is that any integer divides almost all binomial coefficients. More precisely, fix an integer d and let f(N) denote the number of binomial coefficients (nk) with n < N such that d divides (nk). Then

limNf(N)N(N+1)/2=1.

Since the number of binomial coefficients (nk) with n < N is N(N + 1) / 2, this implies that the density of binomial coefficients divisible by d goes to 1.

Binomial coefficients have divisibility properties related to least common multiples of consecutive integers. For example:[13]

(n+kk) divides lcm(n,n+1,,n+k)n.
(n+kk) is a multiple of lcm(n,n+1,,n+k)nlcm((k0),(k1),,(kk)).

Another fact: An integer Template:Math is prime if and only if all the intermediate binomial coefficients

(n1),(n2),,(nn1)

are divisible by n.

Proof: When p is prime, p divides

(pk)=p(p1)(pk+1)k(k1)1 for all Template:Math

because (pk) is a natural number and p divides the numerator but not the denominator. When n is composite, let p be the smallest prime factor of n and let Template:Math. Then Template:Math and

(np)=n(n1)(n2)(np+1)p!=k(n1)(n2)(np+1)(p1)!≢0(modn)

otherwise the numerator Template:Math has to be divisible by Template:Math, this can only be the case when Template:Math is divisible by p. But n is divisible by p, so p does not divide Template:Math and because p is prime, we know that p does not divide Template:Math and so the numerator cannot be divisible by n.

Bounds and asymptotic formulas

The following bounds for (nk) hold for all values of n and k such that Template:Math: nkkk(nk)nkk!<(nek)k. The first inequality follows from the fact that (nk)=nkn1k1n(k1)1 and each of these k terms in this product is nk. A similar argument can be made to show the second inequality. The final strict inequality is equivalent to ek>kk/k!, that is clear since the RHS is a term of the exponential series ek=j=0kj/j!.

From the divisibility properties we can infer that lcm(nk,,n)(nk)lcm((k0),,(kk))(nk)lcm(nk,,n)nk, where both equalities can be achieved.[13]

The following bounds are useful in information theory:[14]Template:Rp 1n+12nH(k/n)(nk)2nH(k/n) where H(p)=plog2(p)(1p)log2(1p) is the binary entropy function. It can be further tightened to n8k(nk)2nH(k/n)(nk)n2πk(nk)2nH(k/n) for all 1kn1.[15]Template:Rp

Both n and k large

Stirling's approximation yields the following approximation, valid when nk,k both tend to infinity: (nk)n2πk(nk)nnkk(nk)nk Because the inequality forms of Stirling's formula also bound the factorials, slight variants on the above asymptotic approximation give exact bounds. In particular, when n is sufficiently large, one has (2nn)22nnπ and n(2nn)22n1. More generally, for Template:Math and Template:Math (again, by applying Stirling's formula to the factorials in the binomial coefficient), n(mnn)mm(n1)+1(m1)(m1)(n1).

If n is large and k is linear in n, various precise asymptotic estimates exist for the binomial coefficient (nk). For example, if |n/2k|=o(n2/3) then (nk)(nn2)ed2/(2n)2n12nπed2/(2n) where d = n − 2k.[16]

Template:Mvar much larger than Template:Mvar

If Template:Mvar is large and Template:Mvar is Template:Math (that is, if Template:Math), then (nk)(nek)k(2πk)1/2exp(k22n(1+o(1))) where again Template:Mvar is the little o notation.[17]

Sums of binomial coefficients

A simple and rough upper bound for the sum of binomial coefficients can be obtained using the binomial theorem: i=0k(ni)i=0kni1ki(1+n)k More precise bounds are given by 18nε(1ε)2H(ε)ni=0k(ni)2H(ε)n, valid for all integers n>k1 with εk/n1/2.[18]

Generalized binomial coefficients

The infinite product formula for the gamma function also gives an expression for binomial coefficients (1)k(zk)=(z+k1k)=1Γ(z)1(k+1)z+1j=k+1(1+1j)z11z+1j which yields the asymptotic formulas (zk)(1)kΓ(z)kz+1and(z+kk)=kzΓ(z+1)(1+z(z+1)2k+𝒪(k2)) as k.

This asymptotic behaviour is contained in the approximation (z+kk)ez(Hkγ)Γ(z+1) as well. (Here Hk is the k-th harmonic number and γ is the Euler–Mascheroni constant.)

Further, the asymptotic formula (z+kj)(kj)(1jk)zand(jjk)(jzjk)(jk)z hold true, whenever k and j/kx for some complex number x.

Generalizations

Generalization to multinomials

Template:Main Binomial coefficients can be generalized to multinomial coefficients defined to be the number:

(nk1,k2,,kr)=n!k1!k2!kr!

where

i=1rki=n.

While the binomial coefficients represent the coefficients of Template:Math, the multinomial coefficients represent the coefficients of the polynomial

(x1+x2++xr)n.

The case r = 2 gives binomial coefficients:

(nk1,k2)=(nk1,nk1)=(nk1)=(nk2).

The combinatorial interpretation of multinomial coefficients is distribution of n distinguishable elements over r (distinguishable) containers, each containing exactly ki elements, where i is the index of the container.

Multinomial coefficients have many properties similar to those of binomial coefficients, for example the recurrence relation:

(nk1,k2,,kr)=(n1k11,k2,,kr)+(n1k1,k21,,kr)++(n1k1,k2,,kr1)

and symmetry:

(nk1,k2,,kr)=(nkσ1,kσ2,,kσr)

where (σi) is a permutation of (1, 2, ..., r).

Taylor series

Using Stirling numbers of the first kind the series expansion around any arbitrarily chosen point z0 is

(zk)=1k!i=0kzisk,i=i=0k(zz0)ij=ik(z0ji)sk+ij,i(k+ij)!=i=0k(zz0)ij=ikz0ji(ji)sk,jk!.

Binomial coefficient with Template:Math

The definition of the binomial coefficients can be extended to the case where n is real and k is integer.

In particular, the following identity holds for any non-negative integer k:

(1/2k)=(2kk)(1)k+122k(2k1).

This shows up when expanding 1+x into a power series using the Newton binomial series :

1+x=k0(1/2k)xk.

Products of binomial coefficients

One can express the product of two binomial coefficients as a linear combination of binomial coefficients:

(zm)(zn)=k=0min(m,n)(m+nkk,mk,nk)(zm+nk),

where the connection coefficients are multinomial coefficients. In terms of labelled combinatorial objects, the connection coefficients represent the number of ways to assign Template:Math labels to a pair of labelled combinatorial objects—of weight m and n respectively—that have had their first k labels identified, or glued together to get a new labelled combinatorial object of weight Template:Math. (That is, to separate the labels into three portions to apply to the glued part, the unglued part of the first object, and the unglued part of the second object.) In this regard, binomial coefficients are to exponential generating series what falling factorials are to ordinary generating series.

The product of all binomial coefficients in the nth row of the Pascal triangle is given by the formula:

k=0n(nk)=k=1nk2kn1.

Partial fraction decomposition

The partial fraction decomposition of the reciprocal is given by

1(zn)=i=0n1(1)n1i(ni)nizi,1(z+nn)=i=1n(1)i1(ni)iz+i.

Newton's binomial series

Template:Main Newton's binomial series, named after Sir Isaac Newton, is a generalization of the binomial theorem to infinite series:

(1+z)α=n=0(αn)zn=1+(α1)z+(α2)z2+.

The identity can be obtained by showing that both sides satisfy the differential equation Template:Math.

The radius of convergence of this series is 1. An alternative expression is

1(1z)α+1=n=0(n+αn)zn

where the identity

(nk)=(1)k(kn1k)

is applied.

Multiset (rising) binomial coefficient

Template:Main Binomial coefficients count subsets of prescribed size from a given set. A related combinatorial problem is to count multisets of prescribed size with elements drawn from a given set, that is, to count the number of ways to select a certain number of elements from a given set with the possibility of selecting the same element repeatedly. The resulting numbers are called multiset coefficients;[19] the number of ways to "multichoose" (i.e., choose with replacement) k items from an n element set is denoted ((nk)).

To avoid ambiguity and confusion with n's main denotation in this article,
let Template:Math and Template:Math.

Multiset coefficients may be expressed in terms of binomial coefficients by the rule (fk)=((rk))=(r+k1k). One possible alternative characterization of this identity is as follows: We may define the falling factorial as (f)k=fk_=(fk+1)(f3)(f2)(f1)f, and the corresponding rising factorial as r(k)=rk=r(r+1)(r+2)(r+3)(r+k1); so, for example, 1718192021=(21)5=215_=175=17(5). Then the binomial coefficients may be written as (fk)=(f)kk!=(fk+1)(f2)(f1)f12345k, while the corresponding multiset coefficient is defined by replacing the falling with the rising factorial: ((rk))=r(k)k!=r(r+1)(r+2)(r+k1)12345k.

Generalization to negative integers n

Template:Pascal triangle extended.svg For any n,

(nk)=n(n+1)(n+k2)(n+k1)k!=(1)kn(n+1)(n+2)(n+k1)k!=(1)k(n+k1k)=(1)k((nk)).

In particular, binomial coefficients evaluated at negative integers n are given by signed multiset coefficients. In the special case n=1, this reduces to (1)k=(1k)=((kk)).

For example, if n = −4 and k = 7, then r = 4 and f = 10:

(47)=109876541234567=(1)7456789101234567=((77))((47))=(17)(107).

Two real or complex valued arguments

The binomial coefficient is generalized to two real or complex valued arguments using the gamma function or beta function via

(xy)=Γ(x+1)Γ(y+1)Γ(xy+1)=1(x+1)B(y+1,xy+1).

This definition inherits these following additional properties from Γ:

(xy)=sin(yπ)sin(xπ)(y1x1)=sin((xy)π)sin(xπ)(yx1y);

moreover,

(xy)(yx)=sin((xy)π)(xy)π.

The resulting function has been little-studied, apparently first being graphed in Template:Harv. Notably, many binomial identities fail: (nm)=(nnm) but (nm)(nnm) for n positive (so n negative). The behavior is quite complex, and markedly different in various octants (that is, with respect to the x and y axes and the line y=x), with the behavior for negative x having singularities at negative integer values and a checkerboard of positive and negative regions:

  • in the octant 0yx it is a smoothly interpolated form of the usual binomial, with a ridge ("Pascal's ridge").
  • in the octant 0xy and in the quadrant x0,y0 the function is close to zero.
  • in the quadrant x0,y0 the function is alternatingly very large positive and negative on the parallelograms with vertices (n,m+1),(n,m),(n1,m1),(n1,m)
  • in the octant 0>x>y the behavior is again alternatingly very large positive and negative, but on a square grid.
  • in the octant 1>y>x+1 it is close to zero, except for near the singularities.

Generalization to q-series

The binomial coefficient has a q-analog generalization known as the Gaussian binomial coefficient.

Generalization to infinite cardinals

The definition of the binomial coefficient can be generalized to infinite cardinals by defining:

(αβ)=|{BA:|B|=β}|

where Template:Mvar is some set with cardinality α. One can show that the generalized binomial coefficient is well-defined, in the sense that no matter what set we choose to represent the cardinal number α, (αβ) will remain the same. For finite cardinals, this definition coincides with the standard definition of the binomial coefficient.

Assuming the Axiom of Choice, one can show that (αα)=2α for any infinite cardinal α.

See also

Template:Div col

Template:Div col end

Notes

Template:Reflist

References

Template:Refbegin

Template:Refend

Template:PlanetMath attribution

  1. Template:Harvtxt
  2. Lilavati Section 6, Chapter 4 (see Template:Harvtxt).
  3. Template:Harvnb
  4. See Template:Harv, which also defines (nk)=0 for k<0. Alternative generalizations, such as to two real or complex valued arguments using the Gamma function assign nonzero values to (nk) for k<0, but this causes most binomial coefficient identities to fail, and thus is not widely used by the majority of definitions. One such choice of nonzero values leads to the aesthetically pleasing "Pascal windmill" in Hilton, Holton and Pedersen, Mathematical reflections: in a room with many mirrors, Springer, 1997, but causes even Pascal's identity to fail (at the origin).
  5. When α=n is a nonnegative integer, (nk)=0 for k>n because the (k=n+1)-th factor of the numerator is n(n+1)+1=0. Thus, the k-th term is a zero product for all kn+1.
  6. Template:Cite journal
  7. This can be seen as a discrete analog of Taylor's theorem. It is closely related to Newton's polynomial. Alternating sums of this form may be expressed as the Nörlund–Rice integral.
  8. Template:Harvtxt.
  9. Template:Citation.
  10. see induction developed in eq (7) p. 1389 in Template:Citation.
  11. Template:Cite journal
  12. Template:Harvnb
  13. 13.0 13.1 Template:Cite journal
  14. Template:Cite book
  15. Template:Cite book
  16. Template:Cite book
  17. Template:Cite book
  18. see e.g. Template:Harvtxt or Template:Harvtxt.
  19. Template:Citation.