hdu 3923 Invoker polya定理

题意：说的是dota里的召唤师的技能组合，玩过，嘿嘿。。有m种元素，现在选n个元素围成一圈组成一个技能(元素选择可以重复),

其中n个元素通过旋转或者翻转得到的一种已知排列，记为重复。问有多少种不重复组合方案。

题解：纯粹的polya定理题，polya定理的详细介绍请百度，这里简要介绍下：

先介绍循环数。举例来说一个圈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。

对于环的旋转来说，顺时针旋转i格，循环数为gcd(i,n)。

而翻转可分奇偶，奇数：对称轴必定通过一个点，所以由n种翻转情况。而循环则是对称轴左边和右边一一对应，共n/2+1个。偶数：也有n条对称轴，其中n/2是通过两个点，同样一一对应,n/2+1个；n/2条是通过两个的之间，一一对应n/2个。

而polya定理说的是，方案数ans=1/|G|*(∑元素数^旋转循环数+∑元素数^翻转循环数)。其中|G|表示变化的总数，旋转n种，翻转n种。

耗时：515MS

#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
#include <iostream>
#include <algorithm>
using namespace std;
typedef __int64 LL;
const LL mod=1e9+7;
LL gcd(LL a,LL b)//求最大公约数
{
    return b==0?a:gcd(b,a%b);
}
LL pow_mod(LL a,LL b)//快速幂求a^b
{
    LL s=1;
    while(b)
    {
        if(b&1)
            s=(s*a)%mod;
        a=(a*a)%mod;
        b=b>>1;
    }
    return s;
}
LL polya(LL m,LL n)//polya定理
{
    LL i,ans=0;
    for(i=1;i<=n;i++)
        ans+=pow_mod(m,gcd(i,n));//旋转的情况，循环数gcd(i,n)，其中i为顺时针旋转i格
    if(n&1)ans+=n*pow_mod(m,n/2+1);//奇数的翻转情况。共n条对称轴，每条的循环数均为n/2+1
    else ans+=n/2*pow_mod(m,n/2)+n/2*pow_mod(m,n/2+1);//偶数的翻转情况。对称轴共n条，n/2条通过2个点，其余n/2条通过两点之间的中心。
    ans=ans%mod*pow_mod(2*n,mod-2)%mod;
    return ans;
}
int main()
{
    LL T,m,n,tt=0;
    cin>>T;
    while(T--)
    {
        cin>>m>>n;
        cout<<"Case #"<<++tt<<": "<<polya(m,n)<<endl;
    }
    return 0;
}