poj 1286 Necklace of Beads 项链旋转+对称群的polya计数

最新推荐文章于 2018-03-27 11:33:39 发布

sepNINE

最新推荐文章于 2018-03-27 11:33:39 发布

阅读量400

点赞数

CC 4.0 BY-SA版权

分类专栏：数学文章标签： poj 算法

本文链接：https://blog.youkuaiyun.com/sepNINE/article/details/49779257

数学专栏收录该内容

46 篇文章

订阅专栏

本文提供了一个解决POJ1286问题的C++实现方案，通过预计算3的幂次和利用辗转相除法求最大公约数来高效解决该数学问题。

//poj 1286
//sep9
#include<iostream>
using namespace std;
typedef __int64 INT;
INT pow3[32];

INT gcd(INT a,INT b)
{
	return a%b==0?b:gcd(b,a%b);
}

int main()
{
	pow3[0]=1;
	for(INT i=1;i<24;++i)
		pow3[i]=pow3[i-1]*3;
	INT n;
	while(scanf("%I64d",&n)==1&&n>=0){
		if(n==0){
			printf("0\n");
			continue;
		}
		INT ans=n&1?n*pow3[n/2+1]:(n/2)*pow3[n/2]+(n/2)*pow3[n/2+1];
		for(INT i=0;i<n;++i)
			ans+=pow3[gcd(i,n)];	
		printf("%I64d\n",ans/(2*n));
	}
	return 0;
}