HDU 3923 Invoker

普通的ploy，套模版就行，只是有一个注意点，就是因为取模了，所以最后不能直接除2*n。

//
//  main.cpp
//  Richard
//
//  Created by 邵金杰 on 16/8/3.
//  Copyright © 2016年 邵金杰. All rights reserved.
//


#include<cstdio>
using namespace std;
typedef long long LL;
const LL maxn=1000000007;
LL gcd(LL a,LL b)
{
    return b==0?a:gcd(b,a%b);
}
LL Pow(LL a,LL b)
{
    LL result=1;
    LL base=a%maxn;
    while(b)
    {
        if(b&1) result=result*base%maxn;
        base=base*base%maxn;
        b>>=1;
    }
    return result;
}
LL polya(LL m,LL n)
{
    LL result=0;
    for(int i=1;i<=n;i++)
        result+=Pow(m,gcd(n,i));
    if(n%2)
        result+=Pow(m,(n+1)/2)*n;
    else
        result+=(Pow(m,n/2)+Pow(m,n/2+1))*n/2;
    return result%maxn*Pow(2*n,maxn-2)%maxn;
}
int main()
{
    LL m,n;
    int t,kase=0;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld%lld",&m,&n);
        printf("Case #%d: ",++kase);
        printf("%lld\n",polya(m,n));
    }
    return 0;
}