Polya计数

将波利亚计数定理整理在这里,作为一个总结和介绍,也方便以后复习

为什么学习Polya计数定理

通过Polya计数定理,我们可以计算等价类的数量,比如下面这个问题:

mm种颜色给一个正方形染色,如果正方形可以自由转动,求染色方案数

让我们从一些概念开始

1 等价关系

1.1 等价关系的定义

假设VV是一个集合,SS是定义在VV上的一个关系,若SS有如下性质:

  • 自反性
  • 传递性
  • 对称性

那么, SS就是一个等价关系

aabb有关系SS,可以记为aSbaSb

假设定义关系SS,图形aa可以旋转得到bb \Leftrightarrow aSbaSb
例如图中的1方块_12方块_2具有关系SS,即他们可以通过旋转得到彼此,而1方块_12方块_2则没有关系SS

1.2 等价类

通过上图,可以看出来:1方块_12方块_2是同一类的,而1方块_12方块_2则是另外两类,于是可以想到集合VV上的等价关系SS将集合的元素划分到不同的类中,我们把它称为等价类

而包含元素aa的等价类则是由满足aSbaSb的所有元素bb组成的(当然也包含元素aa),即C(a)=bV aSbC(a)={b\in V | \ aSb}

仔细想一想,不难发现两个不同等价类是不相交的

2 置换群

2.1 置换群的定义

假设A=1,2,,nA={1,2,\dotsc, n},通过置换,将AA中的元素重新排列,得到另一个排列a1,a2,,ana_1, a_2, \dotsc, a_n,可以把这个过程写成

(12na1a2an) \left ( \begin{matrix} 1 & 2 & \dotsc & n \newline a_1 & a_2 & \dotsc & a_n \end{matrix} \right )

所以,可以将置换看成一个双射函数f:1,2,,n1,2,,nf:{1,2,\dotsc, n} \rightarrow {1,2,\dotsc, n}

置换之间也可进行合成运算

π1=(12344213)  π2=(12342143) \pi_1 = \left ( \begin{matrix} 1 & 2 & 3 & 4 \newline 4 & 2 & 1 & 3 \end{matrix} \right ) \ \ \pi_2 = \left ( \begin{matrix} 1 & 2 & 3 & 4 \newline 2 & 1 & 4 & 3 \end{matrix} \right )
π1π2=(12342431) \pi_1 \circ \pi_2 =\left ( \begin{matrix} 1 & 2 & 3 & 4 \newline 2 & 4 & 3 & 1 \end{matrix} \right ) π1π2(1)=π1(π2(1))=2\pi_1 \circ \pi_2 (1) = \pi_1(\pi_2(1))=2

而在置换集合XX上定义的群则称为置换群,置换群也是群也要满足群的性 质

  • 运算封闭 π1X,π2Xπ1π2X\pi_1 \in X, \pi_2 \in X \Rightarrow \pi_1 \circ \pi_2 \in X
  • 满足结合律
  • 有单位元,置换I(a)=aI(a)=a
  • 有逆元,ππ1=I\pi \circ \pi^{-1} = I

但是置换与计算不同染色模式数量有什么关系?图形的旋转等变换都可以用一个只置换来表示

我们把1-4按顺时针放到正方形的四个方块中


π1=(12341234)\pi_1=\left(\begin{matrix} 1 & 2 & 3 & 4 \newline 1 & 2 & 3 & 4 \end{matrix}\right)

旋转9090^\circ

π2=(12344123)\pi_2=\left(\begin{matrix} 1 & 2 & 3 & 4 \newline 4 & 1 & 2 & 3 \end{matrix}\right)

2.2 置换群衍生的等价关系

假设G=(X,)G=(X,\circ),通过置换πX\pi\in X,元素aa可以被置换为其他元素π(a)=b\pi(a)=b,可以发现这里也存在一种等价关系, 可以定义等价关系SS:

aSbπG,π(a)=b aSb \Leftrightarrow \exist \pi \in G, \pi(a)=b

πG也可以来表示πX\pi \in G 也可以来表示 \pi \in X

而包含aa的等价类C(a)C(a)是由所有满足aSbaSb的元素bb组成的,因此我们也可以说 C(a)=π(a)  πG C(a)={\pi(a) \ | \ \pi \in G}

