Pólya Enumeration

Prerequisites: Permutations and arrangements

Introduction

The Pólya enumeration theorem is usually used to solve counting problems involving “essentially different” objects.

This article may involve group theory

This article may involve group theory. It gives brief explanations of the group-theoretic concepts it uses, so that readers unfamiliar with them can understand and apply the Pólya enumeration theorem. For rigorous statements and discussions of group theory, see the sections on basic concepts of abstract algebra, group theory, and others.

“Symmetry groups”, “symmetric groups”, and “permutation groups”

This article inevitably uses the names of all three kinds of groups. Although they may be easy to confuse, they refer to different concepts. Given a geometric structure, its symmetry operations are the geometric transformations that make it coincide with itself, and its symmetry group is the set of those operations. The symmetric group is the set of all permutations of a given set. A permutation group is a subgroup of a symmetric group, consisting of some, not necessarily all, permutations. Later we explain how to represent the symmetry group of a given geometric structure as a permutation group and use it in counting problems.

Burnside's lemma

Related reading: Burnside's lemma

The Pólya enumeration theorem is an application and generalization of Burnside's lemma. Before introducing it, we briefly review Burnside's lemma.

To extract the general pattern, first consider a simple example.

Necklace coloring

We have a necklace of four beads. Each bead can be red or blue. How many essentially different bead colorings are there? (Two colorings are considered the same if one can be rotated to coincide with the other.)

Solution and analysis

This problem is simple enough to solve by enumeration. There are 44 beads, each with 22 possible colors, so there are 24=162^4=16 possible necklace colorings. Group together colorings obtainable from one another by rotation; there are 66 groups, as shown below. (The code for an individual coloring records the colors clockwise starting from the lower-left bead, with BB for blue and RR for red; colorings in the same group have the same background color.)

Necklace coloring
Necklace coloring

This example shows that to count essentially different colorings, the key is to know how many distinct colorings correspond to each equivalence class. In other words, we need to determine the size of every group in the figure.

Colorings in the same group are precisely those that can be transformed into one another by rotations. There are 44 rotations, namely

G={r0,r1,r2,r3},G=\{r_0,r_1,r_2,r_3\},

representing rotations 0,1,2,30,1,2,3 respectively. A rotation by 00 leaves the necklace unchanged.

First consider the group containing RRBBRRBB. Applying these four operations gives

RRBB,RBBR,BBRR,BRRB.RRBB, RBBR, BBRR, BRRB.

The four colorings are distinct, so this group has 44 elements.

Now consider the group containing BRBRBRBR. Applying the same four operations gives

BRBR,RBRB,BRBR,RBRB.BRBR, RBRB, BRBR, RBRB.

The result of rotating twice is the same as not rotating, and the result of rotating three times is the same as rotating once. Therefore this group has 22 elements.

Finally, consider the groups containing BBBBBBBB and RRRRRRRR. All four operations give each of them back unchanged, so each group has only 11 element.

Let xx denote a coloring, and let GxGx denote the set of color codes obtainable by applying operations to xx. The example suggests that the action of GG on xx has a kind of “periodicity”.

Let ∣G∣|G| be the total number of operations. This “periodicity” means that if mm distinct operations in GG transform xx into itself, then the results for xx under these operations repeat mm times. Thus, for xx, there are ∣G∣/m|G|/m distinct results, which is the size of the group containing xx.

In this example, only the zero rotation r0r_0 fixes RRBBRRBB, so its group has size 4/1=44/1=4; both r0r_0 and the two-step rotation r2r_2 fix BRBRBRBR, so its group has size 4/2=24/2=2; every rotation fixes BBBBBBBB, so its group has size 4/4=14/4=1.

Below, let GxG_x denote the number of operations that transform xx into itself, so ∣Gx∣|G_x| is the mm above. The size of the group containing XX is ∣G∣/∣Gx∣|G|/|G_x|. To count the groups of colorings, enumerate all possible colorings x∈Xx\in X; for a coloring with group size ∣Gx∣|Gx|, denoted xx, assign the weight 1/∣Gx∣1/|Gx|. Then the number of groups is

∣X/G∣=∑x∈X1∣Gx∣=∑x∈X∣Gx∣∣G∣.|X/G|=\sum_{x\in X}\frac{1}{|Gx|}=\sum_{x\in X}\frac{|G_x|}{|G|}.

The form of this expression is not convenient to use. Let gxgx be the result of applying an operation to a coloring x∈Xx\in X using g∈Gg\in G. Then the set GxG_x described above is {g∈G:gx=x}\{g\in G:gx=x\}, so exchanging the order of summation gives

∑x∈X∣Gx∣=∑x∈X∣{g∈G:gx=x}∣=∑x∈X∑g∈G[gx=x]=∑g∈G∑x∈X[gx=x]=∑g∈G∣{x∈X:gx=x}∣=∑g∈G∣Xg∣.\begin{aligned} \sum_{x\in X}|G_x| &=\sum_{x\in X}|\{g\in G:gx=x\}|\\ &=\sum_{x\in X}\sum_{g\in G}[gx=x]\\ &=\sum_{g\in G}\sum_{x\in X}[gx=x]\\ &=\sum_{g\in G}|\{x\in X:gx=x\}|\\ &=\sum_{g\in G}|X^g|. \end{aligned}

