Polya Enumeration Theorem

⚡ Key Takeaways

Pólya Enumeration Theorem (PET): counts distinct configurations under a group of symmetries by substituting a generating-function "figure inventory" into the cycle index of the group. One-line statement: the generating function for inequivalent colorings is ZG( ∑xi, ∑xi2, … )Z_G(\,\sum x_i,\ \sum x_i^2,\ \dots\,). Applies wherever you must count things up to symmetry (necklaces, graphs, chemical isomers, switching functions). Due to George Pólya (1937), with foundations by J. H. Redfield (1927).

💡 The Core Insight

Burnside's Lemma tells you the number of orbits is the average number of fixed colorings. Pólya's leap is to notice that a permutation gg fixes a coloring precisely when the coloring is constant on each cycle of gg. So a cycle of length kk forces kk positions to share one color — contributing the power sum x1k+x2k+⋯x_1^k + x_2^k + \cdots in a weighted count. Replacing each cycle-index variable sks_k by this power sum automatically averages over the group with full weight bookkeeping, yielding the complete generating function for orbits in one stroke.

Formal Statement

DefinitionGroup Action, Colorings, and Weight

Let a finite group GG act on a finite set DD of nn positions (the "domain", e.g. beads of a necklace). Let RR be a set of colors (the "range"). A coloring is a function f:D→Rf : D \to R. Two colorings f1,f2f_1, f_2 are equivalent if there exists g∈Gg \in G with f1(g⋅d)=f2(d)f_1(g \cdot d) = f_2(d) for all d∈Dd \in D; equivalence classes are orbits under the induced action on RDR^D. Assign each color r∈Rr \in R a weight w(r)w(r) (typically a monomial in formal variables). The weight of a coloring is the product over positions:

1
W(f)=∏d∈Dw(f(d)).\quad W(f) = \prod_{d \in D} w\bigl(f(d)\bigr).

Derivation hint: weight is multiplicative across positions; equivalent colorings have equal weight since gg merely permutes the factors of the product.

DefinitionCycle Index of a Permutation Group

For g∈Gg \in G acting on DD with n=∣D∣n = |D|, let jk(g)j_k(g) be the number of cycles of length kk in the cycle decomposition of gg (so ∑kk jk(g)=n\sum_k k\, j_k(g) = n). The cycle index of GG is the polynomial in indeterminates s1,s2,…,sns_1, s_2, \dots, s_n:

2
ZG(s1,s2,…,sn)=1∣G∣∑g∈G∏k=1nsk jk(g).\quad Z_G(s_1, s_2, \dots, s_n) = \frac{1}{|G|} \sum_{g \in G} \prod_{k=1}^{n} s_k^{\,j_k(g)}.

Derivation hint: each gg contributes one monomial recording its cycle type; the prefactor 1/∣G∣1/|G| averages over the group, mirroring Burnside's averaging.

TheoremPólya Enumeration Theorem

Let GG act on DD (∣D∣=n|D| = n) and let RR have a weight function ww. Define the figure-counting series (the "inventory" of colors) by the power sums

3
pk=∑r∈Rw(r)k.\quad p_k = \sum_{r \in R} w(r)^k.

Derivation hint: pkp_k is the value of substituting wkw^k into the simple color inventory ∑rw(r)\sum_r w(r); raising to the kk-th power matches a cycle of length kk forcing kk positions equal. Then the generating function for orbits, weighted by WW, is obtained by substituting sk↦pks_k \mapsto p_k in the cycle index:

4
∑orbits OW(O)=ZG(p1,p2,…,pn)=ZG ⁣(∑rw(r), ∑rw(r)2, … ).\quad \sum_{\text{orbits } O} W(O) = Z_G(p_1, p_2, \dots, p_n) = Z_G\!\left(\sum_r w(r),\ \sum_r w(r)^2,\ \dots\right).

Derivation hint: this is Burnside's Lemma applied to the weighted count, with the fixed-coloring count for gg factoring over its cycles into ∏kpkjk(g)\prod_k p_k^{j_k(g)}.

Unweighted corollary. Setting every w(r)=1w(r) = 1 gives pk=∣R∣=mp_k = |R| = m for all kk, so the number of inequivalent colorings is

5
N=ZG(m,m,…,m).\quad N = Z_G(m, m, \dots, m).

Derivation hint: with all weights 11, W(O)=1W(O)=1 so the left side of [4] simply counts orbits; each pkp_k collapses to mm.

Why It Holds / Derivation