如果A=1,2,3A={1,2,3}
π1=(123123)\pi_1=\left(\begin{matrix} 1 & 2 & 3 \newline 1 & 2 & 3 \end{matrix}\right)
π2=(123321)\pi_2=\left(\begin{matrix} 1 & 2 & 3 \newline 3 & 2 & 1 \end{matrix}\right)
X=π1,π2X={\pi_1, \pi_2},可以验证G=X,G={X,\circ}是一个置换群,1,3{1,3}是一个等价类,2{2}是一个等价类

可以证明自反、对称、传递是满足的

3 伯恩赛德引理

这个定理将给出置换群衍生出的等价关系的(不同)等价关系数量的计数方法

3.1 伯恩赛德引理

假设GG是一个置换群,aaAA中的元素,如果π(a)=a\pi(a)=a,则称AA中的元素aa在置换π\pi下是不变的(Invariant),Inv(π)Inv(\pi)表示不变元素的数量

定理   假设GG是一个集合AA的置换群,设SSGG上衍生出来的等价关系,那么SS中的等价类的数量由下式给出:

1GπGInv(π) \frac{1}{|G|} \sum_{\pi\in G}Inv(\pi)

G|G|是置换群中置换的数量,让我们来用一下这个定理,假设GG下面的置换组成

π1=(12341234),π2=(12342134), \pi_1=\left(\begin{matrix} 1 & 2 & 3 & 4\newline 1 & 2 & 3 & 4 \end{matrix}\right) , \pi_2=\left(\begin{matrix} 1 & 2 & 3 & 4 \newline 2 & 1 & 3 & 4 \end{matrix}\right) ,

π3=(12341243),π4=(12342143), \pi_3=\left(\begin{matrix} 1 & 2 & 3 & 4\newline 1 & 2 & 4 & 3 \end{matrix}\right) , \pi_4=\left(\begin{matrix} 1 & 2 & 3 & 4 \newline 2 & 1 & 4 & 3 \end{matrix}\right) ,

可以看出来,有两个等价类分别是1,2{1,2}3,4{3, 4},等价类的数量为22

可以验证Inv(π1)=4,Inv(π2)=2,Inv(π3)=2,Inv(π4)=0Inv(\pi_1)=4,Inv(\pi_2)=2,Inv(\pi_3)=2,Inv(\pi_4)=0

等价类的数量=14(4+2+2+0)=2 等价类的数量=\frac{1}{4}(4+2+2+0)=2

3.2 伯恩赛德引理的证明

Part 1

这个证明可以跳过,对后面阅读没有什么影响

对于AA中的元素aa,我们可以定义稳定算子(Stabilizer)的概念,St(a)St(a)是保持元素aa不变的置换的集合,St(a)=πGπ(a)=aSt(a)={\pi \in G | \pi(a) = a}

π1=(123123),π2=(123231),π3=(123312), \pi_1=\left(\begin{matrix} 1 & 2 & 3 \newline 1 & 2 & 3 \end{matrix}\right) , \pi_2=\left(\begin{matrix} 1 & 2 & 3 \newline 2 & 3 & 1 \end{matrix}\right) , \pi_3=\left(\begin{matrix} 1 & 2 & 3 \newline 3 & 1 & 2 \end{matrix}\right) , St(2)=π1St(2) = {\pi_1}
包含22的等价类C(2)=π1(2),π2(2),π3(2)=1,2,3C(2)={\pi_1(2),\pi_2(2),\pi_3(2)}={1, 2, 3}

可以发现St(2)×C(2)=G=3|St(2)| \times |C(2)| = |G| = 3,这不是巧合

引理   假设GG是一个置换群,而且aAa \in A,那么, St(a)×C(a)=G |St(a)| \times |C(a)| = |G|

让我们来简单证明这个定理

假设C(a)=b1,b2,,brC(a)={b_1,b_2,\dotsb, b_r},那么存在一个π1\pi_1,满足π1(a)=b1\pi_1(a)=b_1,同样存在π2(a)=b2\pi_2(a)=b_2,\dotsb,设P=π1,π2,,πrP={\pi_1, \pi_2, \dotsc, \pi_r},而且P=C(a)|P|=|C(a)|,于是转化为证明St(a)×P=G|St(a)|\times |P|=|G|

GG中选一个置换π\pi, π\pi一定将aa置换成一个元素,设该元素为bkb_k,即π(a)=bk\pi(a)=b_k

另一个方面,πk(a)=bk\pi_k(a)=b_k,因此πk1π(a)=a\pi^{-1}_k \circ \pi(a)=a,即πk1πSt(a)\pi^{-1}_k \circ \pi \in St(a)

