poj 2409 polya定理

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

你可能感兴趣的:(Polya)