polya定理主要就是解决一类着色问题,或者说是同构计数问题。对Ω的n个目标用m种颜色着色的图象集种类数。
定理的内容:
Pólya定理:设G={p1,p2,…,pg}是Ω上的一个置换群,C(pk)
是置换pk的循环的个数,用M中的颜色对Ω中的元素着色,
着色方案数为 l = 1/|G|*[c1(p1)+c1(p2)+c1(p3)+...c1(p[g])]=1/|G|*[m^c(p1)+m^c(p2)+m^c(p3)+...m^c(p[g])]
定理的内容:
Pólya定理:设G={p1,p2,…,pg}是Ω上的一个置换群,C(pk)
是置换pk的循环的个数,用M中的颜色对Ω中的元素着色,
着色方案数为 l = 1/|G|*[c1(p1)+c1(p2)+c1(p3)+...c1(p[g])]=1/|G|*[m^c(p1)+m^c(p2)+m^c(p3)+...m^c(p[g])]
|G|为置换的总个数,m颜色数c1(pi)指置换pi的不动点的数目(既循环节为1的点数)c(pi)指的是置换pi的循环个数。
题意:一家项链公司生产手镯。n颗珠子形成一个环,用m种颜色给n颗珠子染色,就得到了各种各样的手镯。但是,经过旋转和翻转使之吻合的算同一种方案。
例如,当用2种颜色对5颗珠子进行染色的方案数为8,
解:显然,对于这n个对象,有n种旋转和n种翻转。
1. 对于旋转,有c(gi) = gcd(n,i),i为转动几颗珠子。
2. 对于翻转,当n为奇数时,c(gi) = n/2+1;
当n为偶数时,有n/2个的循环节数为n/2+1,有n/2个的循环节数为n/2。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
typedef __int64 int64;
int gcd(int a, int b)
{
return b==0?a:gcd(b,a%b);
}
//用m种颜色涂均匀分布在圆环上的n颗珠子的方案数
int64 polya_circle(double m, int n)
{
int64 result = 0;
//旋转
for(int i = 1; i <= n; i++)
result += pow(m,gcd(n,i));
//翻转
if(n%2)
result += n*pow(m,n/2+1);
else
result += (pow(m,n/2)+pow(m,n/2+1))*n/2;
return result/n/2;
}
int main()
{
int c,s;
while(true)
{
scanf("%d %d",&c,&s);
if(c==0&&s==0)
break;
printf("%I64d\n",polya_circle(c,s));
}
return 0;
}
转自:
http://yzmduncan.iteye.com/blog/1402942
http://zhulinb123.blog.163.com/blog/static/184414043201182351033633/