码迷,mamicode.com
首页 > 其他好文 > 详细

【Polya定理】poj1286 Necklace of Beads

时间:2017-04-08 00:50:38      阅读:207      评论:0      收藏:0      [点我收藏+]

标签:ret   循环   class   n+1   res   puts   type   algorithm   namespace   

Polya定理:设G={π1,π2,π3........πn}是X={a1,a2,a3.......an}上一个置换群,用m中颜色对X中的元素进行涂色,那么不同的涂色方案数为:1/|G|*(mC(π1)+mC(π2)+mC(π3)+...+mC(πk)). 其中C(πk)为置换πk的循环节的个数。

Polya定理的基础应用。

你得算出旋转和翻转时,每种置换的循环节数。

旋转时,每种置换的循环节数为gcd(n,i);

翻转时,若n为奇数,共有n个循环节数为n+1>>1的置换,

若n为偶数,共有n/2个循环节数为n+2>>1的置换和n/2个循环节数为n>>1的置换。

#include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long ll;
int n;
ll Pow(int x,int p){
	ll res=1;
	for(int i=1;i<=p;++i){
		res*=(ll)x;
	}
	return res;
}
int main(){
	while(1){
		scanf("%d",&n);
		if(n==-1){
			break;
		}
		if(n==0){
			puts("0");
			continue;
		}
		ll sum=0;
		for(int i=1;i<=n;++i){
			sum+=Pow(3,__gcd(n,i));
		}
		if(n&1){
			sum+=(ll)n*Pow(3,n+1>>1);
		}
		else{
			sum+=(ll)(n>>1)*Pow(3,n+2>>1);
			sum+=(ll)(n>>1)*Pow(3,n>>1);
		}
		cout<<sum/(2ll*(ll)n)<<endl;
	}
	return 0;
}

【Polya定理】poj1286 Necklace of Beads

标签:ret   循环   class   n+1   res   puts   type   algorithm   namespace   

原文地址:http://www.cnblogs.com/autsky-jadek/p/6680449.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!