Multiset

Overview

The number of times an element belongs to the multiset is the multiplicity of that member. The total number of elements in a multiset, including repeated memberships, is the cardinality of the multiset. For example, in the multiset {a, a, b, b, b, c} the multiplicities of the members a, b, and c are respectively 2, 3, and 1, and the cardinality of the multiset is 6. To distinguish between sets and multisets, a notation that incorporates square brackets is sometimes used: the multiset {2, 2, 3} can be represented as [2, 2, 3]. In multisets, as in sets and in contrast to tuples, the order of elements is irrelevant: The multisets {a, a, b} and {a, b, a} are equal.

History

Wayne Blizard traced multisets back to the very origin of numbers, arguing that “in ancient times, the number n was often represented by a collection of n strokes, tally marks, or units.” These and similar collections of objects are multisets as strokes, tally marks, or units are considered indistinguishable. This shows that people implicitly used multisets even before mathematics emerged.

This shows that necessity in this structure has been always so urgent that multisets have been several times rediscovered and appeared in literature under different names. For instance, they were referred to as bags by James Lyle Peterson in 1981, who attributed this term to a 1971 paper by Cerf et al. A multiset has been also called an aggregate, heap, bunch, sample, weighted set, occurrence set, and fireset (finitely repeated element set).

Although multisets were implicitly utilized from ancient times, their explicit exploration happened much later. The first known study of multisets is attributed to the Indian mathematician Bhāskarāchārya circa 1150, who described permutations of multisets. The work of Marius Nizolius (1498–1576) contains another early reference to the concept of multisets. Athanasius Kircher found the number of multiset permutations when one element can be repeated. Jean Prestet published a general rule for multiset permutations in 1675. John Wallis explained this rule in more detail in 1685.

In the explicit form, multisets appeared in the work of Richard Dedekind. Other mathematicians formalized multisets and began to study them as a precise mathematical object (structure) in the 20th century.

Example multisets

One of the simplest and most natural examples is the multiset of prime factors of a number n. Here the underlying set of elements is the set of prime divisors of n. For example, the number 120 has the prime factorization

120 = 2 3 3 1 5 1

which gives the multiset {2, 2, 2, 3, 5}.

A related example is the multiset of solutions of an algebraic equation. A quadratic equation, for example, has two solutions. However, in some cases they are both the same number. Thus the multiset of solutions of the equation could be {3, 5}, or it could be {4, 4}. In the latter case it has a solution of multiplicity 2.

A special case of the above are the eigenvalues of a matrix, if these are defined as the multiset of roots of its characteristic polynomial. However a choice is made here: the (usual) definition of eigenvalues does not refer to the characteristic polynomial, and other possibilities give rise to different notions of multiplicity, and therefore to different multisets. The geometric multiplicity of λ as eigenvalue of a matrix A is the dimension of the kernel of A − λI, which is often smaller than its multiplicity as root of the characteristic polynomial (the algebraic multiplicity) when the latter is larger than 1. The set of eigenvalues of A is also the set of roots of its minimal polynomial, but the multiset of those roots need not be the same either as the one defined using algebraic multiplicity, or using the geometric multiplicity.

Formal definition

Within set theory, a multiset may be formally defined as a 2-tuple (A, m) where A is some set and m : A → N ≥ 1 is a function from A to the set N ≥ 1 = { 1 , 2 , 3 , … } of positive integers. The set A is called the underlying set of elements. For each a in A the multiplicity (that is, number of occurrences) of a is the number m(a). If a universe U in which the elements of A must live is specified, the definition can be simplified to just a multiplicity function m U : U → N ≥ 0 from U to the set N ≥ 0 = { 0 , 1 , 2 , 3 , … } of nonnegative integers, obtained by extending m to U with values 0 for each element not in A. This extended multiplicity function is the multiplicity function called 1_A below. Like any function, the function m may be defined as its graph: the set of ordered pairs { ( a , m ( a ) ) : a ∈ A } . With these definitions the multiset written as {a, a, b} is defined as ({a, b}, {(a, 2), (b, 1)}), and the multiset {a, b} is defined as ({a, b}, {(a, 1), (b, 1)}).

