Necklace of Beads

最新推荐文章于 2021-07-28 21:05:56 发布

「已注销」

最新推荐文章于 2021-07-28 21:05:56 发布

阅读量890

点赞数

CC 4.0 BY-SA版权

分类专栏： polya 组合数学*************** 文章标签： polya

本文链接：https://blog.youkuaiyun.com/bo_jwolf/article/details/12417939

组合数学*************** 同时被 2 个专栏收录

45 篇文章

订阅专栏

polya

3 篇文章

订阅专栏

本文提供了一种解决POJ 1286问题的C++实现方案，采用多项式计数原理解决图形染色问题，通过求最大公约数计算不同旋转情况下的颜色组合。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

http://poj.org/problem?id=1286

// File Name: poj1286.cpp
// Author: bo_jwolf
// Created Time: 2013年10月07日 星期一 21:31:08

#include<vector>
#include<list>
#include<map>
#include<set>
#include<deque>
#include<stack>
#include<bitset>
#include<algorithm>
#include<functional>
#include<numeric>
#include<utility>
#include<sstream>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<ctime>

#define INT long long int

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

INT c = 3, s;
INT polya(){
	INT ans = 0;
	for( INT i = 0; i < s; ++i ){
		ans += (INT)pow( 1.0 * c, gcd( s, i ) );
	}
	if( s % 2 ){
		ans += s * ( INT )pow( 1.0 * c, ( s / 2 + 1 ) );
	}
	else{
		ans += s / 2 * ( INT )pow( 1.0 * c, s / 2 );
		ans += s / 2 * ( INT )pow( 1.0 * c, s / 2 + 1 );
	}
	return ans / 2 / s;
}

int main(){
	INT ans;
	while( scanf( "%lld", &s ) != EOF ){
			if( s == -1 )
				break;
			if( s == 0 )
				ans = 0;
			else
				ans = polya();
			printf( "%lld\n", ans );
	}
return 0;
}