poj 1286 Necklace of Beads polya计数

Polya定理

设

是n个对象的一个置换群, 用m种颜色染图这n个对象，则不同的染色方案数为：

其中 ， 为 的循环节数。

首先考虑 翻转：
        如果有偶数个珠子：
                翻转的轴可以是两个珠子的连线（一共n/2个这样的旋转轴)，这样处于轴上的两个珠子都是1阶循环，其他的n-2个珠子可以通过再一次翻转回到原来发样子，一共有（n-2）/2对这样的珠子。
               如果翻转轴在两边的两个珠子的连线的连线上（一共n/2个这样的旋转轴)，那么每个珠子翻转后再经过依稀翻转就能回原位置，所以一共有n/2个循环节。

         如果有奇数个珠子，那么翻转的轴在一个珠子和整个图形的几何中心的连线处。处于翻转轴上的珠子是一阶循环，其他珠子可以再次翻转回来，所以是（n-1)/2个2阶循环。这样的旋转轴一共有n个。

再来考虑 旋转，旋转可分为转过0个珠子，转过一个珠子，转过两个珠子……转过n-1个珠子（转过n个珠子和转过0个效果一样）。
对于转过i个珠子，我们考虑一个珠子a，它每次要转i个单位，那么lcm(n,i)个单位，也就是经过lcm(n,i)/i次旋转才能回到原位置。也就是说对于每个循环是一个lcm(n,i)/i阶的，那么n个数一共有n/lcm(n,i)/i个这样的循环，也就是gcd(n,i)个。

#include <iostream>
#include<stdio.h>
#include<cmath>
using namespace std;

int gcd(int a,int b)
{
    return b == 0?a:gcd(b,a%b);
}

long long polya(int n)
{
    long long sum = 0;
    for(int i = 0;i < n;i++) //旋转
        sum+=pow(3.0,gcd(n,i));
    if(n&1) //奇数个珠子
        sum+=n*pow(3.0,(n-1)/2+1); //翻转轴在几何中心和一个点的连线上
    else
    {
        sum+=n/2*pow(3.0,n/2);//翻转轴在两边连线的中点
        sum+=n/2*pow(3.0,(n-2)/2+2);//翻转轴在两点的连线上
    }
    return sum/(2*n);
}

int main()
{
    int n;
    while(scanf("%d",&n),n!=-1)
    {
        if(n == 0) printf("0\n");
        else printf("%lld\n",polya(n));
    }
    return 0;
}