# Combination

In mathematics, a

**combination**is a selection of items from a collection, such that (unlike permutations) the order of selection does not matter. For example, given three fruits, say an apple, an orange and a pear, there are three combinations of two that can be drawn from this set: an apple and a pear; an apple and an orange; or a pear and an orange. More formally, a

*k*-

**combination**of a set

*S*is a subset of

*k*distinct elements of

*S*. If the set has

*n*elements, the number of

*k*-combinations is equal to the binomial coefficient

:$inom\; nk\; =\; frac\{n(n-1)dotsb(n-k+1)\}\{k(k-1)dotsb1\},$

which can be written using factorials as $extstylefrac\{n!\}\{k!(n-k)!\}$ whenever $kleq\; n$, and which is zero when $k>n$. The set of all

*k*-combinations of a set

*S*is often denoted by $extstyleinom\; Sk$.

Combinations refer to the combination of

*n*things taken

*k*at a time without repetition. To refer to combinations in which repetition is allowed, the terms

*k*-selection,

*k*-multiset, or

*k*-combination with repetition are often used. If, in the above example, it were possible to have two of any one kind of fruit there would be 3 more 2-selections: one with two apples, one with two oranges, and one with two pears.

Although the set of three fruits was small enough to write a complete list of combinations, with large sets this becomes impractical. For example, a poker hand can be described as a 5-combination (

*k*= 5) of cards from a 52 card deck (

*n*= 52). The 5 cards of the hand are all distinct, and the order of cards in the hand does not matter. There are 2,598,960 such combinations, and the chance of drawing any one hand at random is 1 / 2,598,960.

## Number of *k*-combinations

The number of *k*-combinations from a given set

*S*of