Here [⋅][\cdot] is the Iverson bracket. In the result after exchanging the sums, Xg={x∈X:gx=x}X^g=\{x\in X:gx=x\} is the set of colorings that, under operation gg, remain unchanged, with such a coloring denoted xx; in short, it is the set of fixed points of gg.

We can now write the number of groups as

∣X/G∣=1∣G∣∑g∈G∣Xg∣.|X/G|=\frac{1}{|G|}\sum_{g\in G}|X^g|.

In other words, the number of groups is the average number of fixed points of the rotations.

As an application of this result, recompute the number of necklace colorings. The fixed points of the rotations are:

OperationFixed points
r0r_0XX
r1r_1{BBBB,RRRR}\{BBBB,RRRR\}
r2r_2{BBBB,BRBR,RBRB,RRRR}\{BBBB,BRBR,RBRB,RRRR\}
r3r_3{BBBB,RRRR}\{BBBB,RRRR\}

Therefore, the number of groups is

16+2+4+24=6.\frac{16+2+4+2}{4} = 6.

This recovers the result above.

The example suggests the following general result for solving such counting problems. For convenience, we focus on coloring problems here, though the result also applies to other situations; examples appear at the end.

A coloring problem starts with a structure and colors each of its vertices, producing different colorings. The structure has symmetries, so colorings that appear different can be transformed into one another by symmetry operations. Such colorings are called essentially the same. The problem is to determine the number of essentially different colorings.

Following the analysis of the example, first determine the symmetry operations of the given structure. The set of these operations, GG, is called the spatial symmetry group of the structure. In practice, it is usually unnecessary to know the definition of a group; it is enough to enumerate all spatial symmetry operations without omission or duplication. The structures of several common spatial symmetry groups are analyzed later, where the definition of a group is explained.

Let the set of all colorings be XX, and let an individual coloring be xx. The result of applying operation g∈Gg\in G to a coloring x∈Xx\in X is gxgx. Thus all results obtainable by applying an operation to xx form Gx={gx:g∈G}Gx=\{gx:g\in G\}, which is called the orbit of group GG acting on xx. Different colorings in the same orbit are what “essentially the same” means in these problems. Therefore, the number of essentially different colorings is equal to the number of distinct orbits.

The analysis of the example generalizes to the general case.

Burnside's lemma

Given an action of a group GG on a set XX, the number of distinct orbits is

∣X/G∣=1∣G∣∑g∈G∣Xg∣.|X/G|=\frac{1}{|G|}\sum_{g\in G}|X^g|.

Here, Xg={x∈X:gx=x}X^g=\{x\in X:gx=x\} is the set of fixed points under the action of g∈Gg\in G.

Its proof is almost exactly the analysis of the example above. The example used the observation that the results of the action of group GG on one element xx have a kind of “periodicity”, so the number of repetitions equals the number of operations that transform xx into itself. This fact is correct in general, but because the structure of group GG may be complicated, its “periodicity” may not be as direct as in the example. A rigorous statement of this fact requires the orbit-stabilizer theorem from group theory.

In applications, it is enough to enumerate all symmetry operations and give the number of fixed points corresponding to each operation. Here is a slightly more complicated application.

Cube coloring

Color a cube with three colors and find the number of essentially different colorings (two colorings that become the same after a spatial rotation are considered identical).

Solution

A cube has 66 faces, and each face has 33 coloring choices, so there are 363^6 colorings in total, that is, ∣X∣=36|X|=3^6. Let the spatial symmetry group of the cube be GG.

Figure: cube

We now analyze all operations in GG. For convenience, call the six faces front, back, top, bottom, left, and right:

  • Identity: Every coloring is fixed, so ∣Xg∣=36|X^g|=3^6;
  • A 90∘90^\circ rotation about the line through the centers of two opposite faces: there are 33 choices of opposite faces and 22 choices of rotation direction, so this type gives 66 permutations. If the axis passes through the centers of the front and back faces, the top, bottom, left, and right faces must have the same color for the coloring to remain fixed. There are 33 independently colored regions, so ∣Xg∣=33|X^g|=3^3;
  • A 180∘180^\circ rotation about the line through the centers of two opposite faces: there are 33 choices of opposite faces, and the rotation direction does not matter, so this type gives 33 permutations. If the axis passes through the centers of the front and back faces, the top and bottom faces must have the same color, and the left and right faces must have the same color. There are 44 independently colored regions, so ∣Xg∣=34|X^g|=3^4;
  • A 180∘180^\circ rotation about the line through the midpoints of two opposite edges: there are 66 choices of opposite edges, and the rotation direction again does not affect the permutation, so this type gives 66 permutations. If the opposite edges are the boundaries of the front/top pair and the bottom/back pair, the front and top faces must have the same color, the bottom and back faces must have the same color, and the left and right faces must have the same color. There are 33 independently colored regions, so ∣Xg∣=33|X^g|=3^3;
  • A 120∘120^\circ rotation about the line through two opposite vertices: there are 44 choices of opposite vertices and 22 choices of rotation direction, so this type gives 88 permutations. If the opposite vertices are the upper-right vertex of the front face and the lower-left vertex of the back face, the front, top, and right faces must have the same color, and the back, bottom, and left faces must have the same color. There are 22 independently colored regions, so ∣Xg∣=32|X^g|=3^2.