πk(πk1π)=π \pi_k \circ (\pi_k^{-1} \circ \pi) = \pi 这说明GG中的元素π\pi可以由St(a)St(a)C(a)C(a)中的元素合成得到

然后我们再来说明这是唯一的

假设π=πkγ=πlδ\pi=\pi_k\circ \gamma=\pi_l \circ \deltaγ\gammaδ\delta都在St(a)St(a)中,因此πkγ(a)=bk=πlδ(a)=bl\pi_k \circ \gamma(a)=b_k=\pi_l \circ \delta (a)=b_l,因此l=kl=k,这种组合方式是唯一的,因此St(a)×C(a)=G|St(a)| \times |C(a)| = |G|

Part 2

让我们继续证明伯恩赛德引理

假设A=1,2,,nA={1,2,\dotsc,n}, G=π1,π2,,πmG={\pi_1, \pi_2, \dotsc, \pi_m},让我们看看下面这个式子的是否成立

Inv(π1)+Inv(π2)++Inv(πm)=St(1)+St(2)++St(n) |Inv(\pi_1)|+|Inv(\pi_2)|+\dotsb+|Inv(\pi_m)|=|St(1)|+|St(2)|+\dotsb+|St(n)|

可以发现都计数的是满足π(a)=a\pi(a)=a序对<π,a><\pi,a>的个数,因此等式是成立的

将等式两边同时除以|G|,左边就得到了伯恩赛德引理的型式,而右边,根据Part1中的引理可以转化为

1C(1)+1C(2)++1C(n)    (1) \frac{1}{C(1)}+\frac{1}{C(2)}+\dotsb+\frac{1}{C(n)} \ \ \ \ (1)

C(a)C(a)是一个等价类,两个等价类的交集要么是全集要么是空集,假设C(a)=b1,b2,,brC(a)={b_1,b_2,\dotsc, b_r},那么,

1C(b1)+1C(b2)++1C(br)=1r+1r++1r=1 \frac{1}{C(b_1)}+\frac{1}{C(b_2)}+\dotsb+\frac{1}{C(b_r)}=\frac{1}{r}+\frac{1}{r}+\dotsb+\frac{1}{r}=1

于是(1)式计数的就是等价类的个数

4 等价着色

伯恩赛德引理和计算等价着色数量有什么关系呢?其实伯恩赛德引理不能直接帮助我们计算等价着色的数量,但是我们可以适当的定义等价关系,然后再计算等价着色的数量

假设C(R,D)C(R,D)使用RR中的颜色,来给DD中的元素染色的着色集合,

假设RR是黑色和白色,DD是图中的节点,那么C(R,D)C(R,D)是给图中节点染色的集合 C(R,D)=C1,C2,,C16C(R,D)={C_1,C_2,\dotsb, C_{16}}

定义在C(R,D)C(R,D)上的置换π\pi^*,将一种着色变为另一种着色,更形式一点,假定ff是一个着色,那么新着色πf\pi^*f定义为(πf)(a)(\pi^*f)(a)f(π(a))f(\pi(a))

同样,也可以定义C(R,D)C(R,D)上的置换群,G=π  πGG^*={\pi^* \ | \ \pi \in G}GGGG^*的元素数量一样,因为π\piπ对应于一种相同旋转\pi^*对应于一种相同旋转GG^*也可以衍生出等价关系SS^*,而SS^*就是我们关注的。

假设f,gf,g是两个着色,如果有一个πG\pi \in G(πG\pi^*\in G^*),πf=g\pi^*f=g,我们就认为f,gf,g是等价的

我们要求的SS^*划分出来的等价类的数量

从上图可以看出来GG^*的大小是33,设G={π1,π2,π3}G^*=\{\pi_1^*,\pi_2^*,\pi_3^*\},分别对应逆时针旋转0,120,2400^\circ,120^\circ, 240^\circ

π1=(C1C2C3C4C5C6C7C8C9C10C11C12C13C14C15C16C1C2C3C4C5C6C7C8C9C10C11C12C13C14C15C16) \pi_1^*=\left(\begin{matrix} C_1 & C_2 & C_3 & C_4 &C_5 & C_6 &C_7 & C_8 &C_9 & C_{10} & C_{11} & C_{12} & C_{13} & C_{14} & C_{15} & C_{16} \newline C_1 & C_2 & C_3 & C_4 &C_5 & C_6 &C_7 & C_8 &C_9 & C_{10} & C_{11} & C_{12} & C_{13} & C_{14} & C_{15} & C_{16} \end{matrix}\right)

