POJ 1286 polya计数、burnside定理

最新推荐文章于 2020-05-27 00:40:03 发布

HAI__嗨I起来

最新推荐文章于 2020-05-27 00:40:03 发布

阅读量251

点赞数

分类专栏：组合数学

本文链接：https://blog.csdn.net/became_a_wolf/article/details/52292626

版权

组合数学专栏收录该内容

7 篇文章 0 订阅

订阅专栏

题目大意：

n个珠子串成一个圆，用三种颜色去涂色。问一共有多少种不同的涂色方法。

不同的涂色方法被定义为：如果这种涂色情况翻转，旋转不与其他情况相同就为不同。

解题思路：

Polya定理模版题。

burnside定理：对于一个置换f,若一个着色方案s经过置换后不变，称s为f的不动点。将f

的不动点为C(f).则等价类数目为所有C(f)的平均值。

对于顺时针长度为i的旋转，为pow(3,__gcd(n,i)；

对于翻转，当为奇数时，有：n*pow(3.0,(n-1)/2+1);

当为偶数时，有：n/2*pow(3.0,n/2-1+2)+n/2*pow(3.0,n/2);

一共有2*n种情况，最后要除以2*n

#include<iostream>  ///详解见白书147页
#include<cstdio>
#include<algorithm>
#include<cstring>
#define LL long long
using namespace std;
LL ans,n,m;
LL gcd(LL x,LL y)
{
    return y==0?x:gcd(y,x%y);
}
LL pow(LL num,LL x)
{
	LL ans=1;
	while (x){
		if (x&1)  ans=ans*num;
		x>>=1;
		num=num*num;
	}
	return ans;
}
int main()
{
	n=3;
	while (~scanf("%lld",&m),m!=-1){
		if (m==0){
			printf("0\n");
			continue;
		}
		ans=0;
		for (LL i=1;i<=m;i++)
		 ans+=pow(n,gcd(m,i));
		if (m%2)
		ans+=pow(n,(m-1)/2+1)*m;
		else{
            ans+=pow(n,m/2-1+2)*(m/2),
            ans+=pow(n,m/2)*(m/2);
		}

		printf("%lld\n",ans/(m*2));
	}
}