Therefore, the number of essentially different colorings is

1×36+6×33+3×34+6×33+8×321+6+3+6+8=57.\frac{1\times3^6+6\times3^3+3\times3^4+6\times3^3+8\times3^2}{1+6+3+6+8}=57.

The Pólya enumeration theorem

The statement of Burnside's lemma does not use the fact that XX is the set of all colorings of a structure. In fact, Burnside's lemma is not limited to coloring problems. For coloring problems, the Pólya enumeration theorem provides a more precise calculation. It can be viewed as the application of the general Burnside's lemma to coloring problems.

The improvement over Burnside's lemma is that the Pólya enumeration theorem provides a concrete method for calculating the size ∣Xg∣|X^g| of the fixed-point set in coloring problems.

This is intuitive from the cube-coloring example. For each symmetry operation of a cube, the size of its fixed-point set has the form mc(g)m^{c(g)}, where mm is the number of colors and c(g)c(g) is the number of regions that can be colored independently under operation gg. This observation also holds in general, but we need to make clear how to calculate, for a given gg, the value c(g)c(g).

Choosing a coloring of a structure means, mathematically, choosing a map from the set XX of colorable objects in the structure (such as beads in a necklace or faces of a cube) to the set CC of colors, namely f:X→Cf:X\rightarrow C. Thus the set of colorings is CXC^X. The spatial symmetry group GG of the structure acts on the structure and therefore naturally acts on the set XX as well. Each such symmetry operation is a bijection on XX, namely a permutation. 1

Now analyze the structure of the fixed-point set (CX)g(C^X)^g. Given gg as a permutation of XX, the example shows that if a position in XX, denoted xx, can, after finitely many applications of gg, be moved to a position yy, then a fixed point f∈(CX)gf\in (C^X)^g must satisfy f(x)=f(y)f(x)=f(y). In the language of orbits from the previous section, positions xx and yy lie in the same orbit under the action of gg2, so they must be colored the same. In the language of permutations, for permutation gg, the cycle decomposition shows that positions xx and yy lie in the same cycle, so they must be colored the same. Different cycles in the cycle decomposition need not have the same color and can be colored independently. Thus the number of independently colorable regions is c(g)c(g), namely the number of cycles in the cycle decomposition of gg.

It follows that the number of fixed points of operation gg is ∣C∣c(g)|C|^{c(g)}. Substituting this conclusion into Burnside's lemma gives the unweighted version of the Pólya enumeration theorem.

Pólya enumeration theorem (unweighted version)

Given an action of a group GG on a set XX and a color set CC, the number of distinct colorings is

∣CX/G∣=1∣G∣∑g∈Gmc(g),|C^X/G|=\frac{1}{|G|}\sum_{g\in G}m^{c(g)},

where mm is the number of colors and c(g)c(g) is the number of cycles in the cycle decomposition of the permutation representation of an element g∈Gg\in G.

Meaning of the group GG

There is a slight abuse of notation here. If group GG acts on XX, the action on the coloring set CXC^X should be defined separately, but the distinction is omitted here.

As a simple application of the Pólya enumeration theorem, recompute the example from earlier.

Another solution to the necklace-coloring problem

Number the four beads 1∼41\sim 4. The elements of GG in the example have the following permutation representations, all written as cycle decompositions:

  • Zero rotations r0=(1)r_0=(1), with 44 cycles in total (the omitted 11-cycles should be noted);
  • One rotation r1=(1234)r_1=(1234), with 11 cycle;
  • Two rotations r2=(13)(24)r_2=(13)(24), with 22 cycles;
  • Three rotations r3=(1432)r_3=(1432), with 11 cycle.

Therefore, the number of essentially different colorings is

24+21+22+214=6.\frac{2^4+2^1+2^2+2^1}{4}=6.
Another solution to the cube-coloring problem

The preceding analysis effectively gave the cycle representations of each type of permutation, although it did not write them explicitly in numerical notation, so we do not repeat that analysis here. We illustrate only the case of a 180∘180^\circ rotation about the line joining the midpoints of opposite edges. Number the front, back, top, bottom, left, and right faces 1∼61\sim6 in that order. The corresponding permutation is (13)(24)(56)(13)(24)(56), so c(g)=3c(g)=3. Other types of permutation can be analyzed similarly, and the final counting expression is exactly the same as above.

