Polya计数#
将波利亚计数定理整理在这里,作为一个总结和介绍,也方便以后复习
为什么学习Polya计数定理
通过Polya计数定理,我们可以计算等价类的数量,比如下面这个问题:
用m种颜色给一个正方形染色,如果正方形可以自由转动,求染色方案数
让我们从一些概念开始
1 等价关系#
1.1 等价关系的定义#
假设V是一个集合,S是定义在V上的一个关系,若S有如下性质:
那么, S就是一个等价关系
a和b有关系S,可以记为aSb
假设定义关系S,图形a可以旋转得到b ⇔ aSb
例如图中的方块1和方块2具有关系S,即他们可以通过旋转得到彼此,而方块1和方块2则没有关系S

1.2 等价类#
通过上图,可以看出来:方块1和方块2是同一类的,而方块1和方块2则是另外两类,于是可以想到集合V上的等价关系S将集合的元素划分到不同的类中,我们把它称为等价类
而包含元素a的等价类则是由满足aSb的所有元素b组成的(当然也包含元素a),即C(a)=b∈V∣ aSb
仔细想一想,不难发现两个不同等价类是不相交的
2 置换群#
2.1 置换群的定义#
假设A=1,2,…,n,通过置换,将A中的元素重新排列,得到另一个排列a1,a2,…,an,可以把这个过程写成
(1a12a2……nan)
所以,可以将置换看成一个双射函数f:1,2,…,n→1,2,…,n
置换之间也可进行合成运算
π1=(14223143) π2=(12213443)
π1∘π2=(12243341)
π1∘π2(1)=π1(π2(1))=2
而在置换集合X上定义的群则称为置换群,置换群也是群也要满足群的性 质
- 运算封闭 π1∈X,π2∈X⇒π1∘π2∈X
- 满足结合律
- 有单位元,置换I(a)=a
- 有逆元,π∘π−1=I
但是置换与计算不同染色模式数量有什么关系?图形的旋转等变换都可以用一个只置换来表示
我们把1-4按顺时针放到正方形的四个方块中