The concept of a multiset is a generalization of the concept of a set. A multiset corresponds to an ordinary set if the multiplicity of every element is one (as opposed to some larger natural number). An indexed family, (a_i), where i is in some index-set, may define a multiset, sometimes written {a_i}, in which the multiplicity of any element x is the number of indices i such that a i = x . The condition for this to be possible is that no element occurs infinitely many times in the family: even in an infinite multiset, the multiplicities must be finite numbers.

It is possible to extend the definition of a multiset by allowing multiplicities of individual elements to be infinite cardinals instead of natural numbers. Not all properties carry over to this generalization however, and this article uses the definition above, with finite multiplicities.

Multiplicity function

The set indicator function of a normal set is generalized to the multiplicity function for multisets. The set indicator function of a subset A of a set X is the function

1 A : X → { 0 , 1 }

defined by

1 A ( x ) = { 1 if  x ∈ A , 0 if  x ∉ A .

The set indicator function of the intersection of sets is the minimum function of the indicator functions

1 A ∩ B ( x ) = min { 1 A ( x ) , 1 B ( x ) } .

The set indicator function of the union of sets is the maximum function of the indicator functions

1 A ∪ B ( x ) = max { 1 A ( x ) , 1 B ( x ) } .

The set indicator function of the set difference of sets is the zero-truncated subtraction of the indicator functions

1 A ∖ B ( x ) = max ( 0 , 1 A ( x ) − 1 B ( x ) ) .

The set indicator function of a subset is smaller than or equal to that of the superset

A ⊆ B ⇔ ∀ x , 1 A ( x ) ≤ 1 B ( x ) .

The set indicator function of a cartesian product is the product of the indicator functions of the cartesian factors

1 A × B ( x , y ) = 1 A ( x ) ⋅ 1 B ( y ) .

The cardinality of a (finite) set is the sum of the indicator function values

| A | = ∑ x ∈ X 1 A ( x ) .

Now generalize the concept of set indicator function by releasing the constraint that the values are 0 and 1 only and allow the values 2, 3, 4 and so on. The resulting function is called a multiplicity function and such a function defines a multiset. The concepts of intersection, union, difference, subset, cartesian product and cardinality of multisets are defined by the above formulas. There are also the following additional operations.

The multiplicity function of a multiset sum, is the sum of the multiplicity functions

1 A ⊎ B ( x ) = 1 A ( x ) + 1 B ( x ) .

The scalar multiplication of a multiset by a natural number n may be defined as:

1 n ⊗ A ( x ) = n × 1 A ( x ) .

Example multiset operations

These examples demonstrate these operations of multiplicity functions.

{ 1 , 1 , 1 , 3 } ∩ { 1 , 1 , 2 } = { 1 , 1 } { 1 , 1 } ∪ { 1 , 2 } = { 1 , 1 , 2 } { 1 , 1 } ∖ { 1 , 2 } = { 1 } { 1 , 1 } ⊆ { 1 , 1 , 1 , 2 } { 1 , 1 } × { 1 , 2 } = { ( 1 , 1 ) , ( 1 , 1 ) , ( 1 , 2 ) , ( 1 , 2 ) } | { 1 , 1 } | = 2 { 1 , 1 } ⊎ { 1 , 2 } = { 1 , 1 , 1 , 2 } 3 ⊗ { 1 , 1 , 2 } = { 1 , 1 , 1 , 1 , 1 , 1 , 2 , 2 , 2 }

Free commutative monoids

The free commutative monoid on a set X (see free object) can be taken to be the set of finite multisets with elements drawn from X, with the monoid operation being multiset sum and the empty multiset as identity element. Such monoids are also known as (finite) formal sums of elements of X with natural coefficients. The free commutative semigroup is the subset of the free commutative monoid which contains all multisets with elements drawn from X except the empty multiset.

Free abelian groups are formal sums (i.e. linear combinations) of elements of X with integer coefficients. Equivalently, they may be seen as signed finite multisets with elements drawn from X.

Counting multisets

The number of multisets of cardinality k, with elements taken from a finite set of cardinality n, is called the multiset coefficient or multiset number. This number is written by some authors as ( ( n k ) ) , a notation that is meant to resemble that of binomial coefficients; it is used for instance in (Stanley, 1997), and could be pronounced "n multichoose k" to resemble "n choose k" for ( n k ) . Unlike for binomial coefficients, there is no "multiset theorem" in which multiset coefficients would occur, and they should not be confused with the unrelated multinomial coefficients that occur in the multinomial theorem.

The value of multiset coefficients can be given explicitly as

( ( n k ) ) = ( n + k − 1 k ) = ( n + k − 1 ) ! k ! ( n − 1 ) ! = n ( n + 1 ) ( n + 2 ) ⋯ ( n + k − 1 ) k ! ,

where the second expression is as a binomial coefficient; many authors in fact avoid separate notation and just write binomial coefficients. So, the number of such multisets is the same as the number of subsets of cardinality k in a set of cardinality n + k − 1. The analogy with binomial coefficients can be stressed by writing the numerator in the above expression as a rising factorial power

( ( n k ) ) = n k ¯ k ! ,

to match the expression of binomial coefficients using a falling factorial power:

( n k ) = n k _ k ! .

There are for example 4 multisets of cardinality 3 with elements taken from the set {1, 2} of cardinality 2 (n = 2, k = 3), namely {1, 1, 1}, {1, 1, 2}, {1, 2, 2}, {2, 2, 2}. There are also 4 subsets of cardinality 3 in the set {1, 2, 3, 4} of cardinality 4 (n + k − 1), namely {1, 2, 3}, {1, 2, 4}, {1, 3, 4}, {2, 3, 4}.

One simple way to prove the equality of multiset coefficients and binomial coefficients given above, involves representing multisets in the following way. First, consider the notation for multisets that would represent {a, a, a, a, a, a, b, b, c, c, c, d, d, d, d, d, d, d} (6 as, 2 bs, 3 cs, 7 ds) in this form:

∙ ∙ ∙ ∙ ∙ ∙ ∣ ∙ ∙ ∣ ∙ ∙ ∙ ∣ ∙ ∙ ∙ ∙ ∙ ∙ ∙

This is a multiset of cardinality k = 18 made of elements of a set of cardinality n = 4. The number of characters including both dots and vertical lines used in this notation is 18 + 4 − 1. The number of vertical lines is 4 − 1. The number of multisets of cardinality 18 is then the number of ways to arrange the 4 − 1 vertical lines among the 18 + 4 − 1 characters, and is thus the number of subsets of cardinality 4 − 1 in a set of cardinality 18 + 4 − 1. Equivalently, it is the number of ways to arrange the 18 dots among the 18 + 4 − 1 characters, which is the number of subsets of cardinality 18 of a set of cardinality 18 + 4 − 1. This is

( 4 + 18 − 1 4 − 1 ) = ( 4 + 18 − 1 18 ) = 1330 ,

thus is the value of the multiset coefficient and its equivalencies:

( ( 4 18 ) ) = ( 21 18 ) = 21 ! 18 ! 3 ! = ( 21 3 ) = ( ( 19 3 ) ) , = 4 ⋅ 5 ⋅ 6 ⋅ 7 ⋅ 8 ⋅ 9 ⋅ 10 ⋅ 11 ⋅ 12 ⋅ 13 ⋅ 14 ⋅ 15 ⋅ 16 ⋅ 17 ⋅ 18 ⋅ 19 ⋅ 20 ⋅ 21 1 ⋅ 2 ⋅ 3 ⋅ 4 ⋅ 5 ⋅ 6 ⋅ 7 ⋅ 8 ⋅ 9 ⋅ 10 ⋅ 11 ⋅ 12 ⋅ 13 ⋅ 14 ⋅ 15 ⋅ 16 ⋅ 17 ⋅ 18 , = 1 ⋅ 2 ⋅ 3 ⋅ 4 ⋅ 5 ⋯ 16 ⋅ 17 ⋅ 18 ⋅ 19 ⋅ 20 ⋅ 21 1 ⋅ 2 ⋅ 3 ⋅ 4 ⋅ 5 ⋯ 16 ⋅ 17 ⋅ 18 × 1 ⋅ 2 ⋅ 3 , = 19 ⋅ 20 ⋅ 21 1 ⋅ 2 ⋅ 3 .

One may define a generalized binomial coefficient

( n k ) = n ( n − 1 ) ( n − 2 ) ⋯ ( n − k + 1 ) k !

in which n is not required to be a nonnegative integer, but may be negative or a non-integer, or a non-real complex number. (If k = 0, then the value of this coefficient is 1 because it is the empty product.) Then the number of multisets of cardinality k in a set of cardinality n is

( ( n k ) ) = ( − 1 ) k ( − n k ) .

Recurrence relation

A recurrence relation for multiset coefficients may be given as

( ( n k ) ) = ( ( n k − 1 ) ) + ( ( n − 1 k ) ) for  n , k > 0

with

( ( n 0 ) ) = 1 , n ∈ N , and ( ( 0 k ) ) = 0 , k > 0.

The above recurrence may be interpreted as follows. Let [n] :=  { 1 , ⋯ , n } be the source set. There is always exactly one (empty) multiset of size 0, and if n = 0 there are no larger multisets, which gives the initial conditions.

Now, consider the case in which n,k > 0. A multiset of cardinality k with elements from [n] might or might not contain any instance of the final element n. If it does appear, then by removing n once, one is left with a multiset of cardinality k − 1 of elements from [n], and every such multiset can arise, which gives a total of

( ( n k − 1 ) ) possibilities.

If n does not appear, then our original multiset is equal to a multiset of cardinality k with elements from [n − 1], of which there are

( ( n − 1 k ) ) .

Thus,

( ( n k ) ) = ( ( n k − 1 ) ) + ( ( n − 1 k ) ) .

Generalization and connection to the negative binomial series

The multiplicative formula allows the definition of multiset coefficients to be extended by replacing n by an arbitrary number α (negative, real, complex):

( ( α k ) ) = α k ¯ k ! = α ( α + 1 ) ( α + 2 ) ⋯ ( α + k − 1 ) k ( k − 1 ) ( k − 2 ) ⋯ 1 for  k ∈ N  and arbitrary  α .

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

( 1 − X ) − α = ∑ k = 0 ∞ ( ( α k ) ) X k .

This Taylor series 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 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 ) − ( − α β ) ,

