HDU --- 3923 Invoker 【poyla定理】

传送门
//这道题就是一个poyla定理的裸题, 题目描述可以转换成用n种颜色去染一个长度为m的项链, 问一共有多少种不同的染法, 其中可以通过旋转和翻转后形成一样的算一种方法.

// 解法: Poyla定理就不详细说了, 这里简单说说. 定理就一句话 :
方案数 res =1/|G|*(∑元素数^旋转循环数+∑元素数^翻转循环数)。其中|G|表示变化的总数. 在这道题中|G| 是2m种旋转m种，翻转m种 . 分开说.

//旋转 : 举个例子 一个圈: 1->2->3->4->5->6,顺序时针旋转2格,5->6->1->2->3->4。其中1到3的位置，3到5的位置, 5到1的位置，成一个循环. 总的有(1,3,5)(2,4,6)两个循环, 所以循环数为2 .
所以我们可以知道对于环长为m的旋转来说, 顺时针旋转i格, 循环数为gcd(i,m) .

//翻转 : 要分m的奇偶, 当m为奇数时, 恰好有m条对称轴经过某一点, 比如说一个三角形, 三个顶点分别为1，2，3，随便假设一点 1 在对称轴上, 那么绕这个轴翻转, (2,3)(1)就是两个循环节. 所以一条就可以产生m/2+1的循环数. 所以一共是m*(m/2+1).

当m为偶数时, 比如说正方形, 1，2，3，4, 那么对称轴就要两种情况, 第一种对角线, (假设是1，3)那么就会产生(2,4)(1)(3)的循环数, 一共两条, 即m/2*(m/2+1)，如果是以对边的中点连线作为对称轴, 则会产生(1,2)(3,4)的循环数, 还是有2条, 即m/2*(m/2).

//分析到这, 就可以敲啦

AC Code

int cas=1;
ll qpow(ll x,ll y)
{
    ll res = 1;
    while(y){
        if(y&1) res = res*x%mod;
        x = x*x%mod;
        y >>= 1;
    }
    return res;
}
void solve()
{
    int n,m;
    scanf("%d%d",&n,&m);
    ll res = 0;
    for(int i=1;i<=m;i++) res += qpow(n,__gcd(i,m));  //旋转的情况
    if(m%2) res += m*qpow(n,m/2+1);    //奇数翻转的情况
    else res += m/2*qpow(n,m/2) + m/2*qpow(n,m/2+1);  //偶数翻转的情况.
    res = res%mod*qpow(2*m,mod-2)%mod; //用逆元
    printf("Case #%d: ", cas++);
    printf("%lld\n",res);
}