poj 2409 Let it Bead (组合数学 polya计数法)

大致题意：给你c种颜色的珠子，和一个长度为s(c,s<32)的项链，用这c种珠子串成这个项链，项链可以旋转和翻转，经过旋转和翻转所得的项链视为同一种项链，现在告诉你颜色总数c和项链的长度s，求共能组成几条不同的项链。

Time Limit : 2000/1000ms (Java/Other) Memory Limit : 20000/10000K (Java/Other)

　　　　　　　　　　　　　　　　　　　　　　


样例：

Sample Input

1 1
2 1
2 2
5 1
2 5
2 6
6 2
0 0

Sample Output

1
2
3
5
8
13
21

一看也是一个polya计数法的题，跟Problem 1286 Necklace of Beads基本是一个题，详细解释写在那个题了，用最简单的方法就能过，这里不再赘述了。

#include<cstdio>
#include<cmath>
using namespace std;
typedef __int64 ll;
int gcd(int a,int b){
if(b==0){
return a;
}
return gcd(b,a%b);
}
int c,s;
ll polya(){
int i,j;
ll ans=0;
for(i=0;i<s;i++){
ans+=(ll)pow(1.0*c,gcd(s,i));
}
if(s%2){
ans+=s*(ll)pow(1.0*c,s/2+1);
}
else{
ans+=s/2*(ll)pow(1.0*c,s/2);
ans+=s/2*(ll)pow(1.0*c,s/2+1);
}
return ans/2/s;
}
int main(){
while(scanf("%d%d",&c,&s)&&(c||s)){
ll ans=polya();
printf("%I64d\n",ans);
}
return 0;
}