and formulas such as these can be used to prove identities for the multiset coefficients.

If α is a nonpositive integer n, then all terms with k > −n are zero, and the infinite series becomes a finite sum. However, for other values of α, including positive integers and rational numbers, the series is infinite.

Polynomial notation

The set {x} may be represented by the monomial x. Its power set, {{}, {x}}, is represented by the binomial 1 + x.

The set {x, y} may be represented by the monomial xy. Its power set, {{}, {x}, {y}, {x, y}}, is represented by the polynomial

( 1 + x ) ( 1 + y ) = 1 + x + y + x y .

The multiset {x, x} may be represented by the monomial x². The multiset of submultisets of {x, x}, {{}, {x}, {x}, {x, x}}, is represented by the polynomial

( 1 + x ) 2 = 1 + 2 x + x 2 .

The multiset of submultisets of the multiset represented by the polynomial xⁿ is

( 1 + x ) n = ∑ k = 0 n ( n k ) x k .

That is why the binomial coefficient counts the number of k-combinations of an n-set.

The multiset represented by the polynomial x^Ky^{N − K}, containing N elements, K of which are hits, is called a statistical population, and a submultiset is called a statistical sample. The set of samples is

( 1 + x ) K ( 1 + y ) N − K

which by the binomial theorem equals

∑ n = 0 N ∑ k = 0 K ( K k ) ( N − K n − k ) x k y n − k .