Weighted generalization

The unweighted Pólya enumeration theorem can count all essentially different coloring problems, but it is insufficient for more detailed questions. For example, if the number of times each available color may be used is specified in the coloring problem above, the Pólya counting formula cannot be applied directly. In practice, this requires deriving the result again using Burnside's lemma; expressing the resulting conclusions as generating functions gives the weighted version of the Pólya enumeration theorem.

Necklace coloring (with restrictions)

We have a necklace of four beads, each of which can be red or blue. Exactly two red beads and two blue beads are available. How many essentially different bead colorings are there? (Two colorings are considered the same if one can be rotated to coincide with the other.)

Solution and analysis

Apply Burnside's lemma. With two red and two blue beads, there are (42)=6\dbinom{4}{2}=6 colorings. The spatial symmetry group G={r0,r1,r2,r3}G=\{r_0,r_1,r_2,r_3\} consists of rotations by 0∼30\sim3 steps. Their fixed-point sets are:

  • Zero rotations r0=(1)r_0=(1): all 66 colorings are fixed;
  • One rotation r1=(1234)r_1=(1234): a fixed point would require all beads to have the same color, so there are no fixed points;
  • Two rotations r2=(13)(24)r_2=(13)(24): there are two independently colored regions, both of size 22. They must be colored red and blue respectively, so the fixed-point set has size 22;
  • Three rotations r3=(1432)r_3=(1432): this is the same as one rotation, so there are no fixed points.

Therefore, by Burnside's lemma, the number of essentially different colorings is

6+0+2+04=2.\frac{6+0+2+0}{4}=2.

The example gives the following calculation method. For problems with restrictions on the numbers of different colors, we again color the cycles of each permutation in the spatial symmetry group, but the numbers of colors used must be exactly the specified numbers. Such combinatorial problems usually have no closed-form solution. Except for special cases that can be calculated by combinatorial methods, they need to be solved as a knapsack problem.

Generating functions can provide the answers to these counting problems. Given a permutation gg, suppose its cycle type is 1α12α2⋯nαn1^{\alpha_1}2^{\alpha_2}\cdots n^{\alpha_n}, meaning that it has αk\alpha_k cycles of length kk. If each cycle can be colored in one of mm colors, then the generating function

∏k=1n(∑i=1mxik)αk\prod_{k=1}^n\left(\sum_{i=1}^mx_i^k\right)^{\alpha_k}

The coefficient of the monomial x1β1x2β2⋯xmβmx_1^{\beta_1}x_2^{\beta_2}\cdots x_m^{\beta_m} is the count in which the ii-th color is used βi\beta_i times. The combinatorial meaning of the expression ∑i=1mxik\sum_{i=1}^mx_i^k in parentheses is as follows: for a cycle of length kk, the number of colorings that use kk instances of color ii is 11, while the count is 00 in all other cases. This exactly describes the requirement that all positions in the same cycle have the same color.

Given the generating function for colorings under a permutation gg, applying Burnside's lemma to each monomial gives the counts of essentially different colorings for every color combination. Since a generating function is linear in its monomials, the generating function for the number of essentially different colorings is

1∣G∣∑g∈G∏k=1n(∑i=1mxik)αk.\frac1{|G|}\sum_{g\in G}\prod_{k=1}^n\left(\sum_{i=1}^mx_i^k\right)^{\alpha_k}.

Expanding this expression, the coefficient of each monomial gives the number of essentially different colorings for the specified color combination.

In the discussion above, the generating function ∑i=1mxik\sum_{i=1}^mx_i^k for coloring each cycle has no special significance and can be replaced by another generating function. Thus we obtain the following general version of the Pólya enumeration theorem.

Cycle index of a permutation group

Given a permutation group GG, the cycle index of group GG is defined as

ZG(t1,t2,⋯ ,tn)=1∣G∣∑g∈Gt1c1(g)t2c2(g)⋯tncn(g),Z_G(t_1,t_2,\cdots,t_n)=\frac{1}{|G|}\sum_{g\in G}t_1^{c_1(g)}t_2^{c_2(g)}\cdots t_n^{c_n(g)},

where ck(g)c_k(g) is the number of cycles in the cycle decomposition of permutation gg that have length kk; that is, 1c1(g)2c2(g)⋯ncn(g)1^{c_1(g)}2^{c_2(g)}\cdots n^{c_n(g)} is the cycle type of permutation gg.

Pólya enumeration theorem (weighted version)

Given an action of a group GG on a set XX, suppose the generating function for the coloring methods at each point is f(x1,x2,⋯ ,xm)f(x_1,x_2,\cdots,x_m). Then the generating function for the essentially different colorings of XX is

ZG(f(x11,x21,⋯ ,xm1),f(x12,x22,⋯ ,xm2),⋯ ,f(x1n,x2n,⋯ ,xmn)),Z_G(f(x_1^1,x_2^1,\cdots,x_m^1),f(x_1^2,x_2^2,\cdots,x_m^2),\cdots,f(x_1^n,x_2^n,\cdots,x_m^n)),