*n*elements is often denoted in elementary combinatorics texts by $C(n,k)$, or by a variation such as $C^n\_k$, $\{\}\_nC\_k$, $\{\}^nC\_k$, $C\_\{n,k\}$ or even $C\_n^k$ (the latter form was standard in French, Romanian, Russian, Chinese and Polish texts). The same number however occurs in many other mathematical contexts, where it is denoted by $binom\; nk$ (often read as "

*n*choose

*k*"); notably it occurs as a coefficient in the binomial formula, hence its name

**binomial coefficient**. One can define $binom\; nk$ for all natural numbers

*k*at once by the relation

:$(1\; +\; X)^n\; =\; sum\_\{kgeq0\}inom\{n\}\{k\}\; X^k,$

from which it is clear that

:$inom\{n\}\{0\}\; =\; inom\{n\}\{n\}\; =\; 1,$

and further,

:$inom\{n\}\{k\}\; =\; 0$ for

*k*>

*n*.

To see that these coefficients count

*k*-combinations from

*S*, one can first consider a collection of

*n*distinct variables

*X*

_{s}labeled by the elements

*s*of

*S*, and expand the product over all elements of

*S*:

:$prod\_\{sin\; S\}(1+X\_s);$

it has 2 distinct terms corresponding to all the subsets of

*S*, each subset giving the product of the corresponding variables

*X*

_{s}. Now setting all of the

*X*

_{s}equal to the unlabeled variable

*X*, so that the product becomes (1 +

*X*), the term for each

*k*-combination from

*S*becomes

*X*, so that the coefficient of that power in the result equals the number of such

*k*-combinations.

Binomial coefficients can be computed explicitly in various ways. To get all of them for the expansions up to (1 +

*X*), one can use (in addition to the basic cases already given) the recursion relation

:$inom\{n\}\{k\}\; =\; inom\{n\; -\; 1\}\{k\; -\; 1\}\; +\; inom\{n\; -\; 1\}\{k\},$

for 0 <

*k*<

*n*, which follows from (1 +

*X*) = (1 +

*X*)(1 +

*X*); this leads to the construction of Pascal's triangle.

For determining an individual binomial coefficient, it is more practical to use the formula

:$inom\; nk\; =\; frac\{n(n-1)(n-2)cdots(n-k+1)\}\{k!\}$.

The numerator gives the number of

*k*-permutations of

*n*, i.e., of sequences of

*k*distinct elements of

*S*, while the denominator gives the number of such

*k*-permutations that give the same

*k*-combination when the order is ignored.

When

*k*exceeds

*n*/2, the above formula contains factors common to the numerator and the denominator, and canceling them out gives the relation

:$inom\; nk\; =\; inom\; n\{n-k\},$

for 0 ≤

*k*≤

*n*. This expresses a symmetry that is evident from the binomial formula, and can also be understood in terms of

*k*-combinations by taking the complement of such a combination, which is an (

*n*−

*k*)-combination.

Finally there is a formula which exhibits this symmetry directly, and has the merit of being easy to remember:

:$inom\; nk\; =\; frac\{n!\}\{k!(n-k)!\},$

where

*n*

*n*. It is obtained from the previous formula by multiplying denominator and numerator by (

*n*−

*k*)!, so it is certainly inferior as a method of computation to that formula.

The last formula can be understood directly, by considering the

*n*

*S*. Each such permutation gives a

*k*-combination by selecting its first

*k*elements. There are many duplicate selections: any combined permutation of the first

*k*elements among each other, and of the final (

*n*−

*k*) elements among each other produces the same combination; this explains the division in the formula.

From the above formulas follow relations between adjacent numbers in Pascal's triangle in all three directions:

:$inom\; nk\; =\; egin\{cases\}\; inom\; n\{k-1\}\; frac\; \{n-k+1\}k\; quad\; ext\{if\; \}\; k\; 0\; \backslash \; inom\; \{n-1\}k\; frac\; n\{n-k\}\; quad\; ext\{if\; \}\; k\; n\; \backslash \; inom\; \{n-1\}\{k-1\}\; frac\; nk\; quad\; ext\{if\; \}\; n,\; k\; 0\; end\{cases\}$.

Together with the basic cases $binom\; n0=1=\; binom\; nn$, these allow successive computation of respectively all numbers of combinations from the same set (a row in Pascal's triangle), of

*k*-combinations of sets of growing sizes, and of combinations with a complement of fixed size

*n*−

*k*.

### Example of counting combinations

As a specific example, one can compute the number of five-card hands possible from a standard fifty-two card deck as::$\{52\; choose\; 5\}\; =\; frac\{52\; imes51\; imes50\; imes49\; imes48\}\{5\; imes4\; imes3\; imes2\; imes1\}\; =\; frac\{311\{,\}875\{,\}200\}\{120\}\; =\; 2\{,\}598\{,\}960.$

Alternatively one may use the formula in terms of factorials and cancel the factors in the numerator against parts of the factors in the denominator, after which only multiplication of the remaining factors is required: :$egin\{alignat\}\{2\}\; \{52\; choose\; 5\}\; =\; frac\{52!\}\{5!47!\}\; \backslash [5pt]\; =\; frac\{52\; imes51\; imes50\; imes49\; imes48\; imescancel\{47!\}\}\{5\; imes4\; imes3\; imes2\; imescancel\{1\}\; imescancel\{47!\}\}\; \backslash [5pt]\; =\; frac\{52\; imes51\; imes50\; imes49\; imes48\}\{5\; imes4\; imes3\; imes2\}\; \backslash [5pt]\; =\; frac\{(26\; imescancel\{2\})\; imes(17\; imescancel\{3\})\; imes(10\; imescancel\{5\})\; imes49\; imes(12\; imescancel\{4\})\}\{cancel\{5\}\; imescancel\{4\}\; imescancel\{3\}\; imescancel\{2\}\}\; \backslash [5pt]\; =\; \{26\; imes17\; imes10\; imes49\; imes12\}\; \backslash [5pt]\; =\; 2\{,\}598\{,\}960.\; end\{alignat\}$

Another alternative computation, equivalent to the first, is based on writing

:$\{n\; choose\; k\}\; =\; frac\; \{\; (\; n\; -\; 0\; )\; \}1\; imes\; frac\; \{\; (\; n\; -\; 1\; )\; \}2\; imes\; frac\; \{\; (\; n\; -\; 2\; )\; \}3\; imes\; cdots\; imes\; frac\; \{\; (\; n\; -\; (k\; -\; 1)\; )\; \}k,$

which gives

:$\{52\; choose\; 5\}\; =\; frac\{52\}1\; imes\; frac\{51\}2\; imes\; frac\{50\}3\; imes\; frac\{49\}4\; imes\; frac\{48\}5\; =\; 2\{,\}598\{,\}960$. When evaluated in the following order, , this can be computed using only integer arithmetic. The reason is that when each division occurs, the intermediate result that is produced is itself a binomial coefficient, so no remainders ever occur.

Using the symmetric formula in terms of factorials without performing simplifications gives a rather extensive calculation:

:$egin\{align\}\; \{52\; choose\; 5\}\; =\; frac\{n!\}\{k!(n-k)!\}\; =\; frac\{52!\}\{5!(52-5)!\}\; =\; frac\{52!\}\{5!47!\}\; \backslash [6pt]\; =\; frac\{80,658,175,170,943,878,571,660,636,856,403,766,975,289,505,440,883,277,824,000,000,000,000\}\{120\; imes258,623,241,511,168,180,642,964,355,153,611,979,969,197,632,389,120,000,000,000\}\; \backslash [6pt]\; =\; 2\{,\}598\{,\}960.\; end\{align\}$

### Enumerating *k*-combinations

One can enumerate all

*k*-combinations of a given set

*S*of

*n*elements in some fixed order, which establishes a bijection from an interval of $binom\; nk$ integers with the set of those

*k*-combinations. Assuming

*S*is itself ordered, for instance

*S*= { 1, 2, …,

*n*}, there are two natural possibilities for ordering its

*k*-combinations: by comparing their smallest elements first (as in the illustrations above) or by comparing their largest elements first. The latter option has the advantage that adding a new largest element to

*S*will not change the initial part of the enumeration, but just add the new

*k*-combinations of the larger set after the previous ones. Repeating this process, the enumeration can be extended indefinitely with

*k*-combinations of ever larger sets. If moreover the intervals of the integers are taken to start at 0, then the

*k*-combination at a given place

*i*in the enumeration can be computed easily from

*i*, and the bijection so obtained is known as the combinatorial number system. It is also known as "rank"/"ranking" and "unranking" in computational mathematics.

There are many ways to enumerate

*k*combinations. One way is to visit all the binary numbers less than 2. Choose those numbers having

*k*nonzero bits, although this is very inefficient even for small

*n*(e.g.

*n*= 20 would require visiting about one million numbers while the maximum number of allowed

*k*combinations is about 186 thousand for

*k*= 10). The positions of these 1 bits in such a number is a specific

*k*-combination of the set { 1, …,

*n*}. Another simple, faster way is to track

*k*index numbers of the elements selected, starting with {0 ..

*k*−1} (zero-based) or {1 ..

*k*} (one-based) as the first allowed

*k*-combination and then repeatedly moving to the next allowed

*k*-combination by incrementing the last index number if it is lower than

*n*-1 (zero-based) or

*n*(one-based) or the last index number

*x*that is less than the index number following it minus one if such an index exists and resetting the index numbers after

*x*to {

*x*+1,

*x*+2, …}.

## Number of combinations with repetition

A*k*-

**combination with repetitions**, or

*k*-

**multicombination**, or

**multisubset**of size

*k*from a set

*S*is given by a sequence of

*k*not necessarily distinct elements of

*S*, where order is not taken into account: two sequences define the same multiset if one can be obtained from the other by permuting the terms. In other words, the number of ways to sample

*k*elements from a set of

*n*elements allowing for duplicates (i.e., with replacement) but disregarding different orderings (e.g. {2,1,2} = {1,2,2}). Associate an index to each element of

*S*and think of the elements of

*S*as

*types*of objects, then we can let $x\_i$ denote the number of elements of type

*i*in a multisubset. The number of multisubsets of size

*k*is then the number of nonnegative integer solutions of the Diophantine equation: :$x\_1\; +\; x\_2\; +\; ldots\; +\; x\_n\; =\; k.$ If

*S*has

*n*elements, the number of such

*k*-multisubsets is denoted by,

::$left(!!inom\{n\}\{k\}!!\; ight),$ a notation that is analogous to the binomial coefficient which counts

*k*-subsets. This expression,

*n*multichoose

*k*, can also be given in terms of binomial coefficients:

:$left(!!inom\{n\}\{k\}!!\; ight)=inom\{n+k-1\}\{k\}.$

This relationship can be easily proved using a representation known as stars and bars.

A solution of the above Diophantine equation can be represented by $x\_1$

*stars*, a separator (a

*bar*), then $x\_2$ more stars, another separator, and so on. The total number of stars in this representation is

*k*and the number of bars is

*n*- 1 (since no separator is needed at the very end). Thus, a string of

*k*+

*n*- 1 symbols (stars and bars) corresponds to a solution if there are

*k*stars in the string. Any solution can be represented by choosing

*k*out of positions to place stars and filling the remaining positions with bars. For example, the solution $x\_1\; =\; 3,\; x\_2\; =\; 2,\; x\_3\; =\; 0,\; x\_4\; =\; 5$ of the equation $x\_1\; +\; x\_2\; +\; x\_3\; +\; x\_4\; =\; 10$ can be represented by

:$left(!!inom\{n\}\{k\}!!\; ight)=left(!!inom\{k+1\}\{n-1\}!!\; ight).$ This identity follows from interchanging the stars and bars in the above representation.

### Example of counting multisubsets

For example, if you have four types of donuts (*n*= 4) on a menu to choose from and you want three donuts (

*k*= 3), the number of ways to choose the donuts with repetition can be calculated as

:$left(!!inom\{4\}\{3\}!!\; ight)\; =\; inom\{4+3-1\}3\; =\; inom\{6\}\{3\}\; =\; frac\{6\; imes5\; imes4\}\{3\; imes2\; imes1\}\; =\; 20.$

This result can be verified by listing all the 3-multisubsets of the set

*S*= {1,2,3,4}. This is displayed in the following table. The second column shows the nonnegative integer solutions $[x\_1,x\_2,x\_3,x\_4]$ of the equation $x\_1\; +\; x\_2\; +\; x\_3\; +\; x\_4\; =\; 3$ and the last column gives the stars and bars representation of the solutions.

## Number of *k*-combinations for all *k*

The number of *k*-combinations for all

*k*is the number of subsets of a set of

*n*elements. There are several ways to see that this number is 2. In terms of combinations, $sum\_\{0leq\{k\}leq\{n\}\}inom\; nk\; =\; 2^n$, which is the sum of the

*n*th row (counting from 0) of the binomial coefficients in Pascal's triangle. These combinations (subsets) are enumerated by the 1 digits of the set of base 2 numbers counting from 0 to 2 − 1, where each digit position is an item from the set of

*n*.

Given 3 cards numbered 1 to 3, there are 8 distinct combinations (subsets), including the empty set:

:$|\; \{\; \{\}\; ;\; \{1\}\; ;\; \{2\}\; ;\; \{1,\; 2\}\; ;\; \{3\}\; ;\; \{1,\; 3\}\; ;\; \{2,\; 3\}\; ;\; \{1,\; 2,\; 3\}\; \}|\; =\; 2^3\; =\; 8$

Representing these subsets (in the same order) as base 2 numerals: ::0 – 000 ::1 – 001 ::2 – 010 ::3 – 011 ::4 – 100 ::5 – 101 ::6 – 110 ::7 – 111

## Probability: sampling a random combination

There are various algorithms to pick out a random combination from a given set or list. Rejection sampling is extremely slow for large sample sizes. One way to select a

*k*-combination efficiently from a population of size

*n*is to iterate across each element of the population, and at each step pick that element with a dynamically changing probability of $frac\{k-\#\; ext\{samples\; chosen\}\}\{n-\; \#\; ext\{samples\; visited\}\}$. (see reservoir sampling).