poj 1286 Necklace of Beads【polya定理+burnside引理】

和poj 2409差不多，就是k变成3了，详见
还有不一样的地方是记得特判n==0的情况不然会RE

#include<iostream>
#include<cstdio>
using namespace std;
long long n,ans;
long long ksm(long long a,long long b)
{
long long r=1;
while(b)
{
if(b&1)
r=r*a;
a=a*a;
b>>=1;
}
return r;
}
long long gcd(long long a,long long b)
{
return !b?a:gcd(b,a%b);
}
int main()
{
while(scanf("%lld",&n)&&n!=-1)
{
if(!n)
{
puts("0");
continue;
}
ans=(n&1)?ksm(3,n/2+1)*n:ksm(3,n/2+1)*n/2+ksm(3,n/2)*n/2;
for(int i=1;i<=n;i++)
ans+=ksm(3,gcd(i,n));
printf("%lld\n",ans/2/n);
}
return 0;
}