poj 2409 polya定理

#include 
#include 
#include 
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种