where ZG(t1,t2,⋯ ,tn)Z_G(t_1,t_2,\cdots,t_n) is the cycle index of group GG.

Here, if the generating function for coloring a single position is f(x1,x2,⋯ ,xm)f(x_1,x_2,\cdots,x_m), then the generating function for coloring a cycle of length kk is f(x1k,x2k,⋯ ,xmk)f(x_1^k,x_2^k,\cdots,x_m^k). This reflects the fact that if a coloring is a fixed point of a given permutation, all positions in the same cycle must have the same color. Evaluating the generating function at xi=1x_i=1 gives the unweighted version of the Pólya enumeration theorem above.

The statement of the theorem uses the concept of the cycle index of a permutation group. It is independent of any particular coloring problem and describes the structure of a permutation group.

Another solution to the restricted necklace-coloring problem

The cycle index of the rotation symmetry group is 14(t14+t22+2t4)\dfrac14\left(t_1^4+t_2^2+2t_4\right), and the generating function for coloring a single point is r+br+b. Thus the generating function for all colorings is

F(r,b)=14((r+b)4+(r2+b2)2+2(r4+b4))=r4+r3b+2r2b2+rb3+b4.\begin{aligned} F(r,b)&=\frac14\left((r+b)^4+(r^2+b^2)^2+2(r^4+b^4)\right)\\ &=r^4+r^3b+2r^2b^2+rb^3+b^4. \end{aligned}

The desired count is the coefficient of r2b2r^2b^2, namely 22 essentially different colorings. Incidentally, this expression also gives the counts under other restrictions.

Applications

The weighted version of the Pólya enumeration theorem plays an important role in combinatorial counting problems. Here we briefly discuss its applications; for a more general discussion, see formal methods for combinatorial problems.

Diamond necklace

We have a necklace with four identical beads, and any number of diamonds can be set into each bead. If there are four diamonds, how many essentially different ways are there to set them? (Two settings are considered the same if one can be rotated to coincide with the other.)

Solution and analysis

The spatial symmetry group of the necklace is still the same as above. Without restricting the total number of diamonds, the generating function for the setting at one position is

f(x)=1+x+x2+⋯=∑i=1∞xi=11−x.f(x)=1+x+x^2+\cdots=\sum_{i=1}^\infty x^i=\frac{1}{1-x}.

Applying the weighted Pólya enumeration theorem, the generating function for all diamond settings is

F(x)=14(f(x)4+f(x2)2+2f(x4))=1+x+3x2+5x3+10x4+⋯ .\begin{aligned} F(x)&=\frac14\left(f(x)^4+f(x^2)^2+2f(x^4)\right)\\ &=1+x+3x^2+5x^3+10x^4+\cdots. \end{aligned}

Thus the desired number of settings is the coefficient of x4x^4, namely 1010. As a check, enumeration gives

4000,3100,3010,3001,2200,2020,2110,2101,2011,1111.4000,3100,3010,3001,2200,2020,2110,2101,2011,1111.

Each group of four digits gives the number of diamonds set into each bead.

This example shows that the weighted version of the Pólya enumeration theorem can solve a much broader class of problems than coloring counts. It provides a method for extending the count at a single point to the count of essentially different configurations on an entire structure. Coloring problems are only a special case of this class.

Common spatial symmetry groups

One difficulty in Pólya-counting problems is analyzing the structure of permutation groups. Here we briefly discuss common spatial symmetry groups and describe them using their cycle indices. Note that for the spatial symmetry group of the same structure, different sets of objects being acted on give different group actions, and therefore different permutation representations. For example, the action of the spatial symmetry group of a cube on its vertices, edges, and faces gives the vertex, edge, and face permutation groups of the cube, respectively. The numbers of vertices, edges, and faces differ, so these permutation groups and their corresponding cycle indices naturally differ as well. Thus, when solving a specific problem, the objects on which the group acts must not be overlooked.

Relationship between spatial symmetry groups and permutation groups

Although the two concepts are very similar, they are not the same object. In group-theoretic language, given a spatial symmetry group GG and its action on a set XX, the permutation representation of the action is a homomorphism from group GG to the symmetric group SXS_X, denoted φ\varphi. In combinatorial counting, this representation is often faithful, meaning ker⁡φ={e}\ker\varphi=\{e\}, so the homomorphism φ\varphi is an embedding of group GG into group SXS_X. The permutation group in this article is the image of this embedding, namely φ(G)\varphi(G), and it is isomorphic to the spatial symmetry group GG itself. Therefore, for the spatial symmetry group GG of the same structure, different choices of the group action can give different permutation groups φ(G)\varphi(G) and hence different cycle indices (isomorphic permutation groups need not have the same cycle index).

