多校第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;}