poj 1286 Necklace of Beads(Polya定理)

题意：用3种颜色给正n边形涂色，问有多少不等价的方案，如过一个方案能通过某一个方案旋转或者轴对称得到，那么被认为相同的。

思路：看了好几天Polya，终于明白了一点点。这题不是很难，基本的Polya计数方法，通过Burnside定理，可以求出不等价的着色数，利用Polya计数算出使得着色通过置换群G中置换保持不变的着色的个数，这个需要计算置换的循环节。对于对称置换，可以分奇偶讨论求出公式，对于旋转置换计数，需要dfs搜一下循环节个数。刚开始以为旋转置换也是有规律的，但是算了好几个还是不对，没办法上网看了下题解，貌似这东西没规律。。搜了下就过了……

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
const int maxn=30;
ll p[maxn];
int c[maxn],a[maxn],cnt;
bool vis[maxn];
void dfs(int u)
{
    if(vis[u]) return ;
    vis[u]=true;
    dfs(a[u]);
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    p[0]=1;
    for(int i=1;i<maxn;++i) p[i]=p[i-1]*3;
    int n;
    while(cin>>n)
    {
        if(n==-1) break;
        if(n==0)
        {
            cout<<0<<endl;
            continue;
        }
        memset(c,0,sizeof(c));
        for(int i=1;i<=n;++i) a[i]=i;
        for(int i=1;i<=n;++i)
        {
            memset(vis,0,sizeof(vis));
            cnt=0;
            for(int j=1;j<=n;++j)
            {
                if(!vis[j])
                {
                    cnt++;
                    dfs(j);
                }
            }
            c[cnt]++;
            int tmp=a[1];
            for(int j=1;j<n;++j) a[j]=a[j+1];
            a[n]=tmp;
        }
        if(n&1)
        {
            c[(n-1)/2+1]+=n;
        }
        else
        {
            c[n/2]+=n/2;
            c[n/2+1]+=n/2;
        }
        ll ans=0;
        for(int i=1;i<=n;++i)
            ans+=c[i]*p[i];
        ans/=2*n;
        cout<<ans<<endl;
    }
    return 0;
}