Given a structure, its spatial symmetry group is the set of all operations that can transform it into itself. It necessarily satisfies the following conditions:

  • Applying two symmetry operations successively to a given structure can be viewed as applying another symmetry operation; that is, the set of symmetry operations is closed under composition;
    • composition of symmetry operations is associative;
    • there is an identity symmetry operation; leaving the given structure unchanged is itself regarded as an operation;
    • every operation has an inverse operation that can cancel the effect of the given operation.

Groups are the abstraction of concepts satisfying all these conditions. The study of group structures is the main subject of group theory. The analysis here focuses on spatial symmetry groups and mainly uses geometric viewpoints. The common examples given below should provide readers with general approaches to analyzing such problems.

Cyclic groups

For a regular nn-gon, the spatial symmetry group formed by all rotations is called the cyclic group (cyclic group), denoted CnC_n. Let the operation of counterclockwise rotation by (360/n)∘(360/n)^\circ be rr. The elements of group CnC_n can be written as

Cn={e,r,r2,⋯ ,rn−1}.C_n=\{e,r,r^2,\cdots,r^{n-1}\}.

Here rkr^k means applying operation rr kk times, namely counterclockwise rotation by (360k/n)∘(360k/n)^\circ, while e=r0e=r^0 is the identity transformation.

Whether we consider the action of the cyclic group on all vertices or all edges of a regular nn-gon, its permutation representation is the same. Analyze the permutation representation of the action on the set of all vertices as an example. Its cycle index is

Z(Cn)=1n∑d∣nφ(d)tdn/d.Z(C_n)=\frac1n\sum_{d\mid n}\varphi(d)t_{d}^{n/d}.

Here φ(⋅)\varphi(\cdot) is the Euler totient function from number theory.

If we count only rotations, the spatial symmetry group of a necklace of length nn is CnC_n.

Analysis

Number the vertices in counterclockwise order as {0,1,⋯ ,n−1}\{0,1,\cdots,n-1\}. Then rk(i)=i+k mod nr^k(i)=i+k\bmod n. The set of vertices in the cycle containing vertex ii is

{i+ℓk mod n:ℓ∈Z}.\{i+\ell k\bmod n:\ell\in\mathbf Z\}.

Clearly, i≡i+ℓk(modn)i\equiv i+\ell k\pmod n if and only if

ngcd⁡(k,n)∣ℓ.\frac{n}{\gcd(k,n)}\mid\ell.

This means that every cycle containing a vertex ii has length ngcd⁡(k,n)\dfrac{n}{\gcd(k,n)}. Therefore permutation rkr^k has gcd⁡(k,n)\gcd(k,n) cycles of equal length. Combining like terms in the cycle-index expression, for a given d∣nd\mid n, the condition gcd⁡(k,n)=n/d\gcd(k,n)=n/d holds for kk values, and their number is φ(d)\varphi(d); all of them correspond to monomials of the form tdn/dt_d^{n/d}. Thus we obtain the cycle-index expression above.

Dihedral groups

For a regular nn-gon, the operations consisting of all rotations and reflections across symmetry axes also form a spatial symmetry group, called the dihedral group (dihedral group), denoted D2nD_{2n}. Let the operation of counterclockwise rotation by (360/n)∘(360/n)^\circ be rr, and let the operation of reflection across a given symmetry axis (for example, the line from the center to a vertex) be ss. Then the operations of group D2nD_{2n} can be written as

D2n={e,r,⋯ ,rn−1,s,sr,⋯ ,srn−1}.D_{2n}=\{e,r,\cdots,r^{n-1},s,sr,\cdots,sr^{n-1}\}.

Here rkr^k is still a rotation. Although srksr^k means first performing kk rotations and then reflecting across the given symmetry axis, it can equivalently be viewed as reflection across another symmetry axis. Therefore, group D2nD_{2n} has 11 identity transformation, (n−1)(n-1) rotations, and nn reflections in total. Its actions on the vertex set and edge set also have the same permutation representation. Its cycle index is

