多校第9场HDU3923Invoker（polya定理）

给n个颜色,m个点的环，求有多少种通过旋转，翻转变化后不同的染色方法

利用polya定理

用大数取模，一直WA，改成数组预处理就AC了。。。

#include <cstdio>

const long long mm=1000000007;

int I=0;
long long bigmod(long long  b,long long  p,long long  m)
{
   long long k;
   for(k=1,b=b%m; p ; p/=2,b=b*b%m) 
    if(p%2) k=k*b%m;
     //printf("%lld\n", k%m);
     return k%m;
}

long long  gcd (long long  a,long long  b)
{
     return b?gcd(b,a%b):a;
}

int main ()
{
    int cas;
    scanf("%d",&cas);
    while (cas--)
    {
          int  n,m;
          long long  ans=0;
          scanf("%d%d",&n,&m);
          long long t=2*m*mm;//
          long long pre[10020];
          pre[0]=1;
          for (int i=0 ; i<m ; ++i)
           pre[i+1]=(pre[i]*n)%t;
          for (int i=1 ; i<=m ; ++i )//旋转 
          {
              //ans=(ans+bigmod(n,gcd(i,m),t))%t;
              ans=(ans+pre[gcd(i,m)])%t;
          }
          //翻转 分奇偶 
          if(m&1)
          {
              //ans=(ans+bigmod(n,m/2+1,t)%t*m%t)%t;
              ans=(ans+pre[m/2+1]*m)%t;
          }
          else 
          {
               //ans=(ans+(bigmod(n,m/2,t)+bigmod(n,m/2+1,t))*(m/2))%t;
               ans=(ans+(pre[m/2]+pre[m/2+1])*m/2)%t;
          }
          //printf("%lld  %lld\n",ans,t);
          printf("Case #%d: %lld\n",++I,(ans/2/m)%mm);
    }
    return 0;
}