#include <stdio.h>
#include <string.h>
#include <math.h>
double s;
int gcd(int a,int b) //最大公约数,用来求旋转循环数
{
if(b==0)
return a;
return gcd(b,a%b);
}
void folya(int m,int n)
{
int i;
for(i=1;i<=n;i++)
s+=pow(1.0*m,gcd(n,i)); //旋转循环数,顺时针旋转i格
if(n%2==1)
s+=n*pow(1.0*m,n/2+1); //翻转,点个数为奇数时,对称线经过一个点
else
s+=n/2*pow(1.0*m,n/2)+n/2*pow(1.0*m,n/2+1); //翻转,点数为偶数,分为由经过点和没经过点
s=s/(2*n);
}
int main()
{
int m,n;
while(scanf("%d%d",&m,&n)!=EOF)
{
if(n==0&&m==0)
break;
s=0;
folya(m,n);
printf("%d\n",(int)s);
}
return 0;
}
// 循环数,以6格点为例,记点为 1,2,3,4,5,6
// 翻转1:以1,2中间点翻转:得到 2,1,6,5,4,3
// 其中1->2,2->1成一个循环,共三个
// 翻转2:经过1,4的对称线翻转: 1,6,5,4,3,2
//循环为(1)(2,6)(3,5)(4)四个
//旋转:最大公约数就是循环数
// 最后本质不同数(即旋转和翻转后一样的都不要)为1/|G|*(∑颜色数^旋转i循环数+∑颜色数^翻转循环数)
//其中|G|为变换的总数,旋转n种,翻转n种
poj 2409 polya定理
最新推荐文章于 2017-11-11 11:50:39 发布