UVA 11255 Necklace（Burnside引理）

题目链接：http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2222

题意：有三种不同颜色的珠子，白（W颗），灰（G颗），黑（B颗）， 其中3<=W+G+B<=40，问用这些珍珠可以串成多少种(W+G+B)长的不同的项链（项链可在平面内转动，也可翻转）。

思路：令n=W+G+B,（1）旋转：想想，对于每一种置换，我们要找出在这种置换下哪些排列在置换后还是这样。旋转有n种情况，分别是旋转0，1,2……，n-1个。比如，旋转2，设n=12,那么所有的分成2组，每组6个。这两组要保持一样，那么WGB在每组中的个数和排列方式是一样的，此时必须WGB都是偶数才行，进而问题转化成W/2,G/2,B/2个在一组中的排列的个数有多少种（另一组跟这组的排列一样无需再计算），这就是一个简单的组合计算。因此，枚举n的约数i,k=n/Gcd(n,i)，WGB都是k的倍数才可以，然后就是W/k、G/k、B/k的排列；（2）翻转：WGB中有0个奇数，绕对边的轴转和绕相对的顶点的轴转；有1个奇数，绕一个顶点一条边中点的轴转；有2个奇数，设WG为奇数，绕两个相对顶点的轴转，两个顶点分别是WG。

void init()
{
int i,j;
C[0][0]=1;
FOR1(i,40)
{
C[i][0]=C[i][i]=1;
FOR1(j,i-1) C[i][j]=C[i-1][j]+C[i-1][j-1];
}
}

int Gcd(int x,int y)
{
return !y?x:Gcd(y,x%y);
}

int x,y,z,n,T;

int main()
{
init();
RD(T);
while(T--)
{
RD(x,y,z);
n=x+y+z;
int i,k;
u64 ans=0;
FOR1(i,n)
{
k=n/Gcd(n,i);
if(x%k==0&&y%k==0&&z%k==0)
{
ans+=C[n/k][x/k]*C[(y+z)/k][y/k];
}
}
int cnt=x%2+y%2+z%2;
if(cnt==1)
{
if(x&1) ans+=C[n/2][(x-1)/2]*C[(y+z)/2][y/2]*n;
else if(y&1) ans+=C[n/2][(y-1)/2]*C[(x+z)/2][x/2]*n;
else ans+=C[n/2][(z-1)/2]*C[(x+y)/2][x/2]*n;
}
else if(cnt==0)
{
ans+=n/2*C[n/2][x/2]*C[(y+z)/2][y/2];
if(x>=2) ans+=n/2*C[n/2-1][x/2-1]*C[(y+z)/2][y/2];
if(y>=2) ans+=n/2*C[n/2-1][y/2-1]*C[(x+z)/2][x/2];
if(z>=2) ans+=n/2*C[n/2-1][z/2-1]*C[(x+y)/2][x/2];
}
else if(cnt==2)
{
if(x%2==0) ans+=n*C[n/2-1][x/2]*C[(y+z-2)/2][(y-1)/2];
else if(y%2==0) ans+=n*C[n/2-1][y/2]*C[(x+z-2)/2][(x-1)/2];
else ans+=n*C[n/2-1][z/2]*C[(x+y-2)/2][(x-1)/2];
}
PR(ans/(2*n));
}
return 0;
}