题目:
题解:
置换群+Polya定理+Burnside引理模板一类喽
在这里说一下这三个布吉岛是什么东西的东西:
置换:
1换成a1,2换成a2,3换成a3,则(a1,a2,a3)为1-3的排列
置换可以看成一个个轮换
[2,5,4,3,1]=(1,2,5)(3,4)(两个轮换)
[2,5,1,3,4]=(1,2,5,4,3)
轮换的数量称为置换的循环节
置换群:
可以用来描述等价关系,方案A通过群中的某个置换能够变成方案B,则A与B等价
Burnside引理:
对于一个置换f,若方案s变换到自己,称它为f的不动点
f的不动点数目记为C(f),等价类数目为所有C(f)的平均值
Polya定理:
给n个元素着k种颜色,设置换的循环节为m(f),则它的不动点数C(f)=k^m(f)
综上所述,这类问题只需要计算置换f的循环节,这需要O(n)
朴素代码:
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
using namespace std;
int a[120][120];
void make(int s)
{
int i,j;
for (i=1;i<=s;i++)
a[1][i]=i;
for (i=2;i<=s;i++)
{
int cnt=0;
for (j=i;j<=s;j++)
a[i][++cnt]=j;
for (j=1;j<=i-1;j++)
a[i][++cnt]=j;
}
a[s+1][1]=1;
int cnt=1;
for (i=s;i>=2;i--)
a[s+1][++cnt]=i;
for (i=2;i<=s;i++)
{
int cnt=0;
for (j=i;j<=s;j++)
a[i+s][++cnt]=a[s+1][j];
for (j=1;j<=i-1;j++)
a[i+s][++cnt]=a[s+1][j];
}
}
int jie(int j,int s)
{
int ans=0,now;
bool p[100];
memset(p,0,sizeof(p));
for (int i=1;i<=s;i++)
if (!p[i])
{
ans++;now=a[1][i];
while (!p[now]) p[now]=true,now=a[j][now];
}
return ans;
}
int answer(int c,int s)
{
int l=s*2;
int ans=0;
for (int i=1;i<=l;i++)
ans+=pow(c,jie(i,s));
return ans/l;
}
int main()
{
int c,s,ans;
scanf("%d%d",&c,&s);
while (c!=0 && s!=0)//c种颜色,s个珠子
{
make(s);
ans=answer(c,s);
printf("%d\n",ans);
scanf("%d%d",&c,&s);
}
}
很强的代码:
#include<iostream>
#include<stdio.h>
using namespace std;
typedef long long LL;
int n, m;
int gcd(int a, int b) {
if (b == 0)return a;
return gcd(b, a % b);
}
LL pow(LL a, LL b) {
LL ans = 1;
while (b--)ans *= a;
return ans;
}
int main() {
int i, j;
while (scanf("%d%d", &m, &n) && n + m) {
LL ans = 0;
for (i = 0; i < n; i++)
ans = ans + pow((LL) m, (LL) gcd(n, i));
if (n & 1) {
ans = ans + n * pow((LL) m, (LL) n / 2 + 1);
} else
ans = ans + n / 2 * (pow((LL) m, (LL) n / 2) + pow((LL) m, (LL) n / 2 + 1));
printf("%lld\n", ans / 2 / n);
}
}
分两种情况讨论:
旋转:
n种旋转方法每种旋转i个格(1<=i<=n)循环结有gcd(i,n)个
翻转:
(1)这种是经过某个顶点i与中心的连线为轴的翻转,由于n为偶数,有对称性,所以此种共n/2种翻转:
(2)这种是以顶点i和i+1的连线的中点与中心的连线为轴的翻转,同样,根据对称性,也有n/2种翻转:
所以给定长度n,共有2n种置换。