Z(D2n)=12Z(Cn)+{12t1t2k,n=2k+1,14(t12t2k−1+t2k),n=2k.Z(D_{2n})=\frac12Z(C_n)+ \begin{cases} \dfrac12t_1t_2^k,&n=2k+1,\\ \dfrac14\left(t_1^2t_2^{k-1}+t_2^k\right),&n=2k. \end{cases}
Analysis

The analysis of the rotations in group D2nD_{2n}, namely rkr^k including the identity transformation, is the same as for cyclic group CnC_n. The key is to analyze the remaining reflections, separating the cases according to the parity of the number of vertices nn.

When n=2k+1n=2k+1, every reflection axis connects a vertex to the midpoint of the opposite edge, and there are nn such axes. After each reflection, the vertex on the axis remains fixed and all other vertices are exchanged in pairs, giving 11 fixed point (a 11-cycle) and kk 22-cycles.

When n=2kn=2k, there are two types of symmetry axes. Half connect opposite vertices. Reflection across such an axis fixes the two vertices on the axis and exchanges all remaining vertices in pairs, giving 22 fixed points (two 11-cycles) and (k−1)(k-1) 22-cycles. The other half connect the midpoints of opposite edges. Reflection across such an axis exchanges all vertices in pairs, giving kk 22-cycles.

This analysis gives the cycle-index expression above.

Symmetric groups

The set of all permutations of nn elements forms a group called the symmetric group of degree nn (symmetric group), denoted SnS_n. It describes all symmetries that these nn vertices can have. It is also the permutation representation of the action of these symmetry operations on the vertex set.

According to the analysis in Permutations and arrangements, its cycle index is

Z(Sn)=∑a1+2α2+⋯+nαn=nt1α1t2α2⋯tnαn1α12α2⋯nαnα1!α2!⋯αn!.Z(S_n)=\sum_{a_1+2\alpha_2+\cdots+n\alpha_n=n}\frac{t_1^{\alpha_1}t_2^{\alpha_2}\cdots t_n^{\alpha_n}}{1^{\alpha_1}2^{\alpha_2}\cdots n^{\alpha_n}\alpha_1!\alpha_2!\cdots\alpha_n!}.

The count of permutations with cycle type 1α12α2⋯nαn1^{\alpha_1}2^{\alpha_2}\cdots n^{\alpha_n} used here is

n!1α12α2⋯nαnα1!α2!⋯αn!.\frac{n!}{1^{\alpha_1}2^{\alpha_2}\cdots n^{\alpha_n}\alpha_1!\alpha_2!\cdots\alpha_n!}.

It satisfies the recurrence relation

Z(Sn)=1n∑k=1ntkZ(Sn−k),Z(S_n)=\frac1n\sum_{k=1}^nt_kZ(S_{n-k}),

The starting point of the recurrence is Z(S0)=1Z(S_0)=1. The combinatorial meaning of this recurrence is that, to construct a permutation of length nn, we can first choose the cycle containing point nn and let its length be kk, and then construct a permutation on the remaining (n−k)(n-k) vertices.

Given a complete graph with nn vertices, its spatial symmetry group is precisely SnS_n. The cycle index of its action on the set of all vertices is given by Z(Sn)Z(S_n) above. However, its permutation representation on the set of all edges is different. For example, the sizes of the sets differ: the total number of edges is n(n−1)/2n(n-1)/2. The edge case requires additional analysis. Here we give a simple example; the general case can be found in the exercises.

Counting undirected simple graphs

Count the number of undirected simple graphs with 44 vertices up to isomorphism.

Solution

This is equivalent to coloring the complete graph with 44 vertices using two colors and counting essentially different colorings. The spatial symmetry group is S4S_4, and we now analyze the cycle index of its edge permutation group S4(2)S_4^{(2)}.

  • Identity (11 operation): the edges also remain fixed, so the corresponding monomial is t16t_1^6;
  • Swapping two vertices (66 operations): suppose aa and bb are swapped. Edges 11 and 33 remain fixed, while edges 22 and 55 are exchanged, and edges 44 and 66 are exchanged, so the corresponding monomial is 6t12t226t_1^2t_2^2;
  • A three-vertex cycle (88 operations): suppose the cycle is (abc)(abc). The edges 1,2,51,2,5 between these vertices cycle accordingly, and the edges connecting them to the fourth vertex dd, namely 4,6,34,6,3, also cycle accordingly, so the corresponding monomial is 8t328t_3^2;
  • Swapping two pairs of vertices (33 operations): suppose vertices aa and bb are swapped, and vertices cc and dd are swapped. Edges 11 and 33 remain fixed, while edges 22 and 44 are exchanged, and edges 55 and 66 are exchanged, so the corresponding monomial is 3t12t223t_1^2t_2^2;
  • A four-vertex cycle (66 operations): suppose the cycle is (abcd)(abcd). The edges 1,2,3,41,2,3,4 between adjacent vertices cycle accordingly, while the edges 5,65,6 between opposite vertices are exchanged, so the corresponding monomial is 6t2t46t_2t_4.

Therefore, the cycle index of the edge permutation group is

Z(S4(2))=124(t16+9t12t22+8t32+6t2t4).Z(S_4^{(2)})=\dfrac{1}{24}(t_1^6+9t_1^2t_2^2+8t_3^2+6t_2t_4).

By the Pólya enumeration theorem, the number of undirected simple graphs with 44 vertices up to isomorphism is

26+9×24+8×22+6×2224=11.\frac{2^6+9\times 2^4+8\times 2^2+6\times 2^2}{24} = 11.

Polyhedral groups

A polyhedral group is the spatial symmetry group of a regular polyhedron. There are only five regular polyhedra: the regular tetrahedron, cube, regular octahedron, regular dodecahedron, and regular icosahedron. If the adjacency relations among vertices, edges, and faces are preserved while vertices and faces are exchanged, we obtain the dual regular polyhedron. The regular tetrahedron is dual to itself, the cube is dual to the regular octahedron, and the regular dodecahedron is dual to the regular icosahedron. The duality relation simplifies the discussion of their spatial symmetry groups.

Counting only rotations that can be performed in three-dimensional space, there are only three spatial symmetry groups.

  • Tetrahedral group, that is, the spatial symmetry group of a regular tetrahedron:

  • Identity;

    • Rotations by 120∘120^\circ and 240∘240^\circ about the line through a vertex and the center of the opposite face;
    • A 180∘180^\circ rotation about the line through the midpoints of opposite edges.

There are 1+2×4+1×3=121+2\times4+1\times3=12 symmetry operations in total.

The cycle indices of its corresponding permutation groups are as follows.

  • Vertex and face permutation groups: 112(t14+8t1t3+3t22)\dfrac1{12}\left(t_1^4+8t_1t_3+3t_2^2\right);

    • Edge permutation group: 112(t16+8t32+3t12t22)\dfrac1{12}\left(t_1^6+8t_3^2+3t_1^2t_2^2\right).
  • Octahedral group, that is, the spatial symmetry group of the cube (and the regular octahedron):

  • Identity;

    • Rotations by 120∘120^\circ and 240∘240^\circ about the line through opposite vertices;
    • A 180∘180^\circ rotation about the line through the midpoints of opposite edges;
    • Rotations by 90∘90^\circ, 180∘180^\circ, and 270∘270^\circ about the line through the centers of opposite faces.

There are 1+2×4+1×6+3×3=241+2\times 4+1\times 6+3\times 3=24 symmetry operations in total.

The cycle indices of the permutation groups of the corresponding cube are as follows.

  • Vertex permutation group: 124(t18+8t12t32+9t24+6t42)\dfrac{1}{24}\left(t_1^8+8t_1^2t_3^2+9t_2^4+6t_4^2\right);
    • Edge permutation group: 124(t112+8t34+6t12t25+6t43+3t26)\dfrac{1}{24}\left(t_1^{12}+8t_3^4+6t_1^2t_2^5+6t_4^3+3t_2^6\right);
    • Face permutation group: 124(t16+8t32+6t23+6t12t4+3t12t22)\dfrac{1}{24}\left(t_1^6+8t_3^2+6t_2^3+6t_1^2t_4+3t_1^2t_2^2\right).

The permutation groups of the regular octahedron are analogous; only the roles of vertices and faces need to be exchanged.

  • Icosahedral group, that is, the spatial symmetry group of the regular dodecahedron (and the regular icosahedron):

  • Identity;

    • Rotations by 120∘120^\circ and 240∘240^\circ about the line through opposite vertices;
    • A 180∘180^\circ rotation about the line through the midpoints of opposite edges;
    • Rotations by 72∘72^\circ, 144∘144^\circ, 216∘216^\circ, and 288∘288^\circ about the line through the centers of opposite faces.

There are 1+2×10+1×15+6×4=601+2\times 10+1\times 15+6\times 4=60 symmetry operations in total.

The cycle indices of the permutation groups of the corresponding regular dodecahedron are as follows.

  • Vertex permutation group: 160(t120+20t12t36+15t210+24t54)\dfrac{1}{60}\left(t_1^{20}+20t_1^2t_3^6+15t_2^{10}+24t_5^4\right);
    • Edge permutation group: 160(t130+20t310+15t12t214+24t56)\dfrac{1}{60}\left(t_1^{30}+20t_3^{10}+15t_1^2t_2^{14}+24t_5^6\right);
    • Face permutation group: 160(t112+20t34+15t26+24t12t52)\dfrac{1}{60}\left(t_1^{12}+20t_3^4+15t_2^6+24t_1^2t_5^2\right).

The permutation groups of the regular icosahedron are analogous; only the roles of vertices and faces need to be exchanged.

The cycle indices given here are for permutation groups acting separately on vertices, edges, faces, and other individual types of objects. If different types of objects are to be colored simultaneously, a joint cycle index must be written.

Exercises

Coloring problems

These problems only require analyzing the structure of permutation groups and applying the Pólya enumeration theorem.

When the allowed color combinations are restricted, the number of ways to color the cycles needs to be found using knapsack DP or combinatorial methods.

Graph-counting problems

The Pólya enumeration theorem can be applied to graph-counting problems. The difficulty in these problems lies in enumerating the edge permutation group of a graph.

Another class of graph-counting problems to which the Pólya enumeration theorem applies requires directly manipulating generating functions.

References and notes

Footnotes

  1. Thus the spatial symmetry group GG can be represented as a permutation group on XX, that is, as a subgroup of the symmetric group SXS_X. ↩

  2. Strictly speaking, this is the action of the subgroup ⟨g⟩≤G\langle g\rangle\le G. ↩

This content comes from OI Wiki: Pólya Enumeration. See the original for authors and revision history. AIOJ has adapted the formatting, in-site links, and resource paths.

Content is shared under CC BY-SA 4.0 and SATA license; additional terms are governed by the original notice. Thanks to the OI Wiki authors and contributors; per the SATA terms, please visit the OI Wiki repository and star it to show your support.