POJ 1286 Necklace of Beads

人生第一道polya(玻璃呀!)定理。

搞了一个上午的群论,唯一看懂的就是伽罗华很会作……

你说你一个搞数学的,怎么就喜欢革命呢。

唯一的遗憾就是过分追求简洁hhhhhhhhhh这评价真是十分中肯啊。

好了不扯了,反正玻璃呀定理就是用来求染色数的,具体的请任意翻开一本组合数学教材,或者跟我一样(安利一下)找学堂在线的课看去。

显然我们可以发现对序列(1,2,3,......n)做i次轮换出来的置换可以被分解为gcd(i,n)个循环,于是这里套一下Polya定理,然后对于翻转来说的话奇偶讨论一下就好了

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
ll mul[30];
int gcd(int a,int b){
	return b?gcd(b,a%b):a;
}
int main(){
	//freopen("a.in","r",stdin);
	mul[0]=1;
	for(int i=1;i<=28;i++)mul[i]=mul[i-1]*3;
	int n;
	while(scanf("%d",&n)&&n>=0){
		if(!n){
			puts("0");
			continue;
		}
		ll sum=0,cnt=n;
		for(int i=1;i<=n;i++){
			sum+=mul[gcd(i,n)];
		}
		if(n&1){
			sum+=mul[(n+1)/2]*n;
			cnt+=n;
		}else{
			sum+=(mul[n/2]+mul[n/2+1])*n/2;
			cnt+=n;
		}
		printf("%lld\n",sum/cnt);
	}
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值