《ybtoj金牌导航》第八部分第五章例题2 彩色项链2

题目

给定 $m$ 种颜色的珠子，每种颜色珠子的个数均不限，将这些珠子做成长度为 $n$ 的项链。

问能做成多少种不重复的项链。

两条项链相同，当且仅当两条项链通过旋转或是沿中轴线翻转后能重合在一起，且对应珠子的颜色相同。

多组数据
n<32,m<32

思路

polya计数定理。
先考虑旋转。
旋转共n种置换，第i种的轮换个数为gcd(n,i)
接下来是对称。
对称分奇偶性讨论
奇数部分有$n$种置换，轮换个数均为$(n+1)/2$。
偶数部分又分2部分
第一部分以中线对称，共$n/2$种，每种轮换均为$n/2+1$。
第二部分以对点连线对称，共$n/2$种，轮换数量均为$n/2$。
所以旋转的公式为
${\sum }_{i=1}^n{m}^{gcd(i,n)}$
奇数公式为
$n\ast m^{n/2+n\mathrm{m}\mathrm{o}\mathrm{d}2}$
偶数公式为
$n/2\ast (m^{n/2+1}+m^{n/2})$
记得最后除$2\ast n$（$|G|$）
code：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
//When I wrote this code,God and I unterstood what was I doing 
inline int read()
{
    int ret,c,f=1;
    while (((c=getchar())> '9'||c< '0')&&c!='-');
    if (c=='-') f=-1,ret=0;
    else ret=c-'0';
    while ((c=getchar())>='0'&&c<='9') ret=(ret<<1)+(ret<<3)+(c^48);
    return ret*f;
}
int n,m;
int ksm(int x,int y)
{
	int ans=1;
	while (y)
	{
		if (y&1) ans=ans*x;
		x=x*x;
		y>>=1;
	}
	return ans;
}
int gcd(int x,int y)
{
	int r=x%y;
	while (r)
	{
		x=y,y=r,r=x%y;
	}
	return y;
}
int main()
{
	n=read(),m=read();
	while (n!=0||m!=0)
	{
		int ans=0;
		for (int i=1;i<=m;i++) ans+=ksm(n,gcd(i,m));
		if (m&1) ans+=m*ksm(n,m/2+m%2);
		else ans+=m/2*(ksm(n,m/2+1)+ksm(n,m/2));
		cout<<ans/2/m<<endl;
		n=read(),m=read();
	}
	return 0;
}
//Now,only God know
//但行好事，莫问前程 

《论提前开金牌导航的感觉》