Polya定理染色问题 （poj 1286）

#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));
}
}