Polya定理
设
是n个对象的一个置换群, 用m种颜色染图这n个对象,则不同的染色方案数为:
![](https://i-blog.csdnimg.cn/blog_migrate/a0941da8171d2f6a2c7af9c526893a21.png)
是n个对象的一个置换群, 用m种颜色染图这n个对象,则不同的染色方案数为:
![](https://i-blog.csdnimg.cn/blog_migrate/05970dbb2787ff50fc36ec0a523449bd.png)
其中
,
为
的循环节数。
首先考虑 翻转:
如果有偶数个珠子:
翻转的轴可以是两个珠子的连线(一共n/2个这样的旋转轴),这样处于轴上的两个珠子都是1阶循环,其他的n-2个珠子可以通过再一次翻转回到原来发样子,一共有(n-2)/2对这样的珠子。
如果翻转轴在两边的两个珠子的连线的连线上(一共n/2个这样的旋转轴),那么每个珠子翻转后再经过依稀翻转就能回原位置,所以一共有n/2个循环节。
如果有奇数个珠子,那么翻转的轴在一个珠子和整个图形的几何中心的连线处。处于翻转轴上的珠子是一阶循环,其他珠子可以再次翻转回来,所以是(n-1)/2个2阶循环。这样的旋转轴一共有n个。
再来考虑 旋转,旋转可分为转过0个珠子,转过一个珠子,转过两个珠子……转过n-1个珠子(转过n个珠子和转过0个效果一样)。
对于转过i个珠子,我们考虑一个珠子a,它每次要转i个单位,那么lcm(n,i)个单位,也就是经过lcm(n,i)/i次旋转才能回到原位置。也就是说对于每个循环是一个lcm(n,i)/i阶的,那么n个数一共有n/lcm(n,i)/i个这样的循环,也就是gcd(n,i)个。
![](https://i-blog.csdnimg.cn/blog_migrate/0b59b89f184414dcbf2f632db2246017.png)
![](https://i-blog.csdnimg.cn/blog_migrate/caa51530324a8fca44e09a0074236b3f.png)
![](https://i-blog.csdnimg.cn/blog_migrate/53025cc35c55b72216b9d02313151bdc.png)
首先考虑 翻转:
如果有偶数个珠子:
翻转的轴可以是两个珠子的连线(一共n/2个这样的旋转轴),这样处于轴上的两个珠子都是1阶循环,其他的n-2个珠子可以通过再一次翻转回到原来发样子,一共有(n-2)/2对这样的珠子。
如果翻转轴在两边的两个珠子的连线的连线上(一共n/2个这样的旋转轴),那么每个珠子翻转后再经过依稀翻转就能回原位置,所以一共有n/2个循环节。
如果有奇数个珠子,那么翻转的轴在一个珠子和整个图形的几何中心的连线处。处于翻转轴上的珠子是一阶循环,其他珠子可以再次翻转回来,所以是(n-1)/2个2阶循环。这样的旋转轴一共有n个。
再来考虑 旋转,旋转可分为转过0个珠子,转过一个珠子,转过两个珠子……转过n-1个珠子(转过n个珠子和转过0个效果一样)。
对于转过i个珠子,我们考虑一个珠子a,它每次要转i个单位,那么lcm(n,i)个单位,也就是经过lcm(n,i)/i次旋转才能回到原位置。也就是说对于每个循环是一个lcm(n,i)/i阶的,那么n个数一共有n/lcm(n,i)/i个这样的循环,也就是gcd(n,i)个。
#include <iostream>
#include<stdio.h>
#include<cmath>
using namespace std;
int gcd(int a,int b)
{
return b == 0?a:gcd(b,a%b);
}
long long polya(int n)
{
long long sum = 0;
for(int i = 0;i < n;i++) //旋转
sum+=pow(3.0,gcd(n,i));
if(n&1) //奇数个珠子
sum+=n*pow(3.0,(n-1)/2+1); //翻转轴在几何中心和一个点的连线上
else
{
sum+=n/2*pow(3.0,n/2);//翻转轴在两边连线的中点
sum+=n/2*pow(3.0,(n-2)/2+2);//翻转轴在两点的连线上
}
return sum/(2*n);
}
int main()
{
int n;
while(scanf("%d",&n),n!=-1)
{
if(n == 0) printf("0\n");
else printf("%lld\n",polya(n));
}
return 0;
}