Burnside引理与Polya定理 学习笔记

原文链接www.cnblogs.com/zhouzhendong/p/Burnside-Polya.html

问题模型

  有一个长度为 $n$ 的序列,序列中的每一个元素有 $m$ 种取值。

  如果两个序列循环同构,那么我们称这两个序列等价。

  求两两不等价的序列个数。

Burnside引理

  假设有若干个置换 $P_1,P_2,\cdots$ ,设由这些置换生成的置换群为 $Q$ 。如果序列 A 可以通过一个 $Q$ 中的置换变成序列 B,那么我们认为 A 和 B 等价。

  对于一个置换 $P$ ,如果序列 A 经过置换 $P$ 之后不变,那么我们称序列 A 为置换 $P$ 的一个“不动点”。令 $C(P)$ 表示置换 $P$ 的不动点数。

  Burnside引理:

    对于一个置换群 $Q$,两两不同的序列个数为

$$\frac 1 {|Q|} \sum_{P\in Q} C(P)$$

Polya定理

  Polya定理是Burnside引理的一个具体化。

  Polya定理:

    假设置换 $P$ 由 $k$ 个轮换组成,那么 

$C(P) = m^{k}$ 

  至此,我们已经可以解决最开始的问题模型了。

$$ans = \frac 1n \sum_{i=0}^{n-1} m ^ {\gcd(i,n)}$$

例题 - POJ2409 Let it Bead

题意

  有一种项链有 $n$ 个珠子,珠子有 $m$ 种。项链是环,但是没有起始点,也没有方向。也就是说,如果两个项链循环同构或者对称翻转后循环同构,那么他们等价。

  求不等价的项链种数。

  $n\cdot m \leq 32$ 。

题解

  这类问题,我们首先考虑一下置换群 $Q$ 中有哪些元素。

  1. 旋转。可以旋转 $0\cdots n-1$ 次,故有 $n$ 种方案。这部分等价于本文开头的问题模型。

  2. 翻转。由于任何形如 “旋转 + 翻转 + 旋转” 的操作都可以由一次翻转得到,所以本质上只有翻转操作。不管 $n$ 是奇数还是偶数,都有 $n$ 种本质不同的翻转操作,而且这些操作与之前的旋转也不同。

  所以,置换群 $Q$ 中共有 $2n$ 种置换。

  现在我们考虑求 $\sum C(P)$ 。

  对于 1. 旋转,之前已经提到,这种情况对答案的贡献就是

$$\sum_{i=0}^{n-1} m ^ {\gcd(i,n)}$$

  对于 2. 翻转,我们分 $n$ 的奇偶性进行讨论:

  n 为奇数:只能将一个顶点与其对边中点所在直线作为对称轴进行翻转:一共有 $n$ 种方案,每种方案都由 $\frac {n+1} 2$ 个轮换构成,故对答案的贡献是 

$$n \cdot m ^ {\frac {n+1} 2}$$

  n 为偶数:有两种情况:

    1. 以两个顶点所在直线为对称轴进行翻转:一共有 $\frac n 2$ 种方案,每种方案都由 $\frac{n+2} 2$ 个轮换构成,故对答案的贡献是

$$\frac n2 \cdot m ^ {\frac {n+2} 2} $$

    2. 以两条边的中点所在直线为对称轴进行翻转:一共有 $\frac n 2$ 种方案,每种方案都由 $\frac{n} 2$ 个轮换构成,故对答案的贡献是

$$\frac n2 \cdot m ^ {\frac {n} 2} $$

  至此,问题解决。

代码

#include <algorithm>
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
using namespace std;
typedef long long LL;
int n,m;
LL Pow(LL x,LL y){
	LL ans=1;
	for (;y;y>>=1,x*=x)
		if (y&1)
			ans*=x;
	return ans;
}
int gcd(int a,int b){
	return b?gcd(b,a%b):a;
}
int main(){
	while (1){
		cin>>m>>n;
		if (!n&&!m)
			break;
		LL ans=0;
		for (int i=0;i<n;i++)
			ans+=Pow(m,gcd(n,i));
		if (n&1)
			ans+=Pow(m,(n+1)>>1)*n;
		else {
			ans+=Pow(m,n>>1)*(n>>1);
			ans+=Pow(m,(n+2)>>1)*(n>>1);
		}
		ans/=n<<1;
		cout<<ans<<endl;
	}
	return 0;
}

  

 

posted @ 2019-04-21 16:57 -zhouzhendong- 阅读(...) 评论(...) 编辑 收藏
希望