So the number of n-samples with k hits is

( K k ) ( N − K n − k ) .

See hypergeometric distribution and inferential statistics for further on the distribution of hits.

The infinite set of finite multisets of elements taken from the set {x}:

{ { } , { x } , { x , x } , { x , x , x } , ⋯ }

may be represented by the formal power series

S = 1 + x + x 2 + x 3 + ⋯ = 1 + x S .

The formal solution,

S = ( 1 − x ) − 1 ,

makes sense as a set of multisets, but the intermediate result, 1 − x, does not make sense as a set of multisets.

The infinite set of finite multisets of elements taken from the set represented by the monomial xy is

( 1 − x ) − 1 ( 1 − y ) − 1 = 1 + ( x + y ) + ( x 2 + x y + y 2 ) + ⋯ .

The infinite multiset of finite multisets of elements taken from the multiset represented by the monomial x² may be taken as a special case where x = y:

( 1 − x ) − 2 = 1 + 2 x + 3 x 2 + ⋯ .

The infinite multiset of finite multisets of elements taken from the multiset represented by the monomial xⁿ is

( 1 − x ) − n = ∑ k = 0 ∞ ( − n k ) ( − x ) k .

This explains why "multisets are negative sets." The negative binomial coefficients count the number of k-multisets of elements taken from an n-set.