ProofFrom Burnside's Lemma to Pólya

We prove [4]. Define the weighted fixed-point sum for each g∈Gg \in G:

6
Fixw(g)=∑f∈RDg⋅f=fW(f).\quad \mathrm{Fix}_w(g) = \sum_{\substack{f \in R^D \\ g \cdot f = f}} W(f).

Derivation hint: this generalizes the count of fixed colorings to a weighted sum; the unweighted version ∣Fix(g)∣|\mathrm{Fix}(g)| recovers ordinary Burnside.

Step 1 — Weighted Burnside. The weighted Cauchy–Frobenius (Burnside) identity states that summing weights over orbits equals the group-average of weighted fixed points:

7
∑orbits OW(O)=1∣G∣∑g∈GFixw(g).\quad \sum_{\text{orbits } O} W(O) = \frac{1}{|G|} \sum_{g \in G} \mathrm{Fix}_w(g).

Derivation hint: all colorings in an orbit share weight W(O)W(O); double-count pairs (g,f)(g, f) with g⋅f=fg\cdot f = f, using that each orbit's stabilizers have total size ∣G∣/∣O∣⋅∣O∣=∣G∣|G|/|O| \cdot |O| = |G| contributions summing per orbit.

Step 2 — A coloring is fixed iff it is constant on cycles. If g⋅f=fg \cdot f = f, then f(gt⋅d)=f(d)f(g^t \cdot d) = f(d) for all tt, so ff is constant along each orbit of ⟨g⟩\langle g \rangle on DD — that is, on each cycle of gg. Conversely any cycle-constant coloring is fixed.

Step 3 — Factor the weighted fixed sum over cycles. A fixed coloring chooses one color per cycle. A cycle of length kk, all kk positions sharing color rr, contributes weight w(r)kw(r)^k. Summing over the color choice gives factor ∑rw(r)k=pk\sum_r w(r)^k = p_k. Since choices on distinct cycles are independent:

8
Fixw(g)=∏k=1n(∑r∈Rw(r)k)jk(g)=∏k=1npk jk(g).\quad \mathrm{Fix}_w(g) = \prod_{k=1}^{n} \left(\sum_{r \in R} w(r)^k\right)^{j_k(g)} = \prod_{k=1}^{n} p_k^{\,j_k(g)}.

Derivation hint: independence of cycles makes the weighted fixed sum a product; each length-kk cycle contributes the kk-th power sum pkp_k, one factor per such cycle.

Step 4 — Assemble. Substitute [8] into [7]:

9
∑orbits OW(O)=1∣G∣∑g∈G∏k=1npk jk(g)=ZG(p1,…,pn).\quad \sum_{\text{orbits } O} W(O) = \frac{1}{|G|}\sum_{g \in G} \prod_{k=1}^{n} p_k^{\,j_k(g)} = Z_G(p_1, \dots, p_n).

Derivation hint: the right side is exactly the cycle index [2] with sks_k replaced by pkp_k — the substitution in [4]. ■\blacksquare

∎
mathematician · George Pólya | 1887–1985

Hungarian-American mathematician of extraordinary breadth, contributing to probability, combinatorics, number theory, complex analysis, and mathematical pedagogy (How to Solve It, 1945). His 1937 paper "Kombinatorische Anzahlbestimmungen für Gruppen, Graphen und chemische Verbindungen" introduced the cycle-index substitution method, motivated explicitly by counting chemical isomers and graphs. Pólya's work unified scattered enumeration techniques into a single mechanical procedure. Decades later it emerged that J. Howard Redfield had published essentially the same theorem in 1927 in an obscure paper; the result is now often credited as the Redfield–Pólya Theorem. Pólya's name endures largely because of the clarity and reach of his exposition and applications.

Building Intuition

🧠 Intuition

Picture a necklace of nn beads you can rotate. A rotation gg permutes the bead positions in cycles. If gg is a rotation by one step, all nn beads form a single nn-cycle — and the only colorings unchanged by it are those where every bead is the same color. That's why a length-nn cycle contributes pn=∑rw(r)np_n = \sum_r w(r)^n: there are mm monochromatic choices, each weighted by w(r)nw(r)^n. The cycle index is a "symmetry fingerprint" of the group: it forgets which positions cycle together and remembers only how many cycles of each length each symmetry produces. Pólya's substitution then injects the color palette (pkp_k) into this fingerprint. The genius is that the same ZGZ_G answers every coloring problem — change the palette, re-substitute, done. The unweighted formula [5] is just the special case "all colors weigh one," collapsing the inventory to mm everywhere.

