Binomial Coefficient

DOWNLOAD Mathematica Notebook EXPLORE THIS TOPIC IN the MathWorld Classroom

The binomial coefficient (n; k) is the number of ways of picking k unordered outcomes from n possibilities, also known as a combination or combinatorial number. The symbols _nC_k and (n; k) are used to denote a binomial coefficient, and are sometimes read as "n choose k."

(n; k) therefore gives the number of k-subsets possible out of a set of n distinct items. For example, The 2-subsets of {1,2,3,4} are the six pairs {1,2}, {1,3}, {1,4}, {2,3}, {2,4}, and {3,4}, so (4; 2)=6. The number of lattice paths from the origin (0,0) to a point (a,b) is the binomial coefficient (a+b; a) (Hilton and Pedersen 1991).

The value of the binomial coefficient for nonnegative n and k is given explicitly by

 _nC_k=(n; k)=(n!)/((n-k)!k!),
(1)

where z! denotes a factorial. Writing the factorial as a gamma function z!=Gamma(z+1) allows the binomial coefficient to be generalized to noninteger arguments (including complex x and y) as

 (x; y)=(Gamma(x+1))/(Gamma(y+1)Gamma(x-y+1)).
(2)

For nonnegative integer arguments, the gamma function reduces to factorials, leading to

 (n; k)={(n!)/(k!(n-k)!)   for 0<=k<n; 0   otherwise,
(3)

which is Pascal's triangle. Using the symmetry formula

 (Gamma(s-a+1))/(Gamma(s-b+1))=(-1)^(b-a)(Gamma(b-s))/(Gamma(a-s))
(4)

for integer a, b and complex s, this definition can be extended to negative integer arguments, making it continuous at all integer arguments as well as continuous for all complex arguments except for negative integer x and noninteger y, in which case it is infinite (Kronenburg 2011). This definition, given by

 (n; k)={(-1)^k(-n+k-1; k)   for k>=0; (-1)^(n-k)(-k-1; n-k)   for k<=n; 0   otherwise
(5)

for negative integer n and k is in agreement with the binomial theorem, and with combinatorial identities with a few special exceptions (Kronenburg 2011).

The binomial coefficient is implemented in the Wolfram Language as Binomial[n, k], which follows the above convention starting in Version 8.

BinomialCoefficient

Plotting the binomial coefficient

 (x; y)=(x!)/(y!(x-y)!)
(6)

in the xy-plane (Fowler 1996) gives the beautiful plot shown above, which has a very complicated graph for negative x and y and is therefore difficult to render using standard plotting programs.

For a positive integer n, the binomial theorem gives

 (x+a)^n=sum_(k=0)^n(n; k)x^ka^(n-k).
(7)

The finite difference analog of this identity is known as the Chu-Vandermonde identity. A similar formula holds for negative integers,

 (x+a)^(-n)=sum_(k=0)^infty(-n; k)x^ka^(-n-k).
(8)

There are a number of elegant binomial sums.

The binomial coefficients satisfy the identities

(n; 0)=(n; n)=1
(9)
(n; k)=(n; n-k)
(10)
=(-1)^k(k-n-1; k)
(11)
(n; k+1)=(n; k)(n-k)/(k+1)
(12)
(n+1; k)=(n; k)+(n; k-1).
(13)

The product of binomial coefficients is given by

 product_(k=0)^n(n; k)=(H^2(n))/((n!)^(n+1)),
(14)

where H(n) is a hyperfactorial and n! is a factorial.

As shown by Kummer in 1852, if p^k is the largest power of a prime p that divides (m+n; m), where m and n are nonnegative integers, then k is the number of carries that occur when m is added to n in base p (Graham et al. 1989, Exercise 5.36, p. 245; Ribenboim 1989; Vardi 1991, p. 68). Kummer's result can also be stated in the form that the exponent of a prime p dividing (n; m) is given by the number of integers j>=0 for which

 frac(m/p^j)>frac(n/p^j),
(15)

where frac(x) denotes the fractional part of x. This inequality may be reduced to the study of the exponential sums sum_(n)Lambda(n)e(x/n), where Lambda(n) is the Mangoldt function. Estimates of these sums are given by Jutila (1973, 1974), but recent improvements have been made by Granville and Ramare (1996).

R. W. Gosper showed that

 f(n)=(n-1; 1/2(n-1))=(-1)^((n-1)/2) (mod n)
(16)

for all primes, and conjectured that it holds only for primes. This was disproved when Skiena (1990) found it also holds for the composite number n=3×11×179. Vardi (1991, p. 63) subsequently showed that n=p^2 is a solution whenever p is a Wieferich prime and that if n=p^k with k>3 is a solution, then so is n=p^(k-1). This allowed him to show that the only solutions for composite n<1.3×10^7 are 5907, 1093^2, and 3511^2, where 1093 and 3511 are Wieferich primes.

Consider the binomial coefficients f(n)=(2n-1; n), the first few of which are 1, 3, 10, 35, 126, ... (OEIS A001700). The generating function is

 1/2[1/(sqrt(1-4x))-1]=x+3x^2+10x^3+35x^4+....
(17)

These numbers are squarefree only for n=2, 3, 4, 6, 9, 10, 12, 36, ... (OEIS A046097), with no others known. It turns out that f(n) is divisible by 4 unless n belongs to a 2-automatic set S_2, which happens to be the set of numbers whose binary representations contain at most two 1s: 1, 2, 3, 4, 5, 6, 8, 9, 10, 12, 16, 17, 18, ... (OEIS A048645). Similarly, f(n) is divisible by 9 unless n belongs to a 3-automatic set S_3, consisting of numbers n for which the representation of 2n in ternary consists entirely of 0s and 2s (except possibly for a pair of adjacent 1s). The initial elements of S_3 are 1, 2, 3, 4, 6, 7, 9, 10, 11, 12, 13, 18, 19, 21, 22, 27, ... (OEIS A051382). If f(n) is squarefree, then n must belong to S=S_2 intersection S_3. It is very probable that S is finite, but no proof is known. Now, squares larger than 4 and 9 might also divide f(n), but by eliminating these two alone, the only possible n for n<=2^(64) are 1, 2, 3, 4, 6, 9, 10, 12, 18, 33, 34, 36, 40, 64, 66, 192, 256, 264, 272, 513, 514, 516, 576 768, 1026, 1056, 2304, 16392, 65664, 81920, 532480, and 545259520. All of these but the last have been checked, establishing that there are no other n such that f(n) is squarefree for n<=545259520.

Erdős showed that the binomial coefficient (n; k) with 3<=k<=n/2 is a power of an integer for the single case (50; 3)=140^2 (Le Lionnais 1983, p. 48). Binomial coefficients T_(n-1)=(n; 2) are squares a^2 when a^2 is a triangular number, which occur for a=1, 6, 35, 204, 1189, 6930, ... (OEIS A001109). These values of a have the corresponding values n=2, 9, 50, 289, 1682, 9801, ... (OEIS A052436).

The binomial coefficients (n; |_n/2_|) are called central binomial coefficients, where |_x_| is the floor function, although the subset of coefficients (2n; n) is sometimes also given this name. Erdős and Graham (1980, p. 71) conjectured that the central binomial coefficient (2n; n) is never squarefree for n>4, and this is sometimes known as the Erdős squarefree conjecture. Sárkőzy's theorem (Sárkőzy 1985) provides a partial solution which states that the binomial coefficient (2n; n) is never squarefree for all sufficiently large n>=n_0 (Vardi 1991). Granville and Ramare (1996) proved that the only squarefree values are n=2 and 4. Sander (1992) subsequently showed that (2n+/-d; n) are also never squarefree for sufficiently large n as long as d is not "too big."

For p, q, and r distinct primes, then the function (◇) satisfies

 f(pqr)f(p)f(q)f(r)=f(pq)f(pr)f(qr) (mod pqr)
(18)

(Vardi 1991, p. 66).

Most binomial coefficients (n; k) with n>=2k have a prime factor p<=n/k, and Lacampagne et al. (1993) conjecture that this inequality is true for all n>17.125k, or more strongly that any such binomial coefficient has least prime factor p<=n/k or p<=17 with the exceptions (62; 6), (959; 56), (474; 66), (284; 28) for which p=19, 19, 23, 29 (Guy 1994, p. 84).

The binomial coefficient (m; n) (mod 2) can be computed using the XOR operation n XOR m, making Pascal's triangle mod 2 very easy to construct.

Sondow (2005) and Sondow and Zudilin (2006) noted the inequality

 1/(4rm)[((r+1)^(r+1))/(r^r)]^m<((r+1)m; m)<[((r+1)^(r+1))/(r^r)]^m
(19)

for m a positive integer and r>=1 a real number.

Wolfram Web Resources

Mathematica »

The #1 tool for creating Demonstrations and anything technical.

Wolfram|Alpha »

Explore anything with the first computational knowledge engine.

Wolfram Demonstrations Project »

Explore thousands of free applications across science, mathematics, engineering, technology, business, art, finance, social sciences, and more.

Computerbasedmath.org »

Join the initiative for modernizing math education.

Online Integral Calculator »

Solve integrals with Wolfram|Alpha.

Step-by-step Solutions »

Walk through homework problems step-by-step from beginning to end. Hints help you try the next step on your own.

Wolfram Problem Generator »

Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.

Wolfram Education Portal »

Collection of teaching and learning tools built by Wolfram education experts: dynamic textbook, lesson plans, widgets, interactive Demonstrations, and more.

Wolfram Language »

Knowledge-based programming for everyone.