poj2409 Let it Bead（置换群+Polya定理+Burnside引理）

题目：

我是超链接

题解：

置换群+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种置换。