Figure
Rotational symmetry of a 6-bead necklace and its cycle decomposition under rotation by 2 steps

```svg

Necklace C6, rotation by 2 positions 0 1 2 3 4 5 Two 3-cycles: (0 2 4)(1 3 5) → cycle index term s₃²

```

Worked Examples

ExampleTwo-color necklaces with 4 beads under rotation

Given. n=4n = 4 beads in a ring; symmetry group G=C4={e,r,r2,r3}G = C_4 = \{e, r, r^2, r^3\} (rotations only). Colors R={black,white}R = \{\text{black}, \text{white}\} with weights w(black)=bw(\text{black}) = b, w(white)=ww(\text{white}) = w.

Step 1 — Cycle types.

  • ee: four fixed points → s14s_1^4
  • rr (rotate 1): one 4-cycle → s4s_4
  • r2r^2 (rotate 2): two 2-cycles → s22s_2^2
  • r3r^3 (rotate 3): one 4-cycle → s4s_4

Step 2 — Cycle index.

10
ZC4=14(s14+s22+2s4).\quad Z_{C_4} = \tfrac{1}{4}\bigl(s_1^4 + s_2^2 + 2 s_4\bigr).

Derivation hint: average the four monomials from each group element by [2].

Step 3 — Substitute sk↦bk+wks_k \mapsto b^k + w^k:

11
ZC4(b+w, b2+w2, b4+w4)=14 ⁣[(b+w)4+(b2+w2)2+2(b4+w4)].\quad Z_{C_4}(b{+}w,\ b^2{+}w^2,\ b^4{+}w^4) = \tfrac14\!\left[(b{+}w)^4 + (b^2{+}w^2)^2 + 2(b^4{+}w^4)\right].

Derivation hint: replace each sks_k in [10] with the power sum pk=bk+wkp_k = b^k + w^k per [3].

Step 4 — Expand. Computing term by term and dividing by 4:

12
b4+b3w+2 b2w2+b w3+w4.\quad b^4 + b^3 w + 2\,b^2 w^2 + b\,w^3 + w^4.

Derivation hint: (b+w)4=b4+4b3w+6b2w2+4bw3+w4(b{+}w)^4 = b^4{+}4b^3w{+}6b^2w^2{+}4bw^3{+}w^4; add (b2+w2)2=b4+2b2w2+w4(b^2{+}w^2)^2 = b^4{+}2b^2w^2{+}w^4 and 2(b4+w4)2(b^4{+}w^4), then divide by 4.

Reading the answer. Coefficients give counts by composition: 1 all-black, 1 with three black, 2 with two of each, 1 with one black, 1 all-white.

Total necklaces=1+1+2+1+1=6 (or ZC4(2,2,2)=14(16+4+4)=6).\boxed{\text{Total necklaces} = 1+1+2+1+1 = 6 \text{ (or } Z_{C_4}(2,2,2)=\tfrac14(16+4+4)=6).}
ExampleColoring the faces of a cube with 3 colors

Given. Color the 6 faces of a cube with R={1,2,3}R = \{1,2,3\} colors, equivalent under the rotation group of the cube, ∣G∣=24|G| = 24. Use unweighted count (w≡1w \equiv 1, so pk=3p_k = 3).

Step 1 — Cycle index of the cube's rotation group acting on faces (standard result):

13
ZG=124(s16+6 s12s4+3 s12s22+8 s32+6 s23).\quad Z_G = \tfrac{1}{24}\bigl(s_1^6 + 6\,s_1^2 s_4 + 3\,s_1^2 s_2^2 + 8\,s_3^2 + 6\,s_2^3\bigr).

Derivation hint: classify the 24 rotations — identity (1), face rotations ±90° (6), face rotations 180° (3), vertex rotations ±120° (8), edge rotations 180° (6) — and record each one's cycle structure on the 6 faces.

Step 2 — Substitute sk=3s_k = 3:

14
N=124(36+6⋅32⋅3+3⋅32⋅32+8⋅32+6⋅33).\quad N = \tfrac{1}{24}\bigl(3^6 + 6\cdot 3^{2}\cdot 3 + 3\cdot 3^{2}\cdot 3^{2} + 8\cdot 3^{2} + 6\cdot 3^{3}\bigr).

Derivation hint: each sk↦3s_k \mapsto 3; e.g. s16→36s_1^6 \to 3^6, s12s4→32⋅3=33s_1^2 s_4 \to 3^2 \cdot 3 = 3^3, s32→32s_3^2 \to 3^2.