Cumulant generating function

A non-negative integer n, can be represented by the monomial xⁿ. A finite multiset of non-negative integers can likewise be represented by a polynomial f(x).

It is convenient to consider the cumulant generating function g(t) = log f(e^t).

The cardinality of the multiset is e^g(0) = f(1).

The derivative of the cumulant generating function is g ′ ( t ) = f ( e t ) − 1 f ′ ( e t ) e t .

The mean value of the multiset is μ = g ′ ( 0 ) = f ( 1 ) − 1 f ′ ( 1 ) .

The variance of the multiset is σ 2 = g ″ ( 0 ) .

For example, the multiset {2, 2, 2, 3, 5} can be represented by the polynomial f(x) = 3x² + x³ + x⁵, with cumulant generating function g(t) = log(3e^2t + e^3t + e^5t), cardinality 3 + 1 + 1 = 5, derivative g ′ ( t ) = ( 3 e 2 t + e 3 t + e 5 t ) − 1 ( 6 e 2 t + 3 e 3 t + 5 e 5 t ) , and mean value μ = (3 + 1 + 1)⁻¹(6 + 3 + 5) = 2.8.

The numbers ( μ , σ 2 , ⋯ ) = ( g ′ ( 0 ) , g ″ ( 0 ) , ⋯ ) are called cumulants.

The infinite set of non-negative integers {0, 1, 2, ⋯} is represented by the formal power series 1 + x + x² + ⋯ = (1 − x)⁻¹. The mean value and standard deviation are undefined. Nevertheless, it has a cumulant-generating function, g(t) = −log(1 − e^t). The derivative of this cumulant-generating function is g ′ ( t ) = ( e − t − 1 ) − 1 .

A finite multiset of real numbers, A = {A_i}, is represented by the cumulant generating function

g A ( t ) = log ⁡ ( ∑ i e t A i ) .

This representation is unique: different multisets have different cumulant generating functions. See partition function for the case where the numbers in question are the energy levels of a physical system.

