Polya 计数

/*==================================================*\
| Polya计数
| c种颜色的珠子, 组成长为s的项链, 项链没有方向和起始位置;
\*==================================================*/
// 注意超long long
int Polya(int c,int s) {
    int i,j,k,p[64], count; p[0]=1; // power of c
    for(k=1; k<=s ;k++) p[k]=p[k-1]*c;
    //reflection part
    count=(s&1)? (s*p[s/2+1]) : ((s/2)*(p[s/2]+p[s/2+1]));
    //rotation part
    for(k=1; k<=s ;k++) count+=p[gcd(k,s)];
    cout/=2*s;
    return count;
}