Polya定理染色问题 （poj 1286）

题目：http://poj.org/problem?id=1286

这个题目感觉是最裸的题目了很简单！

关于polya定理大家可以参考这里，讲的比较详细，我就说了。以下是poj1286的代码：

#include <stdio.h>
#include <iostream>
#include <math.h>
using namespace std;
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
long long rotate(int c,int n)
{
    int i;long long sum=0;
    for(i=1;i<=n;i++)
    {
        sum+=pow(c,gcd(i,n));
    }
    return sum;
}

long long turn(long long c,int n)
{
    long long sum;
    if(n%2)
        sum=n*pow(c,(n+1)/2);
    else    sum=n/2*((pow(c,n/2))+pow(c,(n+2)/2));
    return sum;
}
long long polya(int c,int n)
{
	if(n==0) return 0;
    long long sum=0;
    sum+=rotate(c,n);
    sum+=turn(c,n);
    n*=2;
    return sum/n;
}
int main()
{
    int n;
    while(scanf("%d",&n)&&n!=-1)
    {
        printf("%lld\n",polya(3,n));
    }
}