π1=(11223344)
旋转90∘
π2=(14213243)
2.2 置换群衍生的等价关系#
假设G=(X,∘),通过置换π∈X,元素a可以被置换为其他元素π(a)=b,可以发现这里也存在一种等价关系, 可以定义等价关系S:
aSb⇔∃π∈G,π(a)=b
π∈G也可以来表示π∈X
而包含a的等价类C(a)是由所有满足aSb的元素b组成的,因此我们也可以说
C(a)=π(a) ∣ π∈G
如果A=1,2,3
π1=(112233)
π2=(132231)
X=π1,π2,可以验证G=X,∘是一个置换群,1,3是一个等价类,2是一个等价类
可以证明自反、对称、传递是满足的
3 伯恩赛德引理#
这个定理将给出置换群衍生出的等价关系的(不同)等价关系数量的计数方法
3.1 伯恩赛德引理#
假设G是一个置换群,a是A中的元素,如果π(a)=a,则称A中的元素a在置换π下是不变的(Invariant),Inv(π)表示不变元素的数量
定理 假设G是一个集合A的置换群,设S是G上衍生出来的等价关系,那么S中的等价类的数量由下式给出:
∣G∣1π∈G∑Inv(π)
∣G∣是置换群中置换的数量,让我们来用一下这个定理,假设G下面的置换组成
π1=(11223344),π2=(12213344),
π3=(11223443),π4=(12213443),
可以看出来,有两个等价类分别是1,2,3,4,等价类的数量为2
可以验证Inv(π1)=4,Inv(π2)=2,Inv(π3)=2,Inv(π4)=0
等价类的数量=41(4+2+2+0)=2
3.2 伯恩赛德引理的证明#
Part 1#
这个证明可以跳过,对后面阅读没有什么影响
对于A中的元素a,我们可以定义稳定算子(Stabilizer)的概念,St(a)是保持元素a不变的置换的集合,St(a)=π∈G∣π(a)=a,
π1=(112233),π2=(122331),π3=(132132),
St(2)=π1,
包含2的等价类C(2)=π1(2),π2(2),π3(2)=1,2,3
可以发现∣St(2)∣×∣C(2)∣=∣G∣=3,这不是巧合
引理 假设G是一个置换群,而且a∈A,那么,
∣St(a)∣×∣C(a)∣=∣G∣
让我们来简单证明这个定理
假设C(a)=b1,b2,⋯,br,那么存在一个π1,满足π1(a)=b1,同样存在π2(a)=b2,⋯,设P=π1,π2,…,πr,而且∣P∣=∣C(a)∣,于是转化为证明∣St(a)∣×∣P∣=∣G∣
从G中选一个置换π, π一定将a置换成一个元素,设该元素为bk,即π(a)=bk,
另一个方面,πk(a)=bk,因此πk−1∘π(a)=a,即πk−1∘π∈St(a),
πk∘(πk−1∘π)=π
这说明G中的元素π可以由St(a)和C(a)中的元素合成得到
然后我们再来说明这是唯一的
假设π=πk∘γ=πl∘δ ,γ和δ都在St(a)中,因此πk∘γ(a)=bk=πl∘δ(a)=bl,因此l=k,这种组合方式是唯一的,因此∣St(a)∣×∣C(a)∣=∣G∣
Part 2#
让我们继续证明伯恩赛德引理
假设A=1,2,…,n, G=π1,π2,…,πm,让我们看看下面这个式子的是否成立
∣Inv(π1)∣+∣Inv(π2)∣+⋯+∣Inv(πm)∣=∣St(1)∣+∣St(2)∣+⋯+∣St(n)∣
可以发现都计数的是满足π(a)=a序对<π,a>的个数,因此等式是成立的
将等式两边同时除以|G|,左边就得到了伯恩赛德引理的型式,而右边,根据Part1中的引理可以转化为
C(1)1+C(2)1+⋯+C(n)1 (1)
而C(a)是一个等价类,两个等价类的交集要么是全集要么是空集,假设C(a)=b1,b2,…,br,那么,
C(b1)1+C(b2)1+⋯+C(br)1=r1+r1+⋯+r1=1
于是(1)式计数的就是等价类的个数
4 等价着色#
伯恩赛德引理和计算等价着色数量有什么关系呢?其实伯恩赛德引理不能直接帮助我们计算等价着色的数量,但是我们可以适当的定义等价关系,然后再计算等价着色的数量
假设C(R,D)使用R中的颜色,来给D中的元素染色的着色集合,
假设R是黑色和白色,D是图中的节点,那么C(R,D)是给图中节点染色的集合
C(R,D)=C1,C2,⋯,C16
定义在C(R,D)上的置换π∗,将一种着色变为另一种着色,更形式一点,假定f是一个着色,那么新着色π∗f定义为(π∗f)(a)是f(π(a))
同样,也可以定义C(R,D)上的置换群,G∗=π∗ ∣ π∈G,G和G∗的元素数量一样,因为π和π∗对应于一种相同旋转。G∗也可以衍生出等价关系S∗,而S∗就是我们关注的。
假设f,g是两个着色,如果有一个π∈G(π∗∈G∗),π∗f=g,我们就认为f,g是等价的
我们要求的S∗划分出来的等价类的数量
从上图可以看出来G∗的大小是3,设G∗={π1∗,π2∗,π3∗},分别对应逆时针旋转0∘,120∘,240∘
π1∗=(C1C1C2C2C3C3C4C4C5C5C6C6C7C7C8C8C9C9C10C10C11C11C12C12C13C13C14C14C15C15C16C16)
π2∗=(C1C1C2C3C3C4C4C2C5C5C6C8C7C6C8C7C9C10C10C11C11C9C12C13C13C14C14C12C15C15C16C16)
π3∗=(C1C1C2C4C3C2C4C3C5C5C6C7C7C8C8C6C9C11C10C9C11C10C12C14C13C12C14C13C15C15C16C16)
Inv(π1∗)=16,Inv(π2∗)=4,Inv(π2∗)=4
根据伯恩赛德引理,
等价着色的个数=31(16+4+4)=8
可以从上图中验证这个结果的正确性,为了与S中的等价类区分,我们称S∗中等价类为模式
通过枚举C(R,D)中的元素,再枚举π∗计算Inv(π∗)过于麻烦,可以直接计算Inv(π∗),于是有了下面的定理
5 Polya定理的特殊情况#
5.1 循环分解#
让我们再回到置换,设置换π1=(1223344551)
让我们不断地重复这个置换,那么可以找到循环节1→2→3→4→5→1,那我们就把π简写(12345)
再比如π2=(152136435264)
可以简写成(152)(364),这一过程被称为循环分解
定义cyc(π)是π循环分解中循环节的数量,比如cyc(π1)=1,cyc(π2)=2
5.2 Polya 计数的特殊情况#
定理 假设G是集合D的置换群,且C(R,D)是使用R中颜色给D中元素染色的集合,R中颜色的数量的m,C(R,D)中不同着色数量(S∗中等价类的数量)由下式给出:
∣G∣1(mcyc(π1)+mcyc(π2)+⋯+mcyc(πk))
其中G=π1,π2,⋯,πk
让我们用这个定理再算一次等价着色的个数,假定将中间的点记为1号
0∘旋转,π1=(1)(2)(3)(4)
120∘,π2=(1)(234)
240∘,π3=(1)(234)
等价着色的个数=31(24+22+22)=8
5.3 循环指标PG#
如果可以用一个生成函数来描述一个置换群,那将很有用,这将引出更一般的波利亚计数定理,因此有了循环指标PG的概念
假设π=(1)(23)(45)(6)(789)(10),可以将他编码为x13x22x3
可以看出来编码方式,xi代表长度为i的循环节,xi的幂次b代表有多少个长度为i的循环节
于是可以用每个置换编码的和再除以∣G∣,
PG(x1,x2,…,xk)=G1π∈G∑x1b1x2b2⋯xkbk
来编码一个置换群
对5.2的例子中的置换群进行编码,PG=31(x14+x1x3+x1x3)
而且,波利亚计数的特殊情况就是PG中x1,x2,…,xk都取m的情况
5.4 Polya计数特殊情况的证明#
同样可以跳过这一部分
将Polya计数特殊情况的公式与伯恩赛德引理对比,因此只要证明mcyc(π)=Inv(π∗)
考虑一个具体的情况,假定π=(1)(23)(456)(7),经过置换2,3的颜色交换,4,5,6将颜色给下一个,若f是一个着色,要让π∗(f)=f,那么f在每个循环节中的染色要一致,因此一共有mcyc(π)这么多f
6 Polya计数定理#
为什么还有这个?
如果要求必须要用黑色k次,或者至少一次黑色,上面的定理就不好用了
6.1 给颜色加权#
假设着色f给D染色,使用红色t1次,绿色t2次,蓝色t3次
假定红色的权值w为r,绿色g,蓝色b,
那么着色f的权W(f)=rt1gt2bt3,而如果染色f,g是等价的,那么他们的权是相等的
于是,我们也可以定义一个着色等价类(模式)的权等于不同等价类的权之和,把这个称为着色等价类(模式)的权,或者模式集合的目录,因为他概括了模式的染色方式
我们再来考虑5.2中的例子,模式目录是
w4+2w3b+2w2b2+2wb3+b4
所有系数之和等于不同等价类的个数,
如果要求用了2个白色两个黑的着色方案数,可通过w2b2前的系数得知是2种
下面的定理将给出计算模式目录的方法
6.2 Polya定理内容#
定理 假设G是集合D上的置换群,而C(R,D)是用R中的颜色给D着色的集合,w是颜色的权值,那么C(R,D)的着色的模式目录为:
PG(r∈R∑w(r),r∈R∑w2(r),⋯,r∈R∑wk(r))
其中PG(x1,x2,⋯,xk)是循环指标
再回到之前的例子
我们已经知道PG=31(x14+x1x3+x1x3)
假定我们用白,黑染色,权值分别为w,b
r∈R∑w(r)=w+b,r∈R∑w3(r)=w3+b3
带入PG
31[(w+b)4+2(w+b)(w3+b3)]
通过对颜色权值的简单赋值就能的要一些有用的结论
- 如果想知道以用有多少种染色模式,那么只要把所有的权值都设为1
- 如果不想用某一颜色,就把颜色设为0,其他颜色设为1
- 如果想知道用某个颜色k次的模式数,就把其他颜色设为1,求改颜色权值幂次为k前的系数
- . . .
QAQ, 终于写完了。。。