π2=(C1C2C3C4C5C6C7C8C9C10C11C12C13C14C15C16C1C3C4C2C5C8C6C7C10C11C9C13C14C12C15C16) \pi_2^*=\left(\begin{matrix} C_1 & C_2 & C_3 & C_4 &C_5 & C_6 &C_7 & C_8 &C_9 & C_{10} & C_{11} & C_{12} & C_{13} & C_{14} & C_{15} & C_{16} \newline C_1 & C_3 & C_4 & C_2 &C_5 & C_8 &C_6 & C_7 &C_{10} & C_{11} & C_{9} & C_{13} & C_{14} & C_{12} & C_{15} & C_{16} \end{matrix}\right)

π3=(C1C2C3C4C5C6C7C8C9C10C11C12C13C14C15C16C1C4C2C3C5C7C8C6C11C9C10C14C12C13C15C16) \pi_3^*=\left(\begin{matrix} C_1 & C_2 & C_3 & C_4 &C_5 & C_6 &C_7 & C_8 &C_9 & C_{10} & C_{11} & C_{12} & C_{13} & C_{14} & C_{15} & C_{16} \newline C_1 & C_4 & C_2 & C_3 &C_5 & C_7 &C_8 & C_6 &C_{11} & C_{9} & C_{10} & C_{14} & C_{12} & C_{13} & C_{15} & C_{16} \end{matrix}\right)

Inv(π1)=16,Inv(π2)=4,Inv(π2)=4Inv(\pi_1^*)=16,Inv(\pi_2^*)=4,Inv(\pi_2^*)=4

根据伯恩赛德引理,

等价着色的个数=13(16+4+4)=8 等价着色的个数=\frac{1}{3}(16+4+4)=8 可以从上图中验证这个结果的正确性,为了与SS中的等价类区分,我们称SS^*中等价类为模式

通过枚举C(R,D)C(R,D)中的元素,再枚举π\pi^*计算Inv(π)Inv(\pi^*)过于麻烦,可以直接计算Inv(π)Inv(\pi^*),于是有了下面的定理

5 Polya定理的特殊情况

5.1 循环分解

让我们再回到置换,设置换π1=(1234523451)\pi_1=\left(\begin{matrix} 1 & 2 & 3 & 4 & 5\newline 2 & 3 & 4 & 5 & 1 \end{matrix}\right) 让我们不断地重复这个置换,那么可以找到循环节1234511\rightarrow 2 \rightarrow 3 \rightarrow 4 \rightarrow 5 \rightarrow 1,那我们就把π\pi简写(12345)(12345)

再比如π2=(123456516324)\pi_2=\left(\begin{matrix} 1 & 2 & 3 & 4 & 5 & 6 \newline 5 & 1 & 6 & 3 & 2 & 4 \end{matrix}\right) 可以简写成(152)(364)(152)(364),这一过程被称为循环分解

定义cyc(π)cyc(\pi)π\pi循环分解中循环节的数量,比如cyc(π1)=1cyc(π2)=2cyc(\pi_1)=1,cyc(\pi_2)=2

5.2 Polya 计数的特殊情况

定理   假设GG是集合DD的置换群,且C(R,D)C(R,D)是使用RR中颜色给DD中元素染色的集合,RR中颜色的数量的mm,C(R,D)C(R,D)中不同着色数量(SS^*中等价类的数量)由下式给出:

1G(mcyc(π1)+mcyc(π2)++mcyc(πk)) \frac{1}{|G|} (m^{cyc(\pi_1)}+m^{cyc(\pi_2)}+\dotsb+m^{cyc(\pi_k)}) 其中G=π1,π2,,πkG={\pi_1, \pi_2, \dotsb, \pi_k}

让我们用这个定理再算一次等价着色的个数,假定将中间的点记为11 0旋转,π1=(1)(2)(3)(4)0^\circ旋转,\pi_1=(1)(2)(3)(4)
120π2=(1)(234)120^\circ, \pi_2=(1)(234)
240π3=(1)(234)240^\circ, \pi_3=(1)(234)

等价着色的个数=13(24+22+22)=8 等价着色的个数=\frac{1}{3}(2^4+2^2+2^2)=8

5.3 循环指标PGP_G

