给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;
}