poj2409Let it Bead Polay定理

题意：给出一个长度为m的项链，每个珠子可以用n种颜色涂色。翻转和旋转后相同的算作一种。有多少种不同的项链？

思路：

（1）

对于Burnside引理，G为所有置换集合，|G|为所有置换个数，gi为第i种置换，D(gi)为在第i种置换下保持不动的元素个数。

对于Polay定理，G为所有置换集合，|G|为所有置换个数，gi为第i种置换，n为颜色的种类，c(gi)为第i种置换的循环节个数。

（2）对于Polay，一般解题步骤：确定循环群；计算每种循环下的循环节个数。

对于本题，旋转可以分为旋转i=0,1,2，……m-1个，循环节个数为Gcd(m,i)。对于翻转，若m为奇数，只有通过一顶点和其对边中点的轴翻转，循环节m/2+1；对于m为偶数，两种翻转的轴，第一是两个顶点的连线，循环节m/2+1，第二是对边中点的连线,循环节m/2。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
#define ll long long
int gcd(int a,int b) {return a%b==0?b:gcd(b,a%b);}
ll quickpow(ll m,ll n)
{
    ll res=1;
    while(n)
    {
        if(n&1)
        res=res*m;
        n=n>>1;
        m=m*m;
    }
    return res;
}
int main()
{
    ll c,n;
    while(~scanf("%lld%lld",&c,&n))
    {
        if(c+n==0)
            break;
        ll ans=quickpow(c,n);
        for(int i=1;i<n;i++)
        ans+=quickpow(c,gcd(n,i));
        if(n&1) ans+=n*quickpow(c,n/2+1);
        else ans+=(n/2*quickpow(c,n/2+1)+n/2*(quickpow(c,n/2)));
        printf("%lld\n",ans/(2*n));
    }
    return 0;
}