前言:若您手中有一本《组合数学》,那么请翻到P330看书吧,书上讲的是真的好。本文也基本参考该书讲解】
本篇博客算是后期两周内会写的一篇关于Polya定理的专题的前置知识】
目录
一、置换
1、定义
【描述可能不那么严谨】我们假设有映射函数,对1~n中的每个数都唯一对应1~n中的某个数,并且每个数都被唯一对应,是满射。为了更加清楚地表示这一映射关系,我们用一个2*n的阵型来表示:
这样的映射关系就是置换。举个例子:对于n=3的时候,有以下6种置换:
2、合成运算
设有置换:
,那么两个置换的合成运算就是:
两两上下对应交换位置即可。这就是置换的合成运算。
并且我们可以发现:
这一点要特别注意。
同时,我们也经常用幂次来表示自身置换的次数。比如:置换。
3、置换群【这个跟Polya定理关系比较大
通过前面的讲解,我们已经明白了置换的含义。那么置换群呢?【下文参考《组合数学,其实建议看书会好理解一些】
我们定义的所n!个置换构成的集合为。如果中的置换的非空子集满足一下三条性质,则为的置换群:
1)封闭性:对于G中的所有置换f与g,也必定在集合中。
2)单位元:中的恒等置换属于G。即G包含恒等置换。
3)逆元的封闭性:对于G中的每一个置换f,也属于G。这一点可以根据置换有循环这一点来判定。
这几个性质的判定其实不会很难,往往一个置换所有阶的置换就是一个置换群。
二、置换快速幂运算【参考潘震皓的讲义
1、一般情况
设有置换【这是轮换形式】,那么就有:
对于这个结果,我们可以将它写成轮换形式,也就是从1~6对每个数都走一遍上下对应的轮换:
【2对应5,5对应3,3对应2,这就是一组循环。1,4,6同理.循环内部顺序可以打乱
同理,我们得到T的几个幂次的轮换:【建议手推哦OuO
可以发现一件事情:对于,轮换形式下,循环的个数就是,并且【每个循环】分别是T中的下标的元素组成。
我们继续轮换下去,还可以发现:。所以一定有:【-1和+1是应付下标从0开始这一点,后面置换群还会提到。
2、互质的时候
我们继续看,会发现一些特殊的情况:当的时候,结果一定会是一个大循环对吧?那么如何快速得到这个循环呢?
我们可以再看两个实验【其实从这里开始基本上都是在搬人家讲义的内容了QuQ】
设,则:
其中,是把其奇数项取出来,再接上偶数项,是取出mod3为1,2,0的项……
所以就可以得到一个定理:令数组,并且有则:
看起来有点复杂。但是可以画个环理解理解。【原图来自潘震皓的讲义】
注意,这样一种行走方式是对于轮换形式的,如果我们要求两个置换,那么一定要先转换成轮换形式,这样早做过后再逆变换回来。逆变换的方式可以自己手推一下。
这个东西有什么用呢?放个很经典的例题:
洛谷P2227 [HNOI2001] 洗牌机
题解
很显然,每一次洗牌就是一次置换的合成操作。一共洗s次牌,就相当于给你的置换状态,让你求的置换状态。
首先题目保证了n是奇数,我们的指数又一定是2的幂次,所以保证了两者互质,我们就可以直接套用上面的那个轮换形式下的递推公式啦。
但是明显我们不能直接跑2^s这么多次。所以我们要降低指数。通过对置换群的了解,我们完全可以将指数对n取余。
代码很短的,直接看吧。
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#define maxn 1005
using namespace std;
typedef long long ll;
int read() {
int x = 0, f = 1, ch = getchar();
while(!isdigit(ch)) {if(ch == '-') f = -1; ch = getchar();}
while(isdigit(ch)) x = (x << 1) + (x << 3) + ch - '0', ch = getchar();
return x * f;
}
int n, s, a[maxn], b[maxn];
signed main() {
n = read(), s = read();
for(int i = 1; i <= n; i++) a[i] = read();
for(int i = 1, j = 1; i <= n; i++, j = a[j]) b[i] = a[j];//现在b里面是a的轮换形式
register int k = 1;
for(int i = 1; i <= s; i++) k = k * 2 % n;//看到底要多少次幂,即走多少步
for(int i = 1, j = 1; i <= n; i++, j = (j + k - 1) % n + 1) a[j] = b[i];//这里是在走那个圈圈
for(int i = 2, j = a[1]; i <= n + 1; i++, j = b[j]) b[j] = a[(i - 1) % n + 1];//反轮换形式
for(int i = 1; i <= n; i++) printf("%d ", b[i]);
return 0;
}
呼……群论真的是一个很大的知识点啊。关于单纯置换的题目还有一个:洛谷P4161 [SCOI2009] 游戏
迎评:)
——End——