如果可以用一个生成函数来描述一个置换群,那将很有用,这将引出更一般的波利亚计数定理,因此有了循环指标PGP_G的概念

假设π=(1)(23)(45)(6)(789)(10)\pi=(1)(23)(45)(6)(789)(10),可以将他编码为x13x22x3x_1^3x_2^2x_3

可以看出来编码方式,xix_i代表长度为ii的循环节,xix_i的幂次bb代表有多少个长度为ii的循环节

于是可以用每个置换编码的和再除以G|G|,

PG(x1,x2,,xk)=1GπGx1b1x2b2xkbk P_G(x_1,x_2,\dotsc,x_k)=\frac{1}{G}\sum_{\pi\in G} x_1^{b_1}x_2^{b_2}\dotsb x_k^{b_k} 来编码一个置换群

对5.2的例子中的置换群进行编码,PG=13(x14+x1x3+x1x3)P_G=\frac{1}{3}(x_1^4+x_1x_3+x_1x_3)

而且,波利亚计数的特殊情况就是PGP_Gx1,x2,,xkx_1,x_2,\dotsc,x_k都取mm的情况

5.4 Polya计数特殊情况的证明

同样可以跳过这一部分

将Polya计数特殊情况的公式与伯恩赛德引理对比,因此只要证明mcyc(π)=Inv(π)m^{cyc(\pi)}=Inv(\pi^*)

考虑一个具体的情况,假定π=(1)(23)(456)(7)\pi=(1)(23)(456)(7),经过置换2,3的颜色交换,4,5,64,5,6将颜色给下一个,若ff是一个着色,要让π(f)=f\pi^*(f)=f,那么ff在每个循环节中的染色要一致,因此一共有mcyc(π)m^{cyc(\pi)}这么多ff

6 Polya计数定理

为什么还有这个?
如果要求必须要用黑色kk次,或者至少一次黑色,上面的定理就不好用了

6.1 给颜色加权

假设着色ffDD染色,使用红色t1t_1次,绿色t2t_2次,蓝色t3t_3

假定红色的权值wwrr,绿色gg,蓝色bb

那么着色ff的权W(f)=rt1gt2bt3W(f)=r^{t_1}g^{t_2}b^{t_3},而如果染色f,gf,g是等价的,那么他们的权是相等的

于是,我们也可以定义一个着色等价类(模式)的权等于不同等价类的权之和,把这个称为着色等价类(模式)的权,或者模式集合的目录,因为他概括了模式的染色方式

我们再来考虑5.2中的例子,模式目录是

w4+2w3b+2w2b2+2wb3+b4 w^4+2w^3b+2w^2b^2+2wb^3+b^4

所有系数之和等于不同等价类的个数,

如果要求用了2个白色两个黑的着色方案数,可通过w2b2w^2b^2前的系数得知是2种

下面的定理将给出计算模式目录的方法

6.2 Polya定理内容

定理   假设GG是集合DD上的置换群,而C(R,D)C(R,D)是用RR中的颜色给DD着色的集合,ww是颜色的权值,那么C(R,D)的着色的模式目录为C(R,D)的着色的模式目录为

PG(rRw(r),rRw2(r),,rRwk(r)) P_G\bigg(\sum_{r\in R}w(r),\sum_{r\in R}w^2(r),\dotsb,\sum_{r\in R}w^k(r)\bigg) 其中PG(x1,x2,,xk)P_G(x_1,x_2,\dotsb,x_k)是循环指标

再回到之前的例子
我们已经知道PG=13(x14+x1x3+x1x3)P_G=\frac{1}{3}(x_1^4+x_1x_3+x_1x_3)
假定我们用白,黑染色,权值分别为w,bw,b
rRw(r)=w+b,rRw3(r)=w3+b3 \sum_{r\in R}w(r)=w+b,\sum_{r\in R}w^3(r)=w^3+b^3 带入PGP_G 13[(w+b)4+2(w+b)(w3+b3)] \frac{1}{3}[(w+b)^4+2(w+b)(w^3+b^3)]

通过对颜色权值的简单赋值就能的要一些有用的结论

  • 如果想知道以用有多少种染色模式,那么只要把所有的权值都设为11
  • 如果不想用某一颜色,就把颜色设为00,其他颜色设为11
  • 如果想知道用某个颜色kk次的模式数,就把其他颜色设为11,求改颜色权值幂次为kk前的系数
  • . . .

QAQ, 终于写完了。。。