Step 3 — Arithmetic. 36=7293^6 = 729; 6⋅27=1626\cdot 27 = 162; 3⋅81=2433\cdot 81 = 243; 8⋅9=728\cdot 9 = 72; 6⋅27=1626\cdot 27 = 162. Sum =729+162+243+72+162=1368= 729+162+243+72+162 = 1368.

N=136824=57 distinct 3-colorings of the cube’s faces.\boxed{N = \dfrac{1368}{24} = 57 \text{ distinct 3-colorings of the cube's faces.}}

Common Misconceptions

⚠️ Common mistake

The group acts on the domain (positions), not the colors. The cycle index is built from the action of GG on DD. A frequent error is to permute colors, or to count ∣D∣|D| as the number of group elements. Also: jk(g)j_k(g) counts cycles of length kk, not cycles indexed by kk — fixed points are 1-cycles contributing to s1s_1, and they are easy to miss. Finally, PET requires a group action (closure, inverses, identity); applying it to a set of symmetries that is not closed under composition gives nonsense.

Counterexample

Naïve division by group size fails. One might guess the number of 2-colorings of n=4n=4 beads is "total colorings / group size" =24/4=4= 2^4 / 4 = 4. But the correct answer from [12] is 6. The simple division ignores colorings with nontrivial stabilizers (e.g. the all-black coloring is fixed by all rotations, so it must not be divided away). PET (via Burnside's averaging) correctly handles these symmetric colorings; plain division by ∣G∣|G| only works when the action on colorings is free, which is rarely the case.

Connections

connection

PET sits atop a tower of related counting tools. It is a weighted, refined generalization of <span class="wiki-link-missing">Burnside&#39;s Lemma</span> (the Cauchy–Frobenius orbit-counting lemma), which it reduces to when all weights are set to 1. The cycle index is a central object in Symmetric Functions theory — substituting power-sum symmetric functions pkp_k links PET directly to the power-sum basis of the ring of symmetric functions. The unweighted formula is an instance of the broader orbit-counting problem studied in Group Theory and Group Actions. For counting labeled vs. unlabeled structures, PET interfaces with exponential and ordinary <a class="wiki-link" href="/wiki/generating-functions">Generating Functions</a>, and de Bruijn's generalization extends it to symmetries acting on both domain and range.

Real-World Applications

🔬 Where this appears

Historical Development

⏳ History
- 1845: Cauchy establishes results on the average number of fixed points of permutations, an early form of orbit counting.
- 1887: Frobenius formalizes the orbit-counting lemma; the result is later misattributed to Burnside.
- 1897: Burnside includes the orbit-counting lemma in his text on finite groups; it becomes known (mistakenly) as "Burnside's Lemma."
- 1927: J. Howard Redfield publishes "The Theory of Group-Reduced Distributions," containing the cycle-index enumeration theorem — overlooked for decades.
- 1937: George Pólya publishes his landmark paper applying the cycle-index method to graphs, trees, and chemical compounds, giving the theorem wide visibility.
- 1948: Pólya's results spread through chemistry and combinatorics literature; isomer-counting becomes a showcase application.
- 1959: N. G. de Bruijn generalizes PET to permutation groups acting on both domain and range, vastly broadening its scope.
- 1960: Frank Harary rediscovers Redfield's 1927 paper, prompting the modern dual attribution Redfield–Pólya Theorem.
- 1973–87: Harary, Palmer, and others systematize graphical enumeration via cycle indices in standard combinatorics texts.

Subtleties & Limits

open-problem

PET reduces counting to computing cycle indices, but obtaining ZGZ_G for complex groups can itself be hard. For the symmetric group SnS_n, the cycle index involves a sum over all integer partitions of nn — explosive growth — and counting non-isomorphic graphs, though solved in principle by PET, has no known simple closed form for the count gng_n (its asymptotics involve 2(n2)/n!2^{\binom{n}{2}}/n! as a leading term). PET also assumes finite DD and RR and a fixed group action; problems where symmetries act jointly on positions and colors, or where the structure is dynamic, require de Bruijn's extension or further machinery. Determining cycle indices for automorphism groups of arbitrary structures is closely tied to the graph isomorphism problem, whose precise complexity remains a celebrated open question.

See Also

Burnside's Lemma · Cycle Index · Group Actions · Symmetric Functions · Generating Functions · Graph Enumeration · Necklace Counting · de Bruijn's Generalization