The cumulant-generating function of a multiset of n real numbers having mean μ and standard deviation σ is g(t) = log n + μt + 2⁻¹(σt)² + ⋯ and the derivative is g ′ ( t ) = μ + σ 2 t + ⋯ .

The cumulant-generating function of a set {k} consisting of a single real number k is g(t) = kt. The derivative is the number itself: g ′ ( t ) = k . So the concept of the derivative of the cumulant generating function of a multiset of real numbers is a generalization of the concept of a real number.

The cumulant-generating function of a constant multiset, {k, k, k, k, ⋯, k} of n elements all equal to the same real number k, is g(t) = log n + kt, and the derivative is the number itself: g ′ ( t ) = k , irrespective of n.

The cumulant-generating function of the multiset of sums of elements of two multisets of numbers is the sum of the two cumulant-generating functions:

g A + B ( t ) = log ⁡ ( ∑ i ∑ j e t ( A i + B j ) ) = log ⁡ ( ∑ i ∑ j e t A i e t B j ) = log ⁡ ( ∑ i e t A i ⋅ ∑ j e t B j ) = log ⁡ ( ∑ i e t A i ) + log ⁡ ( ∑ j e t B j ) = g A ( t ) + g B ( t ) .

There is no general formula for computing the cumulant generating function of a product

g A ⋅ B ( t ) = log ⁡ ( ∑ i ∑ j e t A i B j ) ,

but the special case of a constant times a multiset of numbers is:

g k ⋅ A ( t ) = log ⁡ ( ∑ i e t ( k A i ) ) = log ⁡ ( ∑ i e ( t k ) A i ) = g A ( k t ) .

The multiset 2⋅A = {2A_i} is not the same multiset as 2 × A = A + A = {A_i + A_j}. For example, 2 ⋅ { 1 , − 1 } = { 2 , − 2 } while 2 × {1, −1} = {1, −1} + {1, −1} = {1 + 1, 1 − 1, −1 + 1, −1 −1} = {2, 0, 0, −2}. The cumulant generating function of k × A is

g k × A ( t ) = k g A ( t ) .

The standard normal distribution is like a limit of large multisets of numbers.

lim k → ∞ k − 1 ( k 2 × { 1 , − 1 } ) .

This limit does not make sense as a multiset of numbers, but the derivative of the cumulant generating functions of the multisets in question makes sense, and the limit is well defined.

lim k → ∞ g k − 1 ( k 2 × { 1 , − 1 } ) ′ ( t ) = lim k → ∞ d d t ( k 2 log ⁡ ( e + t k − 1 + e − t k − 1 ) ) = lim k → ∞ d d t ( k 2 log ⁡ 2 + 2 − 1 t 2 + ⋯ ) = t .

The constant term k²log(2) vanishes by differentiation. The terms ⋯ vanish in the limit. So for the standard normal distribution, having mean 0 and standard deviation 1, the derivative of the cumulant generating function is g ′ ( t ) = t . For the normal distribution having mean μ and standard deviation σ, the derivative of the cumulant generating function is g ′ ( t ) = μ + σ 2 t .

Applications

Multisets have various applications. They are becoming the main structure of combinatorics because in its search for higher rigorousness, modern combinatorics has been developed not for sets but for multisets. Multisets have become an important tool in databases. For instance, multisets are often used to implement relations in database systems. Multisets also play an important role in computer science.

There are also other applications. For instance, Richard Rado used multisets as a device to investigate the properties of families of sets. He wrote, "The notion of a set takes no account of multiple occurrence of any one of its members, and yet it is just this kind of information which is frequently of importance. We need only think of the set of roots of a polynomial f(x) or the spectrum of a linear operator."

Generalizations

Different generalizations of multisets have been introduced, studied and applied to solving problems.

Fuzzy multisets

Rough multisets

Real-valued multisets (in which multiplicity of an element can be any real number)

Hybrid sets

where multiplicity is any real-valued step function

Soft multisets

Soft fuzzy multisets

Named set (unification of all generalizations of sets)