手镯定理
手镯定理即【Polya定理】就是一个手镯是由n个珠子构成的,现在我们有k种颜色,每个珠子可以染成这k种颜色中的任意一种,让我们来求算一下我们可以染出多少种本质上并不相同的手环来;
例题 给你顶部和底部是正 m 边形的柱体,侧面是边长为 n 的正方形,有 c 种颜色,一个侧面有 nn 个格子可以染色,总共有 mn*n 个格子可以染色。
题目
可以沿着侧面方向旋转,旋转前和旋转后染色样子一样,算同一种方案,求有多少种染色方案。
公式
L = 1 / |G| (k^c[0]+ k^c[1] + k^ c[2] + k^c[3]+k^c[4]+.......k^c[|G-1|]);
L 为方案数,G=m,k=c^(n*n),c[i]=gcd(m,i);
证明 置换群:G;置换群中置换个数:|G|;
我们喜欢用 m 阶循环这种记号来表示 置换。
m 阶循环:相邻元素可以直接到达,第 m 个位置可以到达第 1 个位置。
例如:5 阶循环:(1 5 2 4 3),1->5->2->4->3->1……,3->4->2->5->1->3……
前置技能还有:
共轭类格式:1^(c1) 2^(c2) … (n)^(cn),其中 cn 代表 n 阶置换的个数。
k 不动置换类(Zk):包含 k 的置换阶为 1
等价类(Ek):k 的轨迹
重要定理:|Ek||Zk| = |G|
Burnside 引理:群 G 是在这 n 个对象上用 m 种颜色进行染色后的方案集合上的置换群
不同等价类个数:L = sum{c1(ai)} / |G|, 1 <= i <= |G|; 其中 ai 是第 i 个置换
因为 sum{Zj} = sum{c1(ai)}, 1 <= j <= n, 1 <= i <= |G|; 其中 n 是置换对象个数
不同等价类个数还可以等于:L = sum{Zi} / |G|, 1 <= i <= n;
Polya 定理:群 G 是作用在 n 个对象上的置换群
不同等价类个数:L = sum{m^c(ai)} / |G|, 1 <= i <= |G|; 其中 c(ai) 为置换 ai 的循环节数(有多少个循环),一个置换对象有 m 种不同情况
区别和联系:
Polya 定理的群 G 是作用在 n 个对象上的置换群,相应地,Burnside 定理中的群 G 是在这 n 个对象上用 m 种颜色进行染色后的方案集合上的置换群。
Burnside 的 ai 和 Polya 的 ai 是不一样的,但是都是求在 ai 置换之后不变的染色方案数的和 / |G|,因为每种染色情况都重复出现了 |G| 次
Burnside 的 ai 作用下不变的图像 正好是对应着 Polya 的 ai 的循环节中的对象染以相同的颜色所得到的图像,所以 c1(ai) = m^c(ai)
本题思路:
Polya 定理做法:
对于一个侧面有 k^(nn) 种染色情况,也就是一个置换对象有 k^(nn) 种不同的情况。
因为正 m 边形,所以共有旋转 0 个位置、1个位置、2个位置……旋转 m-1 个位置 等 m 种置换方法。
假设位置 0 到 m-1.
对于旋转 k 个位置这种置换:第 pos 个位置 和 第 (pos+k*t)%m 个位置一样, t >= 0.
pos+gcd(k, m)*t1,其中 0 <= pos+gcd(k, m)*t1 < m,代表我们能到达的所有位置。
每一个轨迹中位置的个数 t = m/gcd(k, m)。由于在本题中不同轨迹中的位置的个数相同,所以轨迹的个数就是 m/t = gcd(k, m)。因此,旋转 k 个位置之后和原来相同的染色方案数就是 (k(n*n))gcd(k, m)
证明链接
代码
#include <bits/stdc++.h>
typedef long long ll;
using namespace std;
const int N=1e6+5;
const int maxn=1e5+5;
const int p=1e9+7;
ll quick(ll a,ll b)
{
ll sum=1;
while(b)
{
if(b&1)
{
sum=sum*a%p;
}
a=a*a%p;
b>>=1;
}
return sum;
}
int main()
{
int n,m,c;
scanf("%d %d %d",&n,&m,&c);
ll k=n*n;
k=quick(c,k);
ll ans=0;
for(int i=0;i<m;i++)
{
int s1=__gcd(m,i);
ans=(ans+quick(k,s1))%p;
}
ans=ans*quick(m,p-2)%p;
printf("%lld\n